A concerted research effort to advance the hydrological understanding of tropical páramos

Páramos, a neotropical alpine grassland-peatland biome of the northern Andes and Central America, play an essential role in regional and global cycles of water, carbon, and nutrients. They act as water towers, delivering water and ecosystem services from the high mountains down to the Pacific, Caribbean, and Amazon regions. Páramos are also widely recognized as a biodiversity and climate change hot spots, yet they are threatened by anthropogenic activities and environmental changes. Despite their importance for water security and carbon storage, and their vulnerability to human activities, only three decades ago, páramos were severely understudied. Increasing awareness of the need for hydrological evidence to guide sustainable management of páramos prompted action for generating data and for filling long-standing knowledge gaps. This has led to a remarkably successful increase in scientific knowledge, induced by a strong interaction between the scientific, policy, and (local) management communities. A combination of well-established and innovative approaches has been applied to data collection, processing, and analysis. In this review, we provide a short overview of the historical development of research and state of knowledge of the hydrometeorology, flux dynamics, anthropogenic impacts, and the influence of extreme events in páramos. We then present emerging technologies for hydrology and water resources research and management applied to páramos. We discuss how converging science and policy efforts have leveraged traditional and new observational techniques to generate an evidence base that can support the sustainable management of páramos. We conclude that this co-evolution of science and policy was able to successfully cover different spatial and temporal scales. Lastly, we outline future research directions to showcase how sustainable long-term data collection can foster the responsible conservation of páramos water towers.     

The information system of the French Peatland Observation Service: Service National d'Observation Tourbières – A valuable tool to assess the impact of global changes on the hydrology and biogeochemistry of temperate peatlands through long term monitoring

Mitigating and adapting to global changes requires a better understanding of the response of the Biosphere to these environmental variations. Human disturbances and their effects act in the long term (decades to centuries) and consequently, a similar time frame is needed to fully understand the hydrological and biogeochemical functioning of a natural system. To this end, the ‘Centre National de la Recherche Scientifique’ (CNRS) promotes and certifies long-term monitoring tools called national observation services or ‘Service National d'Observation’ (SNO) in a large range of hydrological and biogeochemical systems (e.g., cryosphere, catchments, aquifers). The SNO investigating peatlands, the SNO ‘Tourbières’, was certified in 2011 ( https://www.sno-tourbieres.cnrs.fr/ ). Peatlands are mostly found in the high latitudes of the northern hemisphere and French peatlands are located in the southern part of this area. Thus, they are located in environmental conditions that will occur in northern peatlands in coming decades or centuries and can be considered as sentinels. The SNO Tourbières is composed of four peatlands: La Guette (lowland central France), Landemarais (lowland oceanic western France), Frasne (upland continental eastern France) and Bernadouze (upland southern France). Thirty target variables are monitored to study the hydrological and biogeochemical functioning of the sites. They are grouped into four datasets: hydrology, fluvial export of organic matter, greenhouse gas fluxes and meteorology/soil physics. The data from all sites follow a common processing chain from the sensors to the public repository. The raw data are stored on an FTP server. After operator or automatic processing, data are stored in a database, from which a web application extracts the data to make them available ( https://data-snot.cnrs.fr/data-access/ ). Each year at least, an archive of each dataset is stored in Zenodo, with a digital object identifier (DOI) attribution ( https://zenodo.org/communities/sno_tourbieres_data/ ).     

中国高尔夫球场时空演变、影响因素和扩散模式

利用空间分析和数理统计等方法探究了1984—2016年中国高尔夫球场数量时空演变、影响因素和扩散模式。结果表明：1）球场数量变化经历了低速增长期、稳定增长期、加速增长期和紧缩期；2）球场空间分布范围经历了扩张期和收缩期，向北、向西扩散较为明显；3）目前球场主要分布在经济发展水平较高且外向型经济突出的长三角、珠三角和环渤海地区；4）人均GDP、人口规模、区域开放程度、交通优势度、气候类型和土地资源稀缺性与球场空间分布均呈显著正相关，其中交通优势度相关系数最大；5）高尔夫球场空间扩散模式兼具接触性扩散和等级性扩散。     

Coupled groundwater flow and heat transport simulation for estimating transient aquifer–stream exchange at the lowland River Spree (Germany)

