Search NASA⌕ Search

SEARCH · Search NASA

Results for “surface energy balance”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3

Persistent urban heat

Urban surface and near-surface air temperatures are known to be often higher than their rural counterparts, a phenomenon now labeled as the urban heat island effect. However, whether the elevated urban temperatures are more persistent than rural temperatures at timescales commensurate to heat waves has not been addressed despite its importance for human health. Combining numerical simulations by a global climate model with a surface energy balance theory, it is demonstrated here that urban surface and near-surface air temperatures are significantly more persistent than their rural counterparts in cities dominated by impervious materials with large thermal inertia. Further use of these materials will result in even stronger urban temperature persistence, especially for tropical cities. The present findings help pinpoint mitigation strategies that can simultaneously ameliorate the larger magnitude and stronger persistence of urban temperatures.

54 ENVIRONMENTAL SCIENCES↗

The winter central Arctic surface energy budget: A model evaluation using observations from the MOSAiC campaign

This study evaluates the simulation of wintertime (15 October, 2019, to 15 March, 2020) statistics of the central Arctic near-surface atmosphere and surface energy budget observed during the MOSAiC campaign with short-term forecasts from 7 state-of-the-art operational and experimental forecast systems. Five of these systems are fully coupled ocean-sea ice-atmosphere models. Forecast systems need to simultaneously simulate the impact of radiative effects, turbulence, and precipitation processes on the surface energy budget and near-surface atmospheric conditions in order to produce useful forecasts of the Arctic system. This study focuses on processes unique to the Arctic, such as, the representation of liquid-bearing clouds at cold temperatures and the representation of a persistent stable boundary layer. It is found that contemporary models still struggle to maintain liquid water in clouds at cold temperatures. Given the simple balance between net longwave radiation, sensible heat flux, and conductive ground flux in the wintertime Arctic surface energy balance, a bias in one of these components manifests as a compensating bias in other terms. This study highlights the different manifestations of model bias and the potential implications on other terms. Three general types of challenges are found within the models evaluated: representing the radiative impact of clouds, representing the interaction of atmospheric heat fluxes with sub-surface fluxes (i.e., snow and ice properties), and representing the relationship between stability and turbulent heat fluxes.

54 ENVIRONMENTAL SCIENCES↗

Representing Fine‐Scale Topographic Effects on Surface Radiation Balance in Hyper‐Resolution Land Surface Models

Land surface models are increasingly used to simulate land surface processes at hyper-spatial resolutions (e.g., ∼1 km). As model resolution increases, grid-scale topographic effects on surface radiation fluxes and their interactions between adjacent grids become more pronounced. However, current land surface models routinely neglect the fine-scale topographic effects on surface radiation balance. This study developed physically-based and computationally-efficient parameterizations (fineTOP) that explicitly resolve fine-scale topographic effects on downward shortwave and longwave radiation as well as land surface radiative properties. The newly developed parameterizations were implemented and tested in the Energy Exascale Earth System Model (E3SM) Land Model (ELM). Multi-decadal km-resolution ELM simulations over the California Sierra Nevada show that fine-scale topography significantly impacts the surface energy balance and snow processes across seasons. Slope determines the magnitude of topographic effects, while aspect controls their sign. For slopes larger than 30°, topography-induced change in annual surface temperature can be as large as 3.3 K. Regionally, the mean value and standard deviation of topography-induced changes in annual surface temperature are −0.22 ± 0.38 K and +0.25 ± 0.37 K over north-facing and south-facing slopes, respectively. Topography-induced changes in surface radiative properties account for 3.5% ± 13.8% of total topographic effects on annual net radiation. With fineTOP, ELM captures the aspect-dependence of snow cover fraction, snow water equivalent, and land surface temperature found in MODIS satellite observations and a snow reanalysis data set, while the default ELM fails to capture this phenomenon. The enhanced capability to represent fine-scale topographic effects on surface radiation balance can be used to advance understanding of the role of fine-scale topography in land surface processes and land-atmosphere interactions over mountainous regions.

