Search USGSSearch

SEARCH · Search USGS

Results for “Journal of Hydrology”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

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

At least 1,081 records · Page 60Linked to original sources

Responses of biological soil crusts to rehabilitation strategies

Biological soil crusts (biocrusts) are common to dryland ecosystems and can influence a broad suite of soil ecological functions including stability and surface hydrology. Due to long recovery times following disturbance, there is a clear need for rehabilitation strategies to enhance the recovery of biocrust communities. Essential to biocrust recovery are exopolysaccharides (EPS): secretions comprised mainly of high molecular weight polymers that protect cyanobacteria from harsh environmental conditions. We examined whether biocrust rehabilitation strategies (inoculation in combination with surface shading and artificial soil stabilization) promote EPS production. To test if responses varied by soil texture, we measured biocrust recovery on two fine-textured soil types (clay and sandy clay loam) in a cool desert ecosystem. Shade coupled with inoculum addition resulted in the highest biocrust recovery, especially on clay soils. Independent of rehabilitation strategies, natural recovery of biocrusts occurred more rapidly on clay soils, reflected by greater increases in chlorophyll a (chl a). Chl a, a proxy for cyanobacterial biomass, was correlated to EPS amounts, suggesting that cyanobacteria are significant contributors to EPS production in biocrust development. Despite the role of EPS in biocrust establishment, EPS amounts had negligible effects on soil stability due inherent properties of fine soil texture.

Utah

Spatiotemporal variability of modeled watershed scale surface-depression storage and runoff for the conterminous United States

This study uses the explores the viability of a proxy model calibration strategy through assessment of the spatiotemporal variability of surface-depression storage and runoff generated with the U.S. Geological Survey’s National Hydrologic Model (NHM) infrastructure for hydrologic response units (HRUs; n=109,951) across the conterminous United States (CONUS). Simulated values for each HRU of daily surface-depression storage (treated as a decimal fraction of total possible volume) and monthly normalized runoff (0 to 1) values were calculated using Spearman’s rho at monthly and annual aggregations. Locations where values are correlated show where previously-developed proxy calibration strategies are likely to be effective. In addition, differences in the correlation for monthly and annual time scale aggregations show which time scale drives surface-depression storage processes in the NHM. Results show overall long-term (annual) correlation is more common than short-term (monthly) correlation over the CONUS; however, summary statistics for eighty-six ecoregions show five with higher ranges of monthly relative to annual Spearman’s rank coefficient values. This landscape-scale analysis shows simulations aggregated to an annual time scale are generally more dominant for the CONUS; however, simulations aggregated to monthly, short-term time scales are more dominant in focused areas where surface-depression storage processes are investigated.

Journal of the American Water Resources Associatio

Geochemistry and Soils of the Big Smoky Valley Fens, Nevada

Fens are groundwater-fed wetlands that can provide habitat for plants and animals. Due to anthropogenic activities and climate change, many fens around the world are at risk. This paper presents the results of a study of the hydrology and geochemistry of fens in Big Smoky Valley, central Nevada to support the Bureau of Land Management’s activities in the area. A water sample from the largest fen in the study area was analyzed for its water chemistry and compared to a nearby alluvial aquifer and hot spring. The high SiO 2 concentration of the fen sample implies that the fen water may originate from geothermal water. A soil core was taken to analyze radiocarbon age and soil type. A majority of the core was composed of silt and clay interlayered with water-filled voids. Changes in the character of the clay with depth suggest that there may have been changes in the depositional environment over time. Radiocarbon dating of Ruppia seeds showed longevity of the fen, with the minimum 14 C age of the core as 4,375±40 years. This paper provides reconnaissance-level information on the Big Smoky Valley fens, but further information would be needed to better understand the source of water to the fens or how the fen environment has changed over time with climate.

Nevada

Adapting standardized trout monitoring to a changing climate for the upper Yellowstone River, Montana, USA

