Modeling water availability for trees in tropical forests
Highlights
► A daily water balance model for tropical forests. ► A new statistical way to estimate permanent wilting point and field capacity with time domain reflectometer data. ► Explicit modeling of water taken up by tree roots. ► Quantification of the duration and intensity of tree water stress during dry season in rainforests.
Introduction
Despite annual precipitation that always exceeds 1500 mm year−1, most of the Amazon’s neotropical forests experience some annual dry season (less than 100 mm per month), that is variable in both duration and intensity (Malhi and Wright, 2004, Sombroek, 2001, Xiao et al., 2006, Marengo, 1992).
The consequences of annual drought on tropical forest functioning include a decrease in growth primary production and ecosystem respiration (Goulden et al., 2004, Hutyra et al., 2007, Bonal et al., 2008), and a reduction in tropical tree fluxes for both carbon (Bonal et al., 2000, Miranda et al., 2005) and water fluxes (Fisher et al., 2006). Very recently, an analysis of tree responses to the intense 2005 dry season highlighted the vulnerability of neotropical forests to moisture stress, with the potential for positive feedbacks on climate change due to increased tree mortality (Phillips et al., 2009).
Climate modeling scenarios suggest that the dry season in north-eastern Amazonian forests might lengthen during the 21st century (Cox et al., 2000, Cox et al., 2004, Malhi and Wright, 2004, Malhi et al., 2009). The short-term effect of soil water availability deficits on tropical tree growth, mortality and carbon and water fluxes has recently been quantified under experimentally controlled conditions (Fisher et al., 2007, Nepstad et al., 2007). Long-term inventory plots with regular tree censuses (growth, recruitment, mortality) have been set-up widely in the past few decades throughout Amazonia (Phillips et al., 2010, Clark, 2004, Wagner et al., 2010). These plots offer an unexpected opportunity to analyze the impact of soil water availability on tropical forest dynamics on a large temporal and spatial scale (Clark, 2007). However, to the best of our knowledge, no soil water balance model explicitly accounting for tropical soil and climate characteristics, and able to compute available water for the trees on a plot scale, has ever been developed. The relation between amount of rainfall and water availability for trees is not straightforward and determined by various plant characteristics, such as the root distribution, and soil characteristics, such as the permanent wilting point and the field capacity. By contrast, other widely studied climatic variables such as light and temperature give a relatively direct indication on their effect on forest dynamics (Graham et al., 2003, Clark et al., 2010). Soil water availability to the trees, which can be characterized by Relative Extractable Water (REW, i.e. daily available water standardized by maximum available water), depends on soil characteristics such as structure, texture, composition and porosity, as well as on the rate of water uptake by the trees. Different soil water balance models have been used in Amazonian tropical forests to estimate drought implications for forest flammability and tree growth (Nepstad et al., 2004), to reproduce hydrologic processes (Belk et al., 2007), to evaluate soil water controls on evapotranspiration (Fisher et al., 2007), or to evaluate the importance of deep root uptake (Markewitz et al., 2010). However, none of these models aims to estimate REW. The nearest estimate of REW is so-called plant available water (PAW) described by Nepstad et al. (2004). The spatial resolution of PAW, i.e. 8 km, is too large for use in any precise modeling of the impact of soil drought conditions on tree growth, mortality and/or recruitment. Furthermore, the modeling framework never explicitly simulates the amount of water taken up by tree roots. Modeling approaches designed to estimate REW have already been developed for temperate forests and, for instance, were used to assess soil water control on carbon and water dynamics in European forests during the 2003 drought (Granier et al., 2007). Such temperate models are not suitable for tropical forests. For instance, the polynomial rainfall interception submodel is unsuited to the stand characteristics of tropical forests. Another limit is that water extraction is not modeled and field data is needed for root density, meaning that soil pits need to be dug to quantify vertical root distribution, also meaning that the strong assumption must be made that water absorption by roots is proportional to root distribution.
In this paper we introduce a locally parameterized soil water budget model inspired by the BILJOU temperate model (Granier et al., 1999). As performed by BILJOU, this model estimates soil water availability, stand transpiration and rainfall interception in tropical forests with a daily time step and for different soil types (Fig. 1). Model inputs are daily rainfall, annual means of potential evapotranspiration (PET) and solar radiation, and averaged plant area index (PAI). The soil is filled by rainfall water passing through the canopy. The amount of rainfall intercepted by the canopy is computed in a submodel adapted to tropical forests (Gash et al., 1995). In our model, the soil consists of a succession of fine layers, each of which has a unique field capacity and permanent wilting point. We developed a new method using a Bayesian framework to estimate these two parameters using only Time Domain Reflectometer (TDR) measurements. When the water in a given layer exceeds water content at field capacity, drainage occurs and water fills the next layer, etc. Water extraction from soil layers is due to tree transpiration in addition to soil and understorey evapotranspiration. Soil evaporation and understorey transpiration are computed based on equations developed by Granier et al. (1999), and are assumed to be proportional to the energy reaching the understorey; tree transpiration is computed using potential evapotranspiration. Both understorey and tree transpiration are extracted in accordance with estimated root distribution.
This paper has three specific objectives: (i) to present our water balance model and to describe the different submodels it contains; (ii) to present an original statistical method used to estimate permanent wilting point, field capacity, and root distribution based on Time-Domain Reflectometer (TDR) data only; and (iii) to parameterize and validate the model using TDR data collected on a soil topographic gradient in Paracou, French Guiana.
Section snippets
List of symbols and abbreviations
- P
precipitation
- In
rainfall interception
- Th
throughfall
- Tr
tree transpiration
- PET
potential evapotranspiration
- ρ
ratio Tr/ETP
- Eu
understorey and soil evapotranspiration
- Dr
drainage in depth
- θPWP,l
soil water content at permanent wilting point for layer l
- θFC,l
soil water content at field capacity of soil layer l
- λRfd
root density parameter
- EWl
extractable water of soil layer l
maximum extractable water of soil layer l, θFC,l − θPWP,l
- REW
relative extractable water of the soil
- REWcl
REW of soil layer l, critical when
The model
Our daily water balance model was constructed using discrete-time deterministic formalism. The model contains and interconnects four submodels that compute interception, water infiltration and deep drainage, tree transpiration, and soil plus understorey evapotranspiration (Fig. 1). First, part of the precipitation is intercepted by the canopy, then the remaining part reaches the soil surface and fills the soil. If the soil is at field capacity, the excess water is lost by deep drainage. Tree
Site descriptions and experimental setup
The study site used for calibration is located in Paracou, French Guiana (5°18′N, 52°55′W), a lowland tropical rain forest near Sinnamary (Gourlet-Fleury et al., 2004). The forest is typical of Guianan rainforests (ter Steege et al., 2006). More than 550 woody species attaining 2 cm DBH (Diameter at Breast Height) have been described at the site, with an estimated 160–180 species of trees ≥10 cm DBH per hectare. The dominant families at the site include Leguminoseae, Chrysobalanaceae,
Model resolution levels
The best resolution selected with the BICM was the forest level, BICM = − 4148.4 (soil level, BICM = − 2066.7; tube level, BICM = − 1250.0). Extractable water RMSEP at all model resolutions ranged from 1.5 to 4.5% of water content, Table 4. With regard to the running soil level in M2, the RMSEP for Alt, SLD, DhS and UhS corresponded to 2.631, 3.176, 4.184 and 1.570 % of water content, respectively.
Ratio ρ, tree transpiration/PET,
At all model resolution levels, i.e. tube (M1), soil (M2) or forest (M3), forest stand transpiration
Discussion
In this study we constructed a model of soil water balance for tropical forests and parameterized this model using a dataset collected in Paracou, French Guiana. We then validated this model on an independent dataset from the same forest site. Extractable water is accurately predicted at the three levels of model resolution, i.e. TDR tube, soil type and forest stand, Table 4. At the soil level, RMSEP was always less than 4.2%. This is remarkable considering that precipitation is the only daily
Conclusions
In this study we developed, calibrated and validated a daily soil water balance model for use in tropical forests. Precipitation is the only data required with daily precision. The model works with few parameters, most of which are available in the specialized literature. We put forward a novel method to estimate the remaining site-specific parameters, θFC, θPWP and λRfd, using TDR data only, which should facilitate calibration in other tropical forest sites. The use of a Bayesian framework is
Acknowledgements
This manuscript is part of the CLIMFOR project (FRB grant to BH). This study has been possible thanks to the important field work realized at Paracou by B. Burban, J.-Y. Goret, M. Desprez and V. Freycon. We also thank C. Baraloto and two anonymous reviewers who greatly helped improving a previous version of this manuscript.
References (62)
- et al.
Validation of soil moisture simulations from the PAMII model, and an assessment of their sensitivity to uncertainties in soil hydraulic parameters
Agricultural and Forest Meteorology
(2010) - et al.
Interception water-partitioning dynamics for a pristine rainforest in Central Amazonia: marked differences between normal and dry years
Agricultural and Forest Meteorology
(2007) - et al.
Evaluating climatic and soil water controls on evapotranspiration at two Amazonian rainforest sites
Agricultural and Forest Meteorology
(2008) - et al.
Estimating sparse forest rainfall interception with an analytical model
Journal of Hydrology
(1995) - et al.
A lumped water balance model to evaluate duration and intensity of drought constraints in forest stands
Ecological Modelling
(1999) - et al.
Transpiration of natural rain forest and its dependence on climatic factors
Agricultural and Forest Meteorology
(1996) - et al.
Evidence for soil water control on carbon and water dynamics in European forests during the extremely dry year: 2003
Agricultural and Forest Meteorology
(2007) - et al.
The measurement and modelling of rainfall interception by Amazonian rain forest
Agricultural and Forest Meteorology
(1988) - et al.
A review of rainfall interception modelling
Journal of Hydrology
(2009) - et al.
A predictive model of rainfall interception in forests. i. Derivation of the model from observations in a plantation of Corsican pine
Agricultural Meteorology
(1971)
Modelling interception loss for two sparse eucalypt and pine forests in central Portugal using reformulated Rutter and Gash analytical models
Journal of Hydrology
Spatial and temporal variability of canopy structure in a tropical moist forest
Acta Oecologica-International Journal of Ecology
Detecting leaf phenology of seasonally moist tropical forests in South America with multi-temporal MODIS images
Remote Sensing of Environment
Modeling the effects of throughfall reduction on soil water content in a Brazilian Oxisol under a moist tropical forest
Water Resources Research
Water acquisition patterns of two wet tropical canopy tree species of French Guiana as inferred from (H2O)-O-18 extraction profiles
Annals of Forest Science
Impact of severe dry season on net ecosystem exchange in the Neotropical rainforest of French Guiana
Global Change Biology
Temperate forest trees and stands under severe drought: a review of ecophysiological responses, adaptation processes and long-term consequences
Annals of Forest Science
Sources or sinks? The responses of tropical forests to current and future climate and atmospheric composition
Philosophical Transactions of the Royal Society of London Series B-Biological Sciences
Detecting tropical forests’ responses to global climatic and atmospheric change: current challenges and a way forward
Biotropica
Annual wood production in a tropical rain forest in NE Costa Rica linked to climatic variation but not to increasing CO2
Global Change Biology
Fast determination of light availability and leaf area index in tropical forests
Journal of Tropical Ecology
Amazonian forest dieback under climate-carbon cycle projections for the 21st century
Theoretical and Applied Climatology
Acceleration of global warming due to carbon-cycle feedbacks in a coupled climate model
Nature
Topographic position affects the water regime in a semideciduous tropical forest in Panama
Plant and Soil
A comparison of five indirect methods for characterizing the light environment in a tropical forest
Annals of Forest Science
Higher treefall rates on slopes and waterlogged soils result in lower stand biomass and productivity in a tropical rain forest
Journal of Ecology
Evidence from Amazonian forests is consistent with isohydric control of leaf water potential
Plant Cell and Environment
The response of an Eastern Amazonian rain forest to drought stress: results and modelling analyses from a throughfall exclusion experiment
Global Change Biology
Analytical model of rainfall interception by forests
Quarterly Journal of the Royal Meteorological Society
Cited by (53)
Diverging taxonomic and functional trajectories following disturbance in a Neotropical forest
2020, Science of the Total EnvironmentBeyond species richness and biomass: Impact of selective logging and silvicultural treatments on the functional composition of a neotropical forest
2019, Forest Ecology and ManagementCitation Excerpt :Data were collected at the Paracou experimental site (5°18′N, 52°55′W), a lowland tropical rain forest near Sinnamary, French Guiana. The site receives nearly two-thirds of the annual 3041 mm of precipitation between mid-March and mid-June, and <50 mm per month in September and October (Wagner et al., 2011 Agricultural). The floristic composition is typical of littoral guianan rainforests with dominant families including Leguminoseae, Chrysobalanaceae, Lecythidaceae, Sapotaceae and Burseraceae.
Assessing timber volume recovery after disturbance in tropical forests – A new modelling framework
2018, Ecological ModellingCitation Excerpt :The study is based on data from Paracou research station (5°18′ N, 52°55′ W), a long-term large-scale disturbance experiment located in a lowland tropical forest in French Guiana (Gourlet-Fleury et al., 2004). The climate is affected by the north/south movements of the Inter-Tropical Convergence Zone and the site receives nearly two-thirds of its annual 3041 mm of precipitation between mid-March and mid-June, and less than 50 mm per month in September and October (Wagner et al., 2011). The forest composition is typical of the Guyana Shield rainforests (ter Steege et al., 2013), dominated by Chrysobalanaceae, Fabaceae and Lecythidaceae, and with approximately 180 species of trees ≥10 cm DBH per ha.
What drives long-term variations in carbon flux and balance in a tropical rainforest in French Guiana?
2018, Agricultural and Forest MeteorologyA foliar disease simulation model to assist the design of new control methods against black leaf streak disease of banana
2017, Ecological ModellingCitation Excerpt :The importance of the incubation period has already been demonstrated in wheat through sensitivity analysis of a generic model (Casadebaig et al., 2012). To infer parameters, a Metropolis–Hastings within Gibbs algorithm (Gouache et al., 2012; Wagner et al., 2011) was first validated on simulated data (Appendix C) and was then successfully applied to the data (Appendix D) acquired in an experiment specifically designed to quantitatively describe the dynamics of the epidemic at the leaf compartment scale. The posterior distributions obtained for three parameters allowed us to update our knowledge on the pathogen life cycle parameter of a susceptible plantain variety (Table 3) for which no value is available in the literature, as exemplified in the present study by the Macho por Hembra plantain.
Modelling hydraulic functioning of an adult beech stand under non-limiting soil water and severe drought condition
2017, Ecological ModellingCitation Excerpt :Besides simulating water fluxes, other major outputs are water potentials and internal water reservoirs status, both quantifying trees constraint level. Moreover, in most models the connections between atmospheric demand, soil water availability and tree hydraulic pathway are highly simplified (Chuang et al., 2006; Kuchment et al., 2006; Fisher et al., 2007; Wagner et al., 2011) and stomatal functioning poorly described in detail (Whitley et al., 2009). Some models are mostly focusing on hydraulic function (Tyree 1988), or on the soil-root compartments (Williams et al., 2001; Nadezhdina et al., 2008), or do not distinguish compartments in crowns and soil (Steppe et al., 2005; Zweifel et al., 2007).