Does hillslope trenching enhance groundwater recharge and baseflow in the Peruvian Andes?

As Andean glaciers rapidly retreat due to climate change, the balance of groundwater and glacial meltwater contributions to stream discharge in tropical, proglacial watersheds will change, potentially increasing vulnerability of water resources. The Shullcas River Watershed, near Huancayo, Peru, is fed only partly by the rapidly receding Huaytapallana glaciers (<20% of dry season flow). To potentially increase recharge and therefore increase groundwater derived baseflow, the government and not‐for‐profit organizations have installed trenches along large swaths of hillslope in the Shullcas Watershed. Our study focuses on a nonglacierized subcatchment of the Shullcas River Watershed and has 2 objectives: (a) create a model of the Shullcas groundwater system and assess the controls on stream discharge and (b) investigate the impact of the infiltration trenches on recharge and baseflow. We first collected hydrologic data from the field including a year‐long hydrograph (2015–2016), meteorological data (2011–2016), and infiltration measurements. We use a recharge model to evaluate the impact of trenched hillslopes on infiltration and runoff processes. Finally, we use a 3‐dimensional groundwater model, calibrated to the measured dry season baseflow, to determine the impact of trenching on the catchment. Simulations show that trenched hillslopes receive approximately 3.5% more recharge, relative to precipitation, compared with unaltered hillslopes. The groundwater model indicates that because the groundwater flow system is fast and shallow, incorporating trenched hillslopes (~2% of study subcatchment area) only slightly increases baseflow in the dry season. Furthermore, the location of trenching is an important consideration: Trenching higher in the catchment (further from the river) and in flatter terrain provides more baseflow during the dry season. The results of this study may have important implications for Andean landscape management and water resources.     

A regional model to assess the hydrological sensitivity of medium size catchments to climate variability

Catchment storage capacity is an important factor in the determination of catchment sensitivity to climate variability. Quantification of catchment sensitivity is in turn important in the regional assessment of the effects of possible climate change. In the present paper, an empirical regional model is proposed that quantifies catchment sensitivity as the ratio of present maximum reservoir storage to catchment storage capacity. Catchment storage capacity is defined theoretically using readily available catchment variables. Present maximum reservoir storage in a catchment, as determined from recession analysis, is expressed as a fraction of catchment storage capacity; the fraction defines catchment sensitivity and depends on storage capacity and annual net precipitation. Average annual conditions for present maximum reservoir storage and average annual net precipitation are used to test the developed model. Although the study used data from only 15 catchments in the Upper Loire region in France, the model proved statistically valid. Storage capacity calculated with the model compares favourably with the baseflow index and a storage index defined in previous research. Values of storage capacity are probable with respect to reported water resources in the area. With the model catchment sensitivity can easily be assessed. Flood or drought prone catchments can be identified as well as a catchment's sensitivity to a catchment-type transition (baseflow versus direct flow dominated catchments). © 1998 John Wiley & Sons, Ltd.     

Combined use of isotopic and hydrometric data to conceptualize ecohydrological processes in a high‐elevation tropical ecosystem