Hao, Dalei [Pacific Northwest National Laboratory ↗

Evapotranspiration partitioning estimates from 8 methods from 47 NEON sites, 2019-2021

This dataset provides daily estimates of evapotranspiration (ET) and the transpiration-to-evapotranspiration ratio (T/ET) across 47 terrestrial National Ecological Observatory Network (NEON) sites spanning diverse environmental and biome conditions in the United States across three years of data (2019-2021). Daily ET is reported in both energy units (MJ m⁻² day⁻¹) and equivalent water depth (mm day⁻¹), assuming a constant latent heat of vaporization of 2.45 MJ/kg. The primary method uses a hybrid recurrent neural network–Penman–Monteith framework (RNN-PM), which integrates physically based surface energy balance constraints with data-driven learning to partition ET into transpiration and evaporation components. Model inputs include in situ meteorological observations (air temperature, vapor pressure deficit, wind speed, and radiation) combined with satellite-derived land surface temperature, leaf area index, and soil moisture. For benchmarking and uncertainty assessment, T/ET estimates from seven additional models are included: Priestley-Taylor Jet Propulsion Laboratory (PT-JPL), Penman-Monteith (P-M), Two-Source Energy Balance (TSEB), Support Vector Regression (SVR), and Categorical Boosting (CatBoost), among others—spanning empirical, machine-learning, and process-based approaches (see methods section or linked publication for detailed descriptions). Data Package Contents: The dataset a csv files containing daily ET and T/ET estimates for each site and model, along with associated metadata files these variables. Data can be accessed using common spreadsheet software (e.g., Microsoft Excel, LibreOffice) or programming environments such as R or Python. Together, these data support cross-site comparisons of ecosystem water use, evaluation of ET partitioning methods, and development of improved land–atmosphere exchange models.

EARTH SCIENCE > ATMOSPHERE↗

Evaluating integrated surface/subsurface permafrost thermal hydrology models in ATS (v0.88) against observations from a polygonal tundra site

Numerical simulations are essential tools for understanding the complex hydrologic response of Arctic regions to a warming climate. However, strong coupling among thermal and hydrological processes on the surface and in the subsurface and the significant role that subtle variations in surface topography have in regulating flow direction and surface storage lead to significant uncertainties. Careful model evaluation against field observations is thus important to build confidence.We evaluate the integrated surface/subsurface permafrost thermal hydrology models in the Advanced Terrestrial Simulator (ATS) against field observations from polygonal tundra at the Barrow Environmental Observatory. ATS couples a multiphase, 3D representation of subsurface thermal hydrology with representations of overland nonisothermal flows, snow processes, and surface energy balance. We simulated thermal hydrology of a 3D ice-wedge polygon with geometry that is abstracted but broadly consistent with the surface microtopography at our study site. The simulations were forced by meteorological data and observed water table elevations in ice-wedge polygon troughs. With limited calibration of parameters appearing in the soil evaporation model, the 3-year simulations agreed reasonably well with snow depth, summer water table elevations in the polygon center, and high-frequency soil temperature measurements at several depths in the trough, rim, and center of the polygon. Upscaled evaporation is in good agreement with flux tower observations. The simulations were found to be sensitive to parameters in the bare soil evaporation model, snowpack, and the lateral saturated hydraulic conductivity.Timing of fall freeze-up was found to be sensitive to initial snow density, illustrating the importance of including snow aging effects. The study provides new support for an emerging class of integrated surface/subsurface permafrost simulators.

54 ENVIRONMENTAL SCIENCES↗

Climate Change Impacts on the South American Monsoon System and Its Surface–Atmosphere Processes Through RegCM4 CORDEX-CORE Projections

Here, we evaluate projected changes in surface water and energy balances and surface–atmosphere coupling in the South American Monsoon System (SAMS) for the end of the century (2080–2099). The analyses are based on two ensemble datasets, which follow Representative Concentration Pathway 8.5 in the future period, and cover four subdomains (Northern and Southern—NAMZ and SAMZ—Amazon, La Plata Basin—LPB, and Southern Southeast Brazil—SSB). One ensemble consists of three Global Climate Models (HadGEM2-ES, MPI-ESM-MR and NorESM1-M), while the other consists of their dynamically downscaled version at 25 km horizontal grid spacing using Regional Climate Model version 4 (RegCM4). As both ensembles are able in reproducing the annual cycle of the components of the surface water and energy balances in the present climate, they can be used in the study of future climate. During the wet season (November–March), both ensembles project a decrease in precipitation over NAMZ and SAMZ (an exception is RegCM4 that projects a slight increase in SAMZ), and an increase across the LPB and SSB. These changes do not cause retreat or expansion of the monsoon area over the continent, which is similar to the present climate (1995–2014). For the wet season, the ensembles are in line with the presence of a strong surface–atmosphere coupling in LPB and SSB, weak coupling in SAMZ and very weak coupling in NAMZ. For future climate, the coupling is even weaker in NAMZ, which may be a driver for the negative changes in precipitation. For the other subdomains, while the ensembles project similar signals of precipitation changes, they disagree with the surface-atmosphere coupling highlighting the uncertainties in future climate.

54 ENVIRONMENTAL SCIENCES↗

Constraining Physical Understanding of Aerosol Loading, Biogeochemistry, and Snowmelt Hydrology from Hillslope to Watershed Scale in the East River Scientific Focus Area

The mountain snowpack is a critical component of regional hydrology, ecology, biogeochemistry, and climate in the Western US. This project leveraged the East River Scientific Focus Area (SFA), as an outdoor laboratory to address a significant gap in our understanding of the mountain snowpack; namely, how atmospheric constituent deposition on snowpack influences snow energy balance and nutrient/chemical cycling, and how snowmelt timing and intensity exerts controls on emergent biogeochemical and ecohydrologic behavior. This aligned with East River SFA goals to integrate landscape scale measurements and physical based modeling tools to improve understanding of controls on runoff production, ecohydrology, biogeochemical cycling, and land surface energy partitioning in high mountain watersheds. The project included observations of snowpack at multiple scales, measurements of snow and deposited aerosols properties, in situ time series of surface energy balance, and spatially distributed, process-based snowmelt modeling.

54 ENVIRONMENTAL SCIENCES↗

Strong influence of convective heat transfer efficiency on the cooling benefits of green roof irrigation

Abstract By enhancing evapotranspiration (ET), green roofs provide cooling benefits for the urban environment and are recognized as a promising heat mitigation strategy. The evaporative cooling effects of green roofs strongly depend on the soil moisture conditions and thus irrigation may be needed to sustain the cooling benefits. It has been shown that the magnitude of cooling benefits offered by green roof irrigation varies spatially, but its controlling factors remain elusive. In this study, we combine a surface energy balance (SEB) model with global simulations generated by an improved Earth System Model to illustrate the key factors controlling the cooling benefits of green roof irrigation. We employ a simple irrigation scheme, which is only active when there is no ice in soil layers and when the soil moisture is below field capacity. As a result, most of the irrigation water leaves the green roof system via ET. We find that the magnitude and also the spatial variability of the cooling benefits of green roof irrigation are controlled by the irrigation amount, and a surface energy redistribution factor that encodes the efficiencies of different SEB components in transferring heat. Further analysis indicates that the enhancement of latent heat flux due to irrigation is largely balanced by the reduction of sensible heat flux on green roofs. Therefore, the amount of irrigation needed per unit decrease of green roof surface temperature is mainly controlled by the convective heat transfer efficiency. A lower convective heat transfer efficiency (e.g. under a lower wind speed) helps reduce the amount of irrigation needed per unit decrease of green roof surface temperature. This study highlights the importance of SEB in constraining the cooling benefits of green roof irrigation and provides valuable guidance for urban planners and policy makers in terms of heat mitigation and sustainable water management.

Wang, Linying (ORCID:0000000308769220)↗

Land-Surface Temperature Transitions at the Bankhead National Forest (LASTT-BNF) Field Campaign Report

Land surface types that are common in the Southeast United States (e.g., forests, agriculture, urban) can play key roles in the regulation of surface-atmosphere energy exchange and boundary-layer evolution in the region (Hinkle et al. 2024). Heterogeneity in these land-surface types and associated emissivities can drive strong transitions in land surface temperature (LST), which is an important component of the surface energy balance. Within the forested canopy surrounding the U.S. Department of Energy Atmospheric Radiation Measurement (ARM) User Facility’s field deployment at the Bankhead National Forest (BNF) in Alabama main site (M1) (Kuang et al. 2026), the different components of the vegetation and ground surface can create strong gradients in surface temperature that can influence the canopy thermal environment and vegetative function (e.g., stress, transpiration), both of which can have important, yet poorly characterized, controls on lower atmospheric processes (e.g., thermal turbulence, secondary circulations).

54 ENVIRONMENTAL SCIENCES↗

Modulation of springtime surface sensible heating over the Tibetan Plateau on the interannual variability of East Asian dust cycle

Previous observational evidence and numerical simulations have revealed that the surface sensible heating in spring (March–April–May, MAM) over the Tibetan Plateau (TPSH) can affect the Asian regional hydrological cycle, surface energy balance, and climate through altering atmospheric heat source of the Tibetan Plateau (TP). This study aims to investigate the impacts of MAM TPSH on the interannual variability of East Asian dust cycle by using CAM4-BAM (version 4 of the Community Atmosphere Model coupled to a bulk aerosol model), MERRA-2 (version 2 of the Modern-Era Retrospective analysis for Research and Applications) surface dust concentration, and TPSH measurements. Our simulations show that the surface dust concentrations over the East Asian (EA) dust source region and over the northwestern Pacific (NP) in MAM are significantly positively correlated with TPSH, with regionally averaged correlation coefficients of 0.49 for EA and 0.44 for NP. Similar positive correlations are also shown between the MAM TPSH measurements averaged over the 73 observation sites and the surface dust concentration from MERRA-2. Simulation-based comparisons between strongest and weakest TPSH years reveal that, the MAM surface dust concentration in the strongest TPSH years increases with relative differences of 13.1 % over EA and 36.9 % over NP. These corresponding differences are found in MERRA-2 with 22.9 % and 13.3 % over EA and NP, respectively. Further simulated results show that the processes of whole dust cycles (e.g., dust loading, emission, and transport, as well as dust deposition) are also significantly enhanced during the strongest TPSH years over EA and NP. Through enhancing the TP heat source, stronger TPSH in MAM generates an anticyclonic anomaly in middle and upper troposphere over the TP and over the downstream Pacific region, respectively. These atmospheric circulation anomalies induced by the increased TPSH result in increasing the westerly winds over both EA and NP, which in turn increases dust emissions over the dust source, and dust transport over these two regions, as well as the regional dust cycles. These results suggest that addressing the East Asian dust changes in the future climates requires considering not only increasing greenhouse gas emissions but also the variations of the TP's heat source under global warming.

54 ENVIRONMENTAL SCIENCES↗

An Improved Zhang's Dynamic Water Balance Model Using Budyko-Based Snow Representation for Better Streamflow Predictions

Understanding the water balance of a catchment in relation to its regional climate forcings and catchment characteristics is critical for predicting current and future water resources amid changing climate and land cover. This study intends to improve Zhang's monthly water balance model (a physics-based conceptual hydrologic model) that reflects the physical partitioning process of the hydrological cycle at the basin level based on regional climate and catchment characteristics. The existing model does not include snow process and has confronted evident limitations in snow-affected areas, which is a critical aspect since snowmelt water has been a significant source of water resources for many regions, especially in the temperate and frigid zones. We introduce a snow module based on surface energy balance and Budyko-limits on melting and combine it with the existing water balance equations. Moreover, monthly parameterization is applied to the model to better explain the time-varying hydrological characteristics of a catchment. The proposed model involves five different monthly parameters, which determine the physical partitioning process of the hydrological cycle, and they are regionally calibrated and validated under Budyko-type constraints. The model is applied to 1,210 basins across the continental United States (CONUS), and the simulated streamflow is compared to the observed data. The proposed model significantly outperformed the original model, improving the median NSE by 31% (from 0.51 to 0.67) and increasing the number of catchments with an acceptable NSE by 58%. The spatial variability of the basin characteristics across the CONUS is also investigated based on the calibrated parameters.

54 ENVIRONMENTAL SCIENCES↗

Albedo-induced global warming impact of Conservation Reserve Program grasslands converted to annual and perennial bioenergy crops

Abstract Climate benefit assessments of bioenergy crops often focus on biogeochemical impacts, paying little if any attention to biogeophysical impacts. However, land conversions required for large-scale bioenergy crop production are substantial and may directly affect the climate by altering surface energy balance. In the US, such land conversions are likely to be met in part by converting Conservation Reserve Program (CRP) grassland to bioenergy crops. Here, we converted three 22 year old CRP smooth brome grass fields into no-till corn, switchgrass, or restored prairie bioenergy crops. We assessed the biogeophysical climate impact of the conversions using albedo changes relative to unconverted reference CRP grassland. The corn and perennial fields had higher annual albedo than the grassland they replaced—causing cooling of the local climate. The cooling of the corn field occurred solely during the non-growing season—especially when surfaces were snow-covered, whereas the cooling of the perennial fields was more prominent during the growing season. Compared to biogeochemical impacts with fossil fuel offsets for the same land conversions over eight years, the annual albedo-induced climate benefits add ∼35% and ∼78% to the annual biogeochemical benefits provided from the switchgrass and restored prairie fields, respectively, and offset ∼3.3% of the annual greenhouse gas (GHG) emissions from the corn field. We conclude that albedo-induced climate mitigation from conversion of CRP lands to perennial but not annual bioenergy crops can be substantial, and future climate impact assessments of bioenergy crops should include albedo changes in addition to GHG balances in order to better inform climate policies.

54 ENVIRONMENTAL SCIENCES↗

Ground Heat Flux Reconstruction Using Bayesian Uncertainty Quantification Machinery and Surrogate Modeling

Ground heat flux (G 0 ) 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, G 0 is 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 reconstruct G 0 across 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 available G 0 data (measured or modeled) for snow-free period as a reference. When observed G 0 is 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 developed G 0 reconstruction approach provides novel means for assessing the probabilistic structure of the ground heat flux for regional permafrost change studies.

54 ENVIRONMENTAL SCIENCES↗

Cutting out the middleman: calibrating and validating a dynamic vegetation model (ED2-PROSPECT5) using remotely sensed surface reflectance

Canopy radiative transfer is the primary mechanism by which models relate vegetation composition and state to the surface energy balance, which is important to light- and temperature-sensitive plant processes as well as understanding land–atmosphere feedbacks. In addition, certain parameters (e.g., specific leaf area, SLA) that have an outsized influence on vegetation model behavior can be constrained by observations of shortwave reflectance, thus reducing model predictive uncertainty. Importantly, calibrating against radiative transfer outputs allows models to directly use remote sensing reflectance products without relying on highly derived products (such as MODIS leaf area index) whose assumptions may be incompatible with the target vegetation model and whose uncertainties are usually not well quantified. Here, we created the EDR model by coupling the two-stream representation of canopy radiative transfer in the Ecosystem Demography model version 2 (ED2) with a leaf radiative transfer model (PROSPECT-5) and a simple soil reflectance model to predict full-range, high-spectral-resolution surface reflectance that is dependent on the underlying ED2 model state. We then calibrated this model against estimates of hemispherical reflectance (corrected for directional effects) from the NASA Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) and survey data from 54 temperate forest plots in the northeastern United States. The calibration significantly reduced uncertainty in model parameters related to leaf biochemistry and morphology and canopy structure for five plant functional types.

