Biosphere Model Research Articles

Photosynthetic temperature responses in leaves and canopies: why temperature optima may disagree at different scales.

Understanding how canopy-scale photosynthesis responds to temperature is of paramount importance for realistic prediction of the likely impact of climate change on forest growth. The effects of temperature on leaf-scale photosynthesis have been extensively documented but data demonstrating the temperature response of canopy-scale photosynthesis are relatively rare, and the mechanisms that determine the response are not well quantified. Here, we compared leaf- and canopy-scale photosynthesis responses to temperature measured in a whole-tree chamber experiment and tested mechanisms that could explain the difference between leaf and crown scale temperature optima for photosynthesis. We hypothesised that 1) there is a large contribution of non-light saturated leaves to total crown photosynthesis; 2) photosynthetic component processes vary vertically through the canopy following the gradient in incident light; and 3) seasonal temperature acclimation of photosynthetic biochemistry has a significant role in determining the overall temperature response of canopy photosynthesis. We tested these hypotheses using three models of canopy radiation interception and photosynthesis parameterized with leaf-level physiological data and estimates of canopy leaf area. Our results identified the influence of non-light saturated leaves as a key determinant of the lower temperature optimum of canopy photosynthesis, which reduced the temperature optimum of canopy photosynthesis by 6-8°C compared to the leaf scale. Further, we demonstrate the importance of accounting for within-canopy variation and seasonal temperature acclimation of photosynthetic biochemistry in determining the magnitude of canopy photosynthesis. Overall, our study identifies key processes that need to be incorporated in terrestrial biosphere models to accurately predict temperature responses of whole-tree photosynthesis.

Partitioning of Evapotranspiration and Its Response to Eco‐Environmental Factors Over an Alpine Grassland on the Qinghai–Tibetan Plateau Based on the SiB2 Model

ABSTRACTPartitioning evapotranspiration (ET) is challenging but essential for understanding the exchange of energy, water, and carbon between terrestrial ecosystems and the atmosphere. In this study, we applied the simple biosphere model (SiB2) to partition ET at a typical alpine grassland site on the Qinghai–Tibet Plateau (QTP). In addition, through process‐based model scenario experiments, we quantified the effects of four environmental factors on ET components and predicted their evolution under the two future carbon emission scenarios (ssp126 and ssp585). Our findings are summarized as follows: (1) The original version of SiB2, despite its simple structure, effectively simulates ET and its components. (2) The ratios of annual total transpiration (T), soil evaporation (Es), and canopy interception evaporation (Ei) to ET in the alpine grassland ecosystem were 51%, 43%, and 6%, respectively. (3) Each 100 mm increase in annual precipitation results in a significant increase in soil evaporation (2.77%). A 1°C increase in air temperature leads to a significant increase in vegetation transpiration (5.22%) and canopy interception evaporation (5.63%). Each 100 ppm increase in CO2 concentration causes a significant decrease in T (−5.43%) and ET (−2.97%). An increase in LAI (1 m2 m−2) has the largest effect on canopy interception evaporation (4.67%). (4) Under the high carbon emission scenario (ssp585), all ET components in this ecosystem show a significant growth trend, particularly vegetation transpiration and canopy interception evaporation. These findings will facilitate more precise predictions of the water cycle dynamics, reveal land‐atmosphere interaction mechanisms, and aid in the protection of the ecological environment of the QTP.

Understanding Terrestrial Water and Carbon Cycles and Their Interactions Using Integrated SMAP Soil Moisture and OCO‐2 SIF Observations and Land Surface Models

AbstractRecently, more advanced synchronous global‐scale satellite observations, the Soil Moisture Active Passive enhanced Level 3 (SMAP L3) soil moisture product and the Orbiting Carbon Observatory 2 (OCO‐2) solar‐induced chlorophyll fluorescence (SIF) product, provide an opportunity to improve the predictive understanding of both water and carbon cycles in land surface modeling. The Simplified Simple Biosphere Model version 4 (SSiB4) was coupled with the Top‐down Representation of Interactive Foliage and Flora Including Dynamics Model (TRIFFID) and a mechanistic representation of SIF. Incorporating dynamic vegetation processes reduced global SIF root‐mean‐squared error (RMSE) by 12%. Offline experiments were conducted to understand the water and carbon cycles and their interactions using satellite data as constraints. Results indicate that soil hydraulic properties, the soil hydraulic conductivity at saturation (Ks) and the water retention curve, significantly impact soil moisture and SIF simulation, especially in the semi‐arid regions. The wilting point and maximum Rubisco carboxylation rate (Vmax) affect photosynthesis and transpiration, then soil moisture. However, without atmospheric feedback processes, their effects on soil moisture are undermined due to the compensation between soil evaporation and transpiration. With optimized parameters based on SMAP L3 and OCO‐2 data, the global RMSE of soil moisture and SIF simulations decreased by 15% and 12%, respectively. These findings highlight the importance of integrating advanced satellite data and dynamic vegetation processes to improve land surface models, enhancing understanding of terrestrial water and carbon cycles.

