scholarly journals Constraining ecosystem model with Adaptive Metropolis algorithm using boreal forest site eddy covariance measurements

2016 ◽  
Author(s):  
Jarmo Mäkelä ◽  
Jouni Susiluoto ◽  
Tiina Markkanen ◽  
Mika Aurela ◽  
Ivan Mammarella ◽  
...  
Keyword(s):  
Boreal Forest ◽  
Gross Primary Production ◽  
Co2 Concentration ◽  
Moisture Stress ◽  
Ecosystem Model ◽  
Metropolis Algorithm ◽  
Forest Site ◽  
Diurnal Cycles ◽  
Forest Sites ◽  
Eddy Covariance Measurements

Abstract. We examined parameter optimization in JSBACH ecosystem model, applied for two boreal forest sites in Finland. We identified and tested key parameters in soil hydrology and forest water and carbon exchange related formulations and optimized them using the Adaptive Metropolis algorithm for a five year calibration period (2000–2004) followed by a four year validation period (2005–2008). We were able to improve the modelled seasonal, daily and diurnal cycles of gross primary production and evapotranspiration but unable to enhance the models response to dryness. The improvements are mostly accounted for by parameters related to the ratio of leaf internal CO2 concentration to external CO2, relative humidity, transpiration and soil moisture stress.

scholarly journals Constraining ecosystem model with adaptive Metropolis algorithm using boreal forest site eddy covariance measurements

2016 ◽  
Vol 23 (6) ◽  
pp. 447-465 ◽  
Author(s):  
Jarmo Mäkelä ◽  
Jouni Susiluoto ◽  
Tiina Markkanen ◽  
Mika Aurela ◽  
Heikki Järvinen ◽  
...  
Keyword(s):  
Relative Humidity ◽  
Boreal Forest ◽  
Gross Primary Production ◽  
Ecosystem Model ◽  
Soil Parameters ◽  
Large Contribution ◽  
Forest Site ◽  
Validation Period ◽  
Area Index ◽  
Soil Dryness

Abstract. We examined parameter optimisation in the JSBACH (Kaminski et al., 2013; Knorr and Kattge, 2005; Reick et al., 2013) ecosystem model, applied to two boreal forest sites (Hyytiälä and Sodankylä) in Finland. We identified and tested key parameters in soil hydrology and forest water and carbon-exchange-related formulations, and optimised them using the adaptive Metropolis (AM) algorithm for Hyytiälä with a 5-year calibration period (2000–2004) followed by a 4-year validation period (2005–2008). Sodankylä acted as an independent validation site, where optimisations were not made. The tuning provided estimates for full distribution of possible parameters, along with information about correlation, sensitivity and identifiability. Some parameters were correlated with each other due to a phenomenological connection between carbon uptake and water stress or other connections due to the set-up of the model formulations. The latter holds especially for vegetation phenology parameters. The least identifiable parameters include phenology parameters, parameters connecting relative humidity and soil dryness, and the field capacity of the skin reservoir. These soil parameters were masked by the large contribution from vegetation transpiration. In addition to leaf area index and the maximum carboxylation rate, the most effective parameters adjusting the gross primary production (GPP) and evapotranspiration (ET) fluxes in seasonal tuning were related to soil wilting point, drainage and moisture stress imposed on vegetation. For daily and half-hourly tunings the most important parameters were the ratio of leaf internal CO2 concentration to external CO2 and the parameter connecting relative humidity and soil dryness. Effectively the seasonal tuning transferred water from soil moisture into ET, and daily and half-hourly tunings reversed this process. The seasonal tuning improved the month-to-month development of GPP and ET, and produced the most stable estimates of water use efficiency. When compared to the seasonal tuning, the daily tuning is worse on the seasonal scale. However, daily parametrisation reproduced the observations for average diurnal cycle best, except for the GPP for Sodankylä validation period, where half-hourly tuned parameters were better. In general, the daily tuning provided the largest reduction in model–data mismatch. The models response to drought was unaffected by our parametrisations and further studies are needed into enhancing the dry response in JSBACH.

scholarly journals Validation of 3D-CMCC Forest Ecosystem Model (v.5.1) against eddy covariance data for ten European forest sites

2015 ◽  
Vol 8 (8) ◽  
pp. 6867-6929
Author(s):  
A. Collalti ◽  
S. Marconi ◽  
A. Ibrom ◽  
C. Trotta ◽  
A. Anav ◽  
...  
Keyword(s):  
Eddy Covariance ◽  
Forest Ecosystem ◽  
Carbon Allocation ◽  
Structural Characteristics ◽  
Gross Primary Production ◽  
Ecosystem Model ◽  
Summer Drought ◽  
Drought Event ◽  
Forest Sites ◽  
Species Specific