54 ENVIRONMENTAL SCIENCES↗

Daily evapotranspiration changes during heatwaves at 32 NEON sites, 2019-2021

This dataset provides partitioned evapotranspiration (ET, the combined loss of water from soil and plant surfaces) anomalies during heatwave events—soil evaporation (E) and transpiration (T)—for 268 heatwave events across 32 National Ecological Observatory Network (NEON) flux sites in the contiguous United States from 2019–2021. Using an ensemble of four high-frequency turbulence methods (Flux-variance Similarity, Conditional Eddy Covariance [CEC], CEC with Water-Use Efficiency, and Conditional Eddy Accumulation; see Zahn and Bou-Zeid 2024), half-hourly transpiration-to-evapotranspiration (T/ET) ratios were derived from 20 hertz (Hz, cycles per second) eddy covariance measurements of carbon dioxide (CO₂) and water vapor (H₂O) concentrations. The dataset spans six vegetation types including evergreen and deciduous forests, grasslands, cultivated crops, shrublands, and emergent herbaceous wetlands. Data Package Contents: The dataset includes a single CSV (comma-separated values) file containing daily anomalies (deviations from baseline conditions) for transpiration (Delta_T), evaporation (Delta_E), total evapotranspiration (Delta_ET), and T/ET ratio (Delta_T_ET) during each day of identified heatwave events. The file also includes site codes, dates, heatwave event identifiers, and day-of-heatwave indicators. The CSV file can be opened with spreadsheet software (Microsoft Excel, Google Sheets) or programming environments (Python, R, MATLAB). This resource enables researchers to investigate ecosystem-specific responses to thermal extremes, validate land surface model partitioning of ET fluxes, and examine feedbacks between water cycling and surface energy balance during heatwaves. The dataset is particularly valuable for studies linking vegetation hydraulic strategies to climate resilience, as it captures the divergent responses of shallow-rooted versus deep-rooted ecosystems. Potential applications include improving drought early warning systems, informing irrigation management strategies, and advancing our mechanistic understanding of land-atmosphere interactions under extreme heat conditions.