Different model assumptions about plant hydraulics and photosynthetic temperature acclimation yield diverging implications for tropical forest gross primary production under warming

AbstractTropical forest photosynthesis can decline at high temperatures due to (1) biochemical responses to increasing temperature and (2) stomatal responses to increasing vapor pressure deficit (VPD), which is associated with increasing temperature. It is challenging to disentangle the influence of these two mechanisms on photosynthesis in observations, because temperature and VPD are tightly correlated in tropical forests. Nonetheless, quantifying the relative strength of these two mechanisms is essential for understanding how tropical gross primary production (GPP) will respond to climate change, because increasing atmospheric CO2 concentration may partially offset VPD‐driven stomatal responses, but is not expected to mitigate the effects of temperature‐driven biochemical responses. We used two terrestrial biosphere models to quantify how physiological process assumptions (photosynthetic temperature acclimation and plant hydraulic stress) and functional traits (e.g., maximum xylem conductivity) influence the relative strength of modeled temperature versus VPD effects on light‐saturated GPP at an Amazonian forest site, a seasonally dry tropical forest site, and an experimental tropical forest mesocosm. By simulating idealized climate change scenarios, we quantified the divergence in GPP predictions under model configurations with stronger VPD effects compared with stronger direct temperature effects. Assumptions consistent with stronger direct temperature effects resulted in larger GPP declines under warming, while assumptions consistent with stronger VPD effects resulted in more resilient GPP under warming. Our findings underscore the importance of quantifying the role of direct temperature and indirect VPD effects for projecting the resilience of tropical forests in the future, and demonstrate that the relative strength of temperature versus VPD effects in models is highly sensitive to plant functional parameters and structural assumptions about photosynthetic temperature acclimation and plant hydraulics.

The food-water-climate nexus of green infrastructure: Examining ecosystem services trade-offs of peri-urban agriculture

Emission reduction, heat mitigation, and improved access to water and food provision are increasingly critical challenges for urban areas in the context of global climate change adaptation and mitigation. The revival of local agricultural production is often lauded as a potential nature-based solution. However, an expansion of peri-urban agriculture (peri-UA) may entail significant ecosystem trade-offs. This study explores the impacts on the food-water-climate nexus of different scenarios of peri-urban agricultural expansion in a semi-arid, Mediterranean climate, addressing local food provision, freshwater use, local temperature regulation, global climate change mitigation, and the trade-offs thereof. We estimate food provision and irrigation water requirements based on a georeferenced urban metabolism approach along with atmospheric and biosphere models to examine four land-use scenarios in the Metropolitan Area of Barcelona.Our study reveals that a 31 % (+17.27 km2) and 115 % (+64.25 km2) increase in the current peri-UA in the AMB, results in an increase in local food production of 24 % (+16,503 tons year−1) and 86 % (+58,940 tons year−1), and irrigation water requirements by 10.0 % (+3.2 hm3) and 43.5 % (+14.1 hm3), respectively. The expansion of irrigated peri-UA potentially reduces near-surface temperatures by 0.7 °C, albeit temperature reductions in the densest urban areas are minimal. Since the additional peri-UA is achieved by replacing natural non-forested and forest areas, the simulations predict reductions in the net ecosystem productivity of up to 18.5 % and total carbon stocks by 3.3 %. This integrated approach combining urban metabolism and atmospheric modelling to determine the trade-offs appears to be a promising tool for informing land-use decision-making in the context of urban climate adaptation and mitigation.

A novel AMSR2 retrieval algorithm for global C-band vegetation optical depth and soil moisture (AMSR2 IB): Parameters' calibration, evaluation and inter-comparison