Few high‐elevation tropical catchments worldwide are gauged, and even fewer are studied using combined hydrometric and isotopic data. Consequently, we lack information needed to understand processes governing rainfall–runoff dynamics and to predict their influence on downstream ecosystem functioning. To address this need, we present a combination of hydrometric and water stable isotopic observations in the wet Andean páramo ecosystem of the Zhurucay Ecohydrological Observatory (7.53 km²). The catchment is located in the Andes of south Ecuador between 3400 and 3900 m a.s.l. Water samples for stable isotopic analysis were collected during 2 years (May 2011–May 2013), while rainfall and runoff measurements were continuously recorded since late 2010. The isotopic data reveal that andosol soils predominantly situated on hillslopes drain laterally to histosols (Andean páramo wetlands) mainly located at the valley bottom. Histosols, in turn, feed water to creeks and small rivers throughout the year, establishing hydrologic connectivity between wetlands and the drainage network. Runoff is primarily composed of pre‐event water stored in the histosols, which is replenished by rainfall that infiltrates through the andosols. Contributions from the mineral horizon and the top of the fractured bedrock are small and only seem to influence discharge in small catchments during low flow generation (non‐exceedance flows < Q₃₅). Variations in source contributions are controlled by antecedent soil moisture, rainfall intensity, and duration of rainy periods. Saturated hydraulic conductivity of the soils, higher than the year‐round low precipitation intensity, indicates that Hortonian overland flow rarely occurs during high‐intensity precipitation events. Deep groundwater contributions to discharge seem to be minimal. These results suggest that, in this high‐elevation tropical ecosystem, (1) subsurface flow is a dominant hydrological process and (2) (histosols) wetlands are the major source of stream runoff. Our study highlights that detailed isotopic characterization during short time periods provides valuable information about ecohydrological processes in regions where very few basins are gauged. Copyright © 2016 John Wiley & Sons, Ltd.     

How to construct recursive digital filters for baseflow separation

Recursive digital filtering of hydrographs is a baseflow separation method that can easily be automated and has been recommended for providing reproducible results. In the past, different formulations of the most simple filter type, the so‐called one‐parameter filter, have been proposed. In this paper, a theoretical framework is developed for filter algorithms that were constructed under the assumption that the outflow from an aquifer is linearly proportional to its storage. It is shown that these one‐parameter filters describing an exponential baseflow recession are all special cases of a two‐parameter filter whose equation is specified. Its parameters are the recession constant—which can be objectively determined by a recession analysis—and BFI_max, the maximum value of the baseflow index that can be modelled by the algorithm. This introduces a subjective element into the baseflow calculation, since BFI_max is not measurable. A preliminary analysis based on the results of conventional separation techniques shows that it might be possible to find typical BFI_max values for classes of catchments that can be unequivocally distinguished by their hydrological and hydrogeological characteristics. Copyright © 2004 John Wiley & Sons, Ltd.     

The variation and controls of mean transit times in Australian headwater catchments

Determining mean transit times in headwater catchments is critical for understanding catchment functioning and understanding their responses to changes in landuse or climate. Determining whether mean transit times (MTTs) correlate with drainage density, slope angle, area, or land cover permits a better understanding of the controls on water flow through catchments and allows first-order predictions of MTTs in other catchments to be made. This study assesses whether there are identifiable controls on MTTs determined using ³H in headwater catchments of southeast Australia. Despite MTTs at baseflow varying from a few years to >100 years, it was difficult to predict MTTs using single or groups of readily-measured catchment attributes. The lack of readily-identifiable correlations hampers the prediction of MTTs in adjacent catchments even where these have similar geology, land use, and topography. The long MTTs of the Australian headwater catchments are probably in part due to the catchments having high storage volumes in deeply-weathered regolith, combined with low recharge rates due to high evapotranspiration. However, the difficulty in estimating storage volumes at the catchment scale hampers the use of this parameter to estimate MTTs. The runoff coefficient (the fraction of rainfall exported via the stream) is probably also controlled by evapotranspiration and recharge rates. Correlations between the runoff coefficient and MTTs in individual catchments allow predictions of MTTs in nearby catchments to be made. MTTs are shorter in high rainfall periods as the catchments wet up and shallow water stores are mobilized. Despite the contribution of younger water, the major ion geochemistry in individual catchments commonly does not correlate with MTTs, probably reflecting heterogeneous reactions and varying degrees of evapotranspiration. Documenting MTTs in catchments with high storage volumes and/or low recharge rates elsewhere is important for understanding MTTs in diverse environments.     

Assessment of spatial transferability of process‐based hydrological model parameters in two neighbouring catchments in the Himalayan Region