Abstract. This study evaluates the performances of the new version (v.5.1) of 3D-CMCC Forest Ecosystem Model (FEM) in simulating gross primary production (GPP), against eddy covariance GPP data for ten FLUXNET forest sites across Europe. A new carbon allocation module, coupled with new both phenological and autotrophic respiration schemes, was implemented in this new version. Model ability in reproducing timing and magnitude of daily and monthly GPP fluctuations is validated at intra-annual and inter-annual scale, including extreme anomalous seasons. With the purpose to test the 3D-CMCC FEM applicability over Europe without a site-related calibration, the model has been deliberately parameterized with a single set of species-specific parameterizations for each forest ecosystem. The model consistently reproduces both in timing and in magnitude daily and monthly GPP variability across all sites, with the exception of the two Mediterranean sites. We find that 3D-CMCC FEM tends to better simulate the timing of inter-annual anomalies than their magnitude within measurements uncertainty. In six of eight sites where data were available the model well reproduces the 2003 summer drought event. Finally, for three sites we evaluate if a more accurate representation of forest structural characteristics (i.e. cohorts, forest layers) and species composition can improve model results. In two of the three sites results reveal that model slightly increases its performances, although, statistically speaking, not in a relevant way.

scholarly journals Sensitivity of 21st century simulated ecosystem indicators to model parameters, prescribed climate drivers, RCP scenarios and forest management actions for two Finnish boreal forest sites

2020 ◽  
Vol 17 (10) ◽  
pp. 2681-2700
Author(s):  
Jarmo Mäkelä ◽  
Francesco Minunno ◽  
Tuula Aalto ◽  
Annikki Mäkelä ◽  
Tiina Markkanen ◽  
...  
Keyword(s):  
Forest Management ◽  
Boreal Forest ◽  
Parameter Uncertainty ◽  
Ecosystem Model ◽  
Model Parameters ◽  
Rcp Scenarios ◽  
Management Actions ◽  
Uncertainty Sources ◽  
Forest Sites ◽  
Ecosystem Indicators

Abstract. Forest ecosystems are already responding to changing environmental conditions that are driven by increased atmospheric CO2 concentrations. These developments affect how societies can utilise and benefit from the woodland areas in the future, be it for example climate change mitigation as carbon sinks, lumber for wood industry, or preserved for nature tourism and recreational activities. We assess the effect and the relative magnitude of different uncertainty sources in ecosystem model simulations from the year 1980 to 2100 for two Finnish boreal forest sites. The models used in this study are the land ecosystem model JSBACH and the forest growth model PREBAS. The considered uncertainty sources for both models are model parameters and four prescribed climates with two RCP (representative concentration pathway) scenarios. Usually, model parameter uncertainty is not included in these types of uncertainty studies. PREBAS simulations also include two forest management scenarios. We assess the effect of these sources of variation at four different points in time on several ecosystem indicators, e.g. gross primary production (GPP), ecosystem respiration, soil moisture, recurrence of drought, length of the vegetation active period (VAP), length of the snow melting period and the stand volume. The uncertainty induced by the climate models remains roughly the same throughout the simulations and is overtaken by the RCP scenario impact halfway through the experiment. The management actions are the most dominant uncertainty factors for Hyytiälä and as important as RCP scenarios at the end of the simulations, but they contribute only half as much for Sodankylä. The parameter uncertainty is the least influential of the examined uncertainty sources, but it is also the most elusive to estimate due to non-linear and adverse effects on the simulated ecosystem indicators. Our analysis underlines the importance of carefully considering the implementation of forest use when simulating future ecosystem conditions, as human impact is evident and even increasing in boreal forested regions.

scholarly journals Towards a more harmonized processing of eddy covariance CO<sub>2</sub> fluxes: algorithms and uncertainty estimation

2006 ◽  
Vol 3 (4) ◽  
pp. 961-992 ◽  
Author(s):  
D. Papale ◽  
M. Reichstein ◽  
E. Canfora ◽  
M. Aubinet ◽  
C. Bernhofer ◽  
...  
Keyword(s):  
Eddy Covariance ◽  
Ecosystem Respiration ◽  
Gross Primary Production ◽  
Terrestrial Ecosystem ◽  
Uncertainty Estimation ◽  
Climatic Conditions ◽  
High Temporal Resolution ◽  
Forest Sites ◽  
European Database ◽  
Eddy Covariance Measurements