Subsurface flow and heat transport near Freienbrink, NE Germany, was simulated in order to study groundwater–surface water exchange between a floodplains aquifer and a section of the lowland River Spree and an adjacent oxbow. Groundwater exfiltration was the dominant process, and only fast surface water level rises resulted in temporary infiltration into the aquifer. The main groundwater flow paths are identified based on a 3D groundwater flow model. To estimate mass fluxes across the aquifer–surface water interfaces, a 2D flow and heat transport modelling approach along a transect of 12 piezometers was performed. Results of steady‐state and transient water level simulations show an overall high accuracy with a Spearman coefficient ρ = 0.9996 and root mean square error (RMSE) = 0.008 m. Based on small groundwater flow velocities of about 10^?7 to 10^?6 ms^?1, mean groundwater exfiltration rates of 233 l m^?2 d^?1 are calculated. Short periods of surface water infiltration into the aquifer do not exceed 10 days, and the infiltration rates are in the same range. The heat transport was modelled with slightly less accuracy (ρ = 0.8359 and RMSE = 0.34 °C). In contrast to the predominant groundwater exfiltration, surface water temperatures determine the calculated temperatures in the upper aquifer below both surface water bodies down to 10 m during the whole simulation period. These findings emphasize prevailing of heat conduction over advection in the upper aquifer zones, which seems to be typical for lowland streams with sandy aquifer materials and low hydraulic gradients. Moreover, this study shows the potential of coupled numerical flow and heat transport modelling to understand groundwater–surface water exchange processes in detail. Copyright © 2013 John Wiley & Sons, Ltd.     

Identification of groundwater mean transit times of precipitation and riverbank infiltration by two‐component lumped parameter models

Groundwater transit time is an essential hydrologic metric for groundwater resources management. However, especially in tropical environments, studies on the transit time distribution (TTD) of groundwater infiltration and its corresponding mean transit time (mTT) have been extremely limited due to data sparsity. In this study, we primarily use stable isotopes to examine the TTDs and their mTTs of both vertical and horizontal infiltration at a riverbank infiltration area in the Vietnamese Mekong Delta (VMD), representative of the tropical climate in Asian monsoon regions. Precipitation, river water, groundwater, and local ponding surface water were sampled for 3 to 9 years and analysed for stable isotopes (δ¹⁸O and δ²H), providing a unique data set of stable isotope records for a tropical region. We quantified the contribution that the two sources contributed to the local shallow groundwater by a novel concept of two‐component lumped parameter models (LPMs) that are solved using δ¹⁸O records. The study illustrates that two‐component LPMs, in conjunction with hydrological and isotopic measurements, are able to identify subsurface flow conditions and water mixing at riverbank infiltration systems. However, the predictive skill and the reliability of the models decrease for locations farther from the river, where recharge by precipitation dominates, and a low‐permeable aquitard layer above the highly permeable aquifer is present. This specific setting impairs the identifiability of model parameters. For river infiltration, short mTTs (<40 weeks) were determined for sites closer to the river (<200 m), whereas for the precipitation infiltration, the mTTs were longer (>80 weeks) and independent of the distance to the river. The results not only enhance the understanding of the groundwater recharge dynamics in the VMD but also suggest that the highly complex mechanisms of surface–groundwater interaction can be conceptualized by exploiting two‐component LPMs in general. The model concept could thus be a powerful tool for better understanding both the hydrological functioning of mixing processes and the movement of different water components in riverbank infiltration systems.     

Spatially distributed tracer‐aided modelling to explore water and isotope transport,storage and mixing in a pristine,humid tropical catchment

Rapidly transforming headwater catchments in the humid tropics provide important resources for drinking water, irrigation, hydropower, and ecosystem connectivity. However, such resources for downstream use remain unstudied. To improve understanding of the behaviour and influence of pristine rainforests on water and tracer fluxes, we adapted the relatively parsimonious, spatially distributed tracer‐aided rainfall–runoff (STARR) model using event‐based stable isotope data for the 3.2‐km² San Lorencito catchment in Costa Rica. STARR was used to simulate rainforest interception of water and stable isotopes, which showed a significant isotopic enrichment in throughfall compared with gross rainfall. Acceptable concurrent simulations of discharge (Kling–Gupta efficiency [KGE] ~0.8) and stable isotopes in stream water (KGE ~0.6) at high spatial (10 m) and temporal (hourly) resolution indicated a rapidly responding system. Around 90% of average annual streamflow (2,099 mm) was composed of quick, near‐surface runoff components, whereas only ~10% originated from groundwater in deeper layers. Simulated actual evapotranspiration (ET) from interception and soil storage were low (~420 mm/year) due to high relative humidity (average 96%) and cloud cover limiting radiation inputs. Modelling suggested a highly variable groundwater storage (~10 to 500 mm) in this steep, fractured volcanic catchment that sustains dry season baseflows. This groundwater is concentrated in riparian areas as an alluvial–colluvial aquifer connected to the stream. This was supported by rainfall–runoff isotope simulations, showing a “flashy” stream response to rainfall with only a moderate damping effect and a constant isotope signature from deeper groundwater (~400‐mm additional mixing volume) during baseflow. The work serves as a first attempt to apply a spatially distributed tracer‐aided model to a tropical rainforest environment exploring the hydrological functioning of a steep, fractured‐volcanic catchment. We also highlight limitations and propose a roadmap for future data collection and spatially distributed tracer‐aided model development in tropical headwater catchments.     

High-P granulites of the Songshugou area (Qinling Orogen,east-central China): Petrography,phase relations,and U/Pb zircon geochronology