Objective Long‐term standardized monitoring programs are fundamental to assessing how fish populations respond to anthropogenic stressors. Standardized monitoring programs may need to adopt new methods to adapt to rapid environmental changes that are associated with a changing climate. In the upper Yellowstone River, Montana, biologists have used a standardized, mark–recapture monitoring protocol to annually estimate the abundance of trout since 1978 to assess population status and trends. However, within the past two decades, climate change has caused changes in discharge timing that have prevented standardized monitoring from occurring annually. Methods We investigated the feasibility of using two analytical methods, N‐mixture models and mean capture probability, for estimating the abundance of three trout species in the upper Yellowstone River using the historical long‐term data set; these methods allow abundance to be estimated when a mark–recapture estimate cannot be obtained due to hydrologic conditions. Result When compared with abundance estimates from mark–recapture methods, N‐mixture models most often resulted in negatively biased abundance estimates, whereas mean capture probability analyses resulted in positively biased abundance estimates. Additionally, N‐mixture models produced negatively biased estimates when tested against true abundance values from simulated data sets. The bias in the N‐mixture model estimates was caused by poor model fit and variation in capture probability. The bias in the mean capture probability estimates was caused by heterogeneity in capture probability, likely caused by variable environmental conditions, which were not accounted for in the models. Conclusion N‐mixture models and mean capture probability are not viable alternatives for estimating abundance in the upper Yellowstone River. Thus, exploring additional adaptations to sampling methodologies and analytical approaches, including models that require individually marked fish, will be valuable for this system. Climate change will undoubtedly necessitate changes to standardized sampling methods throughout the world; thus, developing alternative sampling and analytical methods will be important for maintaining the utility of long‐term data sets.

Montana

An empirical assessment of which inland floods can be managed

Riverine flooding is a significant global issue. Although it is well documented that the influence of landscape structure on floods decreases as flood size increases, studies that define a threshold flood-return period, above which landscape features such as topography, land cover and impoundments can curtail floods, are lacking. Further, the relative influences of natural versus built features on floods is poorly understood. Assumptions about the types of floods that can be managed have considerable implications for the cost-effectiveness of decisions to invest in transforming land cover (e.g., reforestation) and in constructing structures (e.g., storm-water ponds) to control floods. This study defines parameters of floods for which changes in landscape structure can have an impact. We compare nine flood-return periods across 31 watersheds with widely varying topography and land cover in the southeastern United States, using long-term hydrologic records (≥20 years). We also assess the effects of built flow-regulating features (best management practices and artificial water bodies) on selected flood metrics across urban watersheds. We show that landscape features affect magnitude and duration of only those floods with return periods ≤10 years, which suggests that larger floods cannot be managed effectively by manipulating landscape structure. Overall, urban watersheds exhibited larger (270 m 3 /s) but quicker (0.41 days) floods than non-urban watersheds (50 m 3 /s and 1.5 days). However, urban watersheds with more flow-regulating features had lower flood magnitudes (154 m 3 /s), but similar flood durations (0.55 days), compared to urban watersheds with fewer flow-regulating features (360 m 3 /s and 0.23 days). Our analysis provides insight into the magnitude, duration and count of floods that can be curtailed by landscape structure and its management. Our findings are relevant to other areas with similar climate, topography, and land use, and can help ensure that investments in flood management are made wisely after considering the limitations of landscape features to regulate floods.

Journal of Environmental Management

The Metzger marsh restoration: A vegetation-centric look after 27 years

We investigated wetland vegetation before, during, and after dike construction at the Metzger Marsh project in western Lake Erie, which was designed to restore a 300-ha wetland that had been degraded following the loss of a protective barrier beach. A dike was constructed in 1995 to replace the function of the eroded barrier beach, but it contained a water-control structure to allow managed hydrologic connection to the lake. The control structure contained a fish passageway to allow movement of fish across the dike, while restricting entry of large common carp. Color-infrared aerial photos from project start in 1994 through 2010 (and 2022) were analyzed to track vegetation changes, and major vegetation types were sampled quantitatively. Drawdown of water levels in 1996 after dike construction elicited a response of mudflat species from the seed bank, as well as tree seedlings. Over half of the marsh was vegetated then and in subsequent years. The water-control structure was opened in 1998, and by 2000, invasive Phragmites australis had gained dominance. Most trees were eventually eliminated by herbicide treatment and flooding, and extent of Phragmites was reduced by management actions. Typha spp. and emergents Sagittaria latifolia and Schoenoplectus tabernaemontani became dominant by 2022. This restoration project increased habitat values for fish and wildlife; it also provided lessons for future projects on lands managed by multiple agencies with differing missions. More importantly, it showed that long-term monitoring data are critical for assessing wetland restoration projects and guiding management decisions.

Ohio

Nutrient dynamics in partially drained arctic thaw lakes

