RESUMO
While droughts predominantly induce immediate reductions in plant carbon uptake, they can also exert long-lasting effects on carbon fluxes through associated changes in leaf area, soil carbon, etc. Among other mechanisms, shifts in carbon allocation due to water stress can contribute to the legacy effects of drought on carbon fluxes. However, the magnitude and impact of these allocation shifts on carbon fluxes and pools remain poorly understood. Using data from a wet tropical flux tower site in French Guiana, we demonstrate that drought-induced carbon allocation shifts can be reliably inferred by assimilating Net Biosphere Exchange (NBE) and other observations within the CARbon DAta MOdel fraMework. This model-data fusion system allows inference of optimized carbon and water cycle parameters and states from multiple observational data streams. We then examined how these inferred shifts affected the duration and magnitude of drought's impact on NBE during and after the extreme event. Compared to a static allocation scheme analogous to those typically implemented in land surface models, dynamic allocation reduced average carbon uptake during drought recovery by a factor of 2.8. Additionally, the dynamic model extended the average recovery time by 5 months. The inferred allocation shifts influenced the post-drought period by altering foliage and fine root pools, which in turn modulated gross primary productivity and heterotrophic respiration for up to a decade. These changes can create a bust-boom cycle where carbon uptake is enhanced some years after a drought, compared to what would have occurred under drought-free conditions. Overall, allocation shifts accounted for 65% [45%-75%] of drought legacy effects in modeled NBE. In summary, drought-induced carbon allocation shifts can play a substantial role in the enduring influence of drought on cumulative land-atmosphere CO2 exchanges and should be accounted for in ecosystem models.
Assuntos
Ciclo do Carbono , Secas , Clima Tropical , Guiana Francesa , Florestas , Carbono/metabolismo , Modelos TeóricosRESUMO
Knowledge of the physiological mechanisms underlying species vulnerability to drought is critical for better understanding patterns of tree mortality. Investigating plant adaptive strategies to drought should thus help to fill this knowledge gap, especially in tropical rainforests exhibiting high functional diversity. In a semi-controlled drought experiment using 12 rainforest tree species, we investigated the diversity in hydraulic strategies and whether they determined the ability of saplings to use stored non-structural carbohydrates during an extreme imposed drought. We further explored the importance of water- and carbon-use strategies in relation to drought survival through a modelling approach. Hydraulic strategies varied considerably across species with a continuum between dehydration tolerance and avoidance. During dehydration leading to hydraulic failure and irrespective of hydraulic strategies, species showed strong declines in whole-plant starch concentrations and maintenance, or even increases in soluble sugar concentrations, potentially favouring osmotic adjustments. Residual water losses mediated the trade-off between time to hydraulic failure and growth, indicating that dehydration avoidance is an effective drought-survival strategy linked to the 'fast-slow' continuum of plant performance at the sapling stage. Further investigations on residual water losses may be key to understanding the response of tropical rainforest tree communities to climate change.
Assuntos
Secas , Floresta Úmida , Árvores , Água , Árvores/fisiologia , Árvores/crescimento & desenvolvimento , Água/metabolismo , Desidratação , Clima TropicalRESUMO
Tropical forests play a pivotal role in regulating the global carbon cycle. However, the response of these forests to changes in absorbed solar energy and water supply under the changing climate is highly uncertain. Three-year (2018-2021) spaceborne high-resolution measurements of solar-induced chlorophyll fluorescence (SIF) from the TROPOspheric Monitoring Instrument (TROPOMI) provide a new opportunity to study the response of gross primary production (GPP) and more broadly tropical forest carbon dynamics to differences in climate. SIF has been shown to be a good proxy for GPP on monthly and regional scales. Combining tropical climate reanalysis records and other contemporary satellite products, we find that on the seasonal timescale, the dependence of GPP on climate variables is highly heterogeneous. Following the principal component analyses and correlation comparisons, two regimes are identified: water limited and energy limited. GPP variations over tropical Africa are more correlated with water-related factors such as vapor pressure deficit (VPD) and soil moisture, while in tropical Southeast Asia, GPP is more correlated with energy-related factors such as photosynthetically active radiation (PAR) and surface temperature. Amazonia is itself heterogeneous: with an energy-limited regime in the north and water-limited regime in the south. The correlations of GPP with climate variables are supported by other observation-based products, such as Orbiting Carbon Observatory-2 (OCO2) SIF and FluxSat GPP. In each tropical continent, the coupling between SIF and VPD increases with the mean VPD. Even on the interannual timescale, the correlation of GPP with VPD is still discernable, but the sensitivity is smaller than the intra-annual correlation. By and large, the dynamic global vegetation models in the TRENDY v8 project do not capture the high GPP seasonal sensitivity to VPD in dry tropics. The complex interactions between carbon and water cycles in the tropics illustrated in this study and the poor representation of this coupling in the current suite of vegetation models suggest that projections of future changes in carbon dynamics based on these models may not be robust.
RESUMO
Water stress can cause declines in plant function that persist after rehydration. Recent work has defined 'resilience' traits characterizing leaf resistance to persistent damage from drought, but whether these traits predict resilience in whole-plant function is unknown. It is also unknown whether the coordination between resilience and 'resistance' - the ability to maintain function during drought - observed globally occurs within ecosystems. For eight rainforest species, we dehydrated and subsequently rehydrated leaves, and measured water stress thresholds for declines in rehydration capacity and maximum quantum yield of photosystem II (Fv /Fm ). We tested correlations with embolism resistance and dry season water potentials (ΨMD ), and calculated safety margins for damage (ΨMD - thresholds) and tested correlations with drought resilience in sap flow and growth. Ψ thresholds for persistent declines in Fv /Fm , indicating resilience, were positively correlated with ΨMD and thresholds for leaf vein embolism. Safety margins for persistent declines in Fv /Fm , but not rehydration capacity, were positively correlated with drought resilience in sap flow. Correlations between resistance and resilience suggest that species' differences in performance during drought are perpetuated after drought, potentially accelerating shifts in forest composition. Resilience to photochemical damage emerged as a promising functional trait to characterize whole-plant drought resilience.
Assuntos
Desidratação , Floresta Úmida , Ecossistema , Secas , Folhas de Planta , ÁrvoresRESUMO
In a global context where water will become a scarce resource under temperate latitudes, managing tree plantations with species associations, i.e., forest mixture or agroforestry, could play a major role in optimizing the sustainable use of this resource. Conceptual frameworks in community ecology suggest that, in mixed plantations, environmental resources such as water may be more efficiently used for carbon acquisition and tree growth thanks to niche complementarity among species. To test the hypotheses behind these conceptual frameworks, we estimated water-use efficiency (WUE) for poplar trees grown in a monoculture, in association with alder trees (forest mixture) and in association with clover leys (agroforestry) in an experimental plantation located in northeastern France. Water-use efficiency was estimated (i) at leaf level through gas exchange measurements and analysis of carbon isotope composition, (ii) at wood level through carbon isotope composition and (iii) at tree level with sap flow sensors and growth increment data. We hypothesized that species interactions would increase WUE of poplars in mixtures due to a reduction in competition and/or facilitation effects due to the presence of the N2-fixing species in mixtures. Poplar trees in both mixture types showed higher WUE than those in the monoculture. The differences we found in WUE between the monoculture and the agroforestry treatment were associated to differences in stomatal conductance and light-saturated net CO2 assimilation rate (at the leaf level) and transpiration (at the tree level), while the differences between the monoculture and the forest mixture were more likely due to differences in stomatal conductance at the leaf level and both transpiration and biomass accumulation at the tree level. Moreover, the more WUE was integrated in time (instantaneous gas exchanges < leaf life span < seasonal wood core < whole tree), the more the differences among treatments were marked.
Assuntos
Populus , Árvores , Água , Transpiração Vegetal , Folhas de Planta , Isótopos de Carbono/análiseRESUMO
Relationships among productivity, leaf phenology, and seasonal variation in moisture and light availability are poorly understood for evergreen broadleaved tropical/subtropical forests, which contribute 25% of terrestrial productivity. On the one hand, as moisture availability declines, trees shed leaves to reduce transpiration and the risk of hydraulic failure. On the other hand, increases in light availability promote the replacement of senescent leaves to increase productivity. Here, we provide a comprehensive framework that relates the seasonality of climate, leaf abscission, and leaf productivity across the evergreen broadleaved tropical/subtropical forest biome. The seasonal correlation between rainfall and light availability varies from strongly negative to strongly positive across the tropics and maps onto the seasonal correlation between litterfall mass and productivity for 68 forests. Where rainfall and light covary positively, litterfall and productivity also covary positively and are always greater in the wetter sunnier season. Where rainfall and light covary negatively, litterfall and productivity are always greater in the drier and sunnier season if moisture supplies remain adequate; otherwise productivity is smaller in the drier sunnier season. This framework will improve the representation of tropical/subtropical forests in Earth system models and suggests how phenology and productivity will change as climate change alters the seasonality of cloud cover and rainfall across tropical/subtropical forests.
RESUMO
The leaf economics spectrum1,2 and the global spectrum of plant forms and functions3 revealed fundamental axes of variation in plant traits, which represent different ecological strategies that are shaped by the evolutionary development of plant species2. Ecosystem functions depend on environmental conditions and the traits of species that comprise the ecological communities4. However, the axes of variation of ecosystem functions are largely unknown, which limits our understanding of how ecosystems respond as a whole to anthropogenic drivers, climate and environmental variability4,5. Here we derive a set of ecosystem functions6 from a dataset of surface gas exchange measurements across major terrestrial biomes. We find that most of the variability within ecosystem functions (71.8%) is captured by three key axes. The first axis reflects maximum ecosystem productivity and is mostly explained by vegetation structure. The second axis reflects ecosystem water-use strategies and is jointly explained by variation in vegetation height and climate. The third axis, which represents ecosystem carbon-use efficiency, features a gradient related to aridity, and is explained primarily by variation in vegetation structure. We show that two state-of-the-art land surface models reproduce the first and most important axis of ecosystem functions. However, the models tend to simulate more strongly correlated functions than those observed, which limits their ability to accurately predict the full range of responses to environmental changes in carbon, water and energy cycling in terrestrial ecosystems7,8.
Assuntos
Ciclo do Carbono , Ecossistema , Plantas/metabolismo , Ciclo Hidrológico , Dióxido de Carbono/metabolismo , Clima , Conjuntos de Dados como Assunto , Umidade , Plantas/classificação , Análise de Componente PrincipalRESUMO
Ecosystem respiration is a major component of the global terrestrial carbon cycle and is strongly influenced by temperature. The global extent of the temperature-ecosystem respiration relationship, however, has not been fully explored. Here, we test linear and threshold models of ecosystem respiration across 210 globally distributed eddy covariance sites over an extensive temperature range. We find thresholds to the global temperature-ecosystem respiration relationship at high and low air temperatures and mid soil temperatures, which represent transitions in the temperature dependence and sensitivity of ecosystem respiration. Annual ecosystem respiration rates show a markedly reduced temperature dependence and sensitivity compared to half-hourly rates, and a single mid-temperature threshold for both air and soil temperature. Our study indicates a distinction in the influence of environmental factors, including temperature, on ecosystem respiration between latitudinal and climate gradients at short (half-hourly) and long (annual) timescales. Such climatological differences in the temperature sensitivity of ecosystem respiration have important consequences for the terrestrial net carbon sink under ongoing climate change.
Assuntos
Ciclo do Carbono , Ecossistema , Respiração , Solo , TemperaturaRESUMO
Despite their low contribution to forest carbon stocks, lianas (woody vines) play an important role in the carbon dynamics of tropical forests. As structural parasites, they hinder tree survival, growth and fecundity; hence, they negatively impact net ecosystem productivity and long-term carbon sequestration.Competition (for water and light) drives various forest processes and depends on the local abundance of resources over time. However, evaluating the relative role of resource availability on the interactions between lianas and trees from empirical observations is particularly challenging. Previous approaches have used labour-intensive and ecosystem-scale manipulation experiments, which are infeasible in most situations.We propose to circumvent this challenge by evaluating the uncertainty of water and light capture processes of a process-based vegetation model (ED2) including the liana growth form. We further developed the liana plant functional type in ED2 to mechanistically simulate water uptake and transport from roots to leaves, and start the model from prescribed initial conditions. We then used the PEcAn bioinformatics platform to constrain liana parameters and run uncertainty analyses.Baseline runs successfully reproduced ecosystem gas exchange fluxes (gross primary productivity and latent heat) and forest structural features (leaf area index, aboveground biomass) in two sites (Barro Colorado Island, Panama and Paracou, French Guiana) characterized by different rainfall regimes and levels of liana abundance.Model uncertainty analyses revealed that water limitation was the factor driving the competition between trees and lianas at the drier site (BCI), and during the relatively short dry season of the wetter site (Paracou). In young patches, light competition dominated in Paracou but alternated with water competition between the wet and the dry season on BCI according to the model simulations.The modelling workflow also identified key liana traits (photosynthetic quantum efficiency, stomatal regulation parameters, allometric relationships) and processes (water use, respiration, climbing) driving the model uncertainty. They should be considered as priorities for future data acquisition and model development to improve predictions of the carbon dynamics of liana-infested forests. Synthesis. Competition for water plays a larger role in the interaction between lianas and trees than previously hypothesized, as demonstrated by simulations from a process-based vegetation model.
RESUMO
Wood density (WD) is often used as a proxy for hydraulic traits such as vulnerability to drought-induced xylem cavitation and maximum water transport capacity, with dense-wooded species generally being more resistant to drought-induced xylem cavitation, having lower rates of maximum water transport and lower sapwood capacitance than light-wooded species. However, relationships between WD and the hydraulic traits that they aim to predict have not been well established in tropical forests, where modeling is necessary to predict drought responses for a high diversity of unmeasured species. We evaluated WD and relationships with stem xylem vulnerability by measuring cavitation curves, sapwood water release curves and minimum seasonal water potential (Ψmin) on upper canopy branches of six tree species and three liana species from a single wet tropical forest site in Panama. The objective was to better understand coordination and trade-offs among hydraulic traits and the potential utility of these relationships for modeling purposes. We found that parameters from sapwood water release curves such as capacitance, saturated water content and sapwood turgor loss point (Ψtlp,x) were related to WD, whereas stem vulnerability curve parameters were not. However, the water potential corresponding to 50% loss of hydraulic conductivity (P50) was related to Ψtlp,x and sapwood osmotic potential at full turgor (πo,x). Furthermore, species with lower Ψmin showed lower P50, Ψtlp,x and πo,x suggesting greater drought resistance. Our results indicate that WD is a good easy-to-measure proxy for some traits related to drought resistance, but not others. The ability of hydraulic traits such as P50 and Ψtlp,x to predict mortality must be carefully examined if WD values are to be used to predict drought responses in species without detailed physiological measurements.
Assuntos
Secas , Árvores , Panamá , Folhas de Planta , Água , Madeira , XilemaRESUMO
The carbon sink capacity of tropical forests is substantially affected by tree mortality. However, the main drivers of tropical tree death remain largely unknown. Here we present a pan-Amazonian assessment of how and why trees die, analysing over 120,000 trees representing > 3800 species from 189 long-term RAINFOR forest plots. While tree mortality rates vary greatly Amazon-wide, on average trees are as likely to die standing as they are broken or uprooted-modes of death with different ecological consequences. Species-level growth rate is the single most important predictor of tree death in Amazonia, with faster-growing species being at higher risk. Within species, however, the slowest-growing trees are at greatest risk while the effect of tree size varies across the basin. In the driest Amazonian region species-level bioclimatic distributional patterns also predict the risk of death, suggesting that these forests are experiencing climatic conditions beyond their adaptative limits. These results provide not only a holistic pan-Amazonian picture of tree death but large-scale evidence for the overarching importance of the growth-survival trade-off in driving tropical tree mortality.
Assuntos
Ecologia , Florestas , Árvores/crescimento & desenvolvimento , Biomassa , Brasil , Dióxido de Carbono , Sequestro de Carbono , Ecossistema , Monitoramento Ambiental , Modelos Biológicos , Modelos de Riscos Proporcionais , Fatores de Risco , Clima TropicalRESUMO
Selective logging, fragmentation, and understory fires directly degrade forest structure and composition. However, studies addressing the effects of forest degradation on carbon, water, and energy cycles are scarce. Here, we integrate field observations and high-resolution remote sensing from airborne lidar to provide realistic initial conditions to the Ecosystem Demography Model (ED-2.2) and investigate how disturbances from forest degradation affect gross primary production (GPP), evapotranspiration (ET), and sensible heat flux (H). We used forest structural information retrieved from airborne lidar samples (13,500 ha) and calibrated with 817 inventory plots (0.25 ha) across precipitation and degradation gradients in the eastern Amazon as initial conditions to ED-2.2 model. Our results show that the magnitude and seasonality of fluxes were modulated by changes in forest structure caused by degradation. During the dry season and under typical conditions, severely degraded forests (biomass loss ≥66%) experienced water stress with declines in ET (up to 34%) and GPP (up to 35%) and increases of H (up to 43%) and daily mean ground temperatures (up to 6.5°C) relative to intact forests. In contrast, the relative impact of forest degradation on energy, water, and carbon cycles markedly diminishes under extreme, multiyear droughts, as a consequence of severe stress experienced by intact forests. Our results highlight that the water and energy cycles in the Amazon are driven by not only climate and deforestation but also the past disturbance and changes of forest structure from degradation, suggesting a much broader influence of human land use activities on the tropical ecosystems.
RESUMO
We apply and compare three widely applicable methods for estimating ecosystem transpiration (T) from eddy covariance (EC) data across 251 FLUXNET sites globally. All three methods are based on the coupled water and carbon relationship, but they differ in assumptions and parameterizations. Intercomparison of the three daily T estimates shows high correlation among methods (R between .89 and .94), but a spread in magnitudes of T/ET (evapotranspiration) from 45% to 77%. When compared at six sites with concurrent EC and sap flow measurements, all three EC-based T estimates show higher correlation to sap flow-based T than EC-based ET. The partitioning methods show expected tendencies of T/ET increasing with dryness (vapor pressure deficit and days since rain) and with leaf area index (LAI). Analysis of 140 sites with high-quality estimates for at least two continuous years shows that T/ET variability was 1.6 times higher across sites than across years. Spatial variability of T/ET was primarily driven by vegetation and soil characteristics (e.g., crop or grass designation, minimum annual LAI, soil coarse fragment volume) rather than climatic variables such as mean/standard deviation of temperature or precipitation. Overall, T and T/ET patterns are plausible and qualitatively consistent among the different water flux partitioning methods implying a significant advance made for estimating and understanding T globally, while the magnitudes remain uncertain. Our results represent the first extensive EC data-based estimates of ecosystem T permitting a data-driven perspective on the role of plants' water use for global water and carbon cycling in a changing climate.
Assuntos
Ecossistema , Transpiração Vegetal , Poaceae , Chuva , Solo , ÁguaRESUMO
The FLUXNET2015 dataset provides ecosystem-scale data on CO2, water, and energy exchange between the biosphere and the atmosphere, and other meteorological and biological measurements, from 212 sites around the globe (over 1500 site-years, up to and including year 2014). These sites, independently managed and operated, voluntarily contributed their data to create global datasets. Data were quality controlled and processed using uniform methods, to improve consistency and intercomparability across sites. The dataset is already being used in a number of applications, including ecophysiology studies, remote sensing studies, and development of ecosystem and Earth system models. FLUXNET2015 includes derived-data products, such as gap-filled time series, ecosystem respiration and photosynthetic uptake estimates, estimation of uncertainties, and metadata about the measurements, presented for the first time in this paper. In addition, 206 of these sites are for the first time distributed under a Creative Commons (CC-BY 4.0) license. This paper details this enhanced dataset and the processing methods, now made available as open-source codes, making the dataset more accessible, transparent, and reproducible.
RESUMO
The sensitivity of tropical forest carbon to climate is a key uncertainty in predicting global climate change. Although short-term drying and warming are known to affect forests, it is unknown if such effects translate into long-term responses. Here, we analyze 590 permanent plots measured across the tropics to derive the equilibrium climate controls on forest carbon. Maximum temperature is the most important predictor of aboveground biomass (-9.1 megagrams of carbon per hectare per degree Celsius), primarily by reducing woody productivity, and has a greater impact per °C in the hottest forests (>32.2°C). Our results nevertheless reveal greater thermal resilience than observations of short-term variation imply. To realize the long-term climate adaptation potential of tropical forests requires both protecting them and stabilizing Earth's climate.
Assuntos
Ciclo do Carbono , Mudança Climática , Florestas , Temperatura Alta , Árvores/metabolismo , Clima Tropical , Aclimatação , Biomassa , Carbono/metabolismo , Planeta Terra , MadeiraRESUMO
Transpiration in humid tropical forests modulates the global water cycle and is a key driver of climate regulation. Yet, our understanding of how tropical trees regulate sap flux in response to climate variability remains elusive. With a progressively warming climate, atmospheric evaporative demand [i.e., vapor pressure deficit (VPD)] will be increasingly important for plant functioning, becoming the major control of plant water use in the twenty-first century. Using measurements in 34 tree species at seven sites across a precipitation gradient in the neotropics, we determined how the maximum sap flux velocity (vmax) and the VPD threshold at which vmax is reached (VPDmax) vary with precipitation regime [mean annual precipitation (MAP); seasonal drought intensity (PDRY)] and two functional traits related to foliar and wood economics spectra [leaf mass per area (LMA); wood specific gravity (WSG)]. We show that, even though vmax is highly variable within sites, it follows a negative trend in response to increasing MAP and PDRY across sites. LMA and WSG exerted little effect on vmax and VPDmax, suggesting that these widely used functional traits provide limited explanatory power of dynamic plant responses to environmental variation within hyper-diverse forests. This study demonstrates that long-term precipitation plays an important role in the sap flux response of humid tropical forests to VPD. Our findings suggest that under higher evaporative demand, trees growing in wetter environments in humid tropical regions may be subjected to reduced water exchange with the atmosphere relative to trees growing in drier climates.
Assuntos
Transpiração Vegetal , Árvores , Secas , Florestas , Pressão de Vapor , ÁguaRESUMO
There is mounting empirical evidence that lianas affect the carbon cycle of tropical forests. However, no single vegetation model takes into account this growth form, although such efforts could greatly improve the predictions of carbon dynamics in tropical forests. In this study, we incorporated a novel mechanistic representation of lianas in a dynamic global vegetation model (the Ecosystem Demography Model). We developed a liana-specific plant functional type and mechanisms representing liana-tree interactions (such as light competition, liana-specific allometries, and attachment to host trees) and parameterized them according to a comprehensive literature meta-analysis. We tested the model for an old-growth forest (Paracou, French Guiana) and a secondary forest (Gigante Peninsula, Panama). The resulting model simulations captured many features of the two forests characterized by different levels of liana infestation as revealed by a systematic comparison of the model outputs with empirical data, including local census data from forest inventories, eddy flux tower data, and terrestrial laser scanner-derived forest vertical structure. The inclusion of lianas in the simulations reduced the secondary forest net productivity by up to 0.46 tC ha-1 year-1 , which corresponds to a limited relative reduction of 2.6% in comparison with a reference simulation without lianas. However, this resulted in significantly reduced accumulated above-ground biomass after 70 years of regrowth by up to 20 tC /ha (19% of the reference simulation). Ultimately, the simulated negative impact of lianas on the total biomass was almost completely cancelled out when the forest reached an old-growth successional stage. Our findings suggest that lianas negatively influence the forest potential carbon sink strength, especially for young, disturbed, liana-rich sites. In light of the critical role that lianas play in the profound changes currently experienced by tropical forests, this new model provides a robust numerical tool to forecast the impact of lianas on tropical forest carbon sinks.
Assuntos
Ecossistema , Clima Tropical , Ciclo do Carbono , Demografia , Florestas , Panamá , ÁrvoresRESUMO
Understanding how climate alters plant-soil water dynamics, and its impact on physiological functions, is critical to improved predictions of vegetation responses to climate change. Here we analyzed how belowground interactions for water shift under warming and drought, and associated impacts on plant functions. In a semi-arid woodland, adult trees (piñon and juniper) and perennial grasses (blue grama) were exposed to warming and precipitation reduction. After 6 years of continuous treatment exposure, soil and plant water isotopic composition was measured to assess plant water uptake depths and community-level water source partitioning. Warming and drought modified plant water uptake depths. Under warming, contrasting changes in water sources between grasses and trees reduced belowground water source partitioning, resulting in higher interspecific competition for water. Under drought, shifts in trees and grass water sources to deeper soil layers resulted in the maintenance of the naturally occurring water source partitioning among species. Trees showed higher water stress, and reduced water use and photosynthesis in response to warming and drought. This case study demonstrates that neighboring plants shift their competitive interactions for water under prolonged warming and drought, but regardless of whether changes in moisture sources will result in increased competition among species or maintained partitioning of water resources, these competitive adaptations may easily be overridden by climate extremes.