Key Points:• ECOSTRESS is a state-of-the-art combination of thermal bands, spatial and temporal resolutions, and measurement accuracy and precision • Data from 82 eddy covariance sites were coalesced concurrently with the first year of ECOSTRESS for Stage 1 validation • Clear-sky ET from ECOSTRESS compared well against a wide range of eddy Abstract The ECOsystem Spaceborne Thermal Radiometer Experiment on Space Station (ECOSTRESS) was launched to the International Space Station on 29 June 2018 by the National Aeronautics and Space Administration (NASA). The primary science focus of ECOSTRESS is centered on evapotranspiration (ET), which is produced as Level-3 (L3) latent heat flux (LE) data products. These data are generated from the Level-2 land surface temperature and emissivity product (L2_LSTE), in conjunction with ancillary surface and atmospheric data. Here, we provide the first validation (Stage 1, preliminary) of the global ECOSTRESS clear-sky ET product (L3_ET_PT-JPL, Version 6.0) against LE measurements at 82 eddy covariance sites around the world. Overall, the ECOSTRESS ET product performs well against the site measurements (clear-sky instantaneous/time of overpass: r 2 = 0.88; overall bias = 8%; normalized root-mean-square error, RMSE = 6%). ET uncertainty was generally consistent across climate zones, biome types, and times of day (ECOSTRESS samples the diurnal cycle), though temperate sites are overrepresented. The 70-m-high spatial resolution of ECOSTRESS improved correlations by 85%, and RMSE by 62%, relative to 1-km pixels. This paper serves as a reference for the ECOSTRESS L3 ET accuracy and Stage 1 validation status for subsequent science that follows using these data.

There are no direct methods to evaluate calculated soil heat flux (SHF) at the surface (G0). Instead, validation and cross evaluation of methods for calculating G0 usually rely on the conventional calorimetric method or the degree of the surface energy balance closure. However, there is uncertainty in the calorimetric method itself, and factors apart from G0 also contribute to nonclosure of the surface energy balance. Here we used a novel approach to evaluate nine different methods for calculating SHF, including the calorimetric method and methods based on analytical solutions of the heat diffusion equation. The SHF (Gz) measured by a self‐calibrating SHF plate at a depth of z = 5 cm below the surface (hereafter Gm_5cm) was deployed as a reference. Each SHF calculation method was assessed by comparing the calculated Gz at the same depth (hereafter Gc_5cm) with Gm_5cm. The calorimetric method and simple measurement method performed best in determining Gc_5cm but still underestimated Gm_5cm by 19% during the daytime. Possible causes for this underestimation include errors and uncertainties in SHF measurements and soil thermal properties, as well as the phase lag between Gc_5cm and Gm_5cm. Our results indicate that the calorimetric method achieves the most accurate SHF estimates if self‐calibrating SHF plates are deployed at two depths (e.g., 5 cm and 10 cm), soil temperature and water content measurements are made in a few depths between the two plates, and soil thermal properties are accurately quantified.

The surface energy balance non-closure problem in eddy covariance (EC) studies has been largely attributed to the influence of large-scale turbulent eddies (hereafter large eddies) on latent and sensible heat fluxes. However, how such large eddies concurrently affect CO 2 fluxes remains less studied and mechanistic links between the energy balance non-closure and CO 2 fluxes are not well understood. Here, using turbulence data collected from an EC tower over a sagebrush ecosystem during two growing seasons, we decomposed the turbulence data into small and large eddies at a cutoff frequency and analyzed their contributions to the fluxes. We found that the magnitude of CO 2 fluxes decreased concurrently with decreased sensible and latent heat fluxes (and thus increased energy balance nonclosure), primarily caused by large turbulent eddies. The contributions of such large eddies to fluxes are dependent not only upon their magnitudes of vertical velocity (w) and scalars (i.e. temperature, water vapor density, and CO 2 concentration), but also upon the phase differences between the large eddies of w and scalars via their covariances. Enlarged phase differences between large eddies of w and these scalars simultaneously led to reductions in the magnitudes of both CO 2 and heat fluxes, linking the lower CO 2 fluxes to energy balance non-closure. Such increased phase differences of large eddies were caused by changes in the structures of large eddies from unstable to near neutral conditions. Given widespread observations in non-closure in the flux community, the processes identified here may bias CO 2 fluxes at many sites and cause upscaled regional and global budgets to be underestimated. More studies are needed to investigate how landscape heterogeneity influences CO 2 fluxes through the influence of associated large eddies on flux exchange.