Thaw lakes are ubiquitous on arctic coastal plains (ACPs). While many thaw lakes have steep banks, stable water levels, and static surface areas, others only partially fill their basins and vary in area over the summer. These partially drained lakes (PDLs) are hydrologically connected to the wetlands immediately surrounding them. Heat and nutrient availability limit aquatic productivity on ACPs, and we hypothesized that shallow shorelines and greater hydrologic connectivity with the landscape should result in greater nutrient concentrations and biogeochemical cycling in PDLs. We tested this by monitoring water chemistry in lakes with varying levels of seasonal drainage in sandy and silty peaty lowland sites on the ACP of Alaska. One highly drained lake (N1) was significantly warmer than minimally drained lakes (minDLs) related to earlier ice off, reaching temperatures as high as 16 °C in June when minDLs still contained ice. Ammonia, total dissolved phosphorus, and dissolved organic carbon and nitrogen concentrations were higher in lakes with greater drainage, and concentrations in N1 rivaled those in the small, biologically productive ponds. Many PDLs displayed a midsummer decrease in nutrients consistent with assimilation by the aquatic ecosystem, and a late‐summer increase most likely related to runoff from drained lake margins following precipitation. N1 exported kilograms of ammonium and total dissolved phosphorus to the stream network over the summer. Given increased warming and drying in the arctic, the proportion of PDLs may be changing, which in turn may affect nutrient and organic matter availability in arctic lakes and export to downstream environments.

Journal of Geophysical Research: Biogeosciences

Improved fish counting method accurately quantifies high‐density fish movement in dual‐frequency identification sonar data files from a coastal wetland environment

There are many ways to quantify fish movement through shallow‐water habitats, but most noninvasive methods (e.g., visual counts) are not effective in turbid coastal wetland waters of the Great Lakes. Dual‐frequency identification sonar (DIDSON) technology (Sound Metrics) offers a noninvasive, hydroacoustic‐based approach to characterize fish movement in wetlands and other habitats by collecting highly detailed fish movement data regardless of light and water quality conditions. High‐resolution data can be analyzed to estimate fish movement in areas where visual observations are difficult. However, enumerating a complex mix of fish sizes by manually counting fish visible in echogram files requires training and is very time consuming. Therefore, four counting techniques were tested to estimate fish abundance from DIDSON echograms that were collected at a hydrologically reconnected coastal wetland in the Great Lakes. Briefly, the four counting methods were (1) manually viewing the entire length of the echogram (full‐hour manual count), (2) manually viewing subsections of the echogram before generating fish estimates by per‐minute average (subsample manual count), (3) using Echoview automated software to generate automated estimates, and (4) using DIDSON viewer software to generate automated estimates. Over 800 echogram‐hours were recorded over a 9‐month period at an open‐flow water control structure connecting a coastal wetland to a tributary to Lake Erie. Commercial fish tracking software (Echoview) and custom software scripts from Milne Technologies were used to semi‐automate fish count estimates for a small subset of data. Semi‐automated software counts were compared to manual counts of identical data files to assess differences in accuracy, cost, processing time, and counter effort. Semi‐automated fish count estimates using Echoview and custom pre‐ and postprocessing software scripts did not differ from baseline manual counts, suggesting that the semi‐automated count process could be a reliable tool to increase efficiency when processing large DIDSON data sets.

Great Lakes

Hydrologic landscapes on the Delmarva Peninsula - Part 2: Estimates of base-flow nitrogen load to Chesapeake Bay

ABSTRACT: The relation between landscape characteristics and water chemistry on the Delmarva Peninsula can be determined through a principal-component analysis of basin characteristics. Two basin types were defined by factor scores: (1) well-drained basins, characterized by combinations of a low percentage of forest cover, a low percentage of poorly drained soil, and elevated channel slope; and (2) poorly drained basins, characterized by a combinations of an elevated percentage of forest cover, an elevated percentage of poorly drained soil, and low channel slopes. Results from base-flow sampling of 29 basins during spring 1991 indicate that water chemistry of the two basin types differ significantly. Concentrations of calcium, magnesium, potassium, alkalinity, chloride, and nitrate are elevated in well-drained basins, and specific conductance is elevated. Concentrations of aluminum, dissolved organic carbon, sodium, and silica are elevated in poorly drained basins whereas specific conductance is low. The chemical patterns found in well-drained basins can be attributed to the application of agricultural chemicals, and those in poorly drained basins can be attributed to ground-water flowpaths. These results indicate that basin types determined by a quantitative analysis of basin characteristics can be related statistically to differences in base-flow chemistry, and that the observed statistical differences can be related to major processes that affect water chemistry.

Journal of the American Water Resources Associatio