Effective monitoring of soil and vegetation properties on a global scale is essential for better understanding climate changes, hydrological dynamics, and ecological processes. Passive microwave remote sensing at C-band radio frequency, with long observation period and relatively high penetration capability, has been widely used to retrieve soil moisture (SM) and vegetation optical depth (C-VOD). The retrieval process is generally achieved by inversion of the τ-ω radiative transfer model, which depends on crucial parameters such as effective scattering albedo (ω) and soil roughness (HR) for accurate retrievals. Current SM/C-VOD retrieval algorithms, such as the Land Parameter Retrieval Model (LPRM), predominantly rely on globally-constant ω and HR values, ignoring the inherent sensitivity of those parameters to varying soil conditions and vegetation types. To evaluate the impact of ω and HR variables on SM and C-VOD retrievals and to improve their accuracy, this study proposed and evaluated a novel retrieval approach from AMSR2 C-band observations during 2017–2020 using the C-band Microwave Emission of the Biosphere (C-MEB) model. We evaluated two new retrieval algorithms, considering either a globally-constant calibration or a land cover-based calibration of ω and HR. As a benchmark for the calibration, we optimized the values of ω and HR by evaluating the retrieved SM against in situ measurements from the International Soil Moisture Network (ISMN) and OzNet hydrological monitoring networks. The main originality compared to previous algorithms is that i) it includes a comprehensive calibration exploring the optimal values of ω and HR, applicable globally or tailored to specific land cover; ii) field SM measurements were leveraged to constrain the calibrated value of ω and HR.For the globally-constant calibration, the optimal values of ω = 0.05 and HR = 0.1 were found to yield the best results. For the land cover-based calibration, an inverse relationship between ω/HR and canopy height was observed, with ω ranging from 0.04 to 0.06 and HR ranging from 0.1 to 0.7 for heights between 0 and 30 m. The algorithm employing a land cover-based calibration (INRAE Bordeaux 2, IB2) exhibited better performance than the one utilizing a globally-constant calibration (INRAE Bordeaux 1, IB1) in evaluating retrieved SM against in situ measurements, as well as in evaluating C-VOD vs various vegetation variables including aboveground biomass (AGB), tree cover, canopy height and several optical vegetation indices. Comparison with LPRM suggested that our IB2 C-VOD retrievals present improved performances in terms of both spatial and temporal results with all considered vegetation variables (spatial correlation (R) between various vegetation variables and C-VOD of 0.76–0.83 for IB2 vs 0.69–0.79 for LPRM), and exhibited lower saturation effects when compared with AGB. In addition, the IB2 SM produced lower root mean squared error (RMSE) (0.147 vs 0.217 m3/m3), bias (−0.03 vs 0.09 m3/m3), and ubRMSE (0.066 vs 0.067 m3/m3) when compared with in situ measurements, although it showed a lower R compared to LPRM SM.

Extending our understanding on the retrievals of surface energy fluxes and surface soil moisture from the “triangle” technique

The present study demonstrates the capability of an inversion modelling scheme so-called the “triangle” to retrieve spatiotemporal estimates of surface energy fluxes and soil surface moisture (SSM) at high resolution using ASTER satellite imagery synergistically with SimSphere land biosphere model. In addition, as a further objective of this study is to examine the use of the technique for retrieving the Evaporative (EF) and the Non-Evaporative (NEF) Fractions as representations of the daytime average fluxes. The applicability of the investigated technique, is demonstrated for sixteen calendar days of year 2011 using in-situ data acquired from nine CarboEurope sites representing a variety of climatic, topographic and environmental conditions. Results indicated a close agreement between all the inverted parameters and the corresponding in-situ data. SSM predicted maps showed a small bias of 0.08 vol vol−1, a scatter of 0.18 vol vol−1 and a RMSD of 0.19 vol vol−1. The predicted LE fluxes showed a relatively low overall agreement (RMSD = 65.10 Wm-2), whereas for H flux reported RMSD was 85.02 Wm-2. The results also confirmed the ability of the investigated technique to provide meaningful estimates of the NEF and EF. All in all, the present study findings were at least comparable, or of higher accuracy, to those reported in other similar verification studies of the “triangle” using both high resolution (airborne) and low resolution (satellite) data. To our knowledge, this study represents the first comprehensive evaluation of the performance of this particular methodological implementation at a European setting combining the SimSphere SVAT model and ASTER EO datasets.

The impacts of modelling prescribed vs. dynamic land cover in a high-CO2 future scenario – greening of the Arctic and Amazonian dieback

