Joel A. Biederman

Recent tree die‐off has little effect on streamflow in contrast to expected increases from historical studies

Recent bark beetle epidemics have caused regional-scale tree mortality in many snowmelt-dominated headwater catchments of western North America. Initial expectations of increased streamflow have not been supported by observations, and the basin-scale response of annual streamflow is largely unknown. Here we quantified annual streamflow responses during the decade following tree die-off in eight infested catchments in the Colorado River headwaters and one nearby control catchment. We employed three alternative empirical methods: (i) double-mass comparison between impacted and control catchments, (ii) runoff ratio comparison before and after die-off, and (iii) time-trend analysis using climate-driven linear models. In contrast to streamflow increases predicted by historical paired catchment studies and recent modeling, we did not detect streamflow changes in most basins following die-off, while one basin consistently showed decreased streamflow. The three analysis methods produced generally consistent results, with time-trend analysis showing precipitation was the strongest predictor of streamflow variability (R2 = 74–96%). Time-trend analysis revealed post-die-off streamflow decreased in three catchments by 11–29%, with no change in the other five catchments. Although counter to initial expectations, these results are consistent with increased transpiration by surviving vegetation and the growing body of literature documenting increased snow sublimation and evaporation from the subcanopy following die-off in water-limited, snow-dominated forests. The observations presented here challenge the widespread expectation that streamflow will increase following beetle-induced forest die-off and highlight the need to better understand the processes driving hydrologic response to forest disturbance.

Terrestrial carbon balance in a drier world: the effects of water availability in southwestern North America

Global modeling efforts indicate semiarid regions dominate the increasing trend and interannual variation of net CO2 exchange with the atmosphere, mainly driven by water availability. Many semiarid regions are expected to undergo climatic drying, but the impacts on net CO2 exchange are poorly understood due to limited semiarid flux observations. Here we evaluated 121 site-years of annual eddy covariance measurements of net and gross CO2 exchange (photosynthesis and respiration), precipitation, and evapotranspiration (ET) in 21 semiarid North American ecosystems with an observed range of 100 - 1000 mm in annual precipitation and records of 4-9 years each. In addition to evaluating spatial relationships among CO2 and water fluxes across sites, we separately quantified site-level temporal relationships, representing sensitivity to interannual variation. Across the climatic and ecological gradient, photosynthesis showed a saturating spatial relationship to precipitation, whereas the photosynthesis-ET relationship was linear, suggesting ET was a better proxy for water available to drive CO2 exchanges after hydrologic losses. Both photosynthesis and respiration showed similar site-level sensitivity to interannual changes in ET among the 21 ecosystems. Furthermore, these temporal relationships were not different from the spatial relationships of long-term mean CO2 exchanges with climatic ET. Consequently, a hypothetical 100-mm change in ET, whether short term or long term, was predicted to alter net ecosystem production (NEP) by 64 gCm(-2) yr(-1). Most of the unexplained NEP variability was related to persistent, site-specific function, suggesting prioritization of research on slow-changing controls. Common temporal and spatial sensitivity to water availability increases our confidence that site-level responses to interannual weather can be extrapolated for prediction of CO2 exchanges over decadal and longer timescales relevant to societal response to climate change.

The carbon balance pivot point of southwestern U.S. semiarid ecosystems: Insights from the 21st century drought

Global-scale studies indicate that semiarid regions strongly regulate the terrestrial carbon sink. However, we lack understanding of how climatic shifts, such as decadal drought, impact carbon sequestration across the wide range of structural diversity in semiarid ecosystems. Therefore, we used eddy covariance measurements to quantify how net ecosystem production of carbon dioxide (NEP) differed with relative grass and woody plant abundance over the last decade of drought in four Southwest U.S. ecosystems. We identified a precipitation “pivot point” in the carbon balance for each ecosystem where annual NEP switched from negative to positive. Ecosystems with grass had pivot points closer to the drought period precipitation than the predrought average, making them more likely to be carbon sinks (and a grass-free shrubland, a carbon source) during the current drought. One reason for this is that the grassland located closest to the shrubland supported higher leaf area and photosynthesis at the same water availability. Higher leaf area was associated with a greater proportion of evapotranspiration being transpiration (T/ET), and therefore with higher ecosystem water use efficiency (gross ecosystem photosynthesis/ET). Our findings strongly show that water availability is a primary driver of both gross and net semiarid productivity and illustrate that structural differences may contribute to the speed at which ecosystem carbon cycling adjusts to climatic shifts.

Increased evaporation following widespread tree mortality limits streamflow response