Semiarid ecosystems play a critical role in determining the interannual variability of the global terrestrial carbon sink. Water availability is a critical driver of productivity in semiarid ecosystems, which often alternate between carbon sink/source functioning during wet/dry years. In this study, we investigate how groundwater availability resulting from groundwater‐river water exchange influences net ecosystem exchange of CO2 (NEE), evapotranspiration (ET), and the surface energy balance at two semiarid ecosystems along the Columbia River in central Washington, USA. We examined 1 year of eddy covariance measurements from an upland sagebrush ecosystem primarily fed by rainfall without groundwater access and a riparian grassland ecosystem with groundwater access during the dry season due to lateral groundwater‐river water exchange. The two sites had distinct seasonal patterns of NEE and ET, driven by differences in water availability. While NEE at the upland sagebrush site was strongly constrained by water availability during the dry months, access to groundwater allowed the riparian site to maintain high NEE magnitude and ET during the same dry months. The riparian site had larger annual gross primary productivity than the upland site (612 vs. 424 gC/m2), which was offset by higher ecosystem respiration (558 vs. 363 gC/m2). Thus, the magnitude of the annual NEE at the upland site was larger than that at the riparian site (−62 vs. −54 gC/m2). Our results demonstrate that groundwater access determined by connectivity between groundwater and surface water can be a critical driver of carbon uptake and ET in semiarid ecosystems.

In this study, we use the Community Land Model Version 5 (CLM5) to investigate how irrigation modulates hydrologic and biogeochemical dynamics in the Upper Columbia‐Priest Rapids (UCPR) watershed, a typical semiarid watershed located in the northwestern United States dominated by cropland. To our knowledge, this constitutes the first application of CLM5 with landscape heterogeneity fully resolved over a watershed. The model is calibrated and evaluated against flux measurements from an AmeriFlux site and the Moderate Resolution Imaging Spectroradiometer (MODIS) products. Two numerical experiments (i.e., irrigated and rainfed) are performed at hyperresolution (~1 km) over the period of 2010–2018, accounting for realistic crop types and management practices. Our results show that irrigation fundamentally alters hydrologic and biogeochemical dynamics of the watershed. By adding 79.6 mm year−1 water in addition to the mean annual precipitation of 204.0 mm year−1, irrigation leads to increases in evapotranspiration and runoff, accompanied by shallower groundwater table depths. Increases in crop productivity in response to irrigation result in more carbon storage in the watershed, and drastically large seasonal fluctuations in soil organic carbon in response to changes in soil temperature and moisture. Irrigation also intensifies the rate of denitrification and mineralization during the growing season, enhancing the interactions between soil mineral nitrogen, the atmosphere, and freshwater systems. Our study demonstrates the potential of CLM5 as an effective tool for understanding hydrological and biogeochemical dynamics in highly managed semiarid watersheds.

Drylands are an important component of Earth's carbon, water, and energy budgets. Changes in hydroclimatic conditions can significantly alter land surface fluxes in these ecosystems. In regions with groundwater‐river water exchange, land surface fluxes are strongly modulated by groundwater table fluctuations in response to river stage variations; however, interactions between fluxes and groundwater remain poorly understood in drylands. Here, we use eddy covariance data collected in paired upland and riparian semiarid ecosystems (AmeriFlux sites US‐Hn1 and US‐Hn2) during 4 years with contrasting river flow conditions to examine interactions among river stage, water table dynamics, meteorological drivers, and ecosystem fluxes. Elevated groundwater availability resulting from groundwater‐surface water exchange at the riparian site enhanced the latent heat flux (LE, the energy form of evapotranspiration), which was sustained longer in the growing season. Interannual variation of dry season LE was small at the upland site, while at the riparian site contrasting river flow conditions resulted in large interannual variability in LE, with June LE of high‐flow years being 2.3 times that of low‐flow years. The 4‐year mean summer LE at the riparian site was 2.1 times larger than that at the upland site. Carbon uptake (NEP) and LE were well correlated at the upland site at the monthly and annual timescales, but at the riparian site enhanced LE due to shallow water table depth did not correspond to enhanced NEP. Our results demonstrate that enhanced groundwater availability modulated by river stage is a critical driver of land surface fluxes in semiarid riparian ecosystems.

No abstract

 The non-closure of the surface energy balance increases with the increasing atmospheric instability. The increased non-closure is not caused by the inadequate sampling of large-scale turbulent motions. The scales of turbulent eddies dominating fluxes change as the atmospheric instability increases.