Abstract. Terrestrial biosphere models are a key tool in investigating the role played by land surface in the global climate system. However, few models simulate the geographic distribution of biomes dynamically, opting instead to prescribe them using remote sensing products. While prescribing land cover still allows for the simulation of the impacts of climate change on vegetation growth and the impacts of land use change, it prevents the simulation of climate-change-driven biome shifts, with implications for the projection of future terrestrial carbon sink. Here, we isolate the impacts of prescribed vs. dynamic land cover implementations in a terrestrial biosphere model. We first introduce a new framework for evaluating dynamic land cover (i.e., the spatial distribution of plant functional types across the land surface), which can be applied across terrestrial biosphere models alongside standard benchmarking of energy, water, and carbon cycle variables in model intercomparison projects. After validating simulated land cover, we then show that the simulated terrestrial carbon sink differs significantly between simulations with dynamic vs. prescribed land cover for a high-CO2 future scenario. This is because of important range shifts that are only simulated when dynamic land cover is implemented: tree expansion into the Arctic and Amazonian transition from forest to grassland. In particular, the projected change in net land–atmosphere CO2 flux at the end of the 21st century is twice as large in simulations with dynamic land cover than in simulations with prescribed land cover. Our results illustrate the importance of climate-change-driven biome shifts for projecting future terrestrial carbon sink.

The Modeled Seasonal Cycles of Surface N2O Fluxes and Atmospheric N2O

AbstractNitrous oxide (N2O) is a greenhouse gas and stratospheric ozone‐depleting substance with large and growing anthropogenic emissions. Previous studies identified the influx of N2O‐depleted air from the stratosphere to partly cause the seasonality in tropospheric N2O (aN2O), but other contributions remain unclear. Here, we combine surface fluxes from eight land and four ocean models from phase 2 of the Nitrogen/N2O Model Intercomparison Project with tropospheric transport modeling to simulate aN2O at eight remote air sampling sites for modern and pre‐industrial periods. Models show general agreement on the seasonal phasing of zonal‐average N2O fluxes for most sites, but seasonal peak‐to‐peak amplitudes differ several‐fold across models. The modeled seasonal amplitude of surface aN2O ranges from 0.25 to 0.80 ppb (interquartile ranges 21%–52% of median) for land, 0.14–0.25 ppb (17%–68%) for ocean, and 0.28–0.77 ppb (23%–52%) for combined flux contributions. The observed seasonal amplitude ranges from 0.34 to 1.08 ppb for these sites. The stratospheric contributions to aN2O, inferred by the difference between the surface‐troposphere model and observations, show 16%–126% larger amplitudes and minima delayed by ∼1 month compared to Northern Hemisphere site observations. Land fluxes and their seasonal amplitude have increased since the pre‐industrial era and are projected to grow further under anthropogenic activities. Our results demonstrate the increasing importance of land fluxes for aN2O seasonality. Considering the large model spread, in situ aN2O observations and atmospheric transport‐chemistry models will provide opportunities for constraining terrestrial and oceanic biosphere models, critical for projecting carbon‐nitrogen cycles under ongoing global warming.

A First Exploration of the Ts/VI “Analytical Triangle” Technique with UAV Imagery for Deriving Key Surface Energy Balance Parameters at Very High Spatial Resolution

Knowledge on the spatiotemporal patterns of surface energy balance parameters is crucial for understanding climate system processes. To this end, the assimilation of Earth Observation data with land biosphere models has shown promising results, but they are still hampered by several limitations related to the spatiotemporal resolution of EO sensors and cloud contamination. With the recent developments on Unmanned Aerial Vehicles (UAVs), there is a great opportunity to overcome these challenges and gain knowledge of surface energy balance parameters at unprecedented resolutions. The present study examines, for the first time, the ability of an inversion-modeling scheme, the so-called “analytical triangle” method, to retrieve estimates of surface energy fluxes and soil surface moisture (SSM) at high spatial resolution using UAV data. A further aim of our study was to examine the representativeness of the SSM estimates for the SM measurements taken at different depths. The selected experimental site is an agricultural site of citrus trees located near the city of Palermo on 30 July 2019. The results of comparisons showed that the sensible and latent heat fluxes from UAV were consistent with those measured from the ground, with absolute differences in comparison to ground measurements being 5.00 Wm−2 for the latent heat (LE) flux and 65.02 Wm−2 for H flux, whereas for the daytime fluxes H/Rn and LE/Rn were 0.161 and 0.012, respectively. When comparing analytical triangle SSM estimates with SM measurements made at different depths, it was found that there was a gradual increase in underestimation with increasing measurement depth. All in all, this study’s results provide a credible demonstration of the significant potential of the technique investigated herein as a cost-effective and rapid solution for estimating key parameters characterizing land surface processes. As those parameters are required by a wide range of disciplines and applications, utilization of the investigated technique in research and practical applications is expected to be seen in the future.

Development and evaluation of the interactive Model for Air Pollution and Land Ecosystems (iMAPLE) version 1.0

