We analyze how biases of meteorological drivers impact the calculation of ecosystem CO2, water and energy fluxes by models. To do so, we drive the same ecosystem model by meteorology from gridded ...products and by meteorology from local observation at eddy-covariance flux sites. The study is focused on six flux tower sites in France spanning across a climate gradient of 7–14 °C annual mean surface air temperature and 600–1040 mm mean annual rainfall, with forest, grassland and cropland ecosystems. We evaluate the results of the ORCHIDEE process-based model driven by meteorology from four different analysis data sets against the same model driven by site-observed meteorology. The evaluation is decomposed into characteristic time scales. The main result is that there are significant differences in meteorology between analysis data sets and local observation. The phase of seasonal cycle of air temperature, humidity and shortwave downward radiation is reproduced correctly by all meteorological models (average R2 = 0.90). At sites located in altitude, the misfit of meteorological drivers from analysis data sets and tower meteorology is the largest. We show that day-to-day variations in weather are not completely well reproduced by meteorological models, with R2 between analysis data sets and measured local meteorology going from 0.35 to 0.70. The bias of meteorological driver impacts the flux simulation by ORCHIDEE, and thus would have an effect on regional and global budgets. The forcing error, defined by the simulated flux difference resulting from prescribing modeled instead of observed local meteorology drivers to ORCHIDEE, is quantified for the six studied sites at different time scales. The magnitude of this forcing error is compared to that of the model error defined as the modeled-minus-observed flux, thus containing uncertain parameterizations, parameter values, and initialization. The forcing error is on average smaller than but still comparable to model error, with the ratio of forcing error to model error being the largest on daily time scale (86%) and annual time scales (80%). The forcing error incurred from using a gridded meteorological data set to drive vegetation models is therefore an important component of the uncertainty budget of regional CO2, water and energy fluxes simulations, and should be taken into consideration in up-scaling studies.
Phenology and associated canopy development exert a strong control over seasonal energy and mass exchanges between the earth's surface and the atmosphere. Satellite measurements are used to assess ...main phenological stages of the vegetation at the global scale. The authors propose a method to derive the start, the maximum, the end, and the length of the vegetation cycle, based on the analysis of temporal series of weekly vegetation index, at a resolution of l degree lat X l degree long for year 1986. Global maps of these characteristics of the vegetation are presented, and their zonal distribution is discussed. The start of the vegetation cycle has been related to temperature sums in the case of temperate deciduous forest and to precipitation in the case of savannahs. It is concluded that satellite measurements offer interesting perspectives for global-scale quantitative phenology modeling
Vegetation reconstructions from pollen data for the Last Glacial Maximum (LGM), 21 ky ago, reveal lanscapes radically different from the modern ones, with, in particular, a massive regression of ...forested areas in both hemispheres. Two main factors have to be taken into account to explain these changes in comparison to today's potential vegetation: a generally cooler and drier climate and a lower level of atmospheric CO2 . In order to assess the relative impact of climate and atmospheric CO2 changes on the global vegetation, we simulate the potential modern vegetation and the glacial vegetation with the dynamical global vegetation model ORCHIDEE, driven by outputs from the IPSL_CM4_v1 atmosphere-ocean general circulation model, under modern or glacial CO2 levels for photosynthesis. ORCHIDEE correctly reproduces the broad features of the glacial vegetation. Our modelling results support the view that the physiological effect of glacial CO2 is a key factor to explain vegetation changes during glacial times. In our simulations, the low atmospheric CO2 is the only driver of the tropical forests regression, and explains half of the response of temperate and boreal forests to glacial conditions. Our study shows that the sensitivity to CO2 changes depends on the background climate over a region, and also depends on the vegetation type, needleleaf trees being much more sensitive than broadleaf trees in our model. This difference of sensitivity leads to a dominance of broadleaf types in the remaining simulated forests, which is not supported by pollen data, but nonetheless suggests a potential impact of CO2 on the glacial vegetation assemblages. It also modifies the competitivity between the trees and makes the amplitude of the response to CO2 dependent on the initial vegetation state.
We modeled the effects of climate change and two forest management scenarios on wood production and forest carbon balance in French forests using process-based models of forest growth. We combined ...data from the national forest inventory and soil network survey, which were aggregated over a 50 x 50-km grid, i.e., the spatial resolution of the climate scenario data. We predicted and analyzed the climate impact on potential forest production over the period 1960-2100. All models predicted a slight increase in potential forest yield until 2030-2050, followed by a plateau or a decline around 2070-2100, with overall, a greater increase in yield in northern France than in the south. Gross and net primary productivities were more negatively affected by soil water and atmospheric water vapor saturation deficits in western France because of a more pronounced shift in seasonal rainfall from summer to winter. The rotation-averaged values of carbon flux and production for different forest management options were estimated during four years (1980, 2015, 2045 and 2080). Predictions were made using a two-dimensional matrix covering the range of local soil and climate conditions. The changes in ecosystem fluxes and forest production were explained by the counterbalancing effect of rising CO2 concentration and increasing water deficit. The effect of climate change decreased with rotation length from short rotations with high production rates and low standing biomasses to long rotations with low productivities and greater standing biomasses. Climate effects on productivity, both negative and positive, were greatest on high fertility sites. Forest productivity in northern France was enhanced by climate change, increasingly from west to east, whereas in the southwestern Atlantic region, productivity was reduced by climate change to an increasing degree from west to east.
The construction of a new forest management module (FMM) within the ORCHIDEE global vegetation model (GVM) allows a realistic simulation of biomass changes during the life cycle of a forest, which ...makes many biomass datasets suitable as validation data for the coupled ORCHIDEE-FM GVM. This study uses three datasets to validate ORCHIDEE-FM at different temporal and spatial scales: permanent monitoring plots, yield tables, and the French national inventory data. The last dataset has sufficient geospatial coverage to allow a novel type of validation: inventory plots can be used to produce continuous maps that can be compared to continuous simulations for regional trends in standing volumes and volume increments. ORCHIDEE-FM performs better than simple statistical models for stand-level variables, which include tree density, basal area, standing volume, average circumference and height, when management intensity and initial conditions are known: model efficiency is improved by an average of 0.11, and its average bias does not exceed 25%. The performance of the model is less satisfying for tree-level variables, including extreme circumferences, tree circumference distribution and competition indices, or when management and initial conditions are unknown. At the regional level, when climate forcing is accurate for precipitation, ORCHIDEE-FM is able to reproduce most productivity patterns in France, such as the local lows of needleleaves in the Parisian basin and of broadleaves in south-central France. The simulation of water stress effects on biomass in the Mediterranean region, however, remains problematic, as does the simulation of the wood increment for coniferous trees. These pitfalls pertain to the general ORCHIDEE model rather than to the FMM. Overall, with an average bias seldom exceeding 40%, the performance of ORCHIDEE-FM is deemed reliable to use it as a new modelling tool in the study of the effects of interactions between forest management and climate on biomass stocks of forests across a range of scales from plot to country.
Agro-land surface models (agro-LSM) have been developed from the integration of specific crop processes into large-scale generic land surface models that allow calculating the spatial distribution ...and variability of energy, water and carbon fluxes within the soil–vegetation–atmosphere continuum. When developing agro-LSM models, particular attention must be given to the effects of crop phenology and management on the turbulent fluxes exchanged with the atmosphere, and the underlying water and carbon pools. A part of the uncertainty of agro-LSM models is related to their usually large number of parameters. In this study, we quantify the parameter-values uncertainty in the simulation of sugarcane biomass production with the agro-LSM ORCHIDEE–STICS, using a multi-regional approach with data from sites in Australia, La Réunion and Brazil. In ORCHIDEE–STICS, two models are chained: STICS, an agronomy model that calculates phenology and management, and ORCHIDEE, a land surface model that calculates biomass and other ecosystem variables forced by STICS phenology. First, the parameters that dominate the uncertainty of simulated biomass at harvest date are determined through a screening of 67 different parameters of both STICS and ORCHIDEE on a multi-site basis. Secondly, the uncertainty of harvested biomass attributable to those most sensitive parameters is quantified and specifically attributed to either STICS (phenology, management) or to ORCHIDEE (other ecosystem variables including biomass) through distinct Monte Carlo runs. The uncertainty on parameter values is constrained using observations by calibrating the model independently at seven sites. In a third step, a sensitivity analysis is carried out by varying the most sensitive parameters to investigate their effects at continental scale. A Monte Carlo sampling method associated with the calculation of partial ranked correlation coefficients is used to quantify the sensitivity of harvested biomass to input parameters on a continental scale across the large regions of intensive sugarcane cultivation in Australia and Brazil. The ten parameters driving most of the uncertainty in the ORCHIDEE–STICS modeled biomass at the 7 sites are identified by the screening procedure. We found that the 10 most sensitive parameters control phenology (maximum rate of increase of LAI) and root uptake of water and nitrogen (root profile and root growth rate, nitrogen stress threshold) in STICS, and photosynthesis (optimal temperature of photosynthesis, optimal carboxylation rate), radiation interception (extinction coefficient), and transpiration and respiration (stomatal conductance, growth and maintenance respiration coefficients) in ORCHIDEE. We find that the optimal carboxylation rate and photosynthesis temperature parameters contribute most to the uncertainty in harvested biomass simulations at site scale. The spatial variation of the ranked correlation between input parameters and modeled biomass at harvest is well explained by rain and temperature drivers, suggesting different climate-mediated sensitivities of modeled sugarcane yield to the model parameters, for Australia and Brazil. This study reveals the spatial and temporal patterns of uncertainty variability for a highly parameterized agro-LSM and calls for more systematic uncertainty analyses of such models.
Three terrestrial biosphere models (LPJ, Orchidee, Biome-BGC) were evaluated with respect to their ability to simulate large-scale climate related trends in gross primary production (GPP) across ...European forests. Simulated GPP and leaf area index (LAI) were compared with GPP estimates based on flux separated eddy covariance measurements of net ecosystem exchange and LAI measurements along a temperature gradient ranging from the boreal to the Mediterranean region. The three models capture qualitatively the pattern suggested by the site data: an increase in GPP from boreal to temperate and a subsequent decline from temperate to Mediterranean climates. The models consistently predict higher GPP for boreal and lower GPP for Mediterranean forests. Based on a decomposition of GPP into absorbed photosynthetic active radiation (APAR) and radiation use efficiency (RUE), the overestimation of GPP for the boreal coniferous forests appears to be primarily related to too high simulated LAI - and thus light absorption (APAR) – rather than too high radiation use efficiency. We cannot attribute the tendency of the models to underestimate GPP in the water limited region to model structural deficiencies with confidence. A likely dry bias of the input meteorological data in southern Europe may create this pattern. On average, the models compare similarly well to the site GPP data (RMSE of ~30% or 420 gC/m2/yr) but differences are apparent for different ecosystem types. In terms of absolute values, we find the agreement between site based GPP estimates and simulations acceptable when we consider uncertainties about the accuracy in model drivers, a potential representation bias of the eddy covariance sites, and uncertainties related to the method of deriving GPP from eddy covariance measurements data. Continental to global data-model comparison studies should be fostered in the future since they are necessary to identify consistent model bias along environmental gradients.
Stand-replacing fires are the dominant fire type in North American boreal forests. They leave a historical legacy of a mosaic landscape of different aged forest cohorts. This forest age dynamics must ...be included in vegetation models to accurately quantify the role of fire in the historical and current regional forest carbon balance. The present study adapted the global process-based vegetation model ORCHIDEE to simulate the CO2 emissions from boreal forest fire and the subsequent recovery after a stand-replacing fire; the model represents postfire new cohort establishment, forest stand structure and the self-thinning process. Simulation results are evaluated against observations of three clusters of postfire forest chronosequences in Canada and Alaska. The variables evaluated include: fire carbon emissions, CO2 fluxes (gross primary production, total ecosystem respiration and net ecosystem exchange), leaf area index, and biometric measurements (aboveground biomass carbon, forest floor carbon, woody debris carbon, stand individual density, stand basal area, and mean diameter at breast height). When forced by local climate and the atmospheric CO2 history at each chronosequence site, the model simulations generally match the observed CO2 fluxes and carbon stock data well, with model-measurement mean square root of deviation comparable with the measurement accuracy (for CO2 flux ~100 g C m−2 yr−1, for biomass carbon ~1000 g C m−2 and for soil carbon ~2000 g C m−2). We find that the current postfire forest carbon sink at the evaluation sites, as observed by chronosequence methods, is mainly due to a combination of historical CO2 increase and forest succession. Climate change and variability during this period offsets some of these expected carbon gains. The negative impacts of climate were a likely consequence of increasing water stress caused by significant temperature increases that were not matched by concurrent increases in precipitation. Our simulation results demonstrate that a global vegetation model such as ORCHIDEE is able to capture the essential ecosystem processes in fire-disturbed boreal forests and produces satisfactory results in terms of both carbon fluxes and carbon-stock evolution after fire. This makes the model suitable for regional simulations in boreal regions where fire regimes play a key role in the ecosystem carbon balance.
Aiming at producing improved estimates of carbon source/sink spatial and interannual patterns across Europe (35% croplands), this work uses the ORCHIDEE‐STICS terrestrial biosphere model including a ...more realistic representation of croplands, described in part 1 (Smith et al., 2010). Crop yield is derived from annual Net Primary Productivity and compared with wheat and grain maize harvest data for five European countries. Over a 34 year period, the best correlation coefficient obtained between observed and simulated yield time series is for irrigated maize in Italy (R = 0.73). In the data as well as in the model, 1976 and 2003 appear as climate anomalies causing a ≈40% yield drop in the most affected regions. Simulated interannual yield anomalies and the spatial pattern of the yield drop in 2003 are found to be more realistic than the results from ORCHIDEE with no representation of croplands. The simulated 2003 anomalous carbon source from European ecosystems to the atmosphere due to the 2003 summer heat wave is in good agreement with atmospheric inversions (0.20GtC, from May to October). The anomaly is twice too large in the ORCHIDEE alone simulation, owing to the unrealistically high exposure of herbaceous plants to the extreme summer conditions. The mechanisms linking abnormally high summer temperatures, the crop productivity drop, and significant carbon source from European ecosystems in 2003 are discussed. Overall, this study highlights the importance of accounting for the specific phenologies of crops sown both in winter and in spring and for irrigation applied to summer crops in regional/global models of the terrestrial carbon cycle.
We are comparing spatially explicit process-model based estimates of the terrestrial carbon balance and its components over Africa and confront them with remote sensing based proxies of vegetation ...productivity and atmospheric inversions of land-atmosphere net carbon exchange. Particular emphasis is on characterizing the patterns of interannual variability of carbon fluxes and analyzing the factors and processes responsible for it. For this purpose simulations with the terrestrial biosphere models ORCHIDEE, LPJ-DGVM, LPJ-Guess and JULES have been performed using a standardized modeling protocol and a uniform set of corrected climate forcing data. While the models differ concerning the absolute magnitude of carbon fluxes, we find several robust patterns of interannual variability among the models. Models exhibit largest interannual variability in southern and eastern Africa, regions which are primarily covered by herbaceous vegetation. Interannual variability of the net carbon balance appears to be more strongly influenced by gross primary production than by ecosystem respiration. A principal component analysis indicates that moisture is the main driving factor of interannual gross primary production variability for those regions. On the contrary in a large part of the inner tropics radiation appears to be limiting in two models. These patterns are partly corroborated by remotely sensed vegetation properties from the SeaWiFS satellite sensor. Inverse atmospheric modeling estimates of surface carbon fluxes are less conclusive at this point, implying the need for a denser network of observation stations over Africa.