Embracing fine-root system complexity in terrestrial ecosystem modeling.

Projecting the dynamics and functioning of the biosphere requires a holistic consideration of whole-ecosystem processes. However, biases toward leaf, canopy, and soil modeling since the 1970s have constantly left fine-root systems being rudimentarily treated. As accelerated empirical advances in the last two decades establish clearly functional differentiation conferred by the hierarchical structure of fine-root orders and associations with mycorrhizal fungi, a need emerges to embrace this complexity to bridge the data-model gap in still extremely uncertain models. Here, we propose a three-pool structure comprising transport and absorptive fine roots with mycorrhizal fungi (TAM) to model vertically resolved fine-root systems across organizational and spatial-temporal scales. Emerging from a conceptual shift away from arbitrary homogenization, TAM builds upon theoretical and empirical foundations as an effective and efficient approximation that balances realism and simplicity. A proof-of-concept demonstration of TAM in a big-leaf model both conservatively and radically shows robust impacts of differentiation within fine-root systems on simulating carbon cycling in temperate forests. Theoretical and quantitative support warrants exploiting its rich potentials across ecosystems and models to confront uncertainties and challenges for a predictive understanding of the biosphere. Echoing a broad trend of embracing ecological complexity in integrative ecosystem modeling, TAM may offer a consistent framework where modelers and empiricists can work together toward this grand goal.

Urban warming advances spring phenology but reduces the response of phenology to temperature in the conterminous United States.

Urbanization has caused environmental changes, such as urban heat islands (UHIs), that affect terrestrial ecosystems. However, how and to what extent urbanization affects plant phenology remains relatively unexplored. Here, we investigated the changes in the satellite-derived start of season (SOS) and the covariation between SOS and temperature (RT ) in 85 large cities across the conterminous United States for the period 2001-2014. We found that 1) the SOS came significantly earlier (6.1 ± 6.3 d) in 74 cities and RT was significantly weaker (0.03 ± 0.07) in 43 cities when compared with their surrounding rural areas (P < 0.05); 2) the decreased magnitude in RT mainly occurred in cities in relatively cold regions with an annual mean temperature <17.3 °C (e.g., Minnesota, Michigan, and Pennsylvania); and 3) the magnitude of urban-rural difference in both SOS and RT was primarily correlated with the intensity of UHI. Simulations of two phenology models further suggested that more and faster heat accumulation contributed to the earlier SOS, while a decrease in required chilling led to a decline in RT magnitude in urban areas. These findings provide observational evidence of a reduced covariation between temperature and SOS in major US cities, implying the response of spring phenology to warming conditions in nonurban environments may decline in the warming future.

Anthropogenically driven climate and landscape change effects on inland water carbon dynamics: What have we learned and where are we going?

Inland waters serve as important hydrological connections between the terrestrial landscape and oceans but are often overlooked in global carbon (C) budgets and Earth System Models. Terrestrially derived C entering inland waters from the watershed can be transported to oceans but over 83% is either buried in sediments or emitted to the atmosphere before reaching oceans. Anthropogenic pressures such as climate and landscape changes are altering the magnitude of these C fluxes in inland waters. Here, we synthesize the most recent estimates of C fluxes and the differential contributions across inland waterbody types (rivers, streams, lakes, reservoirs, and ponds), including recent measurements that incorporate improved sampling methods, small waterbodies, and dried areas. Across all inland waters, we report a global C emission estimate of 4.40 Pg C/year (95% confidence interval: 3.95-4.85 Pg C/year), representing a 13% increase from the most recent estimate. We also review the mechanisms by which the most globally widespread anthropogenically driven climate and landscape changes influence inland water C fluxes. The majority of these drivers are expected to influence terrestrial C inputs to inland waters due to alterations in terrestrial C quality and quantity, hydrological pathways, and biogeochemical processing. We recommend four research priorities for the future study of anthropogenic alterations to inland water C fluxes: (1) before-and-after measurements of C fluxes associated with climate change events and landscape changes, (2) better quantification of C input from land, (3) improved assessment of spatial coverage and contributions of small inland waterbodies to C fluxes, and (4) integration of dried and drawdown areas to global C flux estimates. Improved measurements of inland water C fluxes and quantification of uncertainty in these estimates will be vital to understanding both terrestrial C losses and the "moving target" of inland water C emissions in response to rapid and complex anthropogenic pressures.

Multi-hypothesis comparison of Farquhar and Collatz photosynthesis models reveals the unexpected influence of empirical assumptions at leaf and global scales.