Abstract. Land ecosystems are important sources and sinks of atmospheric components. In turn, air pollutants affect the exchange rates of carbon and water fluxes between ecosystems and the atmosphere. However, these biogeochemical processes are usually not well presented in Earth system models, limiting the explorations of interactions between land ecosystems and air pollutants from regional to global scales. Here, we develop and validate the interactive Model for Air Pollution and Land Ecosystems (iMAPLE) by upgrading the Yale Interactive Terrestrial Biosphere Model with process-based water cycles, fire emissions, wetland methane (CH4) emissions, and trait-based ozone (O3) damage. Within iMAPLE, soil moisture and temperature are dynamically calculated based on the water and energy balance in soil layers. Fire emissions are dependent on dryness, lightning, population, and fuel load. Wetland CH4 is produced but consumed through oxidation, ebullition, diffusion, and plant-mediated transport. The trait-based scheme unifies O3 sensitivity of different plant functional types (PFTs) with the leaf mass per area. Validations show correlation coefficients (R) of 0.59–0.86 for gross primary productivity (GPP) and 0.57–0.84 for evapotranspiration (ET) across the six PFTs at 201 flux tower sites and yield an average R of 0.68 for CH4 emissions at 44 sites. Simulated soil moisture and temperature match reanalysis data with high R above 0.86 and low normalized mean biases (NMBs) within 7 %, leading to reasonable simulations of global GPP (R=0.92, NMB=1.3 %) and ET (R=0.93, NMB=-10.4 %) against satellite-based observations for 2001–2013. The model predicts an annual global area burned of 507.1 Mha, close to the observations of 475.4 Mha with a spatial R of 0.66 for 1997–2016. The wetland CH4 emissions are estimated to be 153.45 Tg [CH4] yr−1 during 2000–2014, close to the multi-model mean of 148 Tg [CH4] yr−1. The model also shows reasonable responses of GPP and ET to the changes in diffuse radiation and yields mean O3 damage of 2.9 % to global GPP. iMAPLE provides an advanced tool for studying the interactions between land ecosystems and air pollutants.

Exploring the interplay between soil thermal and hydrological changes and their impact on carbon fluxes in permafrost ecosystems

Accelerated warming of the Arctic can affect the global climate system by thawing permafrost and exposing organic carbon in soils to decompose and release greenhouse gases into the atmosphere. We used a process-based biosphere model (DVM-DOS-TEM) designed to simulate biophysical and biogeochemical interactions between the soil, vegetation, and atmosphere. We varied soil and environmental parameters to assess the impact on cryohydrological and biogeochemical outputs in the model. We analyzed the responses of ecosystem carbon balances to permafrost thaw by running site-level simulations at two long-term tundra ecological monitoring sites in Alaska: Eight Mile Lake (EML) and Imnavait Creek Watershed (IMN), which are characterized by similar tussock tundra vegetation but differing soil drainage conditions and climate. Model outputs showed agreement with field observations at both sites for soil physical properties and ecosystem CO2 fluxes. Model simulations of Net Ecosystem Exchange (NEE) showed an overestimation during the frozen season (higher CO2 emissions) at EML with a mean NEE of 26.98 ± 4.83 gC/m2/month compared to observational mean of 22.01 ± 5.67 gC/m2/month, and during the fall months at IMN, with a modeled mean of 19.21 ± 7.49 gC/m2/month compared to observation mean of 11.9 ± 4.45 gC/m2/month. Our results underscore the importance of representing the impact of soil drainage conditions on the thawing of permafrost soils, particularly poorly drained soils, which will drive the magnitude of carbon released at sites across the high-latitude tundra. These findings can help improve predictions of net carbon releases from thawing permafrost, ultimately contributing to a better understanding of the impact of Arctic warming on the global climate system.

Carbon Cycle Models Quantify for a Green and Low-Carbon Economy

CO2 emissions from the combustion of fossil fuels represent the largest anthropogenic disruption of the natural global carbon cycle and lead to an undesirable increase in atmospheric CO2 levels. The so-called "carbon-neutral" renewable energy source biomass appears to be a promising replacement for fossil fuels. The "Combined Energy and Biosphere Model" (CEBM), was developed as a method for calculating the atmospheric CO2 content after intensive use of biomass. It includes (i) a biosphere part computing the annual cycle of carbon in the biosphere on 2,433 grid elements covering the continents’ plant mass and soil, as well as (ii) an energy economic part for calculating the fossil CO2 emissions of 119 countries. As the history of understanding the carbon cycle (including quantification of plant growth and decay by formulae, inclusion of deforestation, fossil fuel use and the fertilizer effect) shows, such a global carbon cycle model (as the CEBM) provides the needed methodology to assess the net effects of large-scale biomass energy use on the plant’s atmosphere. Results show that biomass (in a systemic view) is “only half as carbon neutral” as previously thought, even on a principal level. The conclusions indicate that (even if biomass is a valuable mitigation strategy) its global potential is limited. When interpreting CEBM scenarios, as a very first priority, reducing the current annual increase in energy demand is most highly needed for preserving climate.