Estimating the hydrological regime of ungauged catchments in the Himalayan region is challenging due to a lack of sufficient monitoring stations. In this paper, the spatial transferability of the model parameters of the process‐oriented J2000 hydrological model was investigated in 2 glaciated subcatchments of the Koshi river basin in eastern Nepal. The catchments have a high degree of similarity with respect to their static landscape features. The model was first calibrated (1986–1991) and validated (1992–1997) in the Dudh Koshi subcatchment. The calibrated and validated model parameters were then transferred to the nearby Tamor catchment (2001–2009). Sensitivity and uncertainty analyses were carried out for both subcatchments to discover the sensitivity range of the parameters in the two catchments. The model represented the overall hydrograph well in both subcatchments, including baseflow, rising and falling limbs; however, the peak flows were underestimated. The efficiency results according to both Nash–Sutcliffe (ENS) and the coefficient of determination (r²) were above 0.84 in both catchments (1986–1997 in Dudh Koshi and 2001–2009 in Tamor). The ranking of the parameters in respect to their sensitivity matched well for both catchments while taking ENS and log Nash–Sutcliffe (LNS) efficiencies into account. However, there were some differences in sensitivity to ENS and LNS for moderately and less‐sensitive parameters, although the majority (13 out of 16 for ENS and 16 out of 16 for LNS) had a sensitivity response in a similar range. The generalized uncertainty likelihood estimation results suggest that the parameter uncertainty are most of the time within the range and the ensemble mean matches very good (ENS: 0.84) with observed discharge. The results indicate that transfer of the J2000 parameters to a neighbouring catchment in the Himalayan region with similar physiographic landscape characteristics is viable. This indicates the possibility of applying a calibrated process‐based J2000 model to other ungauged catchments in the Himalayan region, which could provide important insights into the hydrological system dynamics and provide much needed information to support water resources planning and management.     

Contribution of occult precipitation to the water balance of páramo ecosystems in the Colombian Andes

Páramos are high‐altitudinal neotropical ecosystems located in the upper regions of the northern Andes. Their hydrology is characterized by an extraordinarily high run‐off ratio. One major contributing mechanism is thought to be fog occurrence, which is common in the páramos and occurs by the cooling of near‐surface moist air, as it is forced to higher elevations by topography. However, field‐based observations and quantification of this flux are rare. We present results of monitoring of occult precipitation, understood as the combination of fog and drizzle inputs, combined with meteorological and soil moisture monitoring for periods between 7 to 17 months in 6 sites distributed over 3 páramos catchments in Colombia: three sites in Romerales (Quindío), two in Chingaza (Cundinamarca), and one in Belmira (Antioquia). Occult precipitation inputs were measured with cylindrical fog gauges with a cover on top. We estimate occult precipitation inputs to add between 7% and 28% to rainfall inputs in the study sites. Our results also show that occult precipitation has a large temporal and spatial variability, both within one site and between sites, which make it difficult to upscale and quantify at a catchment scale. Nevertheless, occult precipitation can be important for downstream water supply given that these inputs are especially concentrated during periods with low rainfall. Lastly, we also find evidence for an increase in soil moisture related to occult precipitation during a dry period in Romerales páramo.     

Nonlinear baseflow recession analysis in watersheds with intermittent streamflow

Abstract

Discharge in most rivers consists mainly of baseflow exfiltrating from shallow groundwater reservoirs, while surface or other direct flows cease soon after rain storms or snowmelt. Analysis of observed baseflow recessions of two rivers in Turkey with intermittent flows and different geographical and climatic characteristics yielded nonlinear storage–outflow relationships of the highly seasonal aquifers. Baseflow separation was carried out using a nonlinear reservoir algorithm. Baseflow seasonality is related to the hydro-climatic conditions influencing groundwater recharge and evapotranspiration of groundwater. As intermittent streams generally have zero flows in the dry season, calibration of recession parameters is in many cases a complicated task.