Day of Heatwave↗

What lies beneath: Vertical temperature heterogeneity in a Mediterranean woodland savanna

As the availability of satellite and airborne thermal infrared remote sensing (TIR-RS) data increases and their spatial, temporal, and spectral resolutions improve, researchers are finding diverse applications for TIR-RS measurements. TIR-RS is now commonly applied in regional- and continental-scale analyses, such as those focused on fire and surface energy balance. However, its application lags in plant physiology and ecology, for which a finer-scale understanding of plant canopy temperatures would be useful to elucidate plant water dynamics, for example. In particular, while methods to disaggregate TIR-RS pixels in horizontal space have advanced, possible vertical stratification of plant canopy temperature and its implications for understanding the correspondence between TIR-RS and finer-scale, field-based thermal measurements (e.g. made with a thermal camera) remain unexplored. Here, we use data from a thermal camera deployed concurrently with the recent ECOSTRESS mission to quantify vertical temperature gradients within tree canopies and temperatures of over- vs. under-story plants in a Mediterranean woodland savanna. We then leverage diverse ancillary data to maximize the geometric comparability of ECOSTRESS and thermal camera measurements, in order to assess the extent to which the two forms of thermal measurements correspond. Specifically, we ask: (1) What are the patterns of intra-canopy and over- vs. under-story vertical temperature in a Mediterranean woodland savanna?, and (2) How can vertically-resolved, but spatially-limited field-based temperature measurements be reconciled with spatially-extensive, but surface-only, temperature measurements of a space-borne remote sensor? Further, we found consistent patterns of vertical thermal heterogeneity both within tree canopies and between ecosystem over- and under-stories. The daytime difference between the top and bottom thirds of blue oak canopies was, on average, 0.48 ° C – and sometimes several times larger. Notably, canopy tops are cooler, likely associated with the under-story grass reaching daytime temperatures often exceeding over-story temperatures by 10° C. Given the consistency of the intra-canopy temperature gradients, we expected the ECOSTRESS sensor would be in better agreement with camera measurements of canopy tops than bulk canopies or canopy bottoms. However, within-canopy gradients were overwhelmed by other sources of disagreement between the measurements, in part associated with upscaling camera measurements across space. Overall, thermal camera and ECOSTRESS measurements were largely in agreement at night (pixel RMSE = 1.1°C), but they were more divergent during times of low (but >0 W/m2) and high incoming solar radiation (daytime pixel RMSE = 3.5°C).

