skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


This content will become publicly available on September 15, 2026

Title: Global Lake Evaporation Estimates by Integrating Penman Method with Equilibrium Temperature Approach
Abstract Modeling evaporationEfrom inland water bodies is challenging largely due to the uncertainties of input data, particularly surface water temperature that plays a key role in the available energy, i.e., net radiationRnminus rate of water heat storage changeG. The equilibrium temperature approach (ETA) for estimating water surface temperature offers an alternative method to calculateRnandGusing standard meteorological data. This study evaluates the global lakeEestimates from the widely used Penman model (PM) coupled with the ETA (PM-ETA) against field observations and model simulations from the Lake, Ice, Snow, and Sediment Simulator (LISSS). Our analysis reveals that the PM-ETA tends to overestimateEby approximately 36% and 24% compared to observations and the LISSS simulations, respectively, despite being driven by the same input data. The biases of the PM-ETAEare more pronounced in the cold and polar regions with distinct seasonality ofRnandG. Furthermore, theEtrends from the PM-ETA deviate from the LISSS simulations over the period of 2001–16 due to the bias trends in the available energy. By incorporating the LISSS-simulatedRnandGinto the PM, the bias inEis reduced to less than ±5% compared to the LISSS results. This study highlights the need to improve the available energy input of the PM to improve the estimates of global lakeEfor better water resource management and planning. Significance StatementThis study addresses a crucial challenge in modeling evaporationEfrom inland water bodies—uncertainties in surface water temperature and available energy inputs, particularly net radiationRnand rate of heat storage changeG. By evaluating the widely used Penman model (PM) coupled with the equilibrium temperature approach (ETA), we reveal a tendency for the PM-ETA to overestimateEglobally, with the largest biases observed in cold and polar regions. Incorporating higher-qualityRnandGestimates from the Lake, Ice, Snow, and Sediment Simulator (LISSS) significantly reduces these biases. These findings highlight the importance of alternative higher-quality data products for available energy inputs for improvingEestimates by the PM.  more » « less
Award ID(s):
2006281 2003076
PAR ID:
10636037
Author(s) / Creator(s):
; ; ; ;
Publisher / Repository:
American Meteorological Society
Date Published:
Journal Name:
Journal of Hydrometeorology
Volume:
26
Issue:
9
ISSN:
1525-755X
Page Range / eLocation ID:
1301 to 1313
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Ground heat flux (G0) is a key component of the land‐surface energy balance of high‐latitude regions. Despite its crucial role in controlling permafrost degradation due to global warming,G0is sparsely measured and not well represented in the outputs of global scale model simulation. In this study, an analytical heat transfer model is tested to reconstructG0across seasons using soil temperature series from field measurements, Global Climate Model, and climate reanalysis outputs. The probability density functions of ground heat flux and of model parameters are inferred using availableG0data (measured or modeled) for snow‐free period as a reference. When observedG0is not available, a numerical model is applied using estimates of surface heat flux (dependent on parameters) as the top boundary condition. These estimates (and thus the corresponding parameters) are verified by comparing the distributions of simulated and measured soil temperature at several depths. Aided by state‐of‐the‐art uncertainty quantification methods, the developedG0reconstruction approach provides novel means for assessing the probabilistic structure of the ground heat flux for regional permafrost change studies. 
    more » « less
  2. The instrumented raft on Sparkling Lake is equipped with a dissolved oxygen and CO2 sensors, a thermistor chain, and meteorological sensors that provide fundamental information on lake thermal structure, weather conditions, evaporation rates, and lake metabolism. Estimating the flux of solutes to and from lakes requires accurate water budgets. Evaporation rates are a critical component of the water budget of lakes. Data from the instrumented raft on Sparkling Lake includes micrometeorological parameters from which evaporation can be calculated. Raft measurements of relative humidity and air temperature (2m height), wind velocity (2m) ,and water temperatures (from thermistors placed throughout the water column at intervals varying from 0.5 to 3m) are combined with measurements of total long-wave and short-wave radiation data from a nearby shore station to determine evaporation by the energy budget technique. Comparable evaporation estimates from mass transfer techniques are calibrated against energy budget estimates to produce a lake-specific mass transfer coefficient for use in estimating evaporation rates. After correcting for flux to or from the atmosphere and vertical mixing within the water column, high frequency measurements of dissolved gases such as carbon dioxide and oxygen can be used to estimate gross primary productivity, respiration, and net ecosystem productivity, the basic components of whole lake metabolism. Other parameters measured include precipitation, wind direction (beginning in 2008), and barometric pressure (beginning in 2008). Sampling Frequency: one minute with hourly and daily averages provided. Number of sites: 1. 
    more » « less
  3. Abstract The southern Lake Michigan region of the United States, home to Chicago, Milwaukee, and other densely populated Midwestern cities, frequently experiences high pollutant episodes with unevenly distributed exposure and health burdens. Using the two‐way coupled Weather Research Forecast and Community Multiscale Air Quality Model (WRF‐CMAQ), we investigate criteria pollutants over a southern Lake Michigan domain using 1.3 and 4 km resolution hindcast simulations. We assess WRF‐CMAQ's performance using data from the National Climatic Data Center and Environmental Protection Agency Air Quality System. Our 1.3 km simulation slightly improves on the 4 km simulation's meteorological and chemical performance while also resolving key details in areas of high exposure and impact, that is, urban environments. At 1.3 km, we find that most air quality‐relevant meteorological components of WRF‐CMAQ perform at or above community benchmarks. WRF‐CMAQ's chemical performance also largely meets community standards, with substantial nuance depending on the performance metric and component assessed. For example, hourly simulated NO2and O3are highly correlated with observations (r > 0.6) while PM2.5is less so (r = 0.4). Similarly, hourly simulated NO2and PM2.5have low biases (<10%), whereas O3biases are larger (>30%). Simulated spatial pollutant patterns show distinct urban‐rural footprints, with urban NO2and PM2.520%–60% higher than rural, and urban O36% lower. We use our 1.3 km simulations to resolve high‐pollution areas within individual urban neighborhoods and characterize seasonal changes in O3regimes across tight spatial gradients. Our findings demonstrate both the benefits and limitations of high‐resolution simulations, particularly over urban settings. 
    more » « less
  4. Abstract Utilizing 956 nights of Na lidar nocturnal mesopause region temperature profiles acquired at Fort Collins, CO (40.6°N, 105.1°W) over a 20‐year period (March 1990–2010), we deduce background nightly mean temperatureand the square of the buoyancy frequencyN2(z) at 2‐km resolution between 83 and 105 km. The temperature climatology reveals the two‐level mesopause structure with clarity and sharp mesopause transitions, resulting in 102 days of summer from Days 121 to 222 of the year. The same data set analyzed at 10‐min and 1‐km resolution gives the gravity wave (GW) temperature perturbationsTi'(z) and the wave varianceVar(T′(z)) and GW potential energyEpm(z) between 85 and 100 km. Seasonal averages of GWVar(T′(z)) andEpm(z) between 90 and 100 km, show thatVar(T′) for spring and autumn are comparable and lower than for summer and winter. Due mainly to the higher background stability, or largerN2(z) in summer,Epm(z) between 85 and 100 km is comparable in spring, summer, and autumn seasons, but ∼30%–45% smaller than the winter values at the same altitude. The uncertainties are about 4% for winter and about 5% for the other three seasons. The values forEpmare (156.0, 176.2, 145.6, and 186.2 J/kg) at 85 km for (spring, summer, autumn, and winter) respectively, (125.4, 120.2, 115.2, and 168.7 J/kg) at 93 km, and (207.5, 180.5, 213.1, and 278.6 J/kg) at 100 km. Going up in altitude, all profiles first decrease and then increase, suggesting that climatologically, GWs break below 85 km. 
    more » « less
  5. Abstract Heat storage change (HSC) is a crucial component of lake's thermal energy budget. Conventional temperature profile based models of HSC require location specific parameters such as lakebed topography. Based on the half‐order time‐derivative formula of heat fluxes, an analytical model was formulated for estimating HSC from water surface temperature and solar radiation without using geography dependent parameters. The proposed model was tested against field measurements at Poyang Lake, a shallow inland lake, which has pronounced seasonal variations in water level and lake area. Our analysis indicates that the model accurately simulates diurnal HSC with a coefficient of determination of 0.94 and a root mean squared error (RMSE) of 77.5 ± 21.6 Wm−2for the study period. Larger nighttime RMSE (75.0 ± 26.8 Wm−2) than the daytime value (55.1 ± 19.7 W m−2) is attributable to larger measurement errors of nighttime turbulent fluxes. The estimation of HSC independent of temperature profile and lake‐specific parameters by the proposed model facilitates remote sensing monitoring the HSC of global water bodies. 
    more » « less