Growing season length (GSL) is a key unifying concept in ecology that can be estimated from eddy covariance-derived estimates of net ecosystem production (NEP). Previous studies disagree on how increasing GSLs may affect NEP in evergreen coniferous forests, potentially due to the variety of methods used to quantify GSL from NEP. We calculated GSL and GSL-NEP regressions at eleven evergreen conifer sites across a broad climatic gradient in western North America using three common approaches: (1) variable length (3-7 days) regressions of day of year versus NEP, (2) a smoothed threshold approach, and (3) the carbon uptake period, followed by a new approach of a method-averaged ensemble. The GSL and the GSL-NEP relationship differed among methods, resulting in linear relationships with variable sign, slope, and statistical significance. For all combinations of sites and methods, the GSL explained between 6% and 82% of NEP with p-values ranging from 0.45 to < 0.01. These results demonstrate the variability among GSL methods and the importance of selecting an appropriate method to accurately project the ecosystem carbon cycling response to longer growing seasons in the future. To encourage this approach in future studies, we outline a series of best practices for GSL method selection depending on research goals and the annual NEP dynamics of the study site(s). These results contribute to understanding growing season dynamics at ecosystem and continental scales and underscore the potential for methodological variability to influence forecasts of the evergreen conifer forest response to climate variability.

Terrestrial production and export of dissolved organic and inorganic carbon (DOC and DIC) to streams depends on water flow and biogeochemical processes in and beneath soils. Yet, understanding of these processes in a rapidly changing climate is limited. Using the watershed-scale reactive-transport model BioRT-HBV and stream data from a snow-dominated catchment in the Rockies, we show deeper groundwater flow averaged about 20% of annual discharge, rising to ∼35% in drier years. DOC and DIC production and export peaked during snowmelt and wet years, driven more by hydrology than temperature. DOC was primarily produced in shallow soils (1.94 ± 1.45 gC/m2/year), stored via sorption, and flushed out during snowmelt. Some DOC was recharged to and further consumed in the deeper subsurface via respiration (−0.27 ± 0.02 gC/m2/year), therefore reducing concentrations in deeper groundwater and stream DOC concentrations at low discharge. Consequently, DOC was primarily exported from the shallow zone (1.62 ± 0.96 gC/m2/year, compared to 0.12 ± 0.02 gC/m2/year from the deeper zone). DIC was produced in both zones but at higher rates in shallow soils (1.34 ± 1.00 gC/m2/year) than in the deep subsurface (0.36 ± 0.02 gC/m2/year). Deep respiration elevated DIC concentrations in the deep zone and stream DIC concentrations at low discharge. In other words, deep respiration is responsible for the commonly-observed increasing DOC concentrations (flushing) and decreasing DIC concentrations (dilution) with increasing discharge. DIC export from the shallow zone was ~66% of annual export but can drop to ∼53% in drier years. Numerical experiments suggest lower carbon production and export in a warmer, drier future, and a higher proportion from deeper flow and respiration processes. These results underscore the often-overlooked but growing importance of deeper processes in a warming climate.

The 2015-2016 El Niño event ranks as one of the most severe on record in terms of the magnitude and extent of sea surface temperature (SST) anomalies generated in the tropical Pacific Ocean. Corresponding global impacts on the climate were expected to rival, or even surpass, those of the 1997-1998 severe El Niño event, which had SST anomalies that were similar in size. However, the 2015-2016 event failed to meet expectations for hydrologic change in many areas, including those expected to receive well above normal precipitation. To better understand how climate anomalies during an El Niño event impact soil moisture, we investigate changes in soil moisture in the humid tropics (between ±25ĝˆ ) during the three most recent super El Niño events of 1982-1983, 1997-1998 and 2015-2016, using data from the Global Land Data Assimilation System (GLDAS). First, we use in situ soil moisture observations obtained from 16 sites across five continents to validate and bias-correct estimates from GLDAS (r2Combining double low line0.54). Next, we apply a k-means cluster analysis to the soil moisture estimates during the El Niño mature phase, resulting in four groups of clustered data. The strongest and most consistent decreases in soil moisture occur in the Amazon basin and maritime southeastern Asia, while the most consistent increases occur over eastern Africa. In addition, we compare changes in soil moisture to both precipitation and evapotranspiration, which showed a lack of agreement in the direction of change between these variables and soil moisture most prominently in the southern Amazon basin, the Sahel and mainland southeastern Asia. Our results can be used to improve estimates of spatiotemporal differences in El Niño impacts on soil moisture in tropical hydrology and ecosystem models at multiple scales.

.