Ostracode δ18O and δ13C evidence of Holocene environmental changes in the sediments of two Minnesota lakes

Stable oxygen and carbon isotope geochemistry of ostracode valves, abundance and assemblages of ostracode species, and sedimentological parameters from cores taken in Williams and Shingobee Lakes in north-central Minnesota show changes in climatic and hydrologic history during the Holocene. Isotopic records are consistent with the following scenario: Before 9800 yr B.P. the two lakes were connected. Increasing evaporation through the jack/red pine period (9800-7700 yr B.P.) led to lower lake levels, leaving small separated basins. The prairie period (7700-4000 yr B.P.) reflects high aridity, and lake levels reached low stands shortly before 6500 yr B.P. Low lake levels are associated with groundwater discharge between 6500 and 6000 yr B.P. The hardwood period (4000-3200 yr B.P.) corresponds to long cold winters and warm to cool summers with lower evaporation rates and slower sedimentation. During the white pine period (<3200 yr B.P.) evaporation increased and/or precipitation shifted to the summer months. These changes can be related to shifting atmospheric circulation patterns. Zonal flow was probably dominant during the early Holocene until the end of the prairie period ( c. 4000 yr B.P.). During the hardwood period a combination of zonal and meridional flow patterns caused long and cold winters and wetter summers. During the white pine period wintners were shorter and the meridional flow pattern more significant. Today meridional flow dominates the circulation pattern.

Minnesota

The consequences of neglecting reservoir storage in national-scale hydrologic models: An appraisal of key streamflow statistics

A better understanding of modeled streamflow errors related to basin reservoir storage is needed for large regions, which normally have many ungaged basins with reservoirs. We quantified the difference between modeled and observed streamflows for one process-based and three statistical-transfer hydrologic models, none of which explicitly accounted for reservoir storage. Streamflow statistics representing low to high flows, seasonality, annual variability, and daily autocorrelation were examined at 1082 study basins across the conterminous USA. All models increasingly overpredict (or decreasingly underpredict) observed annual maximum flows with increasing storage. Correlations between absolute values of errors for low-flow statistics and storage are often larger in magnitude than those for signed errors—additional storage is associated with increases in model errors in both directions even when its overall effect in one direction is weak. The rate of increase in absolute values of model errors was nonlinear for most statistics. For low flows, model errors had a change point to larger errors at 48 days of reservoir storage (relative to long-term mean daily flow); mean and high flows had change points at 147 to 176 days. We present predicted-to-observed errors for nine streamflow statistics over a large range of reservoir storage to help modelers and users of modeled streamflow understand the amount of storage for which explicit reservoir modeling is needed.

Journal of the American Water Resources Associatio

Effect of tidal fluctuations on transient dispersion of simulated contaminant concentrations in coastal aquifers

Variable-density groundwater models require extensive computational resources, particularly for simulations representing short-term hydrologic variability such as tidal fluctuations. Saltwater-intrusion models usually neglect tidal fluctuations and this may introduce errors in simulated concentrations. The effects of tides on simulated concentrations in a coastal aquifer were assessed. Three analyses are reported: in the first, simulations with and without tides were compared for three different dispersivity values. Tides do not significantly affect the transfer of a hypothetical contaminant into the ocean; however, the concentration difference between tidal and non-tidal simulations could be as much as 15%. In the second analysis, the dispersivity value for the model without tides was increased in a zone near the ocean boundary. By slightly increasing dispersivity in this zone, the maximum concentration difference between the simulations with and without tides was reduced to as low as 7%. In the last analysis, an apparent dispersivity value was calculated for each model cell using the simulated velocity variations from the model with tides. Use of apparent dispersivity values in models with a constant ocean boundary seems to provide a reasonable approach for approximating tidal effects in simulations where explicit representation of tidal fluctuations is not feasible.

Hydrogeology Journal

Interactions between soil thermal and hydrological dynamics in the response of Alaska ecosystems to fire disturbance