Abstract. Eddy covariance technique to measure CO2, water and energy fluxes between biosphere and atmosphere is widely spread and used in various regional networks. Currently more that 250 eddy covariance sites are active around the world measuring carbon exchange at high temporal resolution for different biomes and climatic conditions. These data are usually acquired using the same method but they need a set of corrections that are often differently applied to each site and in a subjective way. In this paper a new standardized set of corrections are proposed and the uncertainties introduced by these corrections are assessed for 8 different forest sites in Europe with a total of 12 yearly datasets. The uncertainties introduced on the two components GPP (Gross Primary Production) and TER (Terrestrial Ecosystem Respiration) are also discussed and a quantitative analysis presented . The results show that a standardized data processing is needed for an effective comparison across biomes and for underpinning inter-annual variability. The methodology presented in this paper has also been integrated in the European database of the eddy covariance measurements.

Ecosystem modeling of methane and carbon dioxide fluxes for boreal forest sites

2001 ◽  
Vol 31 (2) ◽  
pp. 208-223 ◽  
Author(s):  
Christopher Potter ◽  
Jill Bubier ◽  
Patrick Crill ◽  
Peter Lafleur
Keyword(s):  
Carbon Dioxide ◽  
Boreal Forest ◽  
Land Surface ◽  
Net Primary Production ◽  
Ground Cover ◽  
Methane Flux ◽  
Growing Season ◽  
Ecosystem Model ◽  
Heat Fluxes ◽  
Water Levels

Predicted daily fluxes from an ecosystem model for water, carbon dioxide, and methane were compared with 1994 and 1996 Boreal Ecosystem–Atmosphere Study (BOREAS) field measurements at sites dominated by old black spruce (Picea mariana (Mill.) BSP) (OBS) and boreal fen vegetation near Thompson, Man. Model settings for simulating daily changes in water table depth (WTD) for both sites were designed to match observed water levels, including predictions for two microtopographic positions (hollow and hummock) within the fen study area. Water run-on to the soil profile from neighboring microtopographic units was calibrated on the basis of daily snowmelt and rainfall inputs to reproduce BOREAS site measurements for timing and magnitude of maximum daily WTD for the growing season. Model predictions for daily evapotranspiration rates closely track measured fluxes for stand water loss in patterns consistent with strong controls over latent heat fluxes by soil temperature during nongrowing season months and by variability in relative humidity and air temperature during the growing season. Predicted annual net primary production (NPP) for the OBS site was 158 g C·m–2 during 1994 and 135 g C·m–2 during 1996, with contributions of 75% from overstory canopy production and 25% from ground cover production. Annual NPP for the wetter fen site was 250 g C·m–2 during 1994 and 270 g C·m–2 during 1996. Predicted seasonal patterns for soil CO2 fluxes and net ecosystem production of carbon both match daily average estimates at the two sites. Model results for methane flux, which also closely match average measured flux levels of –0.5 mg CH4·m–2·day–1 for OBS and 2.8 mg CH4·m–2·day–1 for fen sites, suggest that spruce areas are net annual sinks of about –0.12 g CH4·m–2, whereas fen areas generate net annual emissions on the order of 0.3–0.85 g CH4·m–2, depending mainly on seasonal WTD and microtopographic position. Fen hollow areas are predicted to emit almost three times more methane during a given year than fen hummock areas. The validated model is structured for extrapolation to regional simulations of interannual trace gas fluxes over the entire North America boreal forest, with integration of satellite data to characterize properties of the land surface.

scholarly journals The composition and diversity of natural regeneration of tree species in gaps under different intensities of forest disturbance

2021 ◽  
Author(s):  
Maame Esi Hammond ◽  
Radek Pokorný ◽  
Daniel Okae-Anti ◽  
Augustine Gyedu ◽  
Irene Otwuwa Obeng
Keyword(s):  
Tree Species ◽  
Natural Regeneration ◽  
Forest Disturbance ◽  
Forest Site ◽  
Gap Formation ◽  
Tree Species Composition ◽  
Forest Sites ◽  
Significant Difference ◽  
Shade Tolerant Species ◽  
Disturbed Forest

AbstractThe positive ecological interaction between gap formation and natural regeneration has been examined but little research has been carried out on the effects of gaps on natural regeneration in forests under different intensities of disturbance. This study evaluates the composition, diversity, regeneration density and abundance of natural regeneration of tree species in gaps in undisturbed, intermittently disturbed, and disturbed forest sites. Bia Tano Forest Reserve in Ghana was the study area and three gaps each were selected in the three forest site categories. Ten circular subsampling areas of 1 m2 were delineated at 2 m spacing along north, south, east, and west transects within individual gaps. Data on natural regeneration < 350 cm height were gathered. The results show that the intensity of disturbance was disproportional to gap size. Species diversity differed significantly between undisturbed and disturbed sites and, also between intermittently disturbed and disturbed sites for Simpson’s (1-D), Equitability (J), and Berger–Parker (B–P) indices. However, there was no significant difference among forest sites for Shannon diversity (H) and Margalef richness (MI) indices. Tree species composition on the sites differed. Regeneration density on the disturbed site was significantly higher than on the two other sites. Greater abundance and density of shade-dependent species on all sites identified them as opportunistic replacements of gap-dependent pioneers. Pioneer species giving way to shade tolerant species is a natural process, thus make them worst variant in gap regeneration.