54 ENVIRONMENTAL SCIENCES↗

Physical Controls on the Hydrology of Perennially Ice-Covered Lakes, Taylor Valley, Antarctica (1996–2013)

The McMurdo Dry Valleys, Antarctica, are a polar desert populated with numerous closed-watershed, perennially ice-covered lakes primarily fed by glacial melt. Lake levels have varied by as much as 8 m since 1972 and are currently rising after a decade of decreasing. Precipitation falls as snow, so lake hydrology is dominated by energy available to melt glacier ice and to sublimate lake ice. To understand the energy and hydrologic controls on lake level changes and to explain the variability between neighboring lakes, only a few kilometers apart, we model the hydrology for the three largest lakes in Taylor Valley. Here, we apply a physically based hydrological model that includes a surface energy balance model to estimate glacial melt and lake sublimation to constrain mass fluxes to and from the lakes. Results show that lake levels are very sensitive to small changes in glacier albedo, air temperature, and wind speed. We were able to balance the hydrologic budget in two watersheds using meltwater inflow and sublimation loss from the ice-covered lake alone. A third watershed, closest to the coast, required additional inflow beyond model uncertainties. We hypothesize a shallow groundwater system within the active layer, fed by dispersed snow patches, contributes 23% of the inflow to this watershed. The lakes are out of equilibrium with the current climate. If the climate of our study period (1996–2013) persists into the future, the lakes will reach equilibrium starting in 2300, with levels 2–17 m higher, depending on the lake, relative to the 2020 level.