Developing novel ensemble models for predicting soil hydraulic properties in China’s arid region

Accurately quantifying the mass (water, nutrients, and carbon) and energy exchange processes between the Earth’s atmosphere, biosphere, and lithosphere requires the accurate parameterization of soil hydraulic properties (SHPs) and their spatial heterogeneity. Because direct measurements of SHPs are difficult, time-consuming, and impossible at larger spatial scales, various pedotransfer functions (PTFs) have been developed in the last few decades, providing divergent estimates of SHPs from readily measurable variables. However, existing PTFs are mostly developed for specific regions and may not be suitable for other pedoclimatic conditions. Here, PTFs were developed using multiple machine-learning algorithms to estimate SHPs and examined the weaknesses and strengths of each method in estimating the average and variability of SHPs across China’s arid region. The optimal PTFs were applied to a 1 × 1 km2 regional map of texture and bulk density, thus producing maps of the saturated hydraulic conductivity (Ks), the parameters of the van Genuchten (VG) formulation, field capacity (θfc), wilting point (θwp), plant available water (θpa), and soil macroporosity (ϕm) in the 0–2 m soil profile throughout the region. The results indicate that the ensemble model with the averaging method (EMA) is the most robust for estimating all SHPs. The EMA-PTFs for Ks yielded the best performance for sand textures, followed by sandy loam, loam, silty loam, silty clay loam, loamy sand, and clay loam textures. There were no significant differences in estimating soil water retention curve parameters among the different soil texture classes. Using the same observed data set, we demonstrated that the new EMA-PTFs outperformed those of existing PTFs, such as Rosetta and HYPRES, with RMSE values decreasing by 25–83 % depending on the SHPs. Furthermore, our SHP datasets exhibited significantly deeper soil profiles and higher accuracy than other available regional and global SHP products. The VG retention parameter α shows the greatest variation vertically throughout the 0–2 m soil profile, followed by ln(Ks), θr, θpa, θfc, θwp, ϕm, θs and n. All SHPs exhibit much lower variability in the desert than in other areas, likely due to the homogeneous particle size distribution of desert areas. This study provides more accurate SHP estimates in China’s arid region than the existing SHP products by developing advanced PTFs and highlights the inconsistency in SHPs among different global or regional products; the uncertainties induced by PTFs should thus be considered in future terrestrial biosphere modeling.

The regional climate–chemistry–ecology coupling model RegCM-Chem (v4.6)–YIBs (v1.0): development and application

Abstract. The interactions between the terrestrial biosphere, atmospheric chemistry, and climate involve complex feedbacks that have traditionally been modeled separately. We present a new framework that couples the Yale Interactive terrestrial Biosphere (YIBs) model, a dynamic plant-chemistry model, with the RegCM-Chem model. RegCM-Chem–YIBs integrates meteorological variables and atmospheric chemical composition from RegCM-Chem with land surface parameters from YIBs. The terrestrial carbon flux calculated by YIBs is fed back into RegCM-Chem interactively, thereby representing the interactions between fine particulate matter (PM2.5), ozone (O3), and carbon dioxide (CO2). For testing purposes, we carry out a 1-year simulation (2016) at a 30 km horizontal resolution over East Asia with RegCM-Chem–YIBs. The model accurately captures the spatio-temporal distribution of climate, chemical composition, and ecological parameters. In particular, the estimated O3 and PM2.5 are consistent with ground observations, with correlation coefficients (R) of 0.74 and 0.65, respectively. The simulated CO2 concentration is consistent with observations from six sites (R ranged from 0.89 to 0.97) and exhibits a similar spatial pattern when compared with carbon assimilation products. RegCM-Chem–YIBs produces reasonably good gross primary productivity (GPP) and net primary productivity (NPP), showing seasonal and spatial distributions consistent with satellite observations, and mean biases (MBs) of 0.13 and 0.05 kg C m−2 yr−1. This study illustrates that RegCM-Chem–YIBs is a valuable tool to investigate coupled interactions between the terrestrial carbon cycle, atmospheric chemistry, and climate change at a higher resolution on a regional scale.

Optimizing seasonally variable photosynthetic parameters based on joint carbon and water flux constraints