Mechanistic photosynthesis models are at the heart of terrestrial biosphere models (TBMs) simulating the daily, monthly, annual and decadal rhythms of carbon assimilation (A). These models are founded on robust mathematical hypotheses that describe how A responds to changes in light and atmospheric CO2 concentration. Two predominant photosynthesis models are in common usage: Farquhar (FvCB) and Collatz (CBGB). However, a detailed quantitative comparison of these two models has never been undertaken. In this study, we unify the FvCB and CBGB models to a common parameter set and use novel multi-hypothesis methods (that account for both hypothesis and parameter variability) for process-level sensitivity analysis. These models represent three key biological processes: carboxylation, electron transport, triose phosphate use (TPU) and an additional model process: limiting-rate selection. Each of the four processes comprises 1-3 alternative hypotheses giving 12 possible individual models with a total of 14 parameters. To broaden inference, TBM simulations were run and novel, high-resolution photosynthesis measurements were made. We show that parameters associated with carboxylation are the most influential parameters but also reveal the surprising and marked dominance of the limiting-rate selection process (accounting for 57% of the variation in A vs. 22% for carboxylation). The limiting-rate selection assumption proposed by CBGB smooths the transition between limiting rates and always reduces A below the minimum of all potentially limiting rates, by up to 25%, effectively imposing a fourth limitation on A. Evaluation of the CBGB smoothing function in three TBMs demonstrated a reduction in global A by 4%-10%, equivalent to 50%-160% of current annual fossil fuel emissions. This analysis reveals a surprising and previously unquantified influence of a process that has been integral to many TBMs for decades, highlighting the value of multi-hypothesis methods.

Global vegetation biomass production efficiency constrained by models and observations.

Plants use only a fraction of their photosynthetically derived carbon for biomass production (BP). The biomass production efficiency (BPE), defined as the ratio of BP to photosynthesis, and its variation across and within vegetation types is poorly understood, which hinders our capacity to accurately estimate carbon turnover times and carbon sinks. Here, we present a new global estimation of BPE obtained by combining field measurements from 113 sites with 14 carbon cycle models. Our best estimate of global BPE is 0.41 ± 0.05, excluding cropland. The largest BPE is found in boreal forests (0.48 ± 0.06) and the lowest in tropical forests (0.40 ± 0.04). Carbon cycle models overestimate BPE, although models with carbon-nitrogen interactions tend to be more realistic. Using observation-based estimates of global photosynthesis, we quantify the global BP of non-cropland ecosystems of 41 ± 6 Pg C/year. This flux is less than net primary production as it does not contain carbon allocated to symbionts, used for exudates or volatile carbon compound emissions to the atmosphere. Our study reveals a positive bias of 24 ± 11% in the model-estimated BP (10 of 14 models). When correcting models for this bias while leaving modeled carbon turnover times unchanged, we found that the global ecosystem carbon storage change during the last century is decreased by 67% (or 58 Pg C).

Cryptic phenology in plants: Case studies, implications, and recommendations.

Plant phenology-the timing of cyclic or recurrent biological events in plants-offers insight into the ecology, evolution, and seasonality of plant-mediated ecosystem processes. Traditionally studied phenologies are readily apparent, such as flowering events, germination timing, and season-initiating budbreak. However, a broad range of phenologies that are fundamental to the ecology and evolution of plants, and to global biogeochemical cycles and climate change predictions, have been neglected because they are "cryptic"-that is, hidden from view (e.g., root production) or difficult to distinguish and interpret based on common measurements at typical scales of examination (e.g., leaf turnover in evergreen forests). We illustrate how capturing cryptic phenology can advance scientific understanding with two case studies: wood phenology in a deciduous forest of the northeastern USA and leaf phenology in tropical evergreen forests of Amazonia. Drawing on these case studies and other literature, we argue that conceptualizing and characterizing cryptic plant phenology is needed for understanding and accurate prediction at many scales from organisms to ecosystems. We recommend avenues of empirical and modeling research to accelerate discovery of cryptic phenological patterns, to understand their causes and consequences, and to represent these processes in terrestrial biosphere models.

Thermal, water, and land cover factors led to contrasting urban and rural vegetation resilience to extreme hot months.

With continuing global warming and urbanization, it is increasingly important to understand the resilience of urban vegetation to extreme high temperatures, but few studies have examined urban vegetation at large scale or both concurrent and delayed responses. In this study, we performed an urban-rural comparison using the Enhanced Vegetation Index and months that exceed the historical 90th percentile in mean temperature (referred to as "hot months") across 85 major cities in the contiguous United States. We found that hot months initially enhanced vegetation greenness but could cause a decline afterwards, especially for persistent (≥4 months) and intense (≥+2 °C) episodes in summer. The urban responses were more positive than rural in the western United States or in winter, but more negative during spring-autumn in the eastern United States. The east-west difference can be attributed to the higher optimal growth temperatures and lower water stress levels of the western urban vegetation than the rural. The urban responses also had smaller magnitudes than the rural responses, especially in deciduous forest biomes, and least in evergreen forest biomes. Within each biome, analysis at 1 km pixel level showed that impervious fraction and vegetation cover, local urban heat island intensity, and water stress were the key drivers of urban-rural differences. These findings advance our understanding of how prolonged exposure to warm extremes, particularly within urban environments, affects vegetation greenness and vitality. Urban planners and ecosystem managers should prioritize the long and intense events and the key drivers in fostering urban vegetation resilience to heat waves.