scholarly journals Land surface model photosynthesis and parameter calibration for boreal sites with adaptive population importance sampler

2019 ◽  
Author(s):  
Jarmo Mäkelä ◽  
Jürgen Knauer ◽  
Mika Aurela ◽  
Andrew Black ◽  
Martin Heimann ◽  
...  
Keyword(s):  
Stomatal Conductance ◽  
Land Surface ◽  
Carbon Fixation ◽  
Gross Primary Production ◽  
Effect Of Temperature ◽  
Coefficient Of Determination ◽  
Surface Model ◽  
Drought Event ◽  
Forest Site ◽  
Delayed Effect

Abstract. We calibrated the JSBACH model with six different stomatal conductance formulations using measurements from 10 FLUXNET coniferous evergreen sites in the Boreal zone. The parameter posterior distributions were generated by adaptive population importance sampler and the optimal values by a simple stochastic optimisation algorithm. The observations used to constrain the model are evapotranspiration (ET) and gross primary production (GPP). We identified the key parameters in the calibration process. These parameters control the soil moisture stress function and the overall rate of carbon fixation. We were able to improve the coefficient of determination and the model bias with all stomatal conductance formulations. There was no clear candidate for the best stomatal conductance model, although certain versions produced better estimates depending on the examined variable (ET, GPP) and the used metric. We were also able to significantly enhance the model behaviour during a drought event in a Finnish Scots pine forest site. The JSBACH model was also modified to use a delayed effect of temperature for photosynthetic activity. This modification enabled the model to correctly time and replicate the springtime increase in GPP (and ET) for conifers throughout the measurements sites used in this study.

scholarly journals Tiling soil textures for terrestrial ecosystem modelling via clustering analysis: a case study with CLASS-CTEM (version 2.1)

2017 ◽  
Author(s):  
Joe R. Melton ◽  
Reinel Sospedra-Alfonso ◽  
Kelly E. McCusker
Keyword(s):  
Soil Texture ◽  
Land Surface ◽  
Clustering Analysis ◽  
Clustering Algorithms ◽  
Terrestrial Ecosystem ◽  
Grid Cell ◽  
Moisture Stress ◽  
Ecosystem Model ◽  
Global Scale ◽  
Terrestrial Ecosystem Model

Abstract. We investigate the application of clustering algorithms to represent sub-grid scale variability in soil texture for use in a global-scale terrestrial ecosystem model. Our model, the coupled Canadian Land Surface Scheme – Canadian Terrestrial Ecosystem Model (CLASS-CTEM), is typically implemented at a coarse spatial resolution (ca. 2.8° × 2.8°) due to its use as the land surface component of the Canadian Earth System Model (CanESM). CLASS-CTEM can, however, be run with tiling of the land surface as a means to represent sub-grid heterogeneity. We first determined that the model was sensitive to tiling of the soil textures via an idealized test case before attempting to cluster soil textures globally. To cluster a high-resolution soil texture dataset onto our coarse model grid, we use two linked algorithms (OPTICS (Ankerst et al., 1999; Daszykowski et al., 2002) and Sander et al. (2003)) to provide tiles of representative soil textures for use as CLASS-CTEM inputs. The clustering process results in, on average, about three tiles per CLASS-CTEM grid cell with most cells having four or less tiles. Results from CLASS-CTEM simulations conducted with the tiled inputs (Cluster) versus those using a simple grid-mean soil texture (Gridmean) show CLASS-CTEM, at least on a global scale, is relatively insensitive to the tiled soil textures, however differences can be large in arid or peatland regions. The Cluster simulation has generally lower soil moisture and lower overall vegetation productivity than the Gridmean simulation except in arid regions where plant productivity increases. In these dry regions, the influence of the tiling is stronger due to the general state of vegetation moisture stress which allows a single tile, whose soil texture retains more plant available water, to yield much higher productivity. Although the use of clustering analysis appears promising as a means to represent sub-grid heterogeneity, soil textures appear to be reasonably represented for global scale simulations using a simple grid-mean value.

Sign in / Sign up

Export Citation Format

Share Document