Citation Aksoy, H. & Wittenberg, H. (2011) Nonlinear baseflow recession analysis in watersheds with intermittent streamflow. Hydrol. Sci. J. 56(2), 226–237.     

Is there a superior conceptual groundwater model structure for baseflow simulation?

Previous work has shown that streamflow response during baseflow conditions is a function of storage, but also that this functional relationship varies among seasons and catchments. Traditionally, hydrological models incorporate conceptual groundwater models consisting of linear or non‐linear storage–outflow functions. Identification of the right model structure and model parameterization however is challenging. The aim of this paper is to systematically test different model structures in a set of catchments where different aquifer types govern baseflow generation processes. Nine different two‐parameter conceptual groundwater models are applied with multi‐objective calibration to transform two different groundwater recharge series derived from a soil‐atmosphere‐vegetation transfer model into baseflow separated from streamflow data. The relative performance differences of the model structures allow to systematically improve the understanding of baseflow generation processes and to identify most appropriate model structures for different aquifer types. We found more versatile and more aquifer‐specific optimal model structures and elucidate the role of interflow, flow paths, recharge regimes and partially contributing storages. Aquifer‐specific recommendations of storage models were found for fractured and karstic aquifers, whereas large storage capacities blur the identification of superior model structures for complex and porous aquifers. A model performance matrix is presented, which highlights the joint effects of different recharge inputs, calibration criteria, model structures and aquifer types. The matrix is a guidance to improve groundwater model structures towards their representation of the dominant baseflow generation processes of specific aquifer types. Copyright © 2014 John Wiley & Sons, Ltd.     

Potential mechanistic causes of increased baseflow across northern Eurasia catchments underlain by permafrost

Warming in the Arctic is occurring at twice the rate of the global average, resulting in permafrost thaw and a restructuring of the Arctic hydrologic cycle as indicated by increased stream discharge during low-flow periods. In these cold regions, permafrost thaw is postulated to increase low-flow discharge, or baseflow, through either: (a) localized increases in groundwater storage and discharge to streams due to increased aquifer transmissivity from thickening of the freeze–thaw layer above permafrost known as the active layer or (b) long-term increases in regional groundwater circulation via enhancement of groundwater–surface water interactions due to extensive permafrost loss over decades. While increasing baseflow has been observed throughout northern Eurasia, the precise mechanistic causes remain elusive. In this study, we differentiate between where these two subsurface physical mechanisms of baseflow increase are occurring by performing a baseflow recession analysis using daily streamflow records from 1913 to 2003 for 139 stations in northern Eurasia underlain by varying permafrost areal extents. Results indicate that from 1913 to 2003, the majority of catchments underlain by continuous permafrost have an increasing trend in their recession flow intercepts, a proxy for increasing active layer thickness. Alternatively, the majority of catchments underlain by permafrost types that are less spatially extensive (e.g., discontinuous, sporadic, isolated, or no permafrost) have decreasing trends in their recession flow intercepts, indicating that a potential increase in active layer thickness is not the driving factor of baseflow variations in these catchments. This may indicate that in catchments underlain by continuous permafrost, active layer thickening correlates with increases in baseflow, whereas, in other catchments with less extensive permafrost, increases in baseflow may be caused by wholesale permafrost loss and vertical talik expansion that enhances regional groundwater circulation. The results of this work may inform our understanding of the subsurface mechanisms responsible for the changing Arctic hydrologic cycle.     

Alpine Hydrogeology: The Critical Role of Groundwater in Sourcing the Headwaters of the World

