ALBERT

Description: The profitability of hydropower in Costa Rica is affected by soil erosion and sedimentation in dam reservoirs, which are in turn influenced by land use, infiltration and aquifer interactions with surface water. In order to foster the provision and payment for Hydrological Environmental Services (HES), a quantitative assessment of the impact of specific land uses on the functioning of drainage-basins is required. The present paper aims to study the water balance partitioning in a volcanic coffee agroforestry micro-basin (1 km2, steep slopes) in Costa Rica, as a first step towards evaluating sediment or contaminant loads. The main hydrological processes were monitored during one year, using flume, eddy-covariance flux tower, soil water profiles and piezometers. A new Hydro-SVAT lumped model is proposed, that balances SVAT (Soil Vegetation Atmosphere Transfer) and basin-reservoir routines. The purpose of such a coupling was to achieve a trade-off between the expected performance of ecophysiological and hydrological models, which are often employed separately and at different spatial scales, either the plot or the basin. The calibration of the model to perform streamflow yielded a Nash-Sutcliffe (NS) coefficient equal to 0.89 for the year 2009, while the validation of the water balance partitioning was consistent with the independent measurements of actual evapotranspiration (R2 = 0.79, energy balance closed independently), soil water content (R2 = 0.35) and water table level (R2 = 0.84). Eight months of data from 2010 were used to validate modelled streamflow, resulting in a NS = 0.75. An uncertainty analysis showed that the streamflow modelling was precise for nearly every time step, while a sensitivity analysis revealed which parameters mostly affected model precision, depending on the season. It was observed that 64% of the incident rainfall R flowed out of the basin as streamflow and 25% as evapotranspiration, while the remaining 11% is probably explained by deep percolation, measurement errors and/or inter-annual changes in soil and aquifer water stocks. The model indicated an interception loss equal to 4% of R, a surface runoff of 4% and an infiltration component of 92%. The modelled streamflow was constituted by 87% of baseflow originating from the aquifer, 7% of subsurface non-saturated runoff and 6% of surface runoff. Given the low surface runoff observed under the current physical conditions (andisol) and management practices (no tillage, planted trees, bare soil kept by weeding), this agroforestry system on a volcanic soil demonstrated potential to provide valuable HES, such as a reduced superficial displacement-capacity for fertilizers, pesticides and sediments, as well as a streamflow regulation function provided by the highly efficient mechanisms of aquifer recharge and discharge. The proposed combination of experimentation and modelling across ecophysiological and hydrological approaches proved to be useful to account for the behaviour of a given basin, so that it can be applied to compare HES provision for different regions or management alternatives.

Description: The profitability of hydropower in Costa Rica is affected by soil erosion and sedimentation in dam reservoirs, which are in turn influenced by land use, infiltration and aquifer interactions with surface water. In order to foster the provision and payment of Hydrological Environmental Services (HES), a quantitative assessment of the impact of specific land uses on the functioning of drainage-basins is required. The present paper aims to study the water balance partitioning in a volcanic coffee agroforestry micro-basin (1 km2, steep slopes) in Costa Rica, as a first step towards evaluating sediment or contaminant loads. The main hydrological processes were monitored during one year, using flume, eddy-covariance flux tower, soil water profiles and piezometers. A new Hydro-SVAT lumped model is proposed, that balances SVAT (Soil Vegetation Atmosphere Transfer) and basin-reservoir routines. The purpose of such a coupling was to achieve a trade-off between the expected performance of ecophysiological and hydrological models, which are often employed separately and at different spatial scales, either the plot or the basin. The calibration of the model to perform streamflow yielded a NS coefficient equal to 0.80, while the validation of the water balance partitioning was consistent with the independent measurements of actual evapotranspiration (R2=0.79, energy balance closed independently), soil water content (R2=0.49) and water table level (R2=0.90). An uncertainty analysis showed that the streamflow modelling was precise for nearly every time step, while a sensitivity analysis revealed which parameters mostly affected model precision, depending on the season. It was observed that 64% of the incident rainfall R flowed out of the basin as streamflow, 25% as evapotranspiration and the remaining 11% was attributed to deep percolation. The model indicated an interception loss equal to 4% of R, a surface runoff of 5% and an infiltration component of 91%. The modelled streamflow was constituted by 63% of baseflow originating from the aquifer, 29% of subsurface non-saturated runoff and 8% of surface runoff. Given the low surface runoff observed under the current physical conditions (andisol) and management practices (no tillage, planted trees, bare soil kept by weeding), this agroforestry system on a volcanic soil demonstrated potential to provide valuable HES, such as a reduced superficial displacement-capacity for fertilizers, pesticides and sediments, as well as a streamflow regulation function provided by the highly efficient mechanisms of aquifer recharge and discharge. The proposed combination of experimentation and modelling across ecophysiological and hydrological approaches proved to be useful to account for the behaviour of a given basin, so that it can be applied to compare HES provision for different regions or management alternatives.