Soil temperature and moisture are important factors that control many ecosystem processes. However, interactions between soil thermal and hydrological processes are not adequately understood in cold regions, where the frozen soil, fire disturbance, and soil drainage play important roles in controlling interactions among these processes. These interactions were investigated with a new ecosystem model framework, the dynamic organic soil version of the Terrestrial Ecosystem Model, that incorporates an efficient and stable numerical scheme for simulating soil thermal and hydrological dynamics within soil profiles that contain a live moss horizon, fibrous and amorphous organic horizons, and mineral soil horizons. The performance of the model was evaluated for a tundra burn site that had both preburn and postburn measurements, two black spruce fire chronosequences (representing space-for-time substitutions in well and intermediately drained conditions), and a poorly drained black spruce site. Although space-for-time substitutions present challenges in model-data comparison, the model demonstrates substantial ability in simulating the dynamics of evapotranspiration, soil temperature, active layer depth, soil moisture, and water table depth in response to both climate variability and fire disturbance. Several differences between model simulations and field measurements identified key challenges for evaluating/improving model performance that include (1) proper representation of discrepancies between air temperature and ground surface temperature; (2) minimization of precipitation biases in the driving data sets; (3) improvement of the measurement accuracy of soil moisture in surface organic horizons; and (4) proper specification of organic horizon depth/properties, and soil thermal conductivity.

Alaska

Hydrologic resilience from summertime fog and recharge: A case study for coho salmon recovery planning

Fog and low cloud cover (FLCC) and late summer recharge increase stream baseflow and decrease stream temperature during arid Mediterranean climate summers, which benefits salmon especially under climate warming conditions. The potential to discharge cool water to streams during the late summer (hydrologic capacity; HC) furnished by FLCC and recharge were mapped for the 299 subwatersheds ranked Core, Phase 1, or Phase 2 under the National Marine Fisheries Service Recovery Plan that prioritized restoration and threat abatement action for endangered Central California Coast Coho Salmon evolutionarily significant unit. Two spatially continuous gridded datasets were merged to compare HC: average hrs/day FLCC, a new dataset derived from a decade of hourly National Weather Satellite data, and annual average mm recharge from the USGS Basin Characterization Model. Two use‐case scenarios provide examples of incorporating FLCC‐driven HC indices into long‐term recovery planning. The first, a thermal analysis under future climate, projected 65% of the watershed area for 8–19 coho population units as thermally inhospitable under two global climate models and identified several units with high resilience (high HC under the range of projected warming conditions). The second use case investigated HC by subwatershed rank and coho population, and identified three population units with high HC in areas ranked Phase 1 and 2 and low HC in Core. Recovery planning for cold‐water fish species would benefit by including FLCC in vulnerability analyses.

California

Groundwater recharge and agricultural contamination

Agriculture has had direct and indirect effects on the rates and compositions of groundwater recharge and aquifer biogeochemistry. Direct effects include dissolution and transport of excess quantities of fertilizers and associated materials and hydrologic alterations related to irrigation and drainage. Some indirect effects include changes in water–rock reactions in soils and aquifers caused by increased concentrations of dissolved oxidants, protons, and major ions. Agricultural activities have directly or indirectly affected the concentrations of a large number of inorganic chemicals in groundwater, for example NO 3 – , N 2 , Cl, SO 4 2– , H + , P, C, K, Mg, Ca, Sr, Ba, Ra, and As, as well as a wide variety of pesticides and other organic compounds. For reactive contaminants like NO 3 – , a combination of chemical, isotopic, and environmental-tracer analytical approaches might be required to resolve changing inputs from subsequent alterations as causes of concentration gradients in groundwater. Groundwater records derived from multi-component hydrostratigraphic data can be used to quantify recharge rates and residence times of water and dissolved contaminants, document past variations in recharging contaminant loads, and identify natural contaminant-remediation processes. These data indicate that many of the world's surficial aquifers contain transient records of changing agricultural contamination from the last half of the 20th century. The transient agricultural groundwater signal has important implications for long-term trends and spatial heterogeneity in discharge.

Hydrogeology Journal

Reconstruction of past methane availability in an Arctic Alaska wetland indicates climate influenced methane release during the past ~12,000 years

Atmospheric contributions of methane from Arctic wetlands during the Holocene are dynamic and linked to climate oscillations. However, long-term records linking climate variability to methane availability in Arctic wetlands are lacking. We present a multi-proxy ~12,000 year paleoecological reconstruction of intermittent methane availability from a radiocarbon-dated sediment core (LQ-West) taken from a shallow tundra lake (Qalluuraq Lake) in Arctic Alaska. Specifically, stable carbon isotopic values of photosynthetic biomarkers and methane are utilized to estimate the proportional contribution of methane-derived carbon to lake-sediment-preserved benthic (chironomids) and pelagic (cladocerans) components over the last ~12,000 years. These results were compared to temperature, hydrologic, and habitat reconstructions from the same site using chironomid assemblage data, oxygen isotopes of chironomid head capsules, and radiocarbon ages of plant macrofossils. Cladoceran ephippia from ~4,000 cal year BP sediments have δ13C values that range from ~−39 to −31‰, suggesting peak methane carbon assimilation at that time. These low δ13C values coincide with an apparent decrease in effective moisture and development of a wetland that included Sphagnum subsecundum. Incorporation of methane-derived carbon by chironomids and cladocerans decreased from ~2,500 to 1,500 cal year BP, coinciding with a temperature decrease. Live-collected chironomids with a radiocarbon age of 1,640 cal year BP, and fossil chironomids from 1,500 cal year BP in the core illustrate that ‘old’ carbon has also contributed to the development of the aquatic ecosystem since ~1,500 cal year BP. The relatively low δ13C values of aquatic invertebrates (as low as −40.5‰) provide evidence of methane incorporation by lake invertebrates, and suggest intermittent climate-linked methane release from the lake throughout the Holocene.