Groundwater discharge in alpine headwaters sustains baseflow in rivers originating in mountain ranges of the world, which is critically important for aquatic habitats, run-of-river hydropower generation, and downstream water supply. Groundwater storage in alpine watersheds was long considered negligible, but recent field-based studies have shown that aquifers are ubiquitous in the alpine zone with no soil and vegetation. Talus, moraine, and rock glacier aquifers are common in many alpine regions of the world, although bedrock aquifers occur in some geological settings. Alpine aquifers consisting of coarse sediments have a fast recession of discharge after the recharge season (e.g., snowmelt) or rainfall events, followed by a slow recession that sustains discharge over a long period. The two-phase recession is likely controlled by the internal structure of the aquifers. Spatial extent and distribution of individual aquifers determine the groundwater storage-discharge characteristics in first- and second-order watersheds in the alpine zone, which in turn govern baseflow characteristics in major rivers. Similar alpine landforms appear to have similar hydrogeological characteristics in many mountain ranges across the world, suggesting that a common conceptual framework can be used to understand alpine aquifers based on geological and geomorphological settings. Such a framework will be useful for parameterizing storage-discharge characteristics in large river hydrological models.     

Quantifying flows along hydrological pathways by applying a new filtering algorithm in conjunction with master recession curve analysis

Quantifying the proportion of the river hydrograph derived from the different hydrological pathways is essential for understanding the behaviour of a catchment. This paper describes a new approach using the output from master recession curve analysis to inform a new algorithm based on the Lyne and Hollick ‘one‐parameter’ signal analysis filtering algorithm. This approach was applied to six catchments (including two subcatchments of these) in Ireland. The conceptual model for each catchment consists of four main flow pathways: overland flow, interflow, shallow groundwater and deep groundwater. The results were compared with those of the master recession curve analysis, a recharge coefficient approach developed in Ireland and the semi‐distributed, lumped and deterministic hydrological model Nedbør‐Afstrømings‐Model. The new algorithm removes the ‘free variable’ aspect that is typically associated with filtering algorithms and provides a means of estimating the contribution of each pathway that is consistent with the results of hydrograph separation in catchments that are dominated by quick response pathways. These types of catchments are underlain by poorly productive aquifers that are not capable of providing large baseflows in the river. Such aquifers underlie over 73% of Ireland, ensuring that this new algorithm is applicable in the majority of catchments in Ireland and potentially in those catchments internationally that are strongly influenced by the quick‐responding hydrological pathways. Copyright © 2013 John Wiley & Sons, Ltd.     

Runoff and sediment generation on bench‐terraced hillsides: measurements and up‐scaling of a field‐based model

Despite widespread bench‐terracing, stream sediment yields from agricultural hillsides in upland West Java remain high. We studied the causes of this lack of effect by combining measurements at different spatial scales using an erosion process model. Event runoff and sediment yield from two 4‐ha terraced hillside subcatchments were measured and field surveys of land use, bench‐terrace geometry and storage of sediment in the drainage network were conducted for two consecutive years. Runoff was 3·0–3·9% of rainfall and sediment yield was 11–30 t ha⁻¹ yr⁻¹ for different years, subcatchments and calculation techniques. Sediment storage changes in the subcatchment drainage network were less than 2 t ha⁻¹, whereas an additional 0·3–1·5 t ha⁻¹ was stored in the gully between the subcatchment flumes and the main stream. This suggests mean annual sediment delivery ratios of 86–125%, or 80–104% if this additional storage is included. The Terrace Erosion and Sediment Transport (TEST) model developed and validated for the studied environment was parameterized using erosion plot studies, land use surveys and digital terrain analysis to simulate runoff and sediment generation on the terraced hillsides. This resulted in over‐estimates of runoff and under‐estimates of runoff sediment concentration. Relatively poor model performance was attributed to sample bias in the six erosion plots used for model calibration and unaccounted covariance between important terrain attributes such as slope, infiltration capacity, soil conservation works and vegetation cover. Copyright © 2005 John Wiley & Sons, Ltd.     

Geomorphology and vegetation drive hydrochemistry changes in two Northeast Greenland streams