58 GEOSCIENCES↗

Disentangling the Relative Drivers of Seasonal Evapotranspiration Across a Continental-Scale Aridity Gradient

Evapotranspiration (ET) is a significant ecosystem flux, governing the partitioning of energy at the land surface. Understanding the seasonal pattern and magnitude of ET is critical for anticipating a range of ecosystem impacts, including drought, heat-wave events, and plant mortality. Here, in this study, we identified the relative controls of seasonal variability in ET, and how these controls vary among ecosystems. We used overlapping AmeriFlux and PhenoCam time series at a daily timestep from 20 sites to explore these linkages (# site-years >100), and our study area covered a broad climatological aridity gradient in the U.S. and Canada. We focused on disentangling the most important controls of bulk surface conductance (G s ) and evaporative fraction (EF = LE/[H + LE]), where LE and H represent latent and sensible heat fluxes, respectively. Specifically, we investigated how vegetation phenology varied in importance relative to meteorological variables (vapor pressure deficit and antecedent precipitation) as a driver of G s and EF using path analysis, a framework for quantifying and comparing the causal linkages among multiple response and explanatory variables. Our results revealed that the drivers of G s and EF seasonality varied significantly between energy- and water-limited ecosystems. Specifically, precipitation had a much higher effect in water-limited ecosystems, while seasonal patterns in canopy greenness emerged as a stronger control in energy-limited ecosystems. Given that phenology is expected to shift under future climate, our findings provide key information for understanding and predicting how phenology may impact 21st-century hydroclimate regimes and the surface-energy balance.

54 ENVIRONMENTAL SCIENCES↗