High-pressure (HP) granulites provide telling records of mineral reactions at upper mantle to lower crustal levels and key information on the fate of material in subduction systems. The latter especially applies when they abut eclogite and mantle dunite because such rock associations are crucial for understanding the incompletely known processes at the interface of converging plates. A continental arc, active c. 520–395 Ma ago, formed an enigmatic example of such a rock association in the Songshugou area, Qinling Orogen. To unravel the juxtaposition of the distinct rocks, this study combines petrography, phase equilibria modelling, conventional thermobarometry, and zircon U–Th–Pb–Ti–REE analysis. Two mafic HP granulites, which contain the mineral assemblages garnet–clinopyroxene–plagioclase–rutile–mesoperthite–quartz and garnet–clinopyroxene–plagioclase–rutile, experienced peak metamorphic conditions of ≤1.4 GPa, 860°C and ~1.3 GPa, ≥910°C, respectively. During decompression and cooling, at 489 ± 4 Ma, amphibole lamellae unmixed from a clinopyroxene solid solution and orthopyroxene in part replaced garnet. A felsic HP granulite shows equilibration of garnet, perthite, antiperthite, kyanite, quartz, and rutile at 810–860°C, ~1.2 GPa, sillimanite growth during decompression, and upper amphibolite facies cooling at 510 ± 4 Ma. Though the thermobarometric data are just within the methodological errors, the U/Pb zircon ages imply the HP granulites did not evolve coherently. The HP granulites either represent foundered lower arc crust or originated from subduction erosion because their geochemistry is indistinguishable from that of the hanging-wall plate. Published and new pressure–temperature–time–deformation paths converge at ~710°C, ~0.9 GPa, and ≲470 Ma, implying exhumation tectonics juxtaposed the HP granulites with a mélange of eclogite and mantle dunite at lower crustal levels. This study highlights that lower arc crust can comprise material of diverse evolution.     

A rainfall‐runoff model parameterized from GIS and runoff data

A guiding principle in hydrological modelling should be to keep the number of calibration parameters to a minimum. A reduced number of parameters to be calibrated, while maintaining the accuracy and detail required by modern hydrological models, will reduce parameter and model structure uncertainty and improve model diagnostics. In this study, the dynamics of runoff are derived from the distribution of distances from points in the catchments to the nearest stream. This distribution is unique for each catchment and can be determined from a geographical information system. The distribution of distances, will, when a celerity of (subsurface) flow is introduced, provide a distribution of travel times, or a unit hydrograph (UH). For spatially varying levels of saturation deficit, we have different celerities and, hence, different UHs. Runoff is derived from the superposition of the different UHs. This study shows how celerities can be estimated if we assume that recession events represent the combined UHs for different levels of saturation deficit. A new soil moisture routine which estimates saturated and unsaturated volumes of subsurface water and with only one parameter to calibrate is included in the new model. The performance of the new model is compared with that of the Swedish HBV model and is found to perform equally well for eight Norwegian catchments although the number of parameters to be calibrated in the module concerning soil moisture and runoff dynamics is reduced from seven in the HBV model to one in the new model. It is also shown that the new model has a more realistic representation of the subsurface hydrology. Copyright © 2013 John Wiley & Sons, Ltd.     

1990-2020年间念青唐古拉山中段北坡边坝地区冰川变化及气候响应

念青唐古拉山作为青藏高原东南缘重要山古冰川分布区,受季风影响,各区域冰川变化特征差异明显。论文通过Landsat TM/ETM+/OLI资料、ASRTMGDEM与气象数据,采用比值—阈值法、目视解译和VOLTA模型,结合实地考察,对1990—2020年间念青唐古拉山中段北坡边坝地区现代冰川进退状况、面积变化、冰储量变化以及冰川变化对气候变化响应特征进行研究。结果表明：① 1990—2020年5条冰川(玉贡拉冰川、玛拉波冰川、祥格拉冰川、孔嘎冰川、贡日—庚东冰川)末端高程逐渐升高,面积和冰储量分别减少30.38 km²和4.64 km³,总体缩减并呈现加速趋势。② 冰川冰储量减少0.14~1.92 km³,总体变化率为0.40%·a^-1。2020年上述5条冰川储量占1990年冰川储量的比例分别为0.70、0.99、0.98、0.91和0.82,显示出冰川规模越大,在短时间尺度的变化量越小。③ 气象数据分析显示,1990—2020年研究区冰川变化受气温升高主导,平均气温变化率为0.51 ℃。水热组合呈现温度升高—降水减少,且在最后10 a日益显著,预测未来冰川变化仍受气温控制并呈加速退缩趋势。④ 区域对比研究表明,念青唐古拉山冰川面积变化总体呈退缩状态,但各区域冰川变化特征差异明显。同时,不同研究方法对同一冰川区冰储量模拟结果相差较大,相对误差范围为34.45%~115.49%,精确的冰储量可对比研究方法仍有待进一步研究。