Complementary Relationship of evaporation (CR) theory has gained attention in recent years, in part because it relies solely on standard weather data for estimating evaporation, eliminating the need for land surface moisture and resistance information. CR models show varied skill across diverse landscapes and climates, however, differences in model parameterization, calibration, station data, and study locations have precluded the ability to consistently and systematically compare model skill. This study intercompares four widely used CR models, uncalibrated and calibrated, with eddy covariance data from 82 sites worldwide, and assesses how environmental conditions affect model skill and calibration. Eddy covariance data from over 200 FLUXNET and AmeriFLUX stations were quality assured and controlled, filtered, and energy balance closure corrected using open-source software and benchmarking procedures for reproducibility. Systematic intercomparison showed that, overall, the Rescaled Power function (Szilagyi et al., 2022) had the highest skill followed by Sigmoid (Han & Tian, 2018), while Polynomial model (Brutsaert, 2015) generally biased high, and Advection-Aridity (Brutsaert & Stricker, 1979) often produced unrealistic negative values when observed evaporation was low. Calibration greatly improved model skill, and single parameter models yielded similar skill, highlighting the potential to reduce model non-uniqueness. Calibration and model skill rankings were consistent using energy balance closure corrected and uncorrected datasets. Results indicate that environmental conditions such as aridity index, relative humidity, vapor pressure deficit and the ratio between radiative and apparent potential evaporation can explain spatial patterns of CR model parameters, revealing opportunities for future development of improved, calibration-free CR evaporation modeling.