Alaska

An international code comparison study on coupled thermal, hydrologic and geomechanical processes of natural gas hydrate-bearing sediments

Geologic reservoirs containing gas hydrate occur beneath permafrost environments and within marine continental slope sediments, representing a potentially vast natural gas source. Numerical simulators provide scientists and engineers with tools for understanding how production efficiency depends on the numerous, interdependent (coupled) processes associated with potential production strategies for these gas hydrate reservoirs. Confidence in the modeling and forecasting abilities of these gas hydrate reservoir simulators (GHRSs) grows with successful comparisons against laboratory and field test results, but such results are rare, particularly in natural settings. The hydrate community recognized another approach to building confidence in the GHRS: comparing simulation results between independently developed and executed computer codes on structured problems specifically tailored to the interdependent processes relevant for gas hydrate-bearing systems. The United States Department of Energy, National Energy Technology Laboratory (DOE/NETL), sponsored the first international gas hydrate code comparison study, IGHCCS1, in the early 2000s. IGHCCS1 focused on coupled thermal and hydrologic processes associated with producing gas hydrates from geologic reservoirs via depressurization and thermal stimulation. Subsequently, GHRSs have advanced to model more complex production technologies and incorporate geomechanical processes into the existing framework of coupled thermal and hydrologic modeling. This paper contributes to the validation of these recent GHRS developments by providing results from a second GHRS code comparison study, IGHCCS2, also sponsored by DOE/NETL. IGHCCS2 includes participants from an international collection of universities, research institutes, industry, national laboratories, and national geologic surveys. Study participants developed a series of five benchmark problems principally involving gas hydrate processes with geomechanical components. The five problems range from simple geometries with analytical solutions to a representation of the world’s first offshore production test of methane hydrates, which was conducted with the depressurization method off the coast of Japan. To identify strengths and limitations in the various GHRSs, study participants submitted solutions for the benchmark problems and discussed differing results via teleconferences. The GHRSs evolved over the course of IGHCCS2 as researchers modified their simulators to reflect new insights, lessons learned, and suggested performance enhancements. The five benchmark problems, final sample solutions, and lessons learned that are presented here document the study outcomes and serve as a reference guide for developing and testing gas hydrate reservoir simulators.

Journal of Marine and Petroleum Geology

Responses of marginal and intrinsic water-use efficiency to changing aridity using FLUXNET observations

According to classic stomatal optimization theory, plant stomata are regulated to maximize carbon assimilation for a given water loss. A key component of stomatal optimization models is marginal water-use efficiency (mWUE), the ratio of the change of transpiration to the change in carbon assimilation. Although the mWUE is often assumed to be constant, variability of mWUE under changing hydrologic conditions has been reported. However, there has yet to be a consensus on the patterns of mWUE variabilities and their relations with atmospheric aridity. We investigate the dynamics of mWUE in response to vapor pressure deficit (VPD) and aridity index using carbon and water fluxes from 115 eddy covariance towers available from the global database FLUXNET. We demonstrate a non-linear mWUE-VPD relationship at a sub-daily scale in general; mWUE varies substantially at both low and high VPD levels. However, mWUE remains relatively constant within the mid-range of VPD. Despite the highly non-linear relationship between mWUE and VPD, the relationship can be informed by the strong linear relationship between ecosystem-level inherent water-use efficiency (IWUE) and mWUE using the slope, m *. We further identify site-specific m * and its variability with changing site-level aridity across six vegetation types. We suggest accurately representing the relationship between IWUE and VPD using Michaelis–Menten or quadratic functions to ensure precise estimation of mWUE variability for individual sites.

Journal of Geophysical Research Biogeosciences