Climate change is causing drastic landscape changes in the Arctic, but how these changes modify stream biogeochemistry is not clear yet. We examined how catchment properties influence stream nitrogen (N) and dissolved organic carbon concentrations (DOC) in a high-Arctic environment. We sampled two contrasting headwater streams (10–15 stations over 4.8 and 6.8 km, respectively) in Northeast Greenland (74°N). We characterized the geomorphology (i.e., bedrock, solifluction and alluvial types) and the vegetation (i.e., barren, fell field, grassland and tundra types) cover of each subcatchment area draining into each sampling station and collected water samples for hydrochemistry characterization. The two sampled streams differed in geomorphology and vegetation cover in the catchment. Aucellaelv catchment was mostly covered by a ‘bedrock’ geomorphology (71%) and ‘fellfield’ vegetation (51%), whereas Kæerelv was mostly covered by ‘alluvial’ geomorphology (65%) and ‘grassland’ and ‘tundra’ vegetation (42% and 41% respectively). Hydrochemistry also differed between the two study streams, with higher concentrations of inorganic N forms in Aucellaelv and lower DOC concentrations, compared to Kærelv. The results from the linear mixed model selection showed that vegetation and geomorphology had contrasting effects on stream hydrochemistry. Subcatchments with higher solifluction sheets and limited vegetation had higher nitrate concentrations but lower DOC concentrations. Interestingly, we also found high variability on the production and removal of nitrate across subcatchments. These results indicate landscape controls to nutrient and organic matter exports via flow paths, soil organic matter stocks and nutrient retention via terrestrial vegetation. Moreover, the results suggest that climate change induced alterations to vegetation cover and soil physical disturbance in high-Arctic catchments will affect stream hydrochemistry, with potential effects in stream productivity, trophic relations as well as change of solute export to downstream coastal areas.     

Contributing sources to baseflow in pre‐alpine headwaters using spatial snapshot sampling

Mountainous headwaters consist of different landscape units including forests, meadows and wetlands. In these headwaters it is unclear which landscape units contribute what percentage to baseflow. In this study, we analysed spatiotemporal differences in baseflow isotope and hydrochemistry to identify catchment‐scale runoff contribution. Three baseflow snapshot sampling campaigns were performed in the Swiss pre‐alpine headwater catchment of the Zwäckentobel (4.25 km²) and six of its adjacent subcatchments. The spatial and temporal variability of δ²H, Ca, DOC, AT, pH, SO₄, Mg and H₄SiO₄ of streamflow, groundwater and spring water samples was analysed and related to catchment area and wetland percentage using bivariate and multivariate methods. Our study found that in the six subcatchments, with variable arrangements of landscape units, the inter‐ and intra catchment variability of isotopic and hydrochemical compositions was small and generally not significant. Stream samples were distinctly different from shallow groundwater. An upper spring zone located near the water divide above 1,400 m and a larger wetland were identified by their distinct spatial isotopic and hydrochemical composition. The upstream wetland percentage was not correlated to the hydrochemical streamflow composition, suggesting that wetlands were less connected and act as passive features with a negligible contribution to baseflow runoff. The isotopic and hydrochemical composition of baseflow changed slightly from the upper spring zone towards the subcatchment outlets and corresponded to the signature of deep groundwater. Our results confirm the need and benefits of spatially distributed snapshot sampling to derive process understanding of heterogeneous headwaters during baseflow. Copyright © 2015 John Wiley & Sons, Ltd.     

The impact of wildfire on baseflow recession rates in California

