www.nature.com/scientificreports OPEN Soil water content effects on net ecosystem CO exchange and actual evapotranspiration in a Received: 10 April 2017 Mediterranean semiarid savanna of Accepted: 23 May 2018 Published: xx xx xxxx Central Chile 1,2 3 1 4,5 Francisco J. Meza , Carlo Montes , Felipe Bravo-Martínez , Penélope Serrano-Ortiz & 5,6 Andrew S. Kowalski Biosphere-atmosphere water and carbon fluxes depend on ecosystem structure, and their magnitudes and seasonal behavior are driven by environmental and biological factors. We studied the seasonal behavior of net ecosystem CO exchange (NEE), Gross Primary Productivity (GPP), Ecosystem Respiration (RE), and actual evapotranspiration (ETa) obtained by eddy covariance measurements during two years in a Mediterranean Acacia savanna ecosystem (Acacia caven) in Central Chile. The −2 −2 annual carbon balance was −53 g C m in 2011 and −111 g C m in 2012, showing that the ecosystem acts as a net sink of CO , notwithstanding water limitations on photosynthesis observed in this particularly dry period. Total annual ETa was of 128 mm in 2011 and 139 mm in 2012. Both NEE and ETa exhibited strong seasonality with peak values recorded in the winter season (July to September), as a result of ecosystem phenology, soil water content and rainfall occurrence. Consequently, the maximum carbon assimilation rate occurred in wintertime. Results show that soil water content is a major driver of GPP and RE, defining their seasonal patterns and the annual carbon assimilation capacity of the ecosystem, and also modulating the effect that solar radiation and air temperature have on NEE components at shorter time scales. Arid and semiarid regions of the world cover more than 40% of the total land surface, are characterized by low precipitation amounts and high evaporation rates, are exposed to high hydroclimatic variability and usually exhibit relatively low soil fertility, factors that condition biomass accumulation. These regions are of particular interest for the global carbon budget as it has been estimated that their ecosystems are responsible for up to 20% of terrestrial net primary productivity . Recent studies suggest that current global trends in the carbon sink of the biosphere are dominated by sem- iarid ecosystems . In these regions, NEE is controlled by precipitation, and droughts are identified as major ele- 3,4 ments that reduce GPP by limiting photosynthetic rates, and shortening the length of the growing season . This is the case of the exceptionally positive anomaly in global carbon uptake registered in 2011, which is mainly ascribed to increased productivity of semiarid vegetation in the Southern Hemisphere, and especially to the ee ff ct 5,6 of wet conditions in Australia during a La Niña year . However, since La Niña years are typically associated with negative rainfall anomalies in regions such as Central Chile, this major carbon uptake cannot be generalized. Precipitation is a commonly used climatic variable to explain the behavior of xeric vegetation. However, plant growth and development directly depend on soil water content (SWC) and plant water uptake. Therefore, rainfall records can be misleading as intense precipitation events do not always result in a proportional increase of SWC 1 2 Facultad de Agronomía e Ingeniería Forestal, Pontificia Universidad Católica de Chile, Santiago, Chile. Centro Interdisciplinario de Cambio Global, Pontificia Universidad Católica de Chile, Santiago, Chile. NASA Goddard Institute for Space Studies, New York City, NY, USA. Departamento de Ecología, Universidad de Granada, Granada, 5 6 Spain. Andalusian Institute for Earth System Research (CEAMA-IISTA), Granada, Spain. Departamento de Física Aplicada, Universidad de Granada, Granada, Spain. Correspondence and requests for materials should be addressed to F.J.M. (email: email@example.com) ScientiFic Repo R tS | (2018) 8:8570 | DOI:10.1038/s41598-018-26934-z 1 www.nature.com/scientificreports/ Figure 1. (a) Location and (b) topographic map (bar in m.a.s.l.) of the study area in Central Chile, including the monthly climograph showing the mean annual cycle of rainfall and air temperature (1994–2012). The star in (b) shows the measurements site, the thick lines are administrative divisions and the darker polygon in (b) shows the city of Santiago. We used Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) Global Digital Elevation Model Version 2 (GDEM V2) data retrieved from the online Global Data Explorer, courtesy of the NASA Land Processes Distributed Active Archive Center (LP DAAC), USGS/Earth Resources Observation and Science (EROS) Center, Sioux Falls, South Dakota, https://gdex.cr.usgs.gov/gdex/. either because of low infiltration rates (e.g. because of soil compaction) or a relatively small soil water holding capacity. Additionally, in areas with strong rainfall seasonality, the occurrence of precipitation is not always asso- ciated with the onset of the vegetation growth, exhibiting lags of several weeks and even months, and thus reduc- ing the explanatory power of precipitation records. In Mediterranean environments, the dry season is associated with warm temperature that often result in strong water stress conditions for plants due to the high atmospheric evaporative demand and restricted soil water content . As a direct consequence of climate change, it is expected that Mediterranean-type regions will face a 8 9 systematic decrease in rainfall amounts and a warming-enhanced evapotranspiration , with a resulting decrease 10,11 12 in soil water content that will ae ff ct carbon uptake . Furthermore, other studies indicate that aridity overrides the climatological sensitivity of NEE to temperature in dry areas. This could be the case of semiarid Chile, where increasing dryness in spring and summer has been reported during recent decades , which are consistent with 14 15 climate change projections , suggesting high vulnerability of NEE as climate becomes drier and warmer . 16–18 7,19 Multiple studies have assessed NEE patterns in Mediterranean ecosystems of Europe and California , but Mediterranean ecosystems in Latin America such as those found in Central Chile remain poorly studied. Despite similarities in environmental conditions and vegetation traits, fluctuations in total amount and timing of winter/spring precipitation, the literature points to significant among-site differences in annual NEE, mainly 20,21 due to fluctuations in winter/spring precipitation, both in terms of total amount and timing . Ecosystems respond directly to shortages in water availability by reducing their photosynthetic capacity and respiration, and by showing carry-over effects associated with disturbance in stored soil moisture, organic matter and nutrients, which can lead to an increase in plant mortality . In the same way, reduction of annual carbon uptake by up to 45% has been quantified in Mediterranean cork oak woodlands as a result of a lower photosynthetic rate and leaf 11,23 3 area duration . The same effect was observed on woody species , together with a change from net carbon sink into net source in a grassland, likely due to lower soil exploration by roots. Similarly, a shift to net carbon source was reported in a grassland in California due to a shorter growing season . Given the relevance of carbon uptake in dry ecosystems, the strong seasonality of precipitation, and the increasing occurrence of drought spells in Central Chile, our objective was to assess the seasonal variation of NEE and its components (GPP, RE) as well as ETa in a Mediterranean Acacia caven savanna and their relations to soil water content. The study area is characterized by strong seasonality with more than 85% of total precipitation concentrated in winter, a fact that was particularly evident during the observation period where only few rains fell in autumn. Our underlying hypothesis is that soil water content is an important driver of NEE during winter periods, modifying the response of GPP and RE to other environmental factors such as photosynthetically active radiation (PAR) and temperature. As a consequence, the highest photosynthetic activity of this ecosystem occurs between July and August. This study uses eddy covariance time series of CO and water vapor obtained from a Mediterranean Acacia caven savanna (Fig. 1) during two consecutive hydrological years. ScientiFic Repo R tS | (2018) 8:8570 | DOI:10.1038/s41598-018-26934-z 2 www.nature.com/scientificreports/ Figure 2. (a) Time series of daily values of GPP and EVI for the study site (b) Relationship between EVI and GPP values and (c) Boxplots of annual distribution of the Enhanced Vegetation Index (EVI), gray bars represent the relative distribution of annual precipitation during the year (average monthly values divided by total sum of precipitation). Central mark shows the median and the edges are the 25th and 75th percentiles; dashed lines extend to the most extreme values not considered outliers, and outliers are plotted individually (plus sign). Results Environmental conditions during the experiment. The Enhanced Vegetation Index (EVI) was used as a proxy to relate patterns of ecosystem fluxes to the vegetation cycle. EVI and GPP showed similar temporal variability (Fig. 2a EVI peaked close to 0.4 in August-September). The linear relationship between EVI and GPP (Fig. 2b) allow us to detect ecosystem photosynthetic activity especially when herbaceous species are present, accounting for 75% of GPP variance. The annual cycle of monthly mean EVI values (Fig. 2c) shows strong seasonal variability, closely associated with the seasonal precipitation distribution. Minimum EVI values occur in autumn (April, 0.15 ± 0.01) and max- imum values in late winter (August, 0.40 ± 0.08), before and aer t ft he rainy season, respectively. We observe that the highest photosynthetic activity of the system occurs between late winter and early spring (July, August and September), while the rest of the year is characterized by steady and small NEE values that are interrupted only by short pulses of rain triggering peaks in ecosystem respiration. e s Th easonal course in daily rainfall and soil water content (SWC) (Fig. 3a) shows typical Mediterranean-type features, with most events occurring during winter and early spring (from June/July to Aug/Sep.), coinciding with minimum values of daily mean air (T Fig. 3b) and soil temperature (T , data not shown) of around 5 °C. Although a s SWC responds rapidly to precipitation (Fig. 3a), some rainfall events are of such intensity that maximum soil water holding capacity is rapidly reached and excess water is lost due to runo. Th ff e dry season extends from sum- mer to early autumn, and is characterized by very few and short rain events (Fig. 3a), high average T and T (up to a s 25 °C and 37 °C, respectively), and vapor pressure deficit (VPD) (up to 2–2.5 kPa). For clear-sky days, maximum −2 −1 −2 PAR and average net radiation (R ) values in summer were near 2000 µmol m s and 160 W m , respectively, −2 −1 −2 versus 1000 µmol m s and 25 W m in winter. Seasonal patterns of CO , water vapor fluxes and main drivers. Daily ETa (Fig. 3b) showed an increase from the end of the warm season to a maximum in winter and early spring (15-day moving average of −1 −1 3 −3 0.1 mm d and 1 mm d , respectively) when SWC was between 0.2–0.3 m m . Conversely, reference evapo- transpiration (ET which corresponds to the atmospheric water demand acting on a short grass actively growing that is regarded as a reference vegetation surface) reaches its maximum during the warm season as a response to available energy and atmospheric demand. This pattern is only interrupted by soil water recharge after rainfall events. e a Th nnual course of NEE (Fig. 3c) shows that the Acacia ecosystem varies from a net sink to a carbon source depending on the time of year, with a lower/higher magnitude during the warm/cold season. Greater fluctuations are observed from 2010 through 2011, while a steadier carbon sink behavior remained throughout 2012. Between summer and autumn (January to June), when Acacia almost exclusively composes vegetation, NEE rates were −2 −1 similar (varying from −0.7 to 0.7 g C m d ). In 2011, fluctuations were mainly associated with sporadic rain events (Fig. 3a). However, the effect of rainfalls on NEE depends on the season: in autumn and winter, rainfalls were associated with an increase in carbon uptake, while in the warm/dry season they generated net carbon losses. −2 −1 A maximum carbon assimilation rate (-NEE) of around 2.5 g C m d was observed during winter and early spring (July to September). As shown in Fig. 3c, the first rain events in autumn/winter were followed by positive values of NEE, showing a similar response of GPP and RE to increasing SWC as in precedent months, likely as a response to temperatures allowing microbial respiratory activity to exceed photosynthesis (Fig. 3d). Following previous studies , the effect of SWC on NEE was assessed by analyzing the functional relationship between ETa, GPP and RE. This was applied to 10-day accumulated values, separately for the dry (October-May) and wet (June–September) season. A moderate slope and a positive y-axis intercept for both GPP and RE charac- terizes the relationship between ET and carbon fluxes during the dry period (Fig. 4a). On the other hand, a higher ScientiFic Repo R tS | (2018) 8:8570 | DOI:10.1038/s41598-018-26934-z 3 www.nature.com/scientificreports/ Figure 3. Time series of daily values of (a) rainfall and volumetric soil water content (SWC), (b) actual (ETa) and reference (ET ) evapotranspiration, and air temperature (c) net ecosystem exchange (NEE), (d) gross primary productivity (GPP) and ecosystem respiration (RE). Figure 4. Relationship between 10-day accumulated ETa and carbon fluxes (GPP, RE) during (a) the dry and (b) the wet period of the year. R and p denote the Pearson’s correlation coefficient and p-value of statistical significance, respectively. slope was obtained between ETa and carbon fluxes during the wet period (Fig. 4b), and, unlike the dry period, a negative intercept was found for the ETa:GPP relationship. To further explore how NEE components and environmental variables are related on a seasonal basis, a Principal Component Analysis (PCA) was applied to centered and variance-scaled daily GPP, RE, ETa, T , T , a s ScientiFic Repo R tS | (2018) 8:8570 | DOI:10.1038/s41598-018-26934-z 4 www.nature.com/scientificreports/ Figure 5. Biplot displaying the two first principal components (PC 1 and PC 2) and their corresponding explained variance (in %) derived from a PCA applied to daily gross primary productivity (GPP), ecosystem respiration (RE), actual evapotranspiration (ETa), air temperature (T ), soil temperature (T ), net radiation (R ), a s n photosynthetically active radiation (PAR), soil water content (SWC) and vapor pressure deficit (VPD) (n = 683). Vectors denote the contribution of a single variable to the PC. Color bar denotes the day of the year (DOY) for data pairs. R , PAR, SWC and VPD. Results are displayed in Fig. 5 for the two first principal components (PC; eigenvalues >1, 76.9% of total variance) and for the remaining PC in Table S1. The biplot shows that PC 1 separates carbon ux fl es, ETa and SWC from T , T , R , PAR and VPD, which in turn are poorly correlated (angle between vectors a s n near 90°). On the other hand, the direct relationship between SWC and carbon and water vapor fluxes at the actual time scale and analyzed period is supported by their high correlation (angle close to zero). In addition, the sequential ordering according to the day of the year (DOY) of data pairs accounts for strong seasonality, where dry/hot and wet/cold periods are clearly differentiated. Although explaining less variance (17.3%), PC 2 describes variations in GPP, RE and ETa, given their higher correlation with this component (Table S1). The above results are summarized as the average monthly annual cycle of NEE, GPP and RE (Fig. 6). NEE and its components show little variability from January to April, aer ft which the opposing relationship between GPP and RE generates a maximum of NEE in June following the onset of the rainy season. e Th average maximum assimilation period during winter and early spring (June–September) begins aer t ft he divergence between GPP and RE rates to reach their maxima in August, to then remain close to zero in summer. e eff Th ect of SWC on RE was further evaluated looking at the specific period of winter early spring that pre- sents the maximum NEE magnitude (July–September). Table 1 shows the parameters of the Van’t Hoff model t fi ted to the observed nighttime RE values using linear regression. We used bin widths of 2 K and a minimum bin size of n = 10 to calculate mean values of RE and to reduce undesired effects of heteroscedasticity. In winter (July to September) SWC has an effect on ecosystem respiration (RE), increasing its value at the reference temperature (R ; T = 283.16 K), and almost doubling it when relative soil water content is above eco, ref ref 0.6 in comparison to the lower end. In addition, the rate at which temperature modifies ecosystem respiration (parameter B) increases as SWC increases, this effect is particularly evident when changing from low to medium relative soil water content (RWC, equation 2) values. Because of the strong seasonality in rainfall and SWC, most of the observations in the other periods (spring, summer and autumn) fall within a RWC class of low being dif- ficult to assess the effects of soil water content on RE, only autumn data allow us to compare the parameters of the Van’t Hoff model between RWC classes low (R = 0.38 and B = −0.004) and medium (R = 0.52 and eco,ref eco,ref B = 0.032). e p Th arameters of Table 1 were then used to calculate RE during daytime and combined with NEE records to calculate GPP. We plotted GPP values against RWC (Fig. 7), with values grouped in 0.05 (dimensionless) incre- ments to achieve a relatively homogeneous bin size. e r Th elationship between GPP as a function of RWC indicates that soil moisture has a great effect on the ability of the ecosystem to capture carbon. Discussion −2 −1 −2 −1 Integrating over the year, we found annual NEE values of −53 g C m y and −111 g C m y and ETa values of 128 and 139 mm for 2011 and 2012, respectively. The study was carried out during a severe drought (2010– 2012). In this period, annual rainfall was found to be below normal: 9% deficit in 2010 (223 mm), 49% deficit in 2011 (111 mm), and 36% deficit in 2012 (155 mm), the years represent the initial period of a severe drought that ScientiFic Repo R tS | (2018) 8:8570 | DOI:10.1038/s41598-018-26934-z 5 www.nature.com/scientificreports/ −2 −1 2 −2 −1 RWC class R (μmol m s ) B R RMSE (μmol m s ) eco,ref Low 0.832 0.028 0.46 0.123 Medium 1.073 0.081 0.88 0.156 High 1.340 0.100 0.87 0.163 Table 1. Parameters of the Van’t Hoff model representing the influence of temperature on RE for different soil water content categories for the period from June to September. Figure 6. Annual cycle for monthly-averaged daily values of NEE, GPP and RE. Shaded intervals represent the standard deviation. Figure 7. Relationship between instantaneous GPP (expressed as mass of carbon per unit area per day) and Relative Soil Water content (RWC). ScientiFic Repo R tS | (2018) 8:8570 | DOI:10.1038/s41598-018-26934-z 6 www.nature.com/scientificreports/ extends to 2013 and beyond , comparable only to that registered between 1967 and 1969, where three consec- utive years showed significant deficits. This feature suggests that the Acacia ecosystem acts as a net carbon sink, with large seasonal variability. Since eddy covariance (EC) measurements were taken during a particularly dry period, greater carbon uptake should be expected for close-to-normal or wet years, as suggested by the linear relationship between GPP and ETa that in turn responds to precipitation and SWC. Notwithstanding this, the annual carbon fluxes are within the range observed for similar water-limited Mediterranean ecosystems. For −2 −1 −2 −1 −2 −1 instance, decreases in uptake from 140 g C m y to 28 g C m y (2003–2006) and from 388 g C m y to −2 −1 214 g C m y (2011–2012) were found for a similar Mediterranean oak woodland in Portugal due to the effect 3,11 of drought . A noticeable feature in measured carbon fluxes is the high seasonality of their annual course. Although this can be immediately related to the annual cycle in precipitation and the associated suppression of soil and plant metabolism during the warm season, it can be also explained by the phenological lag between the herbaceous cover and the Acacia activity. Thus, the first rains in autumn give way to germination and subsequent growth of grasses until the beginning of the dry season when Acacia begins its growth . This seasonal behavior gen- erates an increase in biomass that is coupled to the rainy season in winter, which has been also described for similar Mediterranean ecosystems . As a consequence, annual carbon flux variations are dominated by rainfall during winter and early spring months (which maintain relatively high SWC levels), when maximum rates of NEE, GPP, RE and ETa occur, due to the activity of the well-adapted understory developing under unshaded and nitrogen-rich conditions prior to the leafing-out of Acacia . In addition, the net increase in late winter produc- tivity by the herbaceous stratum might be enhanced by the increasing soil water retention capacity and nutrient 27,30 content due to the decomposition of recently shed Acacia leaves during the period of higher water availability . Steady rates of carbon assimilation were observed during the dry season when Acacia photosynthesis is dom- inant and SWC restricts soil and plant respiration while grasses are senescent. This dry-season sink behavior coupled with high atmospheric evaporative demand might be possible thanks to deep soil exploration by the root system of Acacia. In addition, the absence of precipitation in summer suggests a possible water supply to Acacia 31,32 trees originating from sources such as groundwater and its redistribution by roots , which would explain per- sistent photosynthetic rates during the dry season. However, further research is needed to validate this hypoth- esis. Furthermore, a moderate slope characterizes the relationship between ET and carbon fluxes during the dry period (Fig. 4a), indicating a small change in GPP/RE as water availability increases, and suggesting a weak dependence of the seasonal ecosystem functionality on sporadic rains. This can also be inferred from the positive y-axis intercept for both GPP and RE in Fig. 4a, indicating photosynthetic and respiratory capacity when ETa is close to zero. In addition, the slightly higher slope and y-axis intercept of ETa:GPP than ETa:RE would account for a weak and almost invariant carbon sink capacity during this period. It is worth mentioning that although the relationship between ETa and RE is not statistically significant (p = 0.18), there is a trend that should be addressed in future studies that examine, for example, a larger data sample. During the onset of winter rainfalls both GPP and RE experienced substantial increases as a result of SWC increase, contributing to the net carbon uptake and sink capacity in spite of the lower temperature and available energy for photosynthesis, as reported in the literature . Comparing the 2011 and 2012 years, both GPP and RE are higher in the cold season of 2012 which had more rainfall. As shown in Fig. 4b, the existence of an ETa value below which RE exceeds GPP could be related to the exist- ence of a threshold ETa and therefore SWC necessary for photosynthesis activation aer dr ft y periods or for seed −1 germination in the transition from dry/warm to wet/cold season, whose magnitude (ETa ~1 mm 10 d ) might be related to direct soil and/or canopy evaporation when GPP approaches zero. This is relevant considering the actual and projected aridity over the region and the associated probability of a reduced carbon uptake capacity or even conversion to a carbon source , represented by the difference between the two lines in Fig. 4b, as precipi- tation decreases in upcoming decades. The latter represents a major issue also because of the possible southward 34 35 expansion of the Mediterranean-type climate and the invasive nature of Acacia sp .Furthermore, the threshold in available water (the intersection between the lines in Fig. 4b) could be reached later at the beginning of the wet season in a drier climate, shortening the period of higher net carbon assimilation. The slope obtained between ETa and carbon fluxes during the wet period suggests a greater sensitivity to water availability when changes in stomatal conductance or biomass can drive abrupt changes in carbon and water fluxes, as observed at the annual scale . In turn, the steep slope for the case of ETa:GPP accounts for the dependence of carbon absorption capacity on water available to the ecosystem. Unlike the dry period, a negative intercept in the relationship ETa:GPP indi- cates that there is a value of ETa below which little or no growth can be observed. Similar to the dry season, the y-axis intercept in ETa:RE would account for a net carbon loss when ETa approaches zero. Moreover, considering −1 the intercept between the two lines as the point from which GPP is higher than RE, an ETa of 2.5 mm 10 d can be estimated as a threshold above which the ecosystem becomes a carbon sink, whose magnitude depends on available water. e PCA a Th pplied to carbon and water fluxes and their environmental drivers (Fig. 5, Table S1) suggests a sea- sonally varying relationship, with SWC as the main factor controlling the seasonal behavior of carbon and water vapor ux fl es. As such, the great importance of soil water content as a factor controlling RE has been reported in dry lands , when warm-season rain events activate autotrophic and heterotrophic respiration causing short −2 −1 17 RE pulses and carbon release of around 1–1.5 g C m d , whose magnitude depends on rainfall amounts . According to the literature , carbon emissions associated with short precipitation pulses are important in arid regions with sparse vegetation given the low water infiltration of bare soil versus below-canopy soil. In addition, authors estimated at 59 mm the necessary spring rainfall to produce a net CO uptake in a semiarid shrub in Arizona, and that short rainfall pulses triggered 40% of CO emissions from grassland in southeastern Spain for the period 2009–2013 . Moreover, RE pulses were described as resulting in a change from a CO sink to a net 40 41 source for a steppe ecosystem and in a Mediterranean grassland . Also, soil moisture was reported as the main ScientiFic Repo R tS | (2018) 8:8570 | DOI:10.1038/s41598-018-26934-z 7 www.nature.com/scientificreports/ driving factor of growing season NEE in a sagebrush-steppe ecosystem, and that VPD is important when the ecosystem is not limited by SWC . In this sense, the high evaporative demand and dryness in summer prevent the development of the herbaceous stratum, which could compensate this carbon loss, so that this effect should be limited to the duration of surface soil wetness from field capacity to a minimum value (i.e., wilting point). Regarding the warm season in 2012, no precipitation events were observed, so that the steady negative NEE accounts for the moderate CO sink capacity of the ecosystem despite the seasonal dryness and high VPD, which is in agreement with the well-known adaptations to Acacia in terms of deep rooting and soil exploration, while herbaceous species are senescent . The relationship between SWC and RE/GPP suggests an important modulation of the effect that temper - ature and PAR can have on those variables respectively (and therefore on NEE) which might compromise the sink capacity of this system under future climate scenarios. Climate change projections from general circulation models for the region show higher temperatures and decreasing rainfall amounts for this region , which in turn create a scenario where soil water content diminishes (as effective precipitation is reduced and temperature causes higher evaporation rates). Moreover, even for the period of greater activity we can expect that both RE and GPP will be reduced, creating a scenario where NEE changes could substantially compromise the sink behavior of this system. Conclusions Eddy covariance CO fluxes were measured in a Mediterranean Acacia caven savanna ecosystem in Central Chile to examine their seasonal patterns and main drivers during 2011 and 2012. Results showed annual carbon uptake −2 −2 of the Acacia savanna of 53 g C m in 2011 and 111 g C m in 2012 with strong variability due mainly to differ - ences in precipitation and therefore in SWC. Although carbon and water ux fl es correlate with environmental variables during the year, the annual NEE pattern was strongly influenced by SWC and the ecosystem phenology. Results from this study provide an overview about the seasonal behavior of carbon exchange by the Acacia savanna ecosystem of Central Chile during years of water limitation and in face of climate change. Greater sen- sitivity to SWC was observed in GPP and RE during the wet/cold season when grasses are present, versus the slow, steady dry/warm season, which nonetheless strongly determines the carbon sink capacity. This is in direct relation with winter precipitations that determine the herbaceous growth as the SWC determines its develop- ment. The interannual carbon balance indicates that the sink capacity of the Acacia savanna is mainly driven by its productivity during winter, which is thereby relevant from a climate change perspective considering the pro- jected decrease in precipitation over the region. Finally, longer measurements of CO ux fl es are needed in order to improve knowledge about the Mediterranean-type ecosystems of Central Chile and to better understand the sink/source cycles of the Acacia savanna under contrasting precipitation regimes (e.g. El Niño/La Niña years). Materials and Methods Study site description. Measurements were conducted on a 24-ha Mediterranean savanna located at the foot of the Andes in Central Chile (33°02′S 70°44′W, elevation 660 m.a.s.l., terrain slope 3%; Fig. 1). In this region, the destruction and replacement of the native evergreen sclerophyllous forests by croplands and animal hus- bandry have given way to an early successional plant community, locally known as espinal, which covers about 2 million ha of the Mediterranean perarid areas near 30°S (100–150 mm of annual rainfall) to perhumid zones near 27,44 36°S (700–1200 mm of annual rainfall) . This ecosystem represents the most widespread agroforestry system of 29,46 Central Chile, coexisting with herbaceous and sparse woody plants, all exposed to grazing pressure . Vegetation is composed of both woody and herbaceous species dominated by the thorny leguminous tree Acacia caven (Mol.) (Fabaceae; hereafter Acacia). Acacia is a woody plant of about 3 m height form- ing a sparse-open canopy (30% crown cover) arboreal stratum . As a leguminous plant, it is a highly efficient nitrogen-fixing shrub whose semi-deciduous foliage allows frequent nitrogen recycling from lower to upper soil horizons . During the rainy season (austral winter), this phenomenon promotes the development of an herba- ceous stratum composed by annual C3/C4 herbs and grasses including Anoda sp, Erodium moschatum, Trifolium sp, Oxalis sp, Urtica urens and Helenium aromaticum.Other woody plants (<10% of cover) include Porlieria chilensis, Prosopis chilensis or Proustia cuneifolia . The climate of the study region is semiarid Mediterranean-type, with monthly average temperature ranging from 9.4 °C to 21.7 °C (15.8 °C annual mean). Precipitation shows an average value of 234 mm (1994–2012; data presented below) and is highly concentrated in the cold season. Interannual variations are largely regulated by the El Niño-Southern Oscillation, which induces positive/negative precipitation anomalies during the El Niño/ La Niña phases. The soil of the site is classified as a shallow Vertic Calcixerolls (Mollisol), with low organic matter content (2.9%), clay loam texture (30:32:38 sand:silt:clay) and a C:N ratio of 12. Data collection. Eddy covariance measurements and data processing. Turbulent CO and energy fluxes (latent and sensible heat) were measured from October-2010 to October-2012 by an eddy covariance (EC) tower installed 6 m above the soil surface on the Acacia site, consisting of a CSAT3 (Campbell Scientific; Logan, UT, USA; hereaer CS ft I) sonic anemometer and an Li-7500(Li-Cor, Lincoln, NE, USA) open-path infrared gas ana- lyzer (IRGA). Data were sampled at 10 Hz and stored half-hourly as average net fluxes. Data post-processing indicated that the site has a 700-m fetch, and footprint analyses showed that fluxes originated exclusively from the espinal. The Flux-Source Area model allowed verifying that for periods of relative stability (friction velocity −1 −1 >0.1 m s and <0.2 m s , negative sensible heat flux), when measured fluxes are generated farthest from the tower, they originated from within the fetch (data not shown). For such stable conditions, a maximum distance to 80% source area isopleths of 300 m was estimated. 48,49 CO and energy fluxes data were corrected for density perturbations and coordinate rotation and checked 50 51 for quality-control .Data generated by dirty lenses were discarded aer a ft filter based on IRGA diagnostics , and ScientiFic Repo R tS | (2018) 8:8570 | DOI:10.1038/s41598-018-26934-z 8 www.nature.com/scientificreports/ −1 52 averaging periods with insufficient turbulence (friction velocity <0.1 m s ) were rejected . Gaps in the database (~20%) were mainly due to dirty lenses (53%), power failures (28%), anomalous peaks due to morning flush and weak turbulence (19%). An energy balance closure of 68% was obtained, which falls within the range reported in the literature . CO and water vapor balances were calculated after gap filling. The gap-filling method delivers not only an estimate of missing values but also its standard error (σ ). To estimate the uncertainty in annual values, we assumed independence in the errors of half-hourly gap-filled estimates and a normal distribution with a mean equal to the corresponding gap-filling estimate and a standard deviation of σ . From 50 Monte Carlo simulations of the observation period, with missing data replaced by an estimate randomly selected from its corresponding Gaussian distribution, 50 estimates of annual NEE and ET and their corresponding standard deviation were obtained. Gaps in half-hourly net CO and water vapor u fl xes were filled and then NEE was partitioned into GPP and RE by using the “night-time data-base estimate” method . Assuming GPP = 0 at night, RE is obtained as an expo- nential relationship between RE and temperature following the Van´t Hoff equation : BT () ‐Tref NEER == ER e nightnight ecor , ef (1) where R is RE at reference temperature (T , set to 283.16 K) considered as the base respiration, B is a site eco,ref ref specific parameter, T the observed soil temperature (K),. The resulting exponential relationship is extrapolated to daytime and then GPP is obtained as the difference between RE and measured NEE. We used bin widths of 2 K and a minimum bin size of n = 10 to calculate mean values of RE and to reduce undesired effects of heteroscedasticity. Daily and ten day estimates of NEE, GPP, RE and ETa were obtained as the sum of the half hour values over the day or ten days, respectively. Annual estimates correspond to the sum of daily values for the entire year from 2010 to 2012. Environmental variables. Incident PAR and R were obtained from a quantum sensor (LI190SB; Li-Cor) and a NR-Lite sensor (Kipp and Zonen, Delft, Holland), respectively, both mounted at 4 m height. The soil heat flux was measured for both sunlit and (canopy) shaded soil using 4 plates installed at 8 cm depth (HFP01, Hukseflux, Delft, Holland), and soil temperature T by two integrated soil thermocouples (TCAV, CSI) at 2 and 6 cm depth. PAR, R soil heat and temperature data were registered every 10 s, and half-hour averages were stored in a data n, logger (CR5000, CSI). Volumetric SWC was measured daily via time domain reflectometry (TDR) (CS616, CSI) at a depth of 4 cm. TDR sensors were calibrated before starting the experiment using soil samples of the location under study with known water content. Along with these data, daily rainfall and screen air temperature data from 1994 through 2012 were obtained from a weather station located 4.3 km southeast of the study site (33°04′33′′S 70°46′07”W), and long-term rainfall data (1950–2013) from a weather station located in Santiago (33°26′42′′S 70°40′58′′W; Fig. 1b) were used to contextualize drought conditions during the experiment. Additionally, the above-mentioned data were used to compute hourly reference evapotranspiration (ET ) according to the FAO-56 methodology . To study the dependence of GPP and RE on SWC we calculated the relative soil water content (RWC) of the most active period (July to September) and classify the observations into the classes: Low (0 < RWC < 0.3); Medium (0.3 < RWC < 0.6) and High (RWC > 0.6). The equation used to calculate RWC was: SWCS − WC min RWC = SWCS − WC (2) maxmin where SWC (0.04) and SWC (0.49) correspond to the minimum and maximum volumetric SWC values, min max respectively. Enhanced Vegetation Index. Vegetation indices derived from satellite imagery are commonly used to assess ecosystem biomass. In this work, the Enhanced Vegetation Index (EVI) from the National Aeronautics and Space Administration’s (NASA) Moderate Resolution Imaging Spectroradiometer (MODIS) was used as an estimate of the area-average active aboveground biomass so as to investigate its relationship with carbon fluxes and the anomalies induced by drought conditions. Sixteen-day composites at 250 m spatial resolution from 2000 to 2012 (MOD13Q1 product) were obtained for the study site by averaging EVI values for the 13 pixels (3.25 ×3.25 km) around the EC tower with the same land use. Differences in annual and monthly EVI mean values were assessed using a one-way ANOVA multiple comparison test (Fisher’s LSD; p < 0.05). References 1. Jia, B. et al. Biophysical controls on net ecosystem CO exchange over a semiarid shrubland in northwest China. Biogeosciences 11, 4679–4693 (2014). 2. Ahlström, A. et al. The dominant role of semi-arid ecosystems in the trend and variability of the land CO sink. Science3 48, 895–899 (2015). 3. Pereira, J. S. et al. Net ecosystem carbon exchange in three contrasting Mediterraneanecosystems - the effect of drought. Biogeosciences 4, 791–802 (2007). 4. Rambal, S. et al. How drought severity constrains gross primary production(GPP)and its partitioning among carbon pools in a Quercus ilex coppice? Biogesciences 11, 6855–6869 (2014). 5. Le Quéré, C. et al. Global carbon budget 2013. Earth Syst. Sci. Data. 6, 235–263 (2014). 6. Poulter, B. et al. Contribution of semi-arid ecosystems to interannualvariability of the global carbon cycle. Nature 509, 600–603 (2014). ScientiFic Repo R tS | (2018) 8:8570 | DOI:10.1038/s41598-018-26934-z 9 www.nature.com/scientificreports/ 7. Ma, S., Baldocchi, D. D., Xu, L. & Hehn, T. Inter-annual variability in carbon dioxide exchange of an oak/grass savanna and open grass land in California. Agric. For. Meteorol. 147, 157–171 (2007). 8. Polade, S. D., Pierce, D. W., Cayan, D. R., Gershunov, A. & Dettinger, M. D.The key role of dry days in changing regional climate and precipitation regimes. Sci. Rep. 4, https://doi.org/10.1038/srep04364 (2014). 9. Mariotti, A. et al. Mediterranean water cycle changes: transition to drier 21st century conditions in observations and CMIP3 simulations. Environ. Res. Lett. 3, 044001 (2008). 10. Ciais, P. et al. Europe-wide reduction in primary productivity caused by the heat and drought in 2003. Nature4 37, 529–533 (2005). 11. Costa-e-Silva, F. et al. Effects of an extremely dry winter on net ecosystem carbon exchange and tree phenology at a cork oak woodland. Agric. Forest Meteorol. 204, 48–57 (2015). 12. Yi, C. et al. Climate control of terrestrial carbon exchange across biomes and continents. Environ. Res. Lett. 5, 034007 (2010). 13. Meza, F. J. Recent trends and ENSO influence on droughts in Northern Chile: An application of the Standardized Precipitation Evapotranspiration Index. Weather Climate Extremes. 1, 51–58 (2013). 14. Fuenzalida, H. et al. Study on Climate Variability for Chile during the 21st century. Technical Report prepared for the National Environmental Committee (in Spanish). Available on-line at http://www.dgf.uchile.cl/PRECIS (2007). 15. Van Oijen, M. et al. Impact of droughts on the carbon cycle in European vegetation:a probabilistic risk analysis using six vegetation models. Biosgeosciences 11, 6357–6375 (2014). 16. Rambal, S., Joffre, R., Ourcival, J. M., Cavender-Bares, J. & Rocheteau, A. e Th growth respiration component in eddy CO ux fl from a Quercus ilex mediterranean forest. Global Change Biol. 10, 1460–1469 (2004). 17. Jongen, M., Santos Pereira, J., Igreja Aires, L. M. & Pio, C. A. The effects of drought and timing of precipitation on the inter-annual variation in ecosystem-atmosphere exchange in a Mediterranean grassland. Agric. Forest Meteorol. 151, 595–606 (2011). 18. Serrano-Ortiz, P. et al. Seasonality of net carbon exchanges of Mediterranean ecosystemsacross an altitudinal gradient. J. Arid Environ. 115, 1–9 (2015). 19. Vaz, M. et al. Drought-induced photosynthetic inhibition and autumn recovery in twoMediterranean oak species (Quercus ilex and Quercus suber). Tree Physiol. 30, 946–956 (2010). 20. Serrano-Ortiz, P. et al. Variations in daytime net carbon and water exchange in a montane shrubland ecosystem in southeast Spain. Photosynthetica. 45, 30–35 (2007). 21. Scott, R. L., Serrano-Ortiz, P., Domingo, F., Hamerlynck, E. P. & Kowalski, A. S. Commonalities of carbon dioxide exchange in semiarid regions with monsoon and Mediterranean climates. J. Arid Environ. 84, 71–79 (2012). 22. van der Molen, M. K. et al. Drought and ecosystem carbon cycle. Agric. Forest Meteorol. 151, 765–773 (2011). 23. Piayda, A. et al. Drought impact on carbon and water cycling in a Mediterranean Quercus suber L. woodland during theextreme drought event in 2012. Biosgeosciences 11, 7159–7178 (2014). 24. Xu, L. K. & Baldocchi, D. D. Seasonal variation in carbon dioxide exchange over a Mediterranean annual grassland in California. Agric. Forest Meteorol. 123, 79–96 (2004). 25. Scott, R. L., Biederman, J. A., Hamerlynck, E. P. & Barron-Gao ff rd, G. A. The carbon balance pivot point of southwestern U.S. semiarid ecosystems: Insights from the 21st century drought. J. Geophys. Res. Biogeosci. 120, 2612–2624, https://doi.org/10.1002/2015JG003181 (2015). 26. Boisier, J. P., Rondanelli, R. Garreaud, R., & Muñoz, F. Natural and anthropogenic contributions to the Southeast Pacific precipitation decline and recent mega-drought in central Chile. Geophys. Res. Lett. 43, https://doi.org/10.1002/2015GL067265 (2016). 27. Aronson, J., Ovalle, C., Aguilera, L. & Leon, P. Phenology of an ‘immigrant’ savanna tree (Acacia caven, Leguminosae) in the Mediterranean climate zone of Chile. J. Arid Environ. 27, 55–70 (1994). 28. Wang, J. et al. Canopy and climate controls of gross primary production of Mediterranean-type deciduous and evergreen oak savannas. Agric. Forest Meteorol. 226-227, 132–147 (2016). 29. Ovalle, C., del Pozo, A., Casado, M. A. & Acosta, B. & de Miguel, J. M. Consequences of landscape heterogeneity on grassland diversityand productivity in the Espinal agroforestry system of central Chile. Landscape Ecol. 21, 585–594 (2006). 30. Joffre, R. & Rambal, S. How tree cover influences the water-balance of Mediterranean Rangelands. Ecology 74, 570–582 (1993). 31. Canadell, J. et al. Maximum rooting depth of vegetation types at the global scale. Oecologia 108, 583–595 (1996). 32. Scott, R. L. et al. Interannual and seasonal variation influxes of water and carbon dioxide from a riparian woodland ecosystem. Agric. Forest Meteorol. 122, 65–84 (2004). 33. Gerstmann, C., Miranda, M. & Condal, A. Description of space-time variability of thepotential productivity of Acacia caven espinales based on MODIS images and theEnhanced Vegetation Index (EVI). Cienc.Investig.Agrar. 37, 63–73 (2010). 34. Klausmeyer, K. R. & Shaw, M. R. Climate change, habitat loss, protected areas and the climate adaptation potential of species inMediterranean ecosystems worldwide. PLoS ONE 4, e6392 (2009). 35. Stella, J. C., Rodríguez-González, P. M., Dufour, S. & Bendix, J. Riparian vegetation research in Mediterranean-climate regions: common patterns, ecological processes, and considerations for management. Hydrobiol. 719, 291–315 (2013). 36. Domingo, F. et al. Carbon and water exchange in semiarid ecosystems in SE Spain. J. Arid Environ. 75, 1271–1281 (2011). 37. Huxman, T. E. et al. Precipitation pulses and carbon fluxes in semiarid and arid ecosystems. Oecologia 141, 254–268 (2004). 38. Emmerich, W. E. & Verdugo, C. L. Precipitation thresholds for CO uptake in grass and shrub plant communities on Walnut Gulch Experimental Watershed. Water Resour. Res. 44, W05S16, https://doi.org/10.1029/2006WR005690 (2008). 39. López-Ballesteros, A. et al. Enhancement of the net CO release of a semiarid grassland in SE Spain by rain pulses. J. Geophys. Res. Biogeosci. 121, 52–66, https://doi.org/10.1002/2015JG003091 (2015). 40. Wang, Y., Zhou, G. & Wang, Y. Environmental effects on net ecosystem CO exchange athalf-hour and month scales over Stipa krylovii steppein northern China. Agric. Forest Meteorol. 148, 714–722 (2008). 41. Flanagan, L. B., Wever, L. A. & Carlson, P. J. Seasonal andinterannual variation in carbon dioxide exchange andcarbon balance in a northern temperate grassland. Global Change Biol. 8, 599–615 (2002). 42. Kwon, H., Pendall, E., Ewers, B. E., Cleary, M. & Naithani, K. Spring drought regulates summernet ecosystem CO exchange in a sagebrush-steppe ecosystem. Agric. Forest Meteorol. 148, 381–391 (2008). 43. Hunt, J. E., Kelliher, F. M., McSeveny, T. M., Ross, D. J. & Whitehead, D. Long-term carbon exchange in a sparse, seasonally dry tussock grassland. Global Change Biol. 10, 1785–1800 (2004). 44. Ovalle, C., Aronson, J., del Pozo, A. & Avendaño, J. e Th espinal: agroforestry systems of the mediterranean - type climate region of Chile. State of the art and prospects for improvement. Agroforest. Syst. 10, 213–239 (1990). 45. Meza, F. J., Vicuna, S., Jelinek, M., Bustos, E. & Bonelli, S. Assessing water demands and coverage sensitivity to climate change in the urban and rural sectors in Central Chile. Journal of Water and Climate Change. 5(2), 192–203 (2014). 46. Van de Wouw, P., Echeverría, C., Rey-Benayas, J. M. & Holmgren, M. Persistent Acacia savannas replace Mediterranean sclerophyllous forests in South America. Forest Ecol. Manag. 262, 1100–1108 (2011). 47. Schmid, H. P. Source areas for scalars and scalar fluxes. Bound-Lay Meteorol. 67, 293–318 (2004). 48. Webb, E., Pearman, G. & Leuning, R. Correction of flux measurements for density effects due to heat and water vapour transfer. Q. J. R. Meteorol. Soc 106, 85–100 (1980). 49. Kowalski, A. S., Anthoni, P. M., Vong, R. J., Delany, A. C. & Maclean, G. D. Deployment and evaluation of a system for ground-based measurement of cloud liquid water turbulent fluxes. J. Atmos. Ocean Tech. 14, 468–479 (1997). 50. Reverter, B. R. et al. Analyzing the major drivers of NEE in a Mediterranean alpine shrubland. Biogeosciences 7, 2601–2611 (2010). ScientiFic Repo R tS | (2018) 8:8570 | DOI:10.1038/s41598-018-26934-z 10 www.nature.com/scientificreports/ 51. Serrano-Ortiz, P., Kowalski, A. S., Domingo, F., Ruiz, B. & Alados-Arboledas, L. Consequences of uncertainties in CO density for estimating net ecosystem CO exchange by open-path eddy covariance. Bound-Lay Meteorol. 126, 209–218 (2008). 52. Gu, L., Falge, E., Boden, T., Baldocchi, D. D. & Black, T. A. Objective threshold determination for nighttime eddy flux filtering. Agric. Forest Meteorol. 128, 179–197 (2005). 53. Kowalski, A. S. et al. Paired comparisons of carbon exchange between undisturbed and regenerating stands in four managed forests in Europe. Global Change Biol. 10, 1707–1723 (2004). 54. Foken, T. The energy balance closure problem: an overview. Ecol. Appl. 18, 1351–1367 (2008). 55. Reichstein, M. et al. On the separation of net ecosystem exchange into assimilation and ecosystem respiration: review and improved algorithm. Global Change Biology 11, 1424–1439, https://doi.org/10.1111/j.1365-2486.2005.001002.x (2005). 56. Van’t Hoff, J. H. Lectures on theoretical and physical chemistry. Part I. Chemical dynamics (translated by Lehfeldt, R. A.), pp 224–229. Edward Arnold. London (1898). 57. Allen, R. G., Pereira, L. C., Raes, D. & Smith, M. Crop Evapotranspiration. Irrigation and Drainage Paper, 56. United Nations FAO, Rome, pp 300 (1998). Acknowledgements This work was funded by FONDECYT projects 1120713 and 1170429, a grant from the Inter-American Institute for Global Change Research (IAI) [grant number CRN3056], which is supported by the US National Science Foundation [grant number GEO-1128040], and the Spanish Ministry of Economy and Competitiveness project GEI Spain (CGL2014-52838-C2-1-R), including ERDF founds. F. Bravo-Martínez is grateful to CONICYT for the grants “Formación de Capital Humano Avanzado-2009′′, “Beca de Apoyo al término de la tesis doctoral-2012′′, and CORFO INNOVA Grant N° 09CN14-5704. We thank to Enrique Pérez Sanchez-Cañete and Borja Ruíz- Reverter for technical support. We also thank “CODELCO–División Andina” for use of the site. C. Montes acknowledges the NASA Postdoctoral Program and to Universities Space Research Association. Author Contributions F.J. Meza and F. Bravo-Martínez were responsible for the design and implementation of the experiment, and contributed to data processing, analysis and manuscript writing. C. Montes contributed in data processing and analysis, and manuscript writing. P. Serrano-Ortiz and A.S. Kowalski contributed to manuscript writing. All authors reviewed and commented on the manuscript. Additional Information Supplementary information accompanies this paper at https://doi.org/10.1038/s41598-018-26934-z. Competing Interests: The authors declare no competing interests. Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Cre- ative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not per- mitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/. © The Author(s) 2018 ScientiFic Repo R tS | (2018) 8:8570 | DOI:10.1038/s41598-018-26934-z 11
Scientific Reports – Springer Journals
Published: Jun 5, 2018
It’s your single place to instantly
discover and read the research
that matters to you.
Enjoy affordable access to
over 18 million articles from more than
15,000 peer-reviewed journals.
All for just $49/month
Query the DeepDyve database, plus search all of PubMed and Google Scholar seamlessly
Save any article or search result from DeepDyve, PubMed, and Google Scholar... all in one place.
Get unlimited, online access to over 18 million full-text articles from more than 15,000 scientific journals.
Read from thousands of the leading scholarly journals from SpringerNature, Elsevier, Wiley-Blackwell, Oxford University Press and more.
All the latest content is available, no embargo periods.
“Hi guys, I cannot tell you how much I love this resource. Incredible. I really believe you've hit the nail on the head with this site in regards to solving the research-purchase issue.”Daniel C.
“Whoa! It’s like Spotify but for academic articles.”@Phil_Robichaud
“I must say, @deepdyve is a fabulous solution to the independent researcher's problem of #access to #information.”@deepthiw
“My last article couldn't be possible without the platform @deepdyve that makes journal papers cheaper.”@JoseServera