the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
![](https://bg.copernicus.org/articles/21/3735/2024/bg-21-3735-2024-avatar-thumb150.png)
Optimizing the terrestrial ecosystem gross primary productivity using carbonyl sulfide (COS) within a two-leaf modeling framework
Huajie Zhu
Xiuli Xing
Mousong Wu
Weimin Ju
Fei Jiang
Accurately modeling gross primary productivity (GPP) is of great importance for diagnosing terrestrial carbon–climate feedbacks. Process-based terrestrial ecosystem models are often subject to substantial uncertainties, primarily attributed to inadequately calibrated parameters. Recent research has identified carbonyl sulfide (COS) as a promising proxy of GPP due to the close linkage between leaf exchange of COS and carbon dioxide (CO2) through their shared pathway of stomatal diffusion. However, most of the current modeling approaches for COS and CO2 do not explicitly consider the vegetation structural impacts, i.e., the differences between the sunlit and shaded leaves in COS uptake. This study used ecosystem COS fluxes from seven sites to optimize GPP estimation across various ecosystems with the Biosphere-atmosphere Exchange Process Simulator (BEPS), which was further developed to simulate the canopy COS uptake under its state-of-the-art two-leaf framework. Our results demonstrated substantial improvement in GPP simulation across various ecosystems through the data assimilation of COS flux into the two-leaf model, with the ensemble mean of the root mean square error (RMSE) for simulated GPP reduced by 20.16 % to 64.12 %. Notably, we also shed light on the remarkable identifiability of key parameters within the BEPS model, including the maximum carboxylation rate of RuBisCO at 25 °C (Vcmax25), minimum stomatal conductance (), and leaf nitrogen content (Nleaf), despite intricate interactions among COS-related parameters. Furthermore, our global sensitivity analysis delineated both shared and disparate sensitivities of COS and GPP to model parameters and suggested the unique treatment of parameters for each site in COS and GPP modeling. In summary, our study deepened insights into the sensitivity, identifiability, and interactions of parameters related to COS and showcased the efficacy of COS in reducing uncertainty in GPP simulations.
- Article
(6578 KB) - Full-text XML
- BibTeX
- EndNote
Over the past 5 decades, terrestrial ecosystems have been absorbing about 30 % of anthropogenic carbon dioxide (CO2) emissions, playing a crucial role in mitigating climate change (Friedlingstein et al., 2022). Driven by the photosynthesis of terrestrial vegetation, gross primary productivity (GPP) is the largest terrestrial carbon flux and plays an important role in understanding terrestrial carbon–climate feedbacks (Luo, 2007; Wang et al., 2021). However, as direct observations of GPP using atmospheric CO2 observations are confounded by respiration (Hilton et al., 2017), and the modeling of GPP is affected by a range of uncertainties such as poorly calibrated parameters (Macbean et al., 2022), the precise quantification of GPP in terrestrial ecosystems has been a major challenge (Canadell et al., 2000; Yuan et al., 2007).
Carbonyl sulfide (COS) is the most abundant sulfur-containing trace gas in the atmosphere with a lifetime of about 2 years (Montzka et al., 2007; Karu et al., 2023). The tropospheric atmospheric mole fraction of COS is approximately 500 parts per trillion (ppt), exhibiting a typical seasonal amplitude of ∼ 100–200 ppt (Montzka et al., 2007; Ma et al., 2021; Hu et al., 2021; Remaud et al., 2022, 2023; Ma et al., 2023). In the past decade, COS has emerged as a promising tracer for terrestrial photosynthesis (Stimler et al., 2010; Asaf et al., 2013; Launois et al., 2015; Kooijmans et al., 2019) and stomatal conductance (Commane et al., 2015; Wehr et al., 2017; Sun et al., 2022) as the leaf exchange of COS and carbon dioxide (CO2) are tightly coupled through stomata (Sandoval-Soto et al., 2005; Seibt et al., 2010; Wohlfahrt et al., 2012; Whelan et al., 2018). Unlike CO2, which is emitted back into the atmosphere via leaf respiration (Sun et al., 2022), COS is completely destroyed by a hydrolysis reaction catalyzed by carbonic anhydrase (Protoschill-Krebs et al., 1996) without back flux in leaves under normal conditions (Stimler et al., 2010). Consequently, the measurement of COS flux is able to provide a direct and independent way to estimate GPP (Sandoval-Soto et al., 2005; Abadie et al., 2023).
In most early studies, GPP was directly estimated by scaling measurement of plant COS uptake with the empirically derived leaf relative uptake (LRU) approach or its extensions that incorporate the effects of temperature, humidity, light, and CO2 concentration on stomatal conductance (Kohonen et al., 2022a; Sun et al., 2022; Abadie et al., 2023) because of the simplicity of this approach and the sufficiency of it in many cases (Sandoval-Soto et al., 2005; Whelan et al., 2018). In contrast, the process-based model that mechanistically simulates COS plant uptake by incorporating stomatal transport processes has also been developed and widely evaluated (Maignan et al., 2021; Kooijmans et al., 2021). Concurrently, the significance of soil COS exchange has also been recognized, leading to the development of a suite of empirical or mechanistic COS soil exchange models (Kesselmeier et al., 1999; Berry et al., 2013; Launois et al., 2015; Sun et al., 2015; Ogée et al., 2016; Whelan et al., 2022). The process-based COS plant uptake model and soil exchange models have been integrated into land surface models (LSMs) (Berry et al., 2013; Maignan et al., 2021; Kooijmans et al., 2021). Consequently, by constraining the model parameters of LSMs with COS through data assimilation, not only are the model variables like GPP expected to be improved, but our understanding of ecosystem processes is also expected to be enhanced.
Currently, several studies have been conducted to refine the model parameters of LSMs through assimilating the COS fluxes and thereby optimizing the modeling of water–carbon fluxes (Chen et al., 2023; Abadie et al., 2023; Zhu et al., 2023). Within a big-leaf framework, Abadie et al. (2023) demonstrated COS could provide mechanistic constraint on stomatal diffusion, and the joint assimilation of COS and GPP is able to improve the model performance of GPP and latent heat. Ecosystem carbon, water, and energy processes interact and are nonlinear, and the changes in one process could induce variations in the other processes. While COS assimilation has proven to be effective in constraining COS-related model parameters and optimizing GPP, there remains a gap in systematic understanding of the ability of COS to optimize model parameters from different processes. For example, how effective is the assimilation of COS in reducing model prediction uncertainty in GPP and in the relevant ecosystem processes in different ecosystems?
Due to the dissimilar illumination conditions, there is a significant variability of leaf photosynthesis between sunlit and shaded leaves (Chen et al., 1999; Pignon et al., 2017; Wang et al., 2018; Bao et al., 2022). It is now clearly recognized that big-leaf models are conceptually flawed and practically inaccurate, and sunlit–shaded leaf stratification is necessary to make accurate canopy-level photosynthesis estimation (Chen et al., 1999; Chen et al., 2012; Luo et al., 2018). Consequently, in the process-based LSM that simulates COS plant uptake and photosynthesis in a coupled manner (Ball et al., 1987; Berry et al., 2013), the application of the two-leaf model shows promise for providing accurate simulation of plant COS uptake. In this context, we have further explored the capacity of COS to constrain the model parameters of an LSM and to optimize GPP within the two-leaf modeling framework.
Our goal is to address the following questions:
-
Which parameters is the COS simulation sensitive to, and what are the differences in parameter sensitivities between COS and GPP?
-
How effective is COS assimilation in improving model prediction and reducing prediction uncertainty in GPP?
-
Which process parameters can be well identified by the assimilation of COS?
-
How do process parameters interact in COS modeling across diverse ecosystems?
To address these questions, we utilized ecosystem COS flux data to optimize GPP across various ecosystems based on the coupling of COS modeling with the two-leaf-based Biosphere-atmosphere Exchange Process Simulator (BEPS). Through Monte Carlo simulations, we conducted a global parameter sensitivity analysis to explore the sensitivity of COS and GPP simulations to model parameters related not only to photosynthesis but also to water and energy. The interaction and identifiability of these parameters were quantified using Monte Carlo-optimized parameter sets. Additionally, the effectiveness of COS in constraining model uncertainty in simulated COS and GPP is evaluated.
2.1 Model description
2.1.1 BEPS basic model
The BEPS model (Liu et al., 1997; Chen et al., 1999; Chen et al., 2012) used in this study is a process-based diagnostic model driven by remotely sensed vegetation parameters, including leaf area index (LAI), clumping index, and land cover type, as well as meteorological and soil data (Chen et al., 2019). With the coupling among terrestrial carbon, water, and nitrogen cycles (He et al., 2021), it simulates photosynthesis, energy balance, and hydrological and soil biogeochemical processes at hourly time steps (Ju et al., 2006; Liu et al., 2015). For photosynthesis, it stratifies whole canopies into sunlit and shaded leaves and calculates GPP for each group of leaves by scaling Farquhar's leaf biochemical model (Farquhar et al., 1980) up to canopy-level with a temporal and spatial scaling scheme (Chen et al., 1999). In this study, the BEPS model stratifies the soil profile into five layers, and the model implicitly solves the soil water content values for these layers (Ju et al., 2010). Over the last few decades, the BEPS model has been continuously improved and has been used in a wide variety of terrestrial ecosystems (Schwalm et al., 2010; Liu et al., 2015). This study uses the BEPS model that simulates water, carbon, and energy processes at hourly intervals, which enables the detection of diel variations in model variables (Xing et al., 2023).
2.1.2 The two-leaf scheme for GPP and COS modeling
The BEPS model simulates the canopy photosynthesis for the sunlit and shaded leaves separately,
where GPPsunlit and GPPshaded denote the GPP per unit area for sunlit and shaded leaves, and LAIsunlit and LAIshaded represent the LAI values of sunlit and shaded leaves, respectively. LAIsunlit and LAIshaded are calculated as (Chen et al., 1999)
where θ is the solar zenith angle, LAItotal is the total leaf area index of the canopy, and Ω is the clumping index.
GPP values of sunlit and shaded leaves are calculated using Farquhar's model (Farquhar et al., 1980) with consideration of the large difference in incident solar irradiance between these two-leaf groups (Chen et al., 2012; Chen et al., 2019). Stomatal conductances of sunlit and shaded leaves are determined separately according to photosynthesis rates of these leaves, atmospheric CO2 concentration, relative humidity, and soil moisture (Ball et al., 1987; Ju et al., 2010). Detailed descriptions of the photosynthesis and stomatal conductance modeling approach of BEPS are illustrated in Appendix A1.
The ecosystem COS flux includes both plant COS uptake, FCOS,plant, and soil COS flux exchange, FCOS,soil (Whelan et al., 2016). In this work, the canopy-level COS plant uptake, FCOS,plant (), was calculated by upscaling the resistance analog model of COS uptake (Berry et al., 2013) with the two-leaf upscaling scheme (Chen et al., 1999). Considering the different responses of foliage to diffuse and direct solar radiation (Gu et al., 2002), FCOS,plant is calculated as
where FCOS,sunlit and FCOS,shaded denote the leaf-level COS uptake rate () for sunlit and shaded leaves. The leaf-level COS uptake rate, FCOS,leaf, is determined by the following formula (Berry et al., 2013):
where COSa represents the COS mole fraction in the bulk air. gsw and gbw are the stomatal conductance and leaf laminar boundary layer conductance to water vapor (H2O). The factors 1.94 and 1.56 account for the smaller diffusivity of COS with respect to H2O. gCOS indicates the apparent conductance for COS uptake from the intercellular airspaces, which combined the mesophyll conductance (Evans et al., 1994) and the biochemical reaction rate of COS and carbonic anhydrase (Badger and Price, 1994). It can be calculated as
where α is a parameter that is calibrated to observations of simultaneous measurements of COS and CO2 uptake (Stimler et al., 2012). Vcmax is the maximum carboxylation rate of RuBisCO. Analyses of these measurements yield estimates of α of ∼ 1400 for C3 and ∼ 7500 for C4 species. With reference to the COS modeling scheme of the Simple Biosphere Model (version 4.2) (Haynes et al., 2020), gCOS can be calculated as
where denotes the C4 plant flag, taking the value of 1 for C4 plants and 0 otherwise. fw is a soil moisture stress factor describing the sensitivity of gsw to soil water availability (Ju et al., 2006). FAPAR is the scaling factor for leaf radiation (Smith et al., 2008), calculated as
The soil COS fluxes are simulated by considering the abiotic and biotic components separately, as done by Whelan et al. (2016). We took the soil COS modeling scheme including the parameterizations from Whelan et al. (2016) and Whelan et al. (2022) in this study (see Appendix A2), given that our focus is the COS and GPP relationships, and previous studies have verified this approach with measurements over multiple sites.
Table 1Site characteristics. Site identification includes the country initials and a three-letter name for each site. Locations of the sites are provided by the latitude (lat) and longitude (long). PFT stands for plant functional type. ENF and DBF denote evergreen needleleaf forest and deciduous broadleaf forest, respectively.
![](https://bg.copernicus.org/articles/21/3735/2024/bg-21-3735-2024-t01-thumb.png)
2.2 Site description
The model was evaluated on seven sites distributed on the Eurasian and American continents in boreal, temperate, and subtropical regions based on field observations collected from several studies. Those sites were representative of different climate regions and land cover types (in the model represented by plant function types and soil textures, as depicted in Table 1).
2.3 Data
Data used in this study include LAI, land cover type, meteorological and soil data, and CO2 and COS mole fraction data. The CO2 and COS mole fractions in the bulk air were regarded as spatially invariant over the globe but assumed to vary annually. The CO2 mole fraction data in this study are taken from the Global Monitoring Laboratory (https://gml.noaa.gov/ccgg/trends/global.html, last access: July 2022). For the COS mole fraction, we utilized the average of observations from sites SPO (the South Pole) and MLO (Mauna Loa, United States) to drive the model. These data are publicly available online at https://gml.noaa.gov/hats/gases/OCS.html (last access: July 2022).
2.3.1 LAI dataset
The LAI datasets used here are the GLOBMAP global leaf area index product (version 3) (see GLOBMAP global leaf area index since 1981, https://doi.org/10.5281/zenodo.4700264; Liu et al., 2021) and the Global Land Surface Satellite (GLASS) LAI product (version 3) (acquired from https://doi.org/10.12041/geodata.GLASS_LAI_MODIS(0.05D).ver1.db, Xiao et al., 2016). They represent LAI at a spatial resolution of 8 km (Liu et al., 2012) and 1 km (Xiao et al., 2016), respectively, and a temporal resolution of 8 d. With reference to the observed LAI at these sites (Wehr et al., 2017; Rastogi et al., 2018a; Spielmann et al., 2019a; Kohonen et al., 2022a), we used GLOBMAP products to drive the BEPS model at most sites (five out of seven) due to its good agreement with the observed LAI. Specifically, as the GLOBMAP product had considerably underestimated LAI at DK-Sor and was not consistent with the vegetation phenology at ES-Lma during the measurement campaign (Spielmann et al., 2019a), GLASS LAI was used at these two sites. In addition, these LAI products were interpolated into daily values by the nearest-neighbor method for the simulation.
2.3.2 Meteorological dataset
Meteorological data required to force the BEPS model include air temperature, shortwave radiation, precipitation, relative humidity, and wind speed. As the simulations were conducted at the site scale, we utilized FLUXNET2015 data (see https://fluxnet.org, last access: June 2022, for AT-Neu, DK-Sor and ES-Lma, and FI-Hyy and US-Ha1) and AmeriFlux data (see https://ameriflux.lbl.gov, last access: June 2022, for US-Ha1 and US-Wrc). As FLUXNET2015 meteorological data for AT-Neu were only accessible for the period of 2002–2012, we conducted a linear fit between its European Centre for Medium-Range Weather Forecasts (ECMWF) Reanalysis v5 (ERA5) data (https://doi.org/10.24381/cds.adbb2d47, Hersbach et al., 2023) and FLUXNET2015 meteorological data for the corresponding period. Then, we used the fitted parameters to adjust the ERA5 data for 2015, thereby obtaining downscaling information for the meteorological data. In addition, we utilized the FLUXNET data in 2012 and AmeriFlux data and ERA5 shortwave radiation data in 2013 to drive the BEPS model at US-Ha1 due to the absence of FLUXNET data in 2013 and the lack of shortwave radiation data of AmeriFlux.
2.3.3 COS and GPP datasets
The hourly ecosystem COS flux observations were utilized to perform optimization and to evaluate the optimization results. They were derived from existing studies with pre-processing with regard to the data quality check, as listed in Table 1. Specifically, following the recommendations regarding the standardized processing of eddy covariance flux measurements of COS by Kohonen et al. (2020), both the measured and the gap-filled COS flux observations are provided in Vesala et al. (2022), and the latter is utilized in this study. To assess the model performance of GPP, the GPP observations were also collected from FLUXNET (DK-Sor, ES-Lma, FI-Hyy, and US-Ha1 in 2012), AmeriFlux (US-Ha1 in 2013), and existing studies for AT-Neu and IT-Soy (Spielmann et al., 2019a) and for US-Wrc (Rastogi et al., 2018a). Given that only CO2 turbulent flux (FC) or net ecosystem exchange (NEE) data were available at AT-Neu, IT-Soy, and US-Ha1 in 2013, a night flux partitioning model (Reichstein et al., 2005) was employed to derive GPP. This model assumes that nighttime NEE represents ecosystem respiration, Reco, and thus partitions FC or NEE into GPP and Reco based on the semi-empirical models of respiration, which use air temperature as a driver (Lloyd and Taylor, 1994; Lasslop et al., 2012).
2.4 The Monte Carlo-based parameter optimization approach
To evaluate the sensitivity, equifinality, and interaction of model parameters, the Monte Carlo-based parameter optimization approach was employed here (Fig. 1). The methodology calls for rejecting the concept of a unique global optimum parameter set within some particular model structure, instead recognizing the equifinality of parameter sets that exhibit similarly good performance in producing the observed data (Beven and Freer, 2001). In a Monte Carlo simulation framework, a large number of random sets of parameters are derived across specified parameter ranges (Staudt et al., 2010) and employed to drive the model. Subsequently, model realizations are grouped into behavioral and non-behavioral model runs and associated parameter sets based on the values of the single or multiple performance measures and the predefined threshold value (Houska et al., 2014). The former describes acceptable model realizations conditioned on the available observational data (Blasone et al., 2008; Beven and Binley, 2014). The latter describes parameter sets that produce behavior inconsistent with observed behavior. Given the gradual transition of performance measures from behavioral to non-behavioral model runs within the Monte Carlo framework, the threshold value used to distinguish between behavioral and non-behavioral parameter sets was often determined by an acceptable sample rate, i.e., ranking model runs and taking the top X% as behavioral (Beven and Binley, 2014). In the past few decades, this approach has been extensively used in ecosystem modeling with multiple parameters to be calibrated and has shown high ability in constraining multiple ecosystem processes (Tonkin and Doherty, 2009; Houska et al., 2014; He et al., 2016; Wu et al., 2019; Wu et al., 2020; Xing et al., 2023).
![https://bg.copernicus.org/articles/21/3735/2024/bg-21-3735-2024-f01](https://bg.copernicus.org/articles/21/3735/2024/bg-21-3735-2024-f01-thumb.png)
Figure 1Flowchart of the Monte Carlo-based parameter optimization approach used in this study. RMSE signifies root mean square error.
Table 2Descriptions of the nine parameters that were selected to be calibrated. The prior values and prior ranges (in parentheses) of these parameters are given for each plant function type (PFT) or for each soil texture or globally according to the parameter dependency. ENF and DBF denote evergreen needleleaf forest and deciduous broadleaf forest, respectively.
![](https://bg.copernicus.org/articles/21/3735/2024/bg-21-3735-2024-t02-thumb.png)
2.4.1 Parameter selection and sampling strategy
Based on current understanding of COS exchange (Wohlfahrt et al., 2012; Berry et al., 2013; Whelan et al., 2016; Whelan et al., 2018; Cho et al., 2023), photosynthesis (Ball et al., 1987; Raines, 2003; Blankenship, 2021), and related parameter sensitivity studies (Liu et al., 2011; Chen et al., 2012; Chen et al., 2023; Xing et al., 2023; Abadie et al., 2023; Zhu et al., 2023), nine parameters were selected to be calibrated in this study (for details see Table 2). These parameters are related to formulas describing four processes: (1) photosynthesis (Vcmax25, VJslope, Nleaf), (2) soil hydrology (Ksatscalar, bscalar, rdecay), (3) stomatal gas exchange (, ), and (4) energy balance (fleaf). Specifically, Ksatscalar and bscalar are scaling factors designed to optimize the saturated hydraulic conductivity (Ksat) and the Campbell parameter (b) for each soil layer in the BEPS model. The prior values and prior ranges for these parameters (Table 2) were chosen based on the literature (Jackson et al., 1996; Medlyn et al., 1999; Kattge et al., 2009; Miner et al., 2017; Ryu et al., 2018) and default model settings. Uniform distributions were assigned to all parameters, and 20 000 sets of parameters were generated through random sampling.
2.4.2 Selection of behavioral simulations
To measure the agreement between model simulations and observations, a variety of performance metrics have been proposed and utilized in previous studies (Beven and Binley, 1992; Moradkhani et al., 2005; Staudt et al., 2010). In this study, we employed the root mean square error (RMSE) to distinguish between behavioral and non-behavioral simulations:
where N is the total number of observations; obs and sim denote the observations and simulations, respectively; and simi denotes the simulation corresponding to the ith observation, obsi.
Specifically, here we chose an acceptable sample rate of 0.5 %, i.e., the top 100 model runs with the lowest RMSE values for COS as behavioral simulations. Thus, the deterministic model prediction is given by the ensemble mean of the 100 behavioral simulations.
2.5 Uncertainty quantification
The model prediction limits or uncertainty bounds can be determined by forming the cumulative density function (CDF) of the ensemble of simulations (Beven and Binley, 2014), normally chosen at the 5 % and 95 % confidence level in most of the previous studies (Blasone et al., 2008). Similarly, we chose the 5 % and 95 % quantiles of the 20 000 simulations and the 100 behavioral simulations to quantify the model output uncertainty in this study.
2.6 Parameter sensitivity
In order to take full advantage of the Monte Carlo simulations, a density-based global sensitivity analysis approach (Plischke et al., 2013) was used to investigate the sensitivity of COS and GPP simulations to the selected model parameters via the Sensitivity Analysis Library (SALib) (Iwanaga et al., 2022). This approach aims at assessing the influence of the entire input distribution on the entire output distribution without reference to a particular moment of the output (Borgonovo, 2007). According to Borgonovo (2007), the sensitivity index (δ) is always between 0 and 1, it equals 0 if the output is not dependent on the model parameter, and it equals unity if all model parameters are considered.
2.7 Parameter uncertainty
Due to the functional and structural complexity of ecosystems, ecosystem models often require a substantial number of parameters to realize the modeling of various ecosystem processes, and some parameters compensate for each other (Mo et al., 2008). While the parameter interactions related to photosynthesis have been systematically studied (Tang and Zhuang, 2009; Lu et al., 2013; Wu et al., 2019; Xing et al., 2023), the parameter interactions related to COS flux simulation have not been reported. Based on the Monte Carlo-based methodology, the numerous behavioral parameter sets around the “optimum” (Beven and Freer, 2001) provide us with the opportunity to analyze the interactions between the selected parameters. In this study, the Pearson correlation coefficient and the confidence level were employed to identify the parameter interactions.
Parameter identifiability (PI) is the concept of whether uncertain parameters can be correctly estimated from the observed data (Yi et al., 2019). The failure in PI is supposed to be caused by “over-parameterization” and parameter interactions (due to high nonlinearity of model equations) (Gan et al., 2014). Inspired by Yi et al. (2019), who used likelihood confidence interval as a measure of PI, here we used parameter distribution range for the same purpose. Taking into account the influence of the prior distribution of the behavioral parameter sets, the PI is defined as the reduction in the parameter range width. Hence, a large value of PI indicates that the parameter is well identified in the optimization process.
3.1 Parameter sensitivity
The sensitivity indexes of COS and GPP simulations with respect to the model parameters for the seven sites are illustrated in Fig. 2. It can be seen that both COS and GPP simulations exhibit high sensitivity to leaf nitrogen content (Nleaf) and the maximum carboxylation rate of RuBisCO at 25 °C (Vcmax25), while showing low sensitivity to the energy-balance-related parameter fleaf and soil-hydrology-related parameters (including bscalar, Ksatscalar, and rdecay). With the average values of the sensitivity index of 0.11 and 0.10, the photosynthesis-related parameter VJslope and stomata-conductance-related parameter can significantly impact the simulation of GPP. However, those parameters do not exhibit high sensitivity in the modeling of COS. Our results also highlight the crucial role of the intercept of the Ball–Berry model () in the modeling of COS, yet its impact on the simulation of GPP is limited. In summary, our results suggest that the simulated COS and GPP share some similarities in their sensitivities to parameters, but there are also notable differences. Specifically, the parameters and VJslope strongly influence GPP simulations but have minimal impacts on COS simulations. Conversely, the parameter plays a more crucial role in COS simulation.
With mean values of 0.33, 0.29, and 0.09, the sensitivity of COS simulations to Nleaf, Vcmax25, and , respectively, is much larger than that of GPP simulations. However, the patterns of the sensitivity of these parameters for COS and GPP simulations are very similar across these sites. Our results reveal that the simulated COS and GPP are more sensitive to Nleaf and less influenced by Vcmax25. In comparison to other sites, the role of Vcmax25 in the simulation of COS and GPP at IT-Soy is less significant. Additionally, we observed that holds greater importance in the modeling of GPP at US-Wrc than at other sites. Moreover, our results suggest that the modeling of GPP at deciduous broadleaf forest sites (DK-Sor and US-Ha1) is more sensitive to VJslope and less sensitive to compared with other sites.
![https://bg.copernicus.org/articles/21/3735/2024/bg-21-3735-2024-f03](https://bg.copernicus.org/articles/21/3735/2024/bg-21-3735-2024-f03-thumb.png)
Figure 3Cumulative frequency distributions and boxplots for the posterior model parameters obtained by COS assimilation. The grey area represents uniform parameter distributions, while the colored areas denote posterior CDF distributions, with parameters for different sites represented using different colors. The box extends from the first quartile to the third quartile of the parameter values, with a line at the median. The “×” markers denote outliers, and the whiskers represent the lowest or highest parameter values excluding any outliers. The black squares represent the prior parameter values, and the axis ranges denote the prior ranges of the parameters. PI denotes parameter identifiability, defined as the reduction in the parameter range width.
3.2 Posterior parameter distributions
The cumulative frequency distribution and the boxplots of each of the parameters for the 0.5 % best runs are plotted in Fig. 3, with a comparison to the uniform parameter distributions and the prior parameter values. As shown in Fig. 3, the posterior distributions of these parameters differ significantly, indicating that the response of these parameters to the assimilation of COS is quite different. Our results demonstrated that COS fluxes have similar constraining effects on the same parameters in different ecosystems, although the posterior distributions of the same parameter at different sites depicted variations. In general, parameters related to plant growth and stomatal conductance were strongly constrained by the assimilation of COS, while the parameters related to energy balance and soil hydrology were inadequately constrained.
Table 3Mean posterior values for seven study sites in comparison with the default values and the PFT-grouping values (mean ± standard deviation) in Miner et al. (2017). Within the compilation of Miner et al. (2017), FI-Hyy and US-Wrc are classified under the PFT of evergreen gymnosperm tree, while DK-Sor and US-Ha1 fall under the PFT of deciduous angiosperm tree.
![](https://bg.copernicus.org/articles/21/3735/2024/bg-21-3735-2024-t03-thumb.png)
With a distinct shape and remarkably narrow range of the cumulative frequency curves, (the intercept of the Ball–Berry model, representing minimum stomatal conductance) was strongly constrained by the assimilation of COS in this study. For most sites (AT-Neu, DK-Sor, FI-Hyy, IT-Soy, and US-Ha1), the values of were confined to a very limited range of 0 to 0.08 . At these five sites, the average values of the posterior were all located from 0.01 to 0.03 , aligning well with the default value of (0.0175 ) for the BEPS model. In contrast, with posterior values ranging from 0.03 to 0.18 and 0.03 to 0.91 , the default value of for the BEPS model was rejected by the assimilation of COS at ES-Lma and US-Wrc. Despite the broad distribution of posterior at US-Wrc, the cumulative frequency curve still indicates that is well constrained at this site, with 80 % of the posterior located in a narrow range of 0.15 to 0.50 . Overall, our results are reasonable as literature-documented values of are highly variable, and they align well with the compilation provided by Miner et al. (2017), in which more than five-sixths of the values are located between 0 and 0.18 , and about half are located between 0 and 0.04 . Moreover, the mean values of posterior for most (five out of seven) sites are larger than the default value of the BEPS model, suggesting that the current value utilized in BEPS may be underestimated.
Identified as the most sensitive parameters in COS modeling, the plant-growth-related parameters Vcmax25 and Nleaf were generally well constrained in this study. However, unlike , which is strongly constrained at all sites, the posterior cumulative frequency curves of Vcmax25 and Nleaf exhibit considerable variation across sites. Except for US-Ha1 and US-Wrc, the posterior Vcmax25 and Nleaf were mostly distributed in the upper half of the parameter range. Particularly, all of the lower-half values of Vcmax25 and Nleaf were excluded by the behavioral parameter sets at ES-Lma. In contrast, the posterior cumulative frequency curves of Vcmax25 deviated slightly from the original uniform distribution at US-Ha1, indicating that they are not well constrained by the assimilation of COS. As for US-Wrc, both the largest 7 % and the smallest 4 % values of Nleaf are effectively excluded by the assimilation of COS.
Another stomatal-conductance-related parameter, , demonstrated effective constraint through COS assimilation at specific sites (AT-Neu, DK-Sor, ES-Lma, FI-Hyy, and IT-Soy), with parameter range width reductions comparable to Vcmax25 and Nleaf. However, at US-Ha1 and US-Wrc, the posterior cumulative frequency curves of show minimal deviation from the original uniform distribution. Nevertheless, the optimization of is generally achievable through COS assimilation, as supported by our results, in good agreement with the compilation of Miner et al. (2017), in which the average historical values of grouped by PFT (referred to as the PFT-grouping values below) are provided. As indicated in Table 3, the average absolute bias between the default and the PFT-grouping value reached as high as 2.87 for these sites. Through COS assimilation, the mean absolute bias was reduced to 2.59.
The photosynthesis-related parameters VJslope and fleaf also influence COS simulation. However, the posterior distributions of fleaf resemble the original uniform distribution, suggesting that it is not a crucial parameter for COS simulations. The posterior cumulative frequency curve of VJslope also generally deviates slightly from the uniform distribution. Yet, at DK-Sor and US-Ha1, more than two-thirds of the posterior VJslope values are situated in the upper half of the parameter range, indicating that VJslope can also be well constrained by the assimilation of COS in specific cases.
Among these seven sites, the soil-hydrology-related parameters Ksatscalar and bscalar did not exhibit a strong response during the assimilation of COS. However, the posterior cumulative frequency curves of rdecay show notable deviations from the uniform distribution in certain cases. At US-Wrc, higher values of rdecay are more prevalent within the behavioral parameter sets, leading to the posterior mean of rdecay being much greater than the prior mean. Moreover, the largest 14 % values of rdecay are effectively excluded by the assimilation of COS at IT-Soy.
Table 4Comparison of model performance indices for the prior and posterior COS simulations. The RMSEmean values of the prior and posterior simulations are the mean values of the RMSEs of 20 000 prior COS simulations and 100 behavioral COS simulations with COS observations, respectively. The range widths of the prior and posterior COS simulations are defined as the mean values of the difference between the 95th and 5th percentile of the prior and posterior simulations, respectively. The reduction (%) in RMSEmean and range width is calculated as (1−posterior/prior) × 100.
![](https://bg.copernicus.org/articles/21/3735/2024/bg-21-3735-2024-t04-thumb.png)
3.3 The optimization performance in COS fluxes
The posterior simulated COS fluxes were evaluated against the prior simulations and observations. Table 4 lists the mean RMSEs and range widths of the prior and posterior simulated COS fluxes for all the sites. The RMSEmean values of the posterior COS simulations are smaller than those of the prior COS simulations, and the mean RMSE reduction for all sites is 35.51 % ± 13.72 % (mean ± SD). At the same time, the simulation range widths of COS fluxes are also well constrained, with a mean reduction of 81.46 % ± 7.32 % from the prior. The reduction in RMSEmean and range width is particularly significant in US-Wrc, with a value of 61.14 % and 94.70 %, respectively.
![https://bg.copernicus.org/articles/21/3735/2024/bg-21-3735-2024-f04](https://bg.copernicus.org/articles/21/3735/2024/bg-21-3735-2024-f04-thumb.png)
Figure 4Comparison of prior and posterior simulated ecosystem COS fluxes. The ensemble means of the prior (red) and posterior (blue) simulations are plotted around the uncertainty bounds (5th and 95th quantiles). The mean observed COS and its uncertainty (estimated by the standard deviation) are represented by black dots with error bars. The means and uncertainties of these hourly observations and simulations are calculated and plotted on a daily or monthly scale. Error bars are not plotted when more than three-quarters of the observations are missing. The subplot numbers are assigned based on the alphabetical order of the site names.
In Fig. 4, the daily or monthly variations in COS during the observation period at each site are shown. It can be observed that both the prior and the posterior simulations are able to accurately capture the daily variation in or the seasonal cycle of COS across these sites, with the exception of IT-Soy. As IT-Soy is a temporary observatory with no continuous in situ meteorological observations available, the ERA5 meteorological data were used to drive the model for this site, resulting in the simulation not being able to characterize the COS changes very well. Although the simulations perform well in modeling the variations in COS for other sites, our results also suggest that they tend to underestimate the magnitude of COS fluxes at both ends of the growing season (e.g., Fig. 4d). Furthermore, the model markedly underestimates the magnitude of COS during rainy days (DOY 126–134) at ES-Lma (Fig. 4c). These findings suggest substantial deficiencies in modeling the mechanistic process of COS exchange. Nevertheless, it can be stated that the fusion of COS observations with the BEPS model has the capacity to constrain the predictive uncertainty in COS, as evidenced by significantly reduced uncertainty bounds that largely encapsulate observations.
The prior simulations significantly underestimate the COS fluxes at ES-Lma, with the ensemble mean of prior simulations being only about one-third of that of the observations. After optimization, the simulated COS fluxes show a substantial increase and generally align with the observations. However, some observed peaks are still not included in the posterior simulation uncertainty bounds. In contrast, the prior simulations tend to overestimate COS fluxes at forest sites FI-Hyy, US-Ha1, and US-Wrc. At US-Wrc and US-Ha1, the ensemble means of prior simulations are 65.70 % and 64.81 % larger than the observations. The assimilation of COS effectively corrected the overestimation but, at the same time, led to a slight underestimation of the simulated COS for US-Wrc. With the downregulation of COS simulations, the model–observation difference at both ends of the growing season for FI-Hyy further increased. Particularly, significant underestimation is found in the posterior simulations in 2017 for FI-Hyy, despite the posterior simulations showing a remarkable improvement in reproducing COS fluxes over the entire period (2013–2017). As the prior simulations neither noticeably overestimate nor underestimate, there is little difference between the ensemble mean of the prior and posterior simulations at the remaining three sites (AT-Neu, DK-Sor, and IT-Soy). Nevertheless, the assimilation of COS resulted in a remarkable reduction in both RMSEmean and uncertainty bounds for COS simulations at these sites, with mean reductions of 23.93 % and 75.11 %, respectively.
Overall, there are considerable uncertainties in the prior simulations, with the uncertainty bounds comparable to or much larger than the uncertainties in observations, and the ensemble mean strongly deviates from observations in some sites, i.e., ES-Lma. Our results suggest that significant improvement in both the ensemble mean and the predictive uncertainty in COS simulations can be achieved through the addition of the information provided by the COS observations with the Monte Carlo-based parameter optimization approach, especially for evergreen needleleaf forest sites. However, limited by various factors, such as uncertainty in model-driven data and model structure (Cho et al., 2023), currently the model often underestimates the simulation at both ends of the growing season and lacks proficiency in modeling the magnitude of COS during rainy days.
Table 5Comparison of model performance indices for the prior and posterior GPP simulations. The RMSEmean values of the prior and posterior simulations are the mean values of the RMSEs of 20 000 prior GPP simulations and 100 behavioral GPP simulations with GPP observations, respectively. The range widths of the prior and posterior GPP simulations are defined as the mean values of the difference between the 95th and 5th percentiles of the prior and posterior simulations, respectively. The reduction (%) in RMSEmean and range width is calculated as (1−posterior/prior) × 100.
![](https://bg.copernicus.org/articles/21/3735/2024/bg-21-3735-2024-t05-thumb.png)
3.4 The performance of simulated GPP
The mean RMSEs and range widths of both prior and posterior simulated GPP for all sites are presented in Table 5. With reduction ratios of RMSEmean ranging from 20.16 % to 64.12 %, the assimilation of COS effectively enhanced the model performance of GPP to varying degrees. Concurrently, the range widths of GPP simulations were well confined, exhibiting a mean reduction ratio of 65.81 % ± 6.77 %. The maximum reduction in RMSEmean for GPP occurred at US-Wrc, aligning with the substantial improvement observed in the posterior simulated COS at this site. In contrast, a relatively limited impact on improving the prediction of GPP was observed at FI-Hyy, as evidenced by the smaller reduction in both the RMSEmean value and the range width of GPP simulations.
![https://bg.copernicus.org/articles/21/3735/2024/bg-21-3735-2024-f05](https://bg.copernicus.org/articles/21/3735/2024/bg-21-3735-2024-f05-thumb.png)
Figure 5Comparison of prior and posterior simulated GPP. The ensemble means of the prior (red) and posterior (blue) simulations are plotted around the uncertainty bounds (5th and 95th quantiles). The mean observed GPP and its uncertainty (estimated by the standard deviation) are represented by black dots with error bars. The means and uncertainties of these hourly observations and simulations are calculated and plotted on a daily or monthly scale. Error bars are not plotted when more than three-quarters of the observations are missing. The subplot numbers are assigned based on the alphabetical order of the site names.
The BEPS model demonstrated excellent performance in capturing the daily variation in and seasonal cycle of GPP, as illustrated in Fig. 5. However, similar to the COS simulations, the ensemble averages of the prior simulated GPP deviated notably from observations at several sites. For example, at DK-Sor and ES-Lma, the ensemble averages of the prior simulated GPP were only approximately half of those of the observations. After the assimilation of COS, GPP simulations exhibited a significant increase, aligning well with observations at DK-Sor and ES-Lma. Conversely, substantial overestimation of prior GPP simulations was effectively corrected through the assimilation of COS at US-Wrc, resulting in a remarkable enhanced modeling performance in both RMSE and range width. For FI-Hyy and US-Ha1, minimal differences were observed between the ensemble mean of prior and posterior simulations, as the ensemble mean of prior simulated GPP had already consistently fit the observations. Nevertheless, our results highlight notable enhancements in the predictive uncertainty in GPP through COS assimilation at these two sites. In Fig. 5e, it is evident that, likely due to the absence of in situ meteorological data at IT-Soy, GPP trends are not well represented, although the ensemble averages of the GPP simulations are very close to those of the observations in magnitude. However, with a reduction in range width as high as 74.72 %, our finding suggests that the assimilation of COS can significantly reduce the predictive uncertainty in GPP, despite the presence of substantial meteorological data uncertainty.
4.1 Parameter sensitivity
As mentioned before, here we utilize the conductance analog model proposed by Berry et al. (2013) to simulate COS plant uptake. Thus, it is not surprising that both stomatal-conductance-related parameters, and , would impact the modeling of COS flux. Considering the stress of soil moisture on stomatal conductance, the stomatal conductance was calculated by a modified version (Woodward et al., 1995; Ju et al., 2010) of the Ball–Berry model (Ball et al., 1987) based on the close relationship between stomatal conductance and photosynthesis rate. Consequently, both the soil-hydrology-related parameters and the photosynthesis-related parameters can ultimately play a role in the simulation of COS plant uptake by influencing the modeling of the stomatal conductance.
It has been recognized that the photosynthetic capacity simulated by terrestrial ecosystem models is highly sensitive to Vcmax, Jmax, and light conditions (Zaehle et al., 2005; Bonan et al., 2011; Rogers, 2014; Sargsyan et al., 2014; Koffi et al., 2015; Rogers et al., 2017; Xing et al., 2023). Our study corroborates these findings, highlighting the pronounced sensitivity of simulated GPP to Vcmax25, followed by VJslope and fleaf. Moreover, our results reveal that the COS simulations are not notably sensitive to fleaf and VJslope, while Vcmax25 plays a crucial role in the modeling of COS. This is because Vcmax25 not only affects the estimation of stomatal conductance through photosynthesis, but is also used to characterize the apparent conductance for COS uptake from the intercellular airspaces, as both mesophyll conductance and carbonic anhydrase activity tend to scale with Vcmax (Badger and Price, 1994; Evans et al., 1994; Berry et al., 2013). Yet, as the hydrolysis reaction of COS by carbonic anhydrase is not dependent on light, VJslope and fleaf do not play a role in the modeling of apparent conductance and thus have little effect on the simulation of COS.
As the COS plant uptake and photosynthesis are tightly coupled through stomata, one would naturally expect similar sensitivity in simulated COS and GPP to the stomatal-conductance-related parameters and . However, the relationship between COS and stomatal conductance differs significantly from that between GPP and stomatal conductance within the model (e.g., Eq. 5 and the Ball–Berry model). Consequently, a notable difference in sensitivity between simulated GPP and COS to and was identified in this study. Specifically, exhibited more pronounced effects on photosynthesis, while played a crucial role in the simulation of COS.
Given that a significant portion of nitrogen is invested in the photosynthetic machinery (Mu and Chen, 2021), there exists a close association between leaf nitrogen content and leaf photosynthetic capacity (Sage and Pearcy, 1987). Additionally, the well-established relationship between leaf nitrogen content and carboxylation capacity (Kattge et al., 2009; Lu et al., 2022) further emphasizes this connection. Specifically, carboxylation capacity in leaf scale is assumed to be linearly related to leaf nitrogen content in the BEPS model (Medlyn et al., 1999; Chen et al., 2012). Consequently, both Vcmax25 and Nleaf play a crucial role in influencing carboxylation capacity, thus having a substantial impact on the simulation of COS.
The soil-hydrology-related parameters can also affect the simulation of COS plant flux as we take the stress effect of soil moisture on both stomatal conductance and mesophyll conductance into account (Ju et al., 2010; Knauer et al., 2020). These parameters also affect the modeling of COS soil exchange since soil moisture is a significant factor in COS soil biotic flux (Whelan et al., 2016). However, given the smaller magnitude of soil COS exchange compared to plant uptake (Whelan et al., 2018) and the minimal impact of soil moisture stress on photosynthetic capacity (Ma et al., 2022), these soil-hydrology-relevant parameters do not significantly influence the modeling of COS.
![https://bg.copernicus.org/articles/21/3735/2024/bg-21-3735-2024-f06](https://bg.copernicus.org/articles/21/3735/2024/bg-21-3735-2024-f06-thumb.png)
Figure 6Parameter correlation matrix plots with significance levels between the parameters of the behavioral parameter sets. Correlation coefficients are shown in red font when the confidence level is greater than 95 % (p < 0.05), with the superscript “∗” indicating a confidence level greater than 99 % (p < 0.01).
4.2 Parameter interactions
For all seven sites, Pearson correlation coefficients and confidence levels between the selected parameters were calculated, as depicted in Fig. 6. Generally, each site exhibits approximately three to eight parameter combinations with significant correlations (p < 0.05). A total of 8 parameter combinations demonstrate significant correlations at more than one site, while 11 parameter combinations exhibit significant correlations at only one site. Specifically, with a mean correlation coefficient of −0.55 ± 0.14 (negative value representing a negative correlation), the correlations between Vcmax25 and Nleaf are very significant (p < 0.01) at all sites, indicating a robust interaction between them. In addition to Vcmax25 and Nleaf, four parameter combinations show highly significant correlations (p < 0.01) at a minimum of two sites: and , and Vcmax25, and Nleaf, and bscalar and VJslope. Such results indicate the strong interactions among parameters related to stomatal conductance, photosynthesis, and soil hydrology, even if some of them do not significantly impact the modeling of COS.
We observed substantial variations in parameter interactions across different sites. For instance, at AT-Neu, and Nleaf exhibited a highly significant negative correlation with a correlation coefficient as high as −0.45. However, these two parameters seemed irrelevant at DK-Sor, with a correlation coefficient of only −0.03. As for soil-hydrology-related parameters, none of them showed significant correlations with any parameter at IT-Soy, FI-Hyy, and US-Ha1, yet there were four related parameter combinations significantly correlated at DK-Sor (Fig. 6b). Furthermore, while Vcmax25 was highly correlated with Nleaf at all sites, the correlation coefficients varied considerably, ranging from −0.41 to −0.87.
We also observed interactions not only between two parameters but also among several parameters (e.g., in Fig. 6a). Figure B1 showcases the intricate interactions among multiple parameters relevant to COS simulations and illustrates the distribution of behavioral parameter sets. The parameter combinations depicted in Fig. B1 are particularly representative as they originate from diverse sites and include nearly all highly significant correlated combinations. Overall, since the six plant-growth-related parameters used in this study are positively correlated with the simulation of COS, they consistently constrain each other, demonstrating a negative correlation, as shown in Fig. 6a, c, and d. However, due to the nonlinearity of the model, there is not a simple linear relationship between these parameters. For example, at AT-Neu, where the COS observations notably exceed the ensemble mean of prior simulations, Vcmax25 and VJslope exhibit a nonlinear correlation, but both tend to be distributed near their upper limits (Fig. B1).
4.3 Parameter identifiability
As the parameter identifiability is quantified based on the range of the behavioral parameter, its results were presented in Fig. 3 along with the plots of the cumulative likelihood distributions of the behavioral parameters. These results underscore the remarkable ability of COS assimilation to identify , with a mean PI of as high as 0.81 ± 0.28. Identified as the most sensitive parameters for COS modeling, Vcmax25 and Nleaf also exhibit remarkable identifiability, with mean PIs of 0.29 ± 0.19 and 0.26 ± 0.22, respectively. demonstrates varying levels of identifiability, with PIs ranging from 0.01 to 0.28. In contrast, the light-reaction-related parameters VJslope and fleaf are not well identified, with the maximum value of PI of only 0.04. The soil-hydrology-related parameters bscalar, Ksatscalar, and rdecay are also generally unidentifiable. Notably, rdecay is well identified at IT-Soy, in which its PI value (0.14) is approximately 7 times that of the other sites.
In this study, the identifiability of a parameter is closely related to the sensitivity of COS simulations to the parameter, although it is known to be influenced by model over-parameterization and parameter interactions (Gan et al., 2014). For instance, at ES-Lma, where the COS simulations exhibited the greatest sensitivity to Nleaf and Vcmax25, these parameters were also found to have the highest identifiability (Figs. 2a and 3). Given the high sensitivity of COS modeling to Vcmax25, Nleaf, and , it is unsurprising that these parameters can be effectively identified by the assimilation of COS. However, our findings indicate that the sensitivity of Vcmax25 and Nleaf is much greater than that of , yet the latter is much more identifiable. This outcome can be attributed to the highly significant correlation between Vcmax25 and Nleaf, as parameter interaction is a primary contributor to parameter unidentifiability (Gan et al., 2014).
In Sect. 3.1, it is demonstrated that the modeling of COS exhibits a low sensitivity to fleaf, , and VJslope. Consequently, it is reasonable that the assimilation of COS did not effectively identify fleaf, , and VJslope (Fig. 3). However, due to their significant correlations with other plant-growth-related parameters, effective identification is possible in specific cases. Notably, combinations such as and –Vcmax25 are very significantly correlated (Fig. 6b), and both and Vcmax25 are well identified at FI-Hyy. As a result, also attains high identifiability at this site.
It has previously been demonstrated that soil-hydrology-related parameters exert a minimal impact on COS simulations (Fig. 2) and cannot be effectively constrained through COS assimilation in general (Fig. 3). Consequently, these parameters exhibit low identifiability, although significant combinations of correlations associated with soil-hydrology-related parameters were observed at certain sites (e.g., DK-Sor).
4.4 Relationship between COS and GPP simulation performances
In this study, we identified the top 100 parameter sets whose corresponding simulations displayed the smallest RMSE concerning COS observations as the behavioral parameter sets. Subsequently, these behavioral parameter sets were employed to derive the posterior simulated COS and GPP and to estimate prediction uncertainty. Therefore, it is necessary to investigate the distribution of RMSEs for COS simulations and GPP simulations and to understand the relationship between the model performance of COS and that of GPP.
![https://bg.copernicus.org/articles/21/3735/2024/bg-21-3735-2024-f07](https://bg.copernicus.org/articles/21/3735/2024/bg-21-3735-2024-f07-thumb.png)
Figure 7Comparison of RMSE for COS (RMSECOS) and GPP (RMSEGPP) in the Monte Carlo simulations. Each data point represents a parameter set, with color indicating data density. The dashed grey line represents the RMSE threshold for COS simulations, calculated as the mean of the 100th and 101st smallest values of the RMSE.
In Fig. 7, scatterplots of RMSEs for COS and GPP are presented. It can be observed that at most sites, where the scatters are most densely distributed, there tend to be relatively small RMSEs for both COS and GPP. These results indicate that the current model is generally capable of simulating COS and GPP well at the same time. However, given the distinct mechanisms of COS and GPP as well as the uncertainties in model structure and driving data, etc., there are also numerous parameter sets that perform well for either COS or GPP but exhibit significant discrepancies with the observations of the other. For example, the model runs with the 3 % highest RMSE for GPP instead exhibit good performance in terms of COS at ES-Lma, with their mean RMSE values (15.42 ) less than those of the prior (17.10 ). Overall, our results suggest that these behavioral parameter sets, which demonstrate good performance in COS simulation, generally also perform well in modeling GPP. However, the parameter sets with relatively good GPP simulation results exhibit significant variability in the performance of COS modeling.
4.5 Caveats and implication
Compared to the big-leaf model, the two-leaf model has been demonstrated to better describe the canopy radiation distribution, GPP, and stomatal conductance (Luo et al., 2018). In this study, we take advantage of the two-leaf model to simulate COS fluxes from plant and soil based on the BEPS model within the two-leaf framework. Ecosystem COS flux data were used to calibrate the model parameters belonging to BEPS and to optimize GPP simulations among diverse ecosystems within the Monte Carlo-based methodology. Our results demonstrate that COS not only improves the accuracy of GPP simulations but also reduces GPP simulation uncertainty. However, due to the lack of in situ COS concentration and flux data, as well as BEPS model driving data (e.g., meteorological data, LAI data, and clumping index data), we still face challenges in evaluating the performance of the two-leaf model compared to other models in COS simulation. Therefore, there is an urgent need for more in situ meteorological data, vegetation canopy structural parameters, and COS observations.
The spatial and temporal variation in atmospheric COS concentrations has a considerable influence on the COS plant uptake (Ma et al., 2021; Kooijmans et al., 2021) due to the linear relationship between the two (Stimler et al., 2010). With a lack of continuous ground-based COS concentration observations, COS concentrations in the bulk air are regarded as spatially invariant over the globe but assumed to vary annually in this study, which may lead to significant biases in COS simulations. Currently, several recent studies have simulated COS vegetation fluxes based on atmospheric-transport-model-derived COS concentration data within the big-leaf framework (Kooijmans et al., 2021; Maignan et al., 2021; Abadie et al., 2023). These COS fluxes simulated based on big-leaf models were in turn used to drive atmospheric transport models (Remaud et al., 2023; Ma et al., 2023). Within an atmosphere inversion framework, recent studies indicate an underestimation of the biosphere COS sink in high-latitude regions of the Northern Hemisphere (NH) (Ma et al., 2021; Remaud et al., 2023). Larger underestimations of ecosystem COS exchange based on the big-leaf model at high latitudes have also been confirmed at the site scale, and the underestimations of COS are consistent with biases in GPP for some sites (Kooijmans et al., 2021). Interestingly, Luo et al. (2018) demonstrated that the reason for the underestimation of GPP by the big-leaf model is that it fails to accurately describe the instantaneous radiation distribution in the canopy, and the underestimation increases with the increase in LAI. The NH high-latitude regions have relatively high LAI (Fang et al., 2019); therefore the deficiency of the big-leaf model in simulating radiation distribution may contribute to the existence of the missing COS sink in the NH high latitude in summer, and this deficiency is amplified by the larger LAI. In fact, the spatial distribution of LAI (i.e., GLOBMAP LAI) retrieved through remote sensing not only in NH high-latitude regions but also in central Africa aligns with the spatial distribution of the missing sink revealed by the “objective” inversion conducted by Ma et al. (2021) (as illustrated in Fig. 7 in Ma et al., 2021), which further validates the reasonableness of this hypothesis. Therefore, conducting COS simulations under the two-leaf framework at a global scale holds the promise of providing insights into the global COS vegetation sink and benefiting the simulation of the spatial and temporal distribution of COS concentrations. Thus, it is necessary to conduct regional and global COS simulations within the two-leaf model framework in the future.
Taking advantage of the Monte Carlo-based parameter optimization approach, we analyzed the global sensitivity, identifiability, and interactions of COS-related parameters in this study. Furthermore, we quantified the uncertainty in simulated COS and GPP, thereby revealing the capacity of COS to constrain the uncertainty in GPP simulations. However, the Monte Carlo-based parameter optimization approach is subject to controversy (Sambridge and Mosegaard, 2002) due to the numerous subjective decisions involved in its implementation, such as the selection of the parameter range, sample size, and performance metric, etc. Further research is needed to investigate the impact of these factors on the parameter optimization results related to COS and the assessment of model prediction uncertainty.
In this study, ecosystem carbonyl sulfide fluxes were utilized to calibrate the ecosystem model parameters and to optimize GPP simulations among various ecosystems within a Monte Carlo-based approach using COS modeling within BEPS. A global parameter sensitivity analysis was conducted to identify the most sensitive parameters among a set of nine pre-selected parameters. The identifiability and interaction of model parameters were investigated by the behavioral parameter sets. The effectiveness of COS in improving the model performance of GPP was evaluated. The major findings are as follows:
(1) Similar to GPP, we found the modeling of COS is sensitive to parameters Vcmax25 and Nleaf and insensitive to soil-hydrology-related parameters and the energy-related parameter fleaf. Unlike GPP, COS is sensitive to and insensitive to and VJslope.
(2) The assimilation of COS within the Monte Carlo-based approach effectively improved model prediction of GPP and significantly reduced the model predictive uncertainty, with a mean RMSE reduction of 40.56 % ± 13.77 % and a mean range width reduction as high as 65.81 % ± 6.77 %.
(3) Complex and significant two-parameter or multi-parameter interactions exist between the model parameters. Particularly, Vcmax25 and Nleaf show highly significant correlations (p < 0.01) at all sites.
(4) Generally, , Vcmax25, and Nleaf can be well identified through the assimilation of COS, especially , whereas the soil-hydrology-related parameters and the light-reaction-related parameters cannot be identified effectively.
A1 BEPS photosynthesis and stomatal conductance modeling approach
In the BEPS model, the net photosynthesis rate (A) is calculated using the Farquhar model (Farquhar et al., 1980; Chen et al., 1999):
where Ai and Aj are RuBisCO- and RuBP-limited gross photosynthetic rates (), respectively. Rd is the dark leaf respiration (). Vcmax is the maximum carboxylation rate of RuBisCO (); J is the electron transport rate (); Ci and Oi are the intercellular carbon dioxide (CO2) and oxygen (O2) concentrations (mol mol−1), respectively; and Kc and Ko are Michaelis–Menten constants for CO2 and O2 (mol mol−1), respectively.
The electron transport rate, J, is dependent on incident photosynthetic photon flux density (PPFD; ) as
where Jmax is the maximum electron transport rate (), and I is the incident PPFD calculated from the incident shortwave radiation RSW (W m−2):
where β* = 4.55 is the energy–quanta conversion factor (µmol J−1), and fleaf is the ratio of photosynthesis active radiation to the shortwave radiation (unitless).
The maximum carboxylation rate of RuBisCO Vcmax was calculated according to the modified Arrhenius temperature function (Medlyn et al., 2002) and the maximum carboxylation rate of RuBisCO at 25 °C (Vcmax25). Vcmax is generally proportional to the leaf nitrogen content. Considering that both the fractions of sunlit and shaded leaf areas to the total leaf area and the leaf nitrogen content vary with the depth into the canopy, the Vcmax values of sunlit (Vcmax,sunlit) and shaded (Vcmax,shaded) leaves can be obtained through vertical integrations with respect to canopy depth (Chen et al., 2012; De Pury and Farquhar, 1997):
where χn (m2 g−1) is the relative change of Vcmax to leaf nitrogen content; Nleaf (g m−2) is the leaf nitrogen content at the top of the canopy; kn is the leaf nitrogen content decay rate with increasing depth into the canopy, taken as 0.3; and L is the canopy depth described in total LAI. Here, k is calculated as
where G(θ) is the projection coefficient, taken as 0.5.
After Vcmax values for the representative sunlit and shaded leaves are obtained, the maximum electronic transport rates for the sunlit and shaded leaves are obtained from Medlyn et al. (1999):
where VJslope (unitless) is the slope of the relationship of Vcmax and Jmax.
The leaf stomatal conductance to water vapor (gsw in ) is estimated using a modified version of the Ball–Berry (BB) empirical model (Ball et al., 1987) following Woodward et al. (1995):
where is the intercept of the BB model, representing the minimum gsw (); is the empirical slope parameter in the BB model (unitless); Rh is the relative humidity at the leaf surface (unitless); fw is a soil moisture stress factor describing the sensitivity of gsw to soil water availability (Ju et al., 2006); and Ca is the atmospheric CO2 concentration (µmol mol−1).
The soil water availability factor fw,i in each layer i is calculated as
where fi(ψi) is a function of matrix suction ψi (m) (Zierl, 2001), and fi(Ts,i) is a function describing the effect of soil temperature (Ts,i in °C) on soil water uptake (Bonan, 1991).
To consider the variable soil water potential at different depths, the scheme of Ju et al. (2006) was employed to calculate the weight of each layer (wi) to fw:
where n is the number of soil layers (five were used in this study) of the BEPS model, and Ri is the root fraction in layer i, calculated as
where cdi is the cumulative depth (m) of layer i. In this study, each soil layer depth (from top to bottom) of the BEPS model is 0.05, 0.10, 0.20, 0.40, and 1.25 m, respectively.
The overall soil water availability fw is then calculated as
The hydraulic conductivity of each soil layer Ki (ms−1) is expressed as
where Ksati is the saturated hydrological conductivity of soil layer i (m s−1); SWCi is the volumetric liquid soil water content of soil layer i (m s−1); θs,i is the porosity of soil layer i (unitless); and bi is the Campbell parameter for soil layer i, determining the change rate of hydraulic conductivity with SWC (unitless). In this study, Ksati and bi are expressed as
where Ksatdf,i and bdf,i are the default values of Ksati and bi, respectively.
A2 BEPS soil COS modeling approach
The total soil COS flux FCOS,soil is the sum of the abiotic COS flux FCOS,abiotic and biotic COS flux FCOS,biotic:
Here, we take the approach developed in Whelan et al. (2016) for the modeling of FCOS,soil. In this approach, FCOS,abiotic is described as an exponential function of the temperature of soil Tsoil (°C):
where α and β are parameters determined using the least-squares fitting approach. We assigned the values of α and β to BEPS according to the parameterization scheme of Whelan et al. (2016).
FCOS,biotic is described as the product of a power function and an exponential function:
Here, a is the curve shape constant. The maximum biotic COS uptake Fopt and the biotic COS uptake are the COS fluxes () at optimum soil water content SWCopt and a secondary soil water content SWCg and SWCg>SWCopt. A more detailed description of the modeling of FCOS,biotic and the parameterization scheme adopted in this study can be found in Whelan et al. (2022).
![https://bg.copernicus.org/articles/21/3735/2024/bg-21-3735-2024-f08](https://bg.copernicus.org/articles/21/3735/2024/bg-21-3735-2024-f08-thumb.png)
Figure B1Scatterplots showing the behavioral parameter sets in 3D parameter space at AT-Neu (a), FI-Hyy (b) and US-Ha1 (c), and US-Wrc. The scatter colors represent the magnitude of the corresponding parameters using red, green, and blue (RGB) values. The projection of the scatter is illustrated with smaller markers.
Measured eddy covariance carbonyl sulfide flux data can be found at https://doi.org/10.5281/zenodo.3406990 (Spielmann et al., 2019b) for AT-Neu, DK-Sor, ES-Lma, and IT-Soy; at https://doi.org/10.5281/zenodo.6940750 (Kohonen et al., 2022b) for FI-Hyy; and from the Harvard Forest Data Archive under record HF214 (https://doi.org/10.6073/pasta/7ed7b4d1fc7ad303998e76143a3b279a, Commane et al., 2016) for US-Ha1. The raw COS concentration data of US-Wrc can be obtained at https://doi.org/10.5281/zenodo.1422820 (Rastogi et al., 2018b). The meteorological data can be obtained from the FLUXNET database (https://fluxnet.org/, last access: June 2022) for AT-Neu (https://doi.org/10.18140/FLX/1440121, Wohlfahrt et al., 2020), DK-Sor (https://doi.org/10.18140/FLX/1440155, Ibrom and Pilegaard, 2020), ES-Lma (https://doi.org/10.18160/FDSD-GVRS, Migliavacca et al., 2020), FI-Hyy (https://doi.org/10.18140/FLX/1440158, Mammarella et al., 2020), and US-Ha1 (https://doi.org/10.18140/FLX/1440071, Harvard University, 2020); from the AmeriFlux database (https://ameriflux.lbl.gov/, last access: June 2022) for US-Ha1 (https://doi.org/10.17190/AMF/1871137, Munger, 2022) and US-Wrc (https://doi.org/10.17190/AMF/1246114, Wharton, 2016); and from the ERA5 dataset (https://doi.org/10.24381/cds.adbb2d47, Hersbach et al., 2023) for AT-Neu, IT-Soy, and US-Ha1. The GPP data can be obtained from the FLUXNET database for DK-Sor, ES-Lma, FI-Hyy, and US-Ha1; from the AmeriFlux database for US-Ha1; from https://doi.org/10.5281/zenodo.3406990 (Spielmann et al., 2019b) for AT-Neu and IT-Soy; and from https://doi.org/10.5281/zenodo.1422820 (Rastogi et al., 2018b) for US-Wrc. The GLASS LAI is available at https://doi.org/10.12041/geodata.GLASS_LAI_MODIS(0.05D).ver1.db (Xiao et al., 2016), and the GLOBMAP LAI is available at https://doi.org/10.5281/zenodo.4700264 (Liu et al., 2021).
MW designed the experiments and developed the model. XX improved the model and performed the Monte Carlo simulations. HZ carried out the analysis and wrote the original manuscript. All the authors contributed to the writing of the paper.
The contact author has declared that none of the authors has any competing interests.
Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors.
This article is part of the special issue “Land surface–atmosphere interactions – from the microbial to the global scale”. It is not associated with a conference.
The authors thank the reviewers for their constructive and useful comments. We thank everyone that contributed to the collection of data. We acknowledge Tim Moore from McGill University for helping us improve the language in an earlier version of the paper.
This study has been supported by the National Key Research and Development Program of China (2023YFB3907402); the National Natural Science Foundation of China (42371486); and the Research Funds for the Frontiers Science Center for Critical Earth Material Cycling, Nanjing University (grant nos. 0209-14380115 and 0904-14380031).
This paper was edited by Elizabeth Keller and reviewed by two anonymous referees.
Abadie, C., Maignan, F., Remaud, M., Ogée, J., Campbell, J. E., Whelan, M. E., Kitz, F., Spielmann, F. M., Wohlfahrt, G., Wehr, R., Sun, W., Raoult, N., Seibt, U., Hauglustaine, D., Lennartz, S. T., Belviso, S., Montagne, D., and Peylin, P.: Global modelling of soil carbonyl sulfide exchanges, Biogeosciences, 19, 2427–2463, https://doi.org/10.5194/bg-19-2427-2022, 2022.
Abadie, C., Maignan, F., Remaud, M., Kohonen, K.-M., Sun, W., Kooijmans, L., Vesala, T., Seibt, U., Raoult, N., Bastrikov, V., Belviso, S., and Peylin, P.: Carbon and Water Fluxes of the Boreal Evergreen Needleleaf Forest Biome Constrained by Assimilating Ecosystem Carbonyl Sulfide Flux Observations, J. Geophys. Res.-Biogeo., 128, e2023JG007407, https://doi.org/10.1029/2023JG007407, 2023.
Asaf, D., Rotenberg, E., Tatarinov, F., Dicken, U., Montzka, S. A., and Yakir, D.: Ecosystem photosynthesis inferred from measurements of carbonyl sulphide flux, Nat. Geosci., 6, 186–190, 2013.
Badger, M. R. and Price, G. D.: The role of carbonic anhydrase in photosynthesis, Annu Rev. Plant Biol., 45, 369–392, 1994.
Ball, J. T., Woodrow, I. E., and Berry, J. A.: A Model Predicting Stomatal Conductance and its Contribution to the Control of Photosynthesis under Different Environmental Conditions, in: Progress in Photosynthesis Research: Volume 4 Proceedings of the VIIth International Congress on Photosynthesis Providence, Rhode Island, USA, August 10–15, 1986, edited by: Biggins, J., Springer Netherlands, Dordrecht, 221–224, https://doi.org/10.1007/978-94-017-0519-6_48, 1987.
Bao, S., Ibrom, A., Wohlfahrt, G., Koirala, S., Migliavacca, M., Zhang, Q., and Carvalhais, N.: Narrow but robust advantages in two-big-leaf light use efficiency models over big-leaf light use efficiency models at ecosystem level, Agr. Forest Meteorol., 326, 109185, https://doi.org/10.1016/j.agrformet.2022.109185, 2022.
Berry, J., Wolf, A., Campbell, J. E., Baker, I., Blake, N., Blake, D., Denning, A. S., Kawa, S. R., Montzka, S. A., Seibt, U., Stimler, K., Yakir, D., and Zhu, Z.: A coupled model of the global cycles of carbonyl sulfide and CO2: A possible new window on the carbon cycle, J. Geophys Res.-Biogeo., 118, 842–852, https://doi.org/10.1002/jgrg.20068, 2013.
Beven, K. and Binley, A.: The future of distributed models: model calibration and uncertainty prediction, Hydrol. Process., 6, 279–298, 1992.
Beven, K. and Binley, A.: GLUE: 20 years on, Hydrol. Process., 28, 5897–5918, 2014.
Beven, K. and Freer, J.: Equifinality, data assimilation, and uncertainty estimation in mechanistic modelling of complex environmental systems using the GLUE methodology, J. Hydrol., 249, 11–29, 2001.
Blankenship, R. E.: Molecular mechanisms of photosynthesis, John Wiley & Sons, https://doi.org/10.1002/9780470758472, 2021.
Blasone, R.-S., Vrugt, J. A., Madsen, H., Rosbjerg, D., Robinson, B. A., and Zyvoloski, G. A.: Generalized likelihood uncertainty estimation (GLUE) using adaptive Markov Chain Monte Carlo sampling, Adv. Water Resour., 31, 630–648, 2008.
Bonan, G. B.: A biophysical surface energy budget analysis of soil temperature in the boreal forests of interior Alaska, Water Resour. Res., 27, 767–781, 1991.
Bonan, G. B., Lawrence, P. J., Oleson, K. W., Levis, S., Jung, M., Reichstein, M., Lawrence, D. M., and Swenson, S. C.: Improving canopy processes in the Community Land Model version 4 (CLM4) using global flux fields empirically inferred from FLUXNET data, J. Geophys. Res.-Biogeo., 116, G02014, https://doi.org/10.1029/2010JG001593, 2011.
Borgonovo, E.: A new uncertainty importance measure, Reliab. Eng. Syst. Safe., 92, 771–784, 2007.
Canadell, J. G., Mooney, H. A., Baldocchi, D. D., Berry, J. A., Ehleringer, J. R., Field, C. B., Gower, S. T., Hollinger, D. Y., Hunt, J. E., Jackson, R. B., Running, S. W., Shaver, G. R., Steffen, W., Trumbore, S. E., Valentini, R., and Bond, B. Y.: Commentary: Carbon Metabolism of the Terrestrial Biosphere: A Multitechnique Approach for Improved Understanding, Ecosystems, 3, 115–130, https://doi.org/10.1007/s100210000014, 2000.
Chen, B., Wang, P., Wang, S., Ju, W., Liu, Z., and Zhang, Y.: Simulating canopy carbonyl sulfide uptake of two forest stands through an improved ecosystem model and parameter optimization using an ensemble Kalman filter, Ecol. Model., 475, 110212, https://doi.org/10.1016/j.ecolmodel.2022.110212, 2023.
Chen, J., Liu, J., Cihlar, J., and Goulden, M.: Daily canopy photosynthesis model through temporal and spatial scaling for remote sensing applications, Ecol. Model., 124, 99–119, 1999.
Chen, J. M., Mo, G., Pisek, J., Liu, J., Deng, F., Ishizawa, M., and Chan, D.: Effects of foliage clumping on the estimation of global terrestrial gross primary productivity, Global Biogeochem. Cy., 26, GB1019, https://doi.org/10.1029/2010GB003996, 2012.
Chen, J. M., Ju, W., Ciais, P., Viovy, N., Liu, R., Liu, Y., and Lu, X.: Vegetation structural change since 1981 significantly enhanced the terrestrial carbon sink, Nat. Commun., 10, 4259, https://doi.org/10.1038/s41467-019-12257-8, 2019.
Cho, A., Kooijmans, L. M. J., Kohonen, K.-M., Wehr, R., and Krol, M. C.: Optimizing the carbonic anhydrase temperature response and stomatal conductance of carbonyl sulfide leaf uptake in the Simple Biosphere model (SiB4), Biogeosciences, 20, 2573–2594, https://doi.org/10.5194/bg-20-2573-2023, 2023.
Commane, R., Meredith, L. K., Baker, I. T., Berry, J. A., Munger, J. W., Montzka, S. A., Templer, P. H., Juice, S. M., Zahniser, M. S., and Wofsy, S. C.: Seasonal fluxes of carbonyl sulfide in a midlatitude forest, P. Natl. Acad. Sci. USA, 112, 14162–14167, 2015.
Commane, R., Wofsy, S., and Weir, R.: Fluxes of Carbonyl Sulfide at Harvard Forest EMS Tower since 2010 ver 4, Environmental Data Initiative [data set], https://doi.org/10.6073/pasta/7ed7b4d1fc7ad303998e76143a3b279a, 2016.
De Pury, D. and Farquhar, G.: Simple scaling of photosynthesis from leaves to canopies without the errors of big-leaf models, Plant Cell Environ., 20, 537–557, 1997.
Evans, J. R., Caemmerer, S., Setchell, B. A., and Hudson, G. S.: The relationship between CO2 transfer conductance and leaf anatomy in transgenic tobacco with a reduced content of Rubisco, Functional Plant Biology, 21, 475–495, 1994.
Fang, H., Baret, F., Plummer, S., and Schaepman-Strub, G.: An overview of global leaf area index (LAI): Methods, products, validation, and applications, Rev. Geophys., 57, 739–799, 2019.
Farquhar, G. D., von Caemmerer, S. v., and Berry, J. A.: A biochemical model of photosynthetic CO2 assimilation in leaves of C3 species, Planta, 149, 78–90, 1980.
Friedlingstein, P., O'Sullivan, M., Jones, M. W., Andrew, R. M., Gregor, L., Hauck, J., Le Quéré, C., Luijkx, I. T., Olsen, A., Peters, G. P., Peters, W., Pongratz, J., Schwingshackl, C., Sitch, S., Canadell, J. G., Ciais, P., Jackson, R. B., Alin, S. R., Alkama, R., Arneth, A., Arora, V. K., Bates, N. R., Becker, M., Bellouin, N., Bittig, H. C., Bopp, L., Chevallier, F., Chini, L. P., Cronin, M., Evans, W., Falk, S., Feely, R. A., Gasser, T., Gehlen, M., Gkritzalis, T., Gloege, L., Grassi, G., Gruber, N., Gürses, Ö., Harris, I., Hefner, M., Houghton, R. A., Hurtt, G. C., Iida, Y., Ilyina, T., Jain, A. K., Jersild, A., Kadono, K., Kato, E., Kennedy, D., Klein Goldewijk, K., Knauer, J., Korsbakken, J. I., Landschützer, P., Lefèvre, N., Lindsay, K., Liu, J., Liu, Z., Marland, G., Mayot, N., McGrath, M. J., Metzl, N., Monacci, N. M., Munro, D. R., Nakaoka, S.-I., Niwa, Y., O'Brien, K., Ono, T., Palmer, P. I., Pan, N., Pierrot, D., Pocock, K., Poulter, B., Resplandy, L., Robertson, E., Rödenbeck, C., Rodriguez, C., Rosan, T. M., Schwinger, J., Séférian, R., Shutler, J. D., Skjelvan, I., Steinhoff, T., Sun, Q., Sutton, A. J., Sweeney, C., Takao, S., Tanhua, T., Tans, P. P., Tian, X., Tian, H., Tilbrook, B., Tsujino, H., Tubiello, F., van der Werf, G. R., Walker, A. P., Wanninkhof, R., Whitehead, C., Willstrand Wranne, A., Wright, R., Yuan, W., Yue, C., Yue, X., Zaehle, S., Zeng, J., and Zheng, B.: Global Carbon Budget 2022, Earth Syst. Sci. Data, 14, 4811–4900, https://doi.org/10.5194/essd-14-4811-2022, 2022.
Gan, Y., Duan, Q., Gong, W., Tong, C., Sun, Y., Chu, W., Ye, A., Miao, C., and Di, Z.: A comprehensive evaluation of various sensitivity analysis methods: A case study with a hydrological model, Environ. Modell. Softw., 51, 269–285, 2014.
Gu, L., Baldocchi, D., Verma, S. B., Black, T., Vesala, T., Falge, E. M., and Dowty, P. R.: Advantages of diffuse radiation for terrestrial ecosystem productivity, J. Geophys. Res.-Atmos., 107, ACL 2-1–ACL 2-23, 2002.
Harvard University: FLUXNET2015 US-Ha1 Harvard Forest EMS Tower (HFR1), FLUXNET [data set], https://doi.org/10.18140/FLX/1440071, 2020.
Haynes, K., Baker, I., and Denning, S.: Simple biosphere model version 4.2 (SiB4) technical description, Colorado State University, Fort Collins, CO, USA, 2020.
He, H., Jansson, P.-E., Svensson, M., Meyer, A., Klemedtsson, L., and Kasimir, Å.: Factors controlling Nitrous Oxide emission from a spruce forest ecosystem on drained organic soil, derived using the CoupModel, Ecol. Model., 321, 46–63, 2016.
He, Q., Ju, W., Dai, S., He, W., Song, L., Wang, S., Li, X., and Mao, G.: Drought risk of global terrestrial gross primary productivity over the last 40 years detected by a remote sensing-driven process model, J. Geophys. Res.-Biogeo., 126, e2020JG005944, 2021.
Hersbach, H., Bell, B., Berrisford, P., Biavati, G., Horányi, A., Muñoz Sabater, J., Nicolas, J., Peubey, C., Radu, R., Rozum, I., Schepers, D., Simmons, A., Soci, C., Dee, D., and Thépaut, J.-N.: ERA5 hourly data on single levels from 1940 to present, Copernicus Climate Change Service (C3S) Climate Data Store (CDS) [data set], https://doi.org/10.24381/cds.adbb2d47, 2023.
Hilton, T. W., Whelan, M. E., Zumkehr, A., Kulkarni, S., Berry, J. A., Baker, I. T., Montzka, S. A., Sweeney, C., Miller, B. R., and Elliott Campbell, J.: Peak growing season gross uptake of carbon in North America is largest in the Midwest USA, Nat. Clim. Change, 7, 450–454, 2017.
Houska, T., Multsch, S., Kraft, P., Frede, H.-G., and Breuer, L.: Monte Carlo-based calibration and uncertainty analysis of a coupled plant growth and hydrological model, Biogeosciences, 11, 2069–2082, https://doi.org/10.5194/bg-11-2069-2014, 2014.
Hu, L., Montzka, S. A., Kaushik, A., Andrews, A. E., Sweeney, C., Miller, J., Baker, I. T., Denning, S., Campbell, E., Shiga, Y. P., Tans, P., Siso, M. C., Crotwell, M., McKain, K., Thoning, K., Hall, B., Vimont, I., Elkins, J. W., Whelan, M. E., and Suntharalingam, P.: COS-derived GPP relationships with temperature and light help explain high-latitude atmospheric CO2 seasonal cycle amplification, P. Natl. Acad. Sci., 118, e2103423118, https://doi.org/10.1073/pnas.2103423118, 2021.
Ibrom, A. and Pilegaard, K.: FLUXNET2015 DK-Sor Soroe, FLUXNET2015 [data set], https://doi.org/10.18140/FLX/1440155, 2020.
Iwanaga, T., Usher, W., and Herman, J.: Toward SALib 2.0: Advancing the accessibility and interpretability of global sensitivity analyses, Socio-Environmental Systems Modelling, 4, 18155–18155, 2022.
Jackson, R. B., Canadell, J., Ehleringer, J. R., Mooney, H. A., Sala, O. E., and Schulze, E.-D.: A global analysis of root distributions for terrestrial biomes, Oecologia, 108, 389–411, 1996.
Ju, W., Chen, J. M., Black, T. A., Barr, A. G., Liu, J., and Chen, B.: Modelling multi-year coupled carbon and water fluxes in a boreal aspen forest, Agr. Forest Meteorol., 140, 136–151, 2006.
Ju, W., Gao, P., Wang, J., Zhou, Y., and Zhang, X.: Combining an ecological model with remote sensing and GIS techniques to monitor soil water content of croplands with a monsoon climate, Agr. Water Manage., 97, 1221–1231, 2010.
Karu, E., Li, M., Ernle, L., Brenninkmeijer, C. A., Lelieveld, J., and Williams, J.: Carbonyl Sulfide (OCS) in the upper troposphere/Lowermost stratosphere (UT/LMS) region: Estimates of lifetimes and fluxes, Geophys. Res. Lett., 50, e2023GL105826, https://doi.org/10.1029/2023GL105826, 2023.
Kattge, J., Knorr, W., Raddatz, T., and Wirth, C.: Quantifying photosynthetic capacity and its relationship to leaf nitrogen content for global-scale terrestrial biosphere models, Glob. Change Biol., 15, 976–991, 2009.
Kesselmeier, J., Teusch, N., and Kuhn, U.: Controlling variables for the uptake of atmospheric carbonyl sulfide by soil, J. Geophys. Res.-Atmos., 104, 11577–11584, 1999.
Knauer, J., Zaehle, S., De Kauwe, M. G., Haverd, V., Reichstein, M., and Sun, Y.: Mesophyll conductance in land surface models: effects on photosynthesis and transpiration, Plant J., 101, 858–873, 2020.
Koffi, E. N., Rayner, P. J., Norton, A. J., Frankenberg, C., and Scholze, M.: Investigating the usefulness of satellite-derived fluorescence data in inferring gross primary productivity within the carbon cycle data assimilation system, Biogeosciences, 12, 4067–4084, https://doi.org/10.5194/bg-12-4067-2015, 2015.
Kohonen, K.-M., Kolari, P., Kooijmans, L. M. J., Chen, H., Seibt, U., Sun, W., and Mammarella, I.: Towards standardized processing of eddy covariance flux measurements of carbonyl sulfide, Atmos. Meas. Tech., 13, 3957–3975, https://doi.org/10.5194/amt-13-3957-2020, 2020.
Kohonen, K.-M., Dewar, R., Tramontana, G., Mauranen, A., Kolari, P., Kooijmans, L. M. J., Papale, D., Vesala, T., and Mammarella, I.: Intercomparison of methods to estimate gross primary production based on CO2 and COS flux measurements, Biogeosciences, 19, 4067–4088, https://doi.org/10.5194/bg-19-4067-2022, 2022a.
Kohonen, K.-M., Tramontana, G., and Kolari, P.: Dataset for “Intercomparison of methods to estimate gross primary production based on CO2 and COS flux measurements”, Zenodo [data set], https://doi.org/10.5281/zenodo.6940750, 2022b.
Kooijmans, L. M., Sun, W., Aalto, J., Erkkilä, K.-M., Maseyk, K., Seibt, U., Vesala, T., Mammarella, I., and Chen, H.: Influences of light and humidity on carbonyl sulfide-based estimates of photosynthesis, P. Natl. Acad. Sci. USA, 116, 2470–2475, 2019.
Kooijmans, L. M. J., Cho, A., Ma, J., Kaushik, A., Haynes, K. D., Baker, I., Luijkx, I. T., Groenink, M., Peters, W., Miller, J. B., Berry, J. A., Ogée, J., Meredith, L. K., Sun, W., Kohonen, K.-M., Vesala, T., Mammarella, I., Chen, H., Spielmann, F. M., Wohlfahrt, G., Berkelhammer, M., Whelan, M. E., Maseyk, K., Seibt, U., Commane, R., Wehr, R., and Krol, M.: Evaluation of carbonyl sulfide biosphere exchange in the Simple Biosphere Model (SiB4), Biogeosciences, 18, 6547–6565, https://doi.org/10.5194/bg-18-6547-2021, 2021.
Lasslop, G., Migliavacca, M., Bohrer, G., Reichstein, M., Bahn, M., Ibrom, A., Jacobs, C., Kolari, P., Papale, D., Vesala, T., Wohlfahrt, G., and Cescatti, A.: On the choice of the driving temperature for eddy-covariance carbon dioxide flux partitioning, Biogeosciences, 9, 5243–5259, https://doi.org/10.5194/bg-9-5243-2012, 2012.
Launois, T., Peylin, P., Belviso, S., and Poulter, B.: A new model of the global biogeochemical cycle of carbonyl sulfide – Part 2: Use of carbonyl sulfide to constrain gross primary productivity in current vegetation models, Atmos. Chem. Phys., 15, 9285–9312, https://doi.org/10.5194/acp-15-9285-2015, 2015.
Liu, J., Chen, J., Cihlar, J., and Park, W.: A process-based boreal ecosystem productivity simulator using remote sensing inputs, Remote Sens. Environ., 62, 158–175, 1997.
Liu, R., Liu, Y., and Chen, J.: GLOBMAP global Leaf Area Index since 1981 (Version 3.0), Zenodo [data set], https://doi.org/10.5281/zenodo.4700264, 2021.
Liu, Y., Liu, R., and Chen, J. M.: Retrospective retrieval of long-term consistent global leaf area index (1981–2011) from combined AVHRR and MODIS data, J. Geophys. Res.-Biogeo., 117, G04003, https://doi.org/10.1029/2012JG002084, 2012.
Liu, Y., Xiao, J., Ju, W., Zhou, Y., Wang, S., and Wu, X.: Water use efficiency of China's terrestrial ecosystems and responses to drought, Sci. Rep.-UK, 5, 13799, https://doi.org/10.1038/srep13799, 2015.
Liu, Z., Zhou, Y., Ju, W., and Gao, P.: Simulation of soil water content in farm lands with the BEPS ecological model, Transactions of the Chinese Society of Agricultural Engineering, 27, 67–72, 2011.
Lloyd, J. and Taylor, J. A.: On the Temperature Dependence of Soil Respiration, Funct. Ecol., 8, 315–323, https://doi.org/10.2307/2389824, 1994.
Lu, X., Wang, Y.-P., Ziehn, T., and Dai, Y.: An efficient method for global parameter sensitivity analysis and its applications to the Australian community land surface model (CABLE), Agr. Forest Meteorol., 182, 292–303, 2013.
Lu, X., Croft, H., Chen, J. M., Luo, Y., and Ju, W.: Estimating photosynthetic capacity from optimized Rubisco–chlorophyll relationships among vegetation types and under global change, Enviro. Res. Lett., 17, 014028, https://doi.org/10.1088/1748-9326/ac444d, 2022.
Luo, X., Chen, J. M., Liu, J., Black, T. A., Croft, H., Staebler, R., He, L., Arain, M. A., Chen, B., Mo, G., Gonsamo, A., and McCaughey, H.: Comparison of Big-Leaf, Two-Big-Leaf, and Two-Leaf Upscaling Schemes for Evapotranspiration Estimation Using Coupled Carbon-Water Modeling, J. Geophys. Res-Biogeo., 123, 207–225, https://doi.org/10.1002/2017JG003978, 2018.
Luo, Y.: Terrestrial carbon–cycle feedback to climate warming, Annu. Rev. Ecol. Evol. S., 38, 683–712, 2007.
Ma, J., Kooijmans, L. M. J., Cho, A., Montzka, S. A., Glatthor, N., Worden, J. R., Kuai, L., Atlas, E. L., and Krol, M. C.: Inverse modelling of carbonyl sulfide: implementation, evaluation and implications for the global budget, Atmos. Chem. Phys., 21, 3507–3529, https://doi.org/10.5194/acp-21-3507-2021, 2021.
Ma, J., Remaud, M., Peylin, P., Patra, P., Niwa, Y., Rodenbeck, C., Cartwright, M., Harrison, J. J., Chipperfield, M. P., Pope, R. J., Wilson, C., Belviso, S., Montzka, S. A., Vimont, I., Moore, F., Atlas, E. L., Schwartz, E., and Krol, M. C.: Intercomparison of Atmospheric Carbonyl Sulfide (TransCom-COS): 2. Evaluation of Optimized Fluxes Using Ground-Based and Aircraft Observations, J. Geophys. Res.-Atmos., 128, e2023JD039198, https://doi.org/10.1029/2023JD039198, 2023.
Ma, R., Xiao, J., Liang, S., Ma, H., He, T., Guo, D., Liu, X., and Lu, H.: Pixel-level parameter optimization of a terrestrial biosphere model for improving estimation of carbon fluxes with an efficient model–data fusion method and satellite-derived LAI and GPP data, Geosci. Model Dev., 15, 6637–6657, https://doi.org/10.5194/gmd-15-6637-2022, 2022.
MacBean, N., Bacour, C., Raoult, N., Bastrikov, V., Koffi, E. N., Kuppel, S., Maignan, F., Ottlé, C., Peaucelle, M., Santaren, D., and Peylin, P.: Quantifying and Reducing Uncertainty in Global Carbon Cycle Predictions: Lessons and Perspectives From 15 Years of Data Assimilation Studies With the ORCHIDEE Terrestrial Biosphere Model, Global Biogeochem. Cy., 36, e2021GB007177, https://doi.org/10.1029/2021GB007177, 2022.
Maignan, F., Abadie, C., Remaud, M., Kooijmans, L. M. J., Kohonen, K.-M., Commane, R., Wehr, R., Campbell, J. E., Belviso, S., Montzka, S. A., Raoult, N., Seibt, U., Shiga, Y. P., Vuichard, N., Whelan, M. E., and Peylin, P.: Carbonyl sulfide: comparing a mechanistic representation of the vegetation uptake in a land surface model and the leaf relative uptake approach, Biogeosciences, 18, 2917–2955, https://doi.org/10.5194/bg-18-2917-2021, 2021.
Mammarella, I., Keronen, P., Kolari, P., Launiainen, S., Pumpanen, J., Rannik, Ü., Siivola, E., Levula, J., Pohja, T., and Vesala, T.: FLUXNET2015 FI-Hyy Hyytiala, FLUXNET2015 [data set], https://doi.org/10.18140/FLX/1440158, 2020.
Medlyn, B. E., Badeck, F.-W., De Pury, D. G. G., Barton, C. V. M., Broadmeadow, M., Ceulemans, R., De Angelis, P., Forstreuter, M., Jach, M. E., Kellomäki, S., Laitat, E., Marek, M., Philippot, S., Rey, A., Strassemeyer, J., Laitinen, K., Liozon, R., Portier, B., Roberntz, P., Wang, K., and Jstbid, P. G.: Effects of elevated [CO2] on photosynthesis in European forest species: a meta-analysis of model parameters, Plant Cell Environ., 22, 1475–1495, https://doi.org/10.1046/j.1365-3040.1999.00523.x, 1999.
Medlyn, B. E., Dreyer, E., Ellsworth, D., Forstreuter, M., Harley, P. C., Kirschbaum, M. U. F., Le Roux, X., Montpied, P., Strassemeyer, J., Walcroft, A., Wang, K., and Loustau, D.: Temperature response of parameters of a biochemically based model of photosynthesis. II. A review of experimental data, Plant Cell Environ., 25, 1167–1179, https://doi.org/10.1046/j.1365-3040.2002.00891.x, 2002.
Migliavacca, M., El-Madany, T. S., Carrara, A., Reichstein, M., and ICOS Ecosystem Thematic Centre: Drought-2018 ecosystem eddy covariance flux product from Majadas del Tietar North, ICOS [data set], https://doi.org/10.18160/FDSD-GVRS, 2020.
Miner, G. L., Bauerle, W. L., and Baldocchi, D. D.: Estimating the sensitivity of stomatal conductance to photosynthesis: a review, Plant Cell Environ., 40, 1214–1238, 2017.
Mo, X., Chen, J. M., Ju, W., and Black, T. A.: Optimization of ecosystem model parameters through assimilating eddy covariance flux data with an ensemble Kalman filter, Ecol. Model., 217, 157–173, 2008.
Montzka, S., Calvert, P., Hall, B., Elkins, J., Conway, T., Tans, P., and Sweeney, C.: On the global distribution, seasonality, and budget of atmospheric carbonyl sulfide (COS) and some similarities to CO2, J. Geophys. Res.-Atmos., 112, D09302, https://doi.org/10.1029/2006JD007665, 2007.
Moradkhani, H., Hsu, K. L., Gupta, H., and Sorooshian, S.: Uncertainty assessment of hydrologic model states and parameters: Sequential data assimilation using the particle filter, Water Resour. Res., 41, W05012, https://doi.org/10.1029/2004WR003604, 2005.
Mu, X. and Chen, Y.: The physiological response of photosynthesis to nitrogen deficiency, Plant Physiol. Bioch., 158, 76–82, 2021.
Munger, J. W.: AmeriFlux FLUXNET-1F US-Ha1 Harvard Forest EMS Tower (HFR1), AmeriFlux AMP [data set], https://doi.org/10.17190/AMF/1871137, 2022.
Ogée, J., Sauze, J., Kesselmeier, J., Genty, B., Van Diest, H., Launois, T., and Wingate, L.: A new mechanistic framework to predict OCS fluxes from soils, Biogeosciences, 13, 2221–2240, https://doi.org/10.5194/bg-13-2221-2016, 2016.
Pignon, C. P., Jaiswal, D., McGrath, J. M., and Long, S. P.: Loss of photosynthetic efficiency in the shade. An Achilles heel for the dense modern stands of our most productive C4 crops?, J. Exp. Bot., 68, 335–345, 2017.
Plischke, E., Borgonovo, E., and Smith, C. L.: Global sensitivity measures from given data, Eur. J. Oper. Res., 226, 536–550, 2013.
Protoschill-Krebs, G., Wilhelm, C., and Kesselmeier, J.: Consumption of carbonyl sulphide (COS) by higher plant carbonic anhydrase (CA), Atmos. Environ., 30, 3151–3156, 1996.
Raines, C. A.: The Calvin cycle revisited, Photosynth. Res., 75, 1–10, 2003.
Rastogi, B., Berkelhammer, M., Wharton, S., Whelan, M. E., Itter, M. S., Leen, J. B., Gupta, M. X., Noone, D., and Still, C. J.: Large uptake of atmospheric OCS observed at a moist old growth forest: Controls and implications for carbon cycle applications, J. Geophys. Res.-Biogeo., 123, 3424–3438, 2018a.
Rastogi, B., Berkelhammer, M., Wharton, S., Whelan, E. M., Itter, S. M., Leen, B. J., Gupta, X. M., Noone, D., and Still, C. J.: Large uptake of atmospheric OCS observed at a moist old growth forest: Controls and implications for carbon cycle applications (Version 1), Zenodo [data set], https://doi.org/10.5281/zenodo.1422820, 2018b.
Reichstein, M., Falge, E., Baldocchi, D., Papale, D., Aubinet, M., Berbigier, P., Bernhofer, C., Buchmann, N., Gilmanov, T., Granier, A., Grünwald, T., Havránková, K., Ilvesniemi, H., Janous, D., Knohl, A., Laurila, T., Lohila, A., Loustau, D., Matteucci, G., Meyers, T., Miglietta, F., Ourcival, J.-M., Pumpanen, J., Rambal, S., Rotenberg, E., Sanz, M., Tenhunen, J., Seufert, G., Vaccari, F., Vesala, T., Yakir, D., and Valentini, R.: On the separation of net ecosystem exchange into assimilation and ecosystem respiration: review and improved algorithm, Glob. Change Biol., 11, 1424–1439, https://doi.org/10.1111/j.1365-2486.2005.001002.x, 2005.
Remaud, M., Chevallier, F., Maignan, F., Belviso, S., Berchet, A., Parouffe, A., Abadie, C., Bacour, C., Lennartz, S., and Peylin, P.: Plant gross primary production, plant respiration and carbonyl sulfide emissions over the globe inferred by atmospheric inverse modelling, Atmos. Chem. Phys., 22, 2525–2552, https://doi.org/10.5194/acp-22-2525-2022, 2022.
Remaud, M., Ma, J., Krol, M., Abadie, C., Cartwright, M. P., Patra, P., Niwa, Y., Rodenbeck, C., Belviso, S., Kooijmans, L., Lennartz, S., Maignan, F., Chevallier, F., Chipperfield, M. P., Pope, R. J., Harrison, J. J., Vimont, I., Wilson, C., and Peylin, P.: Intercomparison of Atmospheric Carbonyl Sulfide (TransCom-COS; Part One): Evaluating the Impact of Transport and Emissions on Tropospheric Variability Using Ground-Based and Aircraft Data, J. Geophys. Res.-Atmos., 128, e2022JD037817, https://doi.org/10.1029/2022JD037817, 2023.
Rogers, A.: The use and misuse of Vc,max in Earth System Models, Photosynth. Res., 119, 15–29, 2014.
Rogers, A., Medlyn, B. E., Dukes, J. S., Bonan, G., von Caemmerer, S., Dietze, M. C., Kattge, J., Leakey, A. D. B., Mercado, L. M., Niinemets, Ü., Prentice, I. C., Serbin, S. P., Sitch, S., Way, D. A., and Zaehle, S.: A roadmap for improving the representation of photosynthesis in Earth system models, New Phytol., 213, 22–42, https://doi.org/10.1111/nph.14283, 2017.
Ryu, Y., Jiang, C., Kobayashi, H., and Detto, M.: MODIS-derived global land products of shortwave radiation and diffuse and total photosynthetically active radiation at 5 km resolution from 2000, Remote Sens. Environ., 204, 812–825, 2018.
Sage, R. F. and Pearcy, R. W.: The nitrogen use efficiency of C3 and C4 plants: II. Leaf nitrogen effects on the gas exchange characteristics of Chenopodium album (L.) and Amaranthus retroflexus (L.), Plant Physiol., 84, 959–963, 1987.
Sambridge, M. and Mosegaard, K.: Monte Carlo methods in geophysical inverse problems, Rev. Geophys., 40, 3–1-3-29, 2002.
Sandoval-Soto, L., Stanimirov, M., von Hobe, M., Schmitt, V., Valdes, J., Wild, A., and Kesselmeier, J.: Global uptake of carbonyl sulfide (COS) by terrestrial vegetation: Estimates corrected by deposition velocities normalized to the uptake of carbon dioxide (CO2), Biogeosciences, 2, 125–132, https://doi.org/10.5194/bg-2-125-2005, 2005.
Sargsyan, K., Safta, C., Najm, H. N., Debusschere, B. J., Ricciuto, D., and Thornton, P.: Dimensionality reduction for complex models via Bayesian compressive sensing, Int. J. Uncertain. Quan., 4, 63–93, https://doi.org/10.1615/Int.J.UncertaintyQuantification.2013006821, 2014.
Schwalm, C. R., Williams, C. A., Schaefer, K., Anderson, R., Arain, M. A., Baker, I., Barr, A., Black, T. A., Chen, G., Chen, J. M., Ciais, P., Davis, K. J., Desai, A., Dietze, M., Dragoni, D., Fischer, M. L., Flanagan, L. B., Grant, R., Gu, L., Hollinger, D., Izaurralde, R. C., Kucharik, C., Lafleur, P., Law, B. E., Li, L., Li, Z., Liu, S., Lokupitiya, E., Luo, Y., Ma, S., Margolis, H., Matamala, R., McCaughey, H., Monson, R. K., Oechel, W. C., Peng, C., Poulter, B., Price, D. T., Riciutto, D. M., Riley, W., Sahoo, A. K., Sprintsin, M., Sun, J., Tian, H., Tonitto, C., Verbeeck, H., and Verma, S. B.: A model-data intercomparison of CO2 exchange across North America: Results from the North American Carbon Program site synthesis, J. Geophys. Res.-Biogeo., 115, G00H05, https://doi.org/10.1029/2009JG001229, 2010.
Seibt, U., Kesselmeier, J., Sandoval-Soto, L., Kuhn, U., and Berry, J. A.: A kinetic analysis of leaf uptake of COS and its relation to transpiration, photosynthesis and carbon isotope fractionation, Biogeosciences, 7, 333–341, https://doi.org/10.5194/bg-7-333-2010, 2010.
Shaw, D. C., Franklin, J. F., Bible, K., Klopatek, J., Freeman, E., Greene, S., and Parker, G. G.: Ecological setting of the Wind River old-growth forest, Ecosystems, 7, 427–439, 2004.
Smith, B., Knorr, W., Widlowski, J.-L., Pinty, B., and Gobron, N.: Combining remote sensing data with process modelling to monitor boreal conifer forest carbon balances, Forest Ecol. Manag., 255, 3985–3994, 2008.
Spielmann, F. M., Wohlfahrt, G., Hammerle, A., Kitz, F., Migliavacca, M., Alberti, G., Ibrom, A., El-Madany, T. S., Gerdel, K., Moreno, G., Kolle, O., Karl, T., Peressotti, A., and Delle Vedove, G.: Gross Primary Productivity of Four European Ecosystems Constrained by Joint CO2 and COS Flux Measurements, Geophys. Res. Lett., 46, 5284–5293, https://doi.org/10.1029/2019GL082006, 2019a.
Spielmann, F. M., Wohlfahrt, G., Hammerle, A., Kitz, F., Migliavacca, M., Alberti, G., Ibrom, A., El-Madany, T., Gerdel, K., Moreno, G., Kolle, O., Karl, T., Peressotti, A., and Delle Vedove, G.: Dataset for “Gross primary productivity of four European ecosystems constrained by joint CO2 and COS flux measurements” In Geophysical Research Letters (1.0.0), Zenodo [data set], https://doi.org/10.5281/zenodo.3406990, 2019b.
Staudt, K., Falge, E., Pyles, R. D., Paw U, K. T., and Foken, T.: Sensitivity and predictive uncertainty of the ACASA model at a spruce forest site, Biogeosciences, 7, 3685–3705, https://doi.org/10.5194/bg-7-3685-2010, 2010.
Stimler, K., Montzka, S. A., Berry, J. A., Rudich, Y., and Yakir, D.: Relationships between carbonyl sulfide (COS) and CO2 during leaf gas exchange, New Phytol., 186, 869–878, 2010.
Stimler, K., Berry, J. A., and Yakir, D.: Effects of carbonyl sulfide and carbonic anhydrase on stomatal conductance, Plant Physiol., 158, 524–530, 2012.
Sun, W., Maseyk, K., Lett, C., and Seibt, U.: A soil diffusion–reaction model for surface COS flux: COSSM v1, Geosci. Model Dev., 8, 3055–3070, https://doi.org/10.5194/gmd-8-3055-2015, 2015.
Sun, W., Kooijmans, L. M. J., Maseyk, K., Chen, H., Mammarella, I., Vesala, T., Levula, J., Keskinen, H., and Seibt, U.: Soil fluxes of carbonyl sulfide (COS), carbon monoxide, and carbon dioxide in a boreal forest in southern Finland, Atmos. Chem. Phys., 18, 1363–1378, https://doi.org/10.5194/acp-18-1363-2018, 2018.
Sun, W., Berry, J. A., Yakir, D., and Seibt, U.: Leaf relative uptake of carbonyl sulfide to CO2 seen through the lens of stomatal conductance–photosynthesis coupling, New Phytol., 235, 1729–1742, 2022.
Tang, J. and Zhuang, Q.: A global sensitivity analysis and Bayesian inference framework for improving the parameter estimation and prediction of a process-based Terrestrial Ecosystem Model, Journal of Geophys. Res.-Atmos., 114, D15303, https://doi.org/10.1029/2009JD011724, 2009.
Tonkin, M. and Doherty, J.: Calibration-constrained Monte Carlo analysis of highly parameterized models using subspace techniques, Water Resour. Res., 45, W00B10, https://doi.org/10.1029/2007WR006678, 2009.
Vesala, T., Kohonen, K.-M., Kooijmans, L. M. J., Praplan, A. P., Foltýnová, L., Kolari, P., Kulmala, M., Bäck, J., Nelson, D., Yakir, D., Zahniser, M., and Mammarella, I.: Long-term fluxes of carbonyl sulfide and their seasonality and interannual variability in a boreal forest, Atmos. Chem. Phys., 22, 2569–2584, https://doi.org/10.5194/acp-22-2569-2022, 2022.
Wang, J., Jiang, F., Wang, H., Qiu, B., Wu, M., He, W., Ju, W., Zhang, Y., Chen, J. M., and Zhou, Y.: Constraining global terrestrial gross primary productivity in a global carbon assimilation system with OCO-2 chlorophyll fluorescence data, Agr. Forest Meteorol., 304, 108424, https://doi.org/10.1016/j.agrformet.2021.108424, 2021.
Wang, S., Ibrom, A., Bauer-Gottwein, P., and Garcia, M.: Incorporating diffuse radiation into a light use efficiency and evapotranspiration model: An 11-year study in a high latitude deciduous forest, Agr. Forest Meteorol., 248, 479–493, 2018.
Wehr, R., Commane, R., Munger, J. W., McManus, J. B., Nelson, D. D., Zahniser, M. S., Saleska, S. R., and Wofsy, S. C.: Dynamics of canopy stomatal conductance, transpiration, and evaporation in a temperate deciduous forest, validated by carbonyl sulfide uptake, Biogeosciences, 14, 389–401, https://doi.org/10.5194/bg-14-389-2017, 2017.
Wharton, S.: AmeriFlux BASE US-Wrc Wind River Crane Site, AmeriFlux AMP [data set], https://doi.org/10.17190/AMF/1246114, 2016.
Whelan, M. E., Hilton, T. W., Berry, J. A., Berkelhammer, M., Desai, A. R., and Campbell, J. E.: Carbonyl sulfide exchange in soils for better estimates of ecosystem carbon uptake, Atmos. Chem. Phys., 16, 3711–3726, https://doi.org/10.5194/acp-16-3711-2016, 2016.
Whelan, M. E., Lennartz, S. T., Gimeno, T. E., Wehr, R., Wohlfahrt, G., Wang, Y., Kooijmans, L. M. J., Hilton, T. W., Belviso, S., Peylin, P., Commane, R., Sun, W., Chen, H., Kuai, L., Mammarella, I., Maseyk, K., Berkelhammer, M., Li, K.-F., Yakir, D., Zumkehr, A., Katayama, Y., Ogée, J., Spielmann, F. M., Kitz, F., Rastogi, B., Kesselmeier, J., Marshall, J., Erkkilä, K.-M., Wingate, L., Meredith, L. K., He, W., Bunk, R., Launois, T., Vesala, T., Schmidt, J. A., Fichot, C. G., Seibt, U., Saleska, S., Saltzman, E. S., Montzka, S. A., Berry, J. A., and Campbell, J. E.: Reviews and syntheses: Carbonyl sulfide as a multi-scale tracer for carbon and water cycles, Biogeosciences, 15, 3625–3657, https://doi.org/10.5194/bg-15-3625-2018, 2018.
Whelan, M. E., Shi, M., Sun, W., Vries, L. K. d., Seibt, U., and Maseyk, K.: Soil carbonyl sulfide (OCS) fluxes in terrestrial ecosystems: an empirical model, J. Geophys. Res.-Biogeo., 127, e2022JG006858, https://doi.org/10.1029/2022JG006858, 2022.
Wohlfahrt, G., Brilli, F., Hörtnagl, L., Xu, X., Bingemer, H., Hansel, A., and Loreto, F.: Carbonyl sulfide (COS) as a tracer for canopy photosynthesis, transpiration and stomatal conductance: potential and limitations, Plant Cell Environ., 35, 657–667, 2012.
Wohlfahrt, G., Hammerle, A., and Hörtnagl, L.: FLUXNET2015 AT-Neu Neustift, FLUXNET [data set], https://doi.org/10.18140/FLX/1440121, 2020.
Woodward, F. I., Smith, T. M., and Emanuel, W. R.: A global land primary productivity and phytogeography model, Global Biogeochem. Cy., 9, 471–490, 1995.
Wu, M., Ran, Y., Jansson, P.-E., Chen, P., Tan, X., and Zhang, W.: Global parameters sensitivity analysis of modeling water, energy and carbon exchange of an arid agricultural ecosystem, Agr. Forest Meteorol., 271, 295–306, 2019.
Wu, M., Tan, X., Wu, J., Huang, J., Jansson, P.-E., and Zhang, W.: Coupled water transport and heat flux in seasonally frozen soils: uncertainties identification in multi-site calibration, Environ. Earth Sci., 79, 524, https://doi.org/10.1007/s12665-020-09262-2, 2020.
Xiao, Z., Liang, S., Wang, J., Xiang, Y., Zhao, X., and Song, J.: Long-time-series global land surface satellite leaf area index product derived from MODIS and AVHRR surface reflectance, IEEE T. Geosci. Remote, 54, 5301–5318, https://doi.org/10.1109/TGRS.2016.2560522, 2016 (data available at: https://doi.org/10.12041/geodata.GLASS_LAI_MODIS(0.05D).ver1.db).
Xing, X., Wu, M., Zhang, W., Ju, W., Tagesson, T., He, W., Wang, S., Wang, J., Hu, L., Yuan, S., Zhu, T., Wang, X., Ran, Y., Li, S., Wang, C., and Jiang, F.: Modeling China's terrestrial ecosystem gross primary productivity with BEPS model: Parameter sensitivity analysis and model calibration, Agr. Forest Meteorol., 343, 109789, https://doi.org/10.1016/j.agrformet.2023.109789, 2023.
Yi, D. H., Kim, D. W., and Park, C. S.: Parameter identifiability in Bayesian inference for building energy models, Energ. Buildings, 198, 318–328, 2019.
Yuan, W., Liu, S., Zhou, G., Zhou, G., Tieszen, L. L., Baldocchi, D., Bernhofer, C., Gholz, H., Goldstein, A. H., Goulden, M. L., Hollinger, D. Y., Hu, Y., Law, B. E., Stoy, P. C., Vesala, T., and Wofsy, S. C.: Deriving a light use efficiency model from eddy covariance flux data for predicting daily gross primary production across biomes, Agr. Forest Meteorol., 143, 189–207, https://doi.org/10.1016/j.agrformet.2006.12.001, 2007.
Zaehle, S., Sitch, S., Smith, B., and Hatterman, F.: Effects of parameter uncertainties on the modeling of terrestrial biosphere dynamics, Global Biogeochem. Cy., 19, GB3020, https://doi.org/10.1029/2004GB002395, 2005.
Zhu, H., Wu, M., Jiang, F., Vossbeck, M., Kaminski, T., Xing, X., Wang, J., Ju, W., and Chen, J. M.: Assimilation of Carbonyl Sulfide (COS) fluxes within the adjoint-based data assimilation system–Nanjing University Carbon Assimilation System (NUCAS v1.0), EGUsphere [preprint], https://doi.org/10.5194/egusphere-2023-1955, 2023.
Zierl, B.: A water balance model to simulate drought in forested ecosystems and its application to the entire forested area in Switzerland, J. Hydrol., 242, 115–136, 2001.