The effect of wildfire on peak streamflow and annual water yield has been investigated empirically in numerous studies. The effect of wildfire on baseflow recession rates, in contrast, is not well documented. The objective of this paper was to quantify the effect of wildfire on baseflow recession rates in California for both individual watersheds and for all the study watersheds collectively. Two additional variables, antecedent groundwater storage and potential evapotranspiration, were also investigated for their effect on baseflow recession rates and postfire baseflow recession rate response. Differences between prefire and postfire baseflow recession rates were modeled statistically in 8 watersheds using a mixed statistical model that accounted for fixed and random effects. For the all‐watershed model, antecedent groundwater storage, potential evapotranspiration, and wildfire were each found to be significant controls on baseflow recession rates. Wildfire decreased baseflow recession rates 52.5% (37.6% to 66.0%), implying that postfire reductions in above‐ground vegetation (e.g., decreased interception, decreased evapotranspiration) were a stronger control on baseflow recession rate change than hydrophobicity. At an individual watershed scale, baseflow recession rate response to wildfire was found to be sensitive to intraannual differences in antecedent groundwater storage in 2 watersheds, with the effect of wildfire on baseflow recession rates being greater with lower levels of antecedent groundwater storage. Examination of burn severity for a subset of the study watersheds pointed to riparian zone burn severity as a potential primary control on postfire recession rate change. This study demonstrates that wildfire may have a substantial impact on fluxes to and from groundwater storages, altering the rate at which baseflow recedes.     

Transient modelling of lacustrine regressions: two case studies from the Andean Altiplano

A model was developed for estimating the delay between a change in climatic conditions and the corresponding fall of water level in large lakes. The input data include: rainfall, temperature, extraterrestrial radiation and astronomical mid‐month daylight hours. The model uses two empirical coefficients for computing the potential evaporation and one parameter for the soil capacity. The case studies are two subcatchments of the Altiplano (196 000 km²), in which the central low points are Lake Titicaca and a salar corresponding to the desiccation of the Tauca palaeolake. During the Holocene, the two catchments experienced a 100 m fall in water level corresponding to a decrease in water surface area of 3586 km² and 55 000 km², respectively. Under modern climatic conditions with a marked rainy season, the model allows simulation of water levels in good agreement with the observations: 3810 m a.s.l. for Lake Titicaca and lack of permanent wide ponds in the southern subcatchment. Simulations were carried out under different climatic conditions that might explain the Holocene fall in water level. Computed results show quite different behaviour for the two subcatchments. For the northern subcatchment, the time required for the 100 m fall in lake‐level ranges between 200 and 2000 years when, compared with the present conditions, (i) the rainfall is decreased by 15% (640 mm/year), or (ii) the temperature is increased by 5·5 °C, or (iii) rainfall is distributed equally over the year. For the southern subcatchment (Tauca palaeolake), the time required for a 100 m decrease in water level ranges between 50 and 100 years. This decrease requires precipitation values lower than 330 mm/year. Copyright © 2004 John Wiley & Sons, Ltd.     

Effect of land cover and hydro‐meteorological controls on soil water DOC concentrations in a high‐elevation tropical environment

Páramo soils store high amounts of organic carbon. However, the effects of climate change and changes in land cover and use (LC/LU) in this high‐elevation tropical ecosystem may cause a decrease in their carbon storage capacity. Therefore, better understanding of the factors influencing the Páramo soils' carbon storage and export is urgently needed. To fill this knowledge gap, we investigated the differences in dissolved organic carbon (DOC) content in the soil water of four LC/LU types (tussock grass, natural forest, pine plantations, and pasture) and the factors controlling its variability in the Quinuas Ecohydrological Observatory in south Ecuador. Weekly measurements of soil water DOC concentrations, meteorological variables, soil water content, and temperature from various depths and slope positions were monitored within the soils' organic and mineral horizons between October 2014 and January 2017. These data were used to generate regression trees and random forest statistical models to identify the factors controlling soil water DOC concentrations. From high to low concentrations, natural forest depict the highest DOC concentrations followed by pasture, tussock grass, and pine forest. For all LC/LU types, DOC concentrations increase with decreasing soil moisture. Our results also show that LC/LU is the most important predictor of soil water DOC concentrations, followed by sampling depth and soil moisture. Interestingly, atmospheric variables and antecedent evapotranspiration and precipitation conditions show only little influence on DOC concentrations during the monitoring period. Our findings provide unique information that can help improve the management of soil and water resources in the Páramo and other peat dominated ecosystems elsewhere.