A North American epidemic of mountain pine beetle (MPB) has disturbed over 5 million ha of forest containing headwater catchments crucial to water resources. However, there are limited observations of MPB effects on partitioning of precipitation between vapor loss and streamflow, and to our knowledge these fluxes have not been observed simultaneously following disturbance. We combined eddy covariance vapor loss (V), catchment streamflow (Q), and stable isotope indicators of evaporation (E) to quantify hydrologic partitioning over 3 years in MPB-impacted and control sites. Annual control V was conservative, varying only from 573 to 623 mm, while MPB site V varied more widely from 570 to 700 mm. During wet periods, MPB site V was greater than control V in spite of similar above-canopy potential evapotranspiration (PET). During a wet year, annual MPB V was greater and annual Q was lower as compared to an average year, while in a dry year, essentially all water was partitioned to V. Ratios of 2H and 18O in stream and soil water showed no kinetic evaporation at the control site, while MPB isotope ratios fell below the local meteoric water line, indicating greater E and snowpack sublimation (Ss) counteracted reductions in transpiration (T) and sublimation of canopy-intercepted snow (Sc). Increased E was possibly driven by reduced canopy shading of shortwave radiation, which averaged 21 W m−2 during summer under control forest as compared to 66 W m−2 under MPB forest. These results show that abiotic vapor losses may limit widely expected streamflow increases.

Partitioning evapotranspiration using long‐term carbon dioxide and water vapor fluxes

The separate components of evapotranspiration (ET) elucidate the pathways and time scales over which water is returned to the atmosphere, but ecosystem-scale measurements of transpiration (T) and evaporation (E) remain elusive. We propose a novel determination of E and T using multiyear eddy covariance estimates of ET and gross ecosystem photosynthesis (GEP). The method is applicable at water-limited sites over time periods during which a linear regression between GEP (abscissa) and ET (ordinate) yields a positive ET axis intercept, an estimate of E. At four summer-rainfall semiarid sites, T/ET increases to a peak coincident with maximum GEP and remains elevated as the growing season progresses, consistent with previous, direct measurements. The seasonal course of T/ET is related to increasing leaf area index and declining frequency of rainy days—an index of the wet surface conditions that promote E—suggesting both surface and climatic controls on ET partitioning.

Riparian zones attenuate nitrogen loss following bark beetle-induced lodgepole pine mortality

A North American bark beetle infestation has killed billions of trees, increasing soil nitrogen and raising concern for N loss impacts on downstream ecosystems and water resources. There is surprisingly little evidence of stream N response in large basins, which may result from surviving vegetation uptake, gaseous loss, or dilution by streamflow from unimpacted stands. Observations are lacking along hydrologic flow paths connecting soils with streams, challenging our ability to determine where and how attenuation occurs. Here we quantified biogeochemical concentrations and fluxes at a lodgepole pine-dominated site where bark beetle infestation killed 50–60% of trees. We used nested observations along hydrologic flow paths connecting hillslope soils to streams of up to third order. We found soil water NO3 concentrations increased 100-fold compared to prior research at this and nearby southeast Wyoming sites. Nitrogen was lost below the major rooting zone to hillslope groundwater, where dissolved organic nitrogen (DON) increased by 3–10 times (mean 1.65 mg L−1) and NO3-N increased more than 100-fold (3.68 mg L−1) compared to preinfestation concentrations. Most of this N was removed as hillslope groundwater drained through riparian soils, and NO3 remained low in streams. DON entering the stream decreased 50% within 5 km downstream, to concentrations typical of unimpacted subalpine streams (~0.3 mg L−1). Although beetle outbreak caused hillslope N losses similar to other disturbances, up to 5.5 kg ha−1y−1, riparian and in-stream removal limited headwater catchment export to <1 kg ha−1y−1. These observations suggest riparian removal was the dominant mechanism preventing hillslope N loss from impacting streams.

Chlorophyll Fluorescence Better Captures Seasonal and Interannual Gross Primary Productivity Dynamics Across Dryland Ecosystems of Southwestern North America

Satellite remote sensing provides unmatched spatiotemporal information on vegetation gross primary productivity (GPP). Yet understanding of the relationship between GPP and remote sensing observations and how it changes with factors such as scale, biophysical constraint, and vegetation type remains limited. This knowledge gap is especially apparent for dryland ecosystems, which have characteristic high spatiotemporal variability and are under-represented by long-term field measurements. Here we utilize an eddy covariance (EC) data synthesis for southwestern North America in an assessment of how accurately satellite-derived vegetation proxies capture seasonal to interannual GPP dynamics across dryland gradients. We evaluate the enhanced vegetation index, solar-induced fluorescence (SIF), and the photochemical reflectivity index. We find evidence that SIF is more accurately capturing seasonal GPP dynamics particularly for evergreen-dominated EC sites andmore accurately estimating the full magnitude of interannual GPP dynamics for all dryland EC sites. These results suggest that incorporation of SIF could significantly improve satellite-based GPP estimates.