Description: We evaluate three terrestrial biosphere models (LPJ, Orchidee, Biome-BGC) with respect to their capacity to simulate climate related trends in gross primary production (GPP) of forests in Europe. We compare simulated GPP and leaf area index (LAI) with GPP estimates based on flux separated eddy covariance measurements of net ecosystem exchange (NEE) and LAI measurements along a gradient in mean annual temperature from the boreal to the Mediterranean.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 zone appears to be primarily related to too high simulated LAI - and thus light absorption (APAR) – rather than too high radiation use efficiency. 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. Given 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, we find the agreement between site data and simulations acceptable, providing confidence in simulations of GPP for European forests.

Description: 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.

Description: Land surface modelling has been widely used to characterize the two-way interactions between climate and human activities in terrestrial ecosystems such as deforestation, agricultural expansion, and urbanization. Towards an effort to quantify the effects of forests to oil palm conversion occurring in the tropics on land–atmosphere carbon, water and energy fluxes, we introduce a new perennial crop plant functional type (PFT) for oil palm. Due to the modular and sequential nature of oil palm growth (around 40 stacked phytomers) and yield (fruit bunches axillated on each phytomer), we developed a specific sub-canopy structure for simulating palm's growth and yield within the framework of the Community Land Model (CLM4.5). In this structure each phytomer has its own prognostic leaf growth and fruit yield capacity like a PFT but with shared stem and root components among all phytomers. Phenology and carbon and nitrogen allocation operate on the different phytomers in parallel but at unsynchronized steps, so that multiple fruit yields per annum are enabled in terms of carbon and nitrogen outputs. An important phenological phase is identified for the palm PFT – the storage growth period of bud and "spear" leaves which are photosynthetically inactive before expansion. Agricultural practices such as transplanting, fertilization, and leaf pruning are represented. Parameters introduced for the new PFT were calibrated and validated with field measurements of leaf area index (LAI) and yield from Sumatra, Indonesia. In calibration with a mature oil palm plantation, the cumulative yields from 2005 to 2014 matched perfectly between simulation and observation (mean percentage error = 4 %). Simulated inter-annual dynamics of PFT-level and phytomer-level LAI were both within the range of field measurements. Validation from eight independent oil palm sites shows the ability of the model to adequately predict the average leaf growth and fruit yield across sites but also indicates that seasonal dynamics and site-to-site variability of yield are driven by processes not yet implemented in the model. The new sub-canopy structure and phenology and allocation functions now allow exploring the effects of tropical land use change, from natural ecosystems to oil palm plantations, on carbon, water and energy cycles and regional climate.

Description: In order to quantify the effects of forests to oil palm conversion occurring in the tropics on land–atmosphere carbon, water and energy fluxes, we develop a new perennial crop sub-model CLM-Palm for simulating a palm plant functional type (PFT) within the framework of the Community Land Model (CLM4.5). CLM-Palm is tested here on oil palm only but is meant of generic interest for other palm crops (e.g., coconut). The oil palm has monopodial morphology and sequential phenology of around 40 stacked phytomers, each carrying a large leaf and a fruit bunch, forming a multilayer canopy. A sub-canopy phenological and physiological parameterization is thus introduced so that each phytomer has its own prognostic leaf growth and fruit yield capacity but with shared stem and root components. Phenology and carbon and nitrogen allocation operate on the different phytomers in parallel but at unsynchronized steps, separated by a thermal period. An important phenological phase is identified for the oil palm – the storage growth period of bud and "spear" leaves which are photosynthetically inactive before expansion. Agricultural practices such as transplanting, fertilization and leaf pruning are represented. Parameters introduced for the oil palm were calibrated and validated with field measurements of leaf area index (LAI), yield and net primary production (NPP) from Sumatra, Indonesia. In calibration with a mature oil palm plantation, the cumulative yields from 2005 to 2014 matched notably well between simulation and observation (mean percentage error = 3 %). Simulated inter-annual dynamics of PFT-level and phytomer-level LAI were both within the range of field measurements. Validation from eight independent oil palm sites shows the ability of the model to adequately predict the average leaf growth and fruit yield across sites and sufficiently represent the significant nitrogen- and age-related site-to-site variability in NPP and yield. Results also indicate that seasonal dynamics of yield and remaining small-scale site-to-site variability of NPP are driven by processes not yet implemented in the model or reflected in the input data. The new sub-canopy structure and phenology and allocation functions in CLM-Palm allow exploring the effects of tropical land-use change, from natural ecosystems to oil palm plantations, on carbon, water and energy cycles and regional climate.