Terrestrial biosphere models (TBMs) often adopt the Farquhar biochemical model coupled with the Ball-Berry stomatal conductance (gs) model to simulate ecosystem carbon and water fluxes. The parameters m, representing the sensitivity of gs to the photosynthetic rate, and Vcmax25, representing the leaf photosynthetic capacity, are two pivotal parameters but the two main sources of uncertainties in TBM simulations. The temporal variations of m in TBMs are still elusive, due to the lack of direct observations. It also remains unclear how accurate estimates of m and Vcmax25 can improve the simulations of carbon and water fluxes. In this study, we used a Bayesian parameter optimization approach to estimate seasonally varying m and Vcmax25 from eddy covariance observations in a mixed forest stand at the Borden Forest Research Station located in southern Ontario, Canada and used in-situ observations of m and Vcmax25 for validation. Three strategies were tested for optimizing m and Vcmax25, including the carbon, water, and carbon-water coupling scenarios. m and Vcmax25 optimized from carbon-water coupling constraints shows best correlations with the measured m (R2 = 0.70) and Vcmax25 (R2 = 0.70). By incorporating optimized m and Vcmax25with seasonal variations, we found considerable improvements in the estimated gross primary productivity (GPP) and evapotranspiration (ET) compared with constant m and Vcmax25, with R2 increasing from 0.78 to 0.85 for GPP, from 0.65 to 0.71 for ET and RMSE reducing from 2.579 g C m−2d−1 to 2.038 g C m−2d−1 for GPP, from 1.151 mm d−1 to 0.137 mm d−1 for ET. This study proposes an effective approach to retrieve m and Vcmax25 for TBMs and demonstrates the efficacy of incorporating seasonally variable m and Vcmax25 for reducing the uncertainties in GPP and ET simulations, which supports accurate quantifications of land-atmosphere exchanges.

Land-neutral negative emissions through biochar-based fertilization—assessing global potentials under varied management and pyrolysis conditions

Climate stabilization is crucial for restabilizing the Earth system but should not undermine biosphere integrity, a second pillar of Earth system functioning. This is of particular concern if it is to be achieved through biomass-based negative emission (NE) technologies that compete for land with food production and ecosystem protection. We assess the NE contribution of land- and calorie-neutral pyrogenic carbon capture and storage (LCN-PyCCS) facilitated by biochar-based fertilization, which sequesters carbon and reduces land demand by increasing crop yields. Applying the global biosphere model LPJmL with an enhanced representation of fast-growing species for PyCCS feedstock production, we calculated a land-neutral global NE potential of 0.20–1.10 GtCO2 year−1 assuming 74% of the biochar carbon remaining in the soil after 100 years (for + 10% yield increase; no potential for + 5%; 0.61–1.88 GtCO2 year−1 for + 15%). The potential is primarily driven by the achievable yield increase and the management intensity of the biomass producing systems. NE production is estimated to be enhanced by + 200–270% if management intensity increases from a marginal to a moderate level. Furthermore, our results show sensitivity to process-specific biochar yields and carbon contents, producing a difference of + 40–75% between conservative assumptions and an optimized setting. Despite these challenges for making world-wide assumptions on LCN-PyCCS systems in modeling, our findings point to discrepancies between the large NE volumes calculated in demand-driven and economically optimized mitigation scenarios and the potentials from analyses focusing on supply-driven approaches that meet environmental and socioeconomic preconditions as delivered by LCN-PyCCS.

Optimal representation of spring phenology on photosynthetic productivity across the northern hemisphere forests

Accurate simulation of the onset of spring is crucial for predicting the photosynthetic productivity of forest ecosystems. Nonetheless, the potential of phenology simulations on predictions of forest gross primary productivity (GPP) remains poorly understood, with previous studies generally focused on predicting phenology dates themselves or on limited scales. Here, we constructed a framework through the critical process of leaf growth to couple a terrestrial biosphere model—FORCCHN2 with five widely used phenology models, including one-phase (involving heat forcing), two-phase (considering chilling demand), and multi-phase models (integrating photoperiod effects) through the critical process of leaf growth. We evaluated the GPP performance from the multiple forest biomes in 74 eddy covariance (EC) sites and compared the total trends to satellite observations across the northern hemisphere forests. The predictions of the five models could reproduce 62.41–67.24 % variations of the GPP observations in all EC sites. In the deciduous broadleaf forest sites, we found the two-phase and multi-phase coupled models performed better than the one-phase coupled model. Chilling may play an essential role in controlling photosynthetic activity in deciduous broadleaf forests. Moreover, the multi-phase coupled model integrated photoperiod effects performed best at capturing GPP changes of the whole northern hemisphere forests. The predictions also showed the GPP in the northern hemisphere ranged from 21.68 to 26.06 Pg C year−1, and the total GPP showed an increasing trend. This study provides local and regional available evidence for the impact of phenology on photosynthesis with climate warming. The results highlight the necessity for enhancing understanding of the phenology factors on GPP and integrating the appropriate phenology representation to improve predictions of photosynthetic productivity for forests, especially in prediction to large scales and long-term trends.