Quantification of human contribution to soil moisture-based terrestrial aridity.

Current knowledge of the spatiotemporal patterns of changes in soil moisture-based terrestrial aridity has considerable uncertainty. Using Standardized Soil Moisture Index (SSI) calculated from multi-source merged data sets, we find widespread drying in the global midlatitudes, and wetting in the northern subtropics and in spring between 45°N-65°N, during 1971-2016. Formal detection and attribution analysis shows that human forcings, especially greenhouse gases, contribute significantly to the changes in 0-10 cm SSI during August-November, and 0-100 cm during September-April. We further develop and apply an emergent constraint method on the future SSI's signal-to-noise (S/N) ratios and trends under the Shared Socioeconomic Pathway 5-8.5. The results show continued significant presence of human forcings and more rapid drying in 0-10 cm than 0-100 cm. Our findings highlight the predominant human contributions to spatiotemporally heterogenous terrestrial aridification, providing a basis for drought and flood risk management.

Simulated projections of boreal forest peatland ecosystem productivity are sensitive to observed seasonality in leaf physiology†.

We quantified seasonal CO2 assimilation capacities for seven dominant vascular species in a wet boreal forest peatland then applied data to a land surface model parametrized to the site (ELM-SPRUCE) to test if seasonality in photosynthetic parameters results in differences in simulated plant responses to elevated CO2 and temperature. We collected seasonal leaf-level gas exchange, nutrient content and stand allometric data from the field-layer community (i.e., Maianthemum trifolium (L.) Sloboda), understory shrubs (Rhododendron groenlandicum (Oeder) Kron and Judd, Chamaedaphne calyculata (L.) Moench., Kalmia polifolia Wangenh. and Vaccinium angustifolium Alton.) and overstory trees (Picea mariana (Mill.) B.S.P. and Larix laricina (Du Roi) K. Koch). We found significant interspecific seasonal differences in specific leaf area, nitrogen content (by area; Na) and photosynthetic parameters (i.e., maximum rates of Rubisco carboxylation (Vcmax25°C), electron transport (Jmax25°C) and dark respiration (Rd25°C)), but minimal correlation between foliar Na and Vcmax25°C, Jmax25°C or Rd25°C, which illustrates that nitrogen alone is not a good correlate for physiological processes such as Rubisco activity that can change seasonally in this system. ELM-SPRUCE was sensitive to the introduction of observed interspecific seasonality in Vcmax25°C, Jmax25°C and Rd25°C, leading to simulated enhancement of net primary production (NPP) using seasonally dynamic parameters as compared with use of static parameters. This pattern was particularly pronounced under simulations with higher temperature and elevated CO2, suggesting a key hypothesis to address with future empirical or observational studies as climate changes. Inclusion of species-specific seasonal photosynthetic parameters should improve estimates of boreal ecosystem-level NPP, especially if impacts of seasonal physiological ontogeny can be separated from seasonal thermal acclimation.

Uncertainty Quantification of Extratropical Forest Biomass in CMIP5 Models over the Northern Hemisphere.

Simplified representations of processes influencing forest biomass in Earth system models (ESMs) contribute to large uncertainty in projections. We evaluate forest biomass from eight ESMs outputs archived in the Coupled Model Intercomparison Project Phase 5 (CMIP5) using the biomass data synthesized from radar remote sensing and ground-based observations across northern extratropical latitudes. ESMs exhibit large biases in the forest distribution, forest fraction, and mass of carbon pools that contribute to uncertainty in forest total biomass (biases range from -20 Pg C to 135 Pg C). Forest total biomass is primarily positively correlated with precipitation variations, with surface temperature becoming equally important at higher latitudes, in both simulations and observations. Relatively small differences in forest biomass between the pre-industrial period and the contemporary period indicate uncertainties in forest biomass were introduced in the pre-industrial model equilibration (spin-up), suggesting parametric or structural model differences are a larger source of uncertainty than differences in transient responses. Our findings emphasize the importance of improved (1) models of carbon allocation to biomass compartments, (2) distribution of vegetation types in models, and (3) reproduction of pre-industrial vegetation conditions, in order to reduce the uncertainty in forest biomass simulated by ESMs.