Modelled forest ecosystem carbon–nitrogen dynamics with integrated mycorrhizal processes under elevated CO2

Abstract. Almost 95 % of all terrestrial plant species form symbioses with mycorrhizal fungi that mediate plant–soil interactions: mycorrhizae facilitate plant nitrogen (N) acquisition and are, therefore, vital for plant growth, but they also build a pathway for plant-assimilated carbon (C) into the rhizosphere. Therefore, mycorrhizae likely play an important role in shaping the response of ecosystems to environmental changes such as rising atmospheric carbon dioxide (CO2) concentrations, which can increase plant N demand and the transfer of plant C assimilation to the soil. While the importance of mycorrhizal fungi is widely recognised, they are rarely represented in current terrestrial biosphere models (TBMs) explicitly. Here, we present a novel, dynamic plant–mycorrhiza–soil model as part of the QUINCY (QUantifying Interactions between terrestrial Nutrient CYcles and the climate system) TBM. This new model is based on mycorrhizal functional types that either actively mine soil organic matter (SOM) for N or enhance soil microbial activity through increased transfer of labile C into the rhizosphere, thereby (passively) priming SOM decomposition. Using the Duke Free-Air CO2 Enrichment (FACE) experiment, we show that mycorrhizal fungi can have important effects on projected SOM turnover and plant nutrition under ambient as well as elevated-CO2 treatments. Specifically, we find that including enhanced active mining of SOM for N in the model allows one to more closely match the observations with respect to observed decadal responses of plant growth, plant N acquisition and soil C dynamics to elevated CO2, whereas a simple enhancement of SOM turnover by increased below-ground C transfer of mycorrhizae is unable to replicate the observed responses. We provide an extensive parameter uncertainty study to investigate the robustness of our findings with respect to model parameters that cannot readily be constrained by observations. Our study points to the importance of implementing mycorrhizal functionalities in TBMs as well as to further observational needs to better constrain mycorrhizal models and to close the existing major knowledge gaps in actual mycorrhizal functioning.

Global datasets of hourly carbon and water fluxes simulated using a satellite-based process model with dynamic parameterizations

Abstract. Diagnostic terrestrial biosphere models (TBMs) forced by remote sensing observations have been a principal tool for providing benchmarks on global gross primary productivity (GPP) and evapotranspiration (ET). However, these models often estimate GPP and ET at coarse daily or monthly steps, hindering analysis of ecosystem dynamics at the diurnal (hourly) scales, and prescribe some essential parameters (i.e., the Ball–Berry slope (m) and the maximum carboxylation rate at 25 °C (Vcmax25)) as constant, inducing uncertainties in the estimates of GPP and ET. In this study, we present hourly estimations of global GPP and ET datasets at a 0.25° resolution from 2001 to 2020 simulated with a widely used diagnostic TBM – the Biosphere–atmosphere Exchange Process Simulator (BEPS). We employed eddy covariance observations and machine learning approaches to derive and upscale the seasonally varied m and Vcmax25 for carbon and water fluxes. The estimated hourly GPP and ET are validated against flux observations, remote sensing, and machine learning-based estimates across multiple spatial and temporal scales. The correlation coefficients (R2) and slopes between hourly tower-measured and modeled fluxes are R2=0.83, regression slope =0.92 for GPP, and R2=0.72, regression slope =1.04 for ET. At the global scale, we estimated a global mean GPP of 137.78±3.22 Pg C yr−1 (mean ± 1 SD) with a positive trend of 0.53 Pg C yr−2 (p<0.001), and an ET of 89.03±0.82×103 km3 yr−1 with a slight positive trend of 0.10×103 km3 yr−2 (p<0.001) from 2001 to 2020. The spatial pattern of our estimates agrees well with other products, with R2=0.77–0.85 and R2=0.74–0.90 for GPP and ET, respectively. Overall, this new global hourly dataset serves as a “handshake” among process-based models, remote sensing, and the eddy covariance flux network, providing a reliable long-term estimate of global GPP and ET with diurnal patterns and facilitating studies related to ecosystem functional properties, global carbon, and water cycles. The hourly GPP and ET estimates are available at https://doi.org/10.57760/sciencedb.ecodb.00163 (Leng et al., 2023a) and the accumulated daily datasets are available at https://doi.org/10.57760/sciencedb.ecodb.00165 (Leng et al., 2023b).