Refine
Year of publication
- 2019 (74) (remove)
Document Type
- Article (55)
- Postprint (7)
- Doctoral Thesis (5)
- Other (4)
- Part of a Book (1)
- Habilitation Thesis (1)
- Review (1)
Language
- English (74) (remove)
Is part of the Bibliography
- yes (74) (remove)
Keywords
- climate change (3)
- hydrology (3)
- Abandonment (2)
- Biodiversity (2)
- Determinants (2)
- Drivers (2)
- Ecosystem services (2)
- GLOF (2)
- Intensity (2)
- National-Parks (2)
Institute
- Institut für Umweltwissenschaften und Geographie (74) (remove)
Optical flow models as an open benchmark for radar-based precipitation nowcasting (rainymotion v0.1)
(2019)
Quantitative precipitation nowcasting (QPN) has become an essential technique in various application contexts, such as early warning or urban sewage control. A common heuristic prediction approach is to track the motion of precipitation features from a sequence of weather radar images and then to displace the precipitation field to the imminent future (minutes to hours) based on that motion, assuming that the intensity of the features remains constant (“Lagrangian persistence”). In that context, “optical flow” has become one of the most popular tracking techniques. Yet the present landscape of computational QPN models still struggles with producing open software implementations. Focusing on this gap, we have developed and extensively benchmarked a stack of models based on different optical flow algorithms for the tracking step and a set of parsimonious extrapolation procedures based on image warping and advection. We demonstrate that these models provide skillful predictions comparable with or even superior to state-of-the-art operational software. Our software library (“rainymotion”) for precipitation nowcasting is written in the Python programming language and openly available at GitHub (https://github.com/hydrogo/rainymotion, Ayzel et al., 2019). That way, the library may serve as a tool for providing fast, free, and transparent solutions that could serve as a benchmark for further model development and hypothesis testing – a benchmark that is far more advanced than the conventional benchmark of Eulerian persistence commonly used in QPN verification experiments.
Optical flow models as an open benchmark for radar-based precipitation nowcasting (rainymotion v0.1)
(2019)
Quantitative precipitation nowcasting (QPN) has become an essential technique in various application contexts, such as early warning or urban sewage control. A common heuristic prediction approach is to track the motion of precipitation features from a sequence of weather radar images and then to displace the precipitation field to the imminent future (minutes to hours) based on that motion, assuming that the intensity of the features remains constant (“Lagrangian persistence”). In that context, “optical flow” has become one of the most popular tracking techniques. Yet the present landscape of computational QPN models still struggles with producing open software implementations. Focusing on this gap, we have developed and extensively benchmarked a stack of models based on different optical flow algorithms for the tracking step and a set of parsimonious extrapolation procedures based on image warping and advection. We demonstrate that these models provide skillful predictions comparable with or even superior to state-of-the-art operational software. Our software library (“rainymotion”) for precipitation nowcasting is written in the Python programming language and openly available at GitHub (https://github.com/hydrogo/rainymotion, Ayzel et al., 2019). That way, the library may serve as a tool for providing fast, free, and transparent solutions that could serve as a benchmark for further model development and hypothesis testing – a benchmark that is far more advanced than the conventional benchmark of Eulerian persistence commonly used in QPN verification experiments.
During the last few decades, the rapid separation of the Small Aral Sea from the isolated basin has changed its hydrological and ecological conditions tremendously. In the present study, we developed and validated the hybrid model for the Syr Darya River basin based on a combination of state-of-the-art hydrological and machine learning models. Climate change impact on freshwater inflow into the Small Aral Sea for the projection period 2007–2099 has been quantified based on the developed hybrid model and bias corrected and downscaled meteorological projections simulated by four General Circulation Models (GCM) for each of three Representative Concentration Pathway scenarios (RCP). The developed hybrid model reliably simulates freshwater inflow for the historical period with a Nash–Sutcliffe efficiency of 0.72 and a Kling–Gupta efficiency of 0.77. Results of the climate change impact assessment showed that the freshwater inflow projections produced by different GCMs are misleading by providing contradictory results for the projection period. However, we identified that the relative runoff changes are expected to be more pronounced in the case of more aggressive RCP scenarios. The simulated projections of freshwater inflow provide a basis for further assessment of climate change impacts on hydrological and ecological conditions of the Small Aral Sea in the 21st Century.
During the last few decades, the rapid separation of the Small Aral Sea from the isolated basin has changed its hydrological and ecological conditions tremendously. In the present study, we developed and validated the hybrid model for the Syr Darya River basin based on a combination of state-of-the-art hydrological and machine learning models. Climate change impact on freshwater inflow into the Small Aral Sea for the projection period 2007-2099 has been quantified based on the developed hybrid model and bias corrected and downscaled meteorological projections simulated by four General Circulation Models (GCM) for each of three Representative Concentration Pathway scenarios (RCP). The developed hybrid model reliably simulates freshwater inflow for the historical period with a Nash-Sutcliffe efficiency of 0.72 and a Kling-Gupta efficiency of 0.77. Results of the climate change impact assessment showed that the freshwater inflow projections produced by different GCMs are misleading by providing contradictory results for the projection period. However, we identified that the relative runoff changes are expected to be more pronounced in the case of more aggressive RCP scenarios. The simulated projections of freshwater inflow provide a basis for further assessment of climate change impacts on hydrological and ecological conditions of the Small Aral Sea in the 21st Century.
Above and underground hydrological processes depend on soil moisture (SM) variability, driven by different environmental factors that seldom are well-monitored, leading to a misunderstanding of soil water temporal patterns. This study investigated the stability of the SM temporal dynamics to different monitoring temporal resolutions around the border between two soil types in a tropical watershed. Four locations were instrumented in a small-scale watershed (5.84 km(2)) within the tropical coast of Northeast Brazil, encompassing different soil types (Espodossolo Humiluvico or Carbic Podzol, and Argissolo Vermelho-Amarelo or Haplic Acrisol), land covers (Atlantic Forest, bush vegetation, and grassland) and topographies (flat and moderate slope). The SM was monitored at a temporal resolution of one hour along the 2013-2014 hydrological year and then resampled a resolutions of 6 h, 12 h, 1 day, 2 days, 4 days, 7 days, and 15 days. Descriptive statistics, temporal variability, time-stability ranking, and hierarchical clustering revealed uneven associations among SM time components. The results show that the time-invariant component ruled SM temporal variability over the time-varying parcel, either at high or low temporal resolutions. Time-steps longer than 2 days affected the mean statistical metrics of the SM time-variant parcel. Additionally, SM at downstream and upstream sites behaved differently, suggesting that the temporal mean was regulated by steady soil properties (slope, restrictive layer, and soil texture), whereas their temporal anomalies were driven by climate (rainfall) and hydrogeological (groundwater level) factors. Therefore, it is concluded that around the border between tropical soil types, the distinct behaviour of time-variant and time-invariant components of SM time series reflects different combinations of their soil properties.
The Value of Empirical Data for Estimating the Parameters of a Sociohydrological Flood Risk Model
(2019)
In this paper, empirical data are used to estimate the parameters of a sociohydrological flood risk model. The proposed model, which describes the interactions between floods, settlement density, awareness, preparedness, and flood loss, is based on the literature. Data for the case study of Dresden, Germany, over a period of 200years, are used to estimate the model parameters through Bayesian inference. The credibility bounds of their estimates are small, even though the data are rather uncertain. A sensitivity analysis is performed to examine the value of the different data sources in estimating the model parameters. In general, the estimated parameters are less biased when using data at the end of the modeled period. Data about flood awareness are the most important to correctly estimate the parameters of this model and to correctly model the system dynamics. Using more data for other variables cannot compensate for the absence of awareness data. More generally, the absence of data mostly affects the estimation of the parameters that are directly related to the variable for which data are missing. This paper demonstrates that combining sociohydrological modeling and empirical data gives additional insights into the sociohydrological system, such as quantifying the forgetfulness of the society, which would otherwise not be easily achieved by sociohydrological models without data or by standard statistical analysis of empirical data.
The improvement of process representations in hydrological models is often only driven by the modelers' knowledge and data availability. We present a comprehensive comparison between two hydrological models of different complexity that is developed to support (1) the understanding of the differences between model structures and (2) the identification of the observations needed for model assessment and improvement. The comparison is conducted on both space and time and by aggregating the outputs at different spatiotemporal scales. In the present study, mHM, a process‐based hydrological model, and ParFlow‐CLM, an integrated subsurface‐surface hydrological model, are used. The models are applied in a mesoscale catchment in Germany. Both models agree in the simulated river discharge at the outlet and the surface soil moisture dynamics, lending their supports for some model applications (drought monitoring). Different model sensitivities are, however, found when comparing evapotranspiration and soil moisture at different soil depths. The analysis supports the need of observations within the catchment for model assessment, but it indicates that different strategies should be considered for the different variables. Evapotranspiration measurements are needed at daily resolution across several locations, while highly resolved spatially distributed observations with lower temporal frequency are required for soil moisture. Finally, the results show the impact of the shallow groundwater system simulated by ParFlow‐CLM and the need to account for the related soil moisture redistribution. Our comparison strategy can be applied to other models types and environmental conditions to strengthen the dialog between modelers and experimentalists for improving process representations in Earth system models.
Movement ecology aims to provide common terminology and an integrative framework of movement research across all groups of organisms. Yet such work has focused on unitary organisms so far, and thus the important group of filamentous fungi has not been considered in this context. With the exception of spore dispersal, movement in filamentous fungi has not been integrated into the movement ecology field. At the same time, the field of fungal ecology has been advancing research on topics like informed growth, mycelial translocations, or fungal highways using its own terminology and frameworks, overlooking the theoretical developments within movement ecology. We provide a conceptual and terminological framework for interdisciplinary collaboration between these two disciplines, and show how both can benefit from closer links: We show how placing the knowledge from fungal biology and ecology into the framework of movement ecology can inspire both theoretical and empirical developments, eventually leading towards a better understanding of fungal ecology and community assembly. Conversely, by a greater focus on movement specificities of filamentous fungi, movement ecology stands to benefit from the challenge to evolve its concepts and terminology towards even greater universality. We show how our concept can be applied for other modular organisms (such as clonal plants and slime molds), and how this can lead towards comparative studies with the relationship between organismal movement and ecosystems in the focus.
From waste to resource
(2019)
Reservoir networks have been established worldwide to ensure water supply, but water availability is endangered quantitatively and qualitatively by sedimentation. Reuse of sediment silted in reservoirs as fertilizer has been proposed, thus transforming nutrient-enriched sediments from waste into resource. The aim of this study is to assess the potential of reusing sediment as a nutrient source for agriculture a semiarid basin in Brazil. where 1029 reservoirs were identified. Sedimentation was modelled for the entire reservoir network, accounting for 7 x 10(5) tons of y(-1)sediment deposition. Nutrients contents in reservoir sediments was analysed and com- pared to nutrients contents of agricultural soils in the catchment. The potential of reusing sediment as fertilizer was assessed for maize crops (Zea mays L) and the sediment mass required to fertilize the soil was computed considering that the crop nitrogen requirement would be fully provided by the sediment. Economic feasibility was analysed by comparing the costs of the proposed practice to those obtained if the area was fertilized by traditional means. Results showed that, where reservoirs fall dry frequently and sediments can be removed by excavation, soil fertilization with sediment presents lower costs than those observed for application of commercial chemical fertilizers. Compared to conventional fertilization, when using sediments with high nutrient content, 25% of costs could be saved, while when using sediments with low nutrient content costs are 9% higher. According to the local conditions, sediments with nitrogen content above 1.5 g kg(-1) are cost efficient as nitrogen source. However, physical and chemical analyses are recommended to define the sediment mass to be used and to identify any constraint to the application of the practice, like the high sodium adsorption ratio observed in one of the studied reservoirs, which can contribute to soil salinization. (C) 2019 Elsevier B.V. All rights reserved.
Remote Sensing technologies allow to map biophysical, biochemical, and earth surface parameters of the land surface. Of especial interest for various applications in environmental and urban sciences is the combination of spectral and 3D elevation information. However, those two data streams are provided separately by different instruments, namely airborne laser scanner (ALS) for elevation and a hyperspectral imager (HSI) for high spectral resolution data. The fusion of ALS and HSI data can thus lead to a single data entity consistently featuring rich structural and spectral information. In this study, we present the application of fusing the first pulse return information from ALS data at a sub-decimeter spatial resolution with the lower-spatial resolution hyperspectral information available from the HSI into a hyperspectral point cloud (HSPC). During the processing, a plausible hyperspectral spectrum is assigned to every first-return ALS point. We show that the complementary implementation of spectral and 3D information at the point-cloud scale improves object-based classification and information extraction schemes. This improvements have great potential for numerous land cover mapping and environmental applications.
The ability to reflect is considered an essential element of Education for Sustainable Development (ESD) and a key competence for learners and educators in ESD (UNECE Strategy for ESD, 2012). In contrast to its high importance, little is known about how reflective thinking can be identified, influenced or increased in the classroom. Therefore, the objective of this study is to address this need by developing an empirical multi-stage model designed to help educators diagnose different levels of reflective thinking and to identify factors that influence students’ reflective thinking about sustainability. Based on a 4–8-week project with grade 10 and 11 students studying sustainability, reflective thinking performance using weblogs as reflective journals was analysed. In addition, qualitative semi-structured interviews were conducted with the teachers to comprehend the learning environment and the personal value they assigned to ESD in their geography class. To determine the levels of reflective thinking achieved by the students, the study built on the work of Dewey (1933) and pre-existing multi-stage models of reflective thinking (Bain, Ballantyne, & Packer, 1999; Chen, Wei, Wu, & Uden, 2009). Using a qualitative, iterative data analysis, the study adapted the stage models to be applicable in ESD and found great differences in the students’ reflection levels. Furthermore, the study identified eight factors that influence students’ reflective thinking about sustainability. The outcomes of this study may be valuable for educators in high school and higher education, who seek to diagnose their students’ reflective thinking performance and facilitate reflection about sustainability.
For effective disaster risk management and adaptation planning, a good understanding of current and projected flood risk is required. Recent advances in quantifying flood risk at the regional and global scale have largely neglected critical infrastructure, or addressed this important sector with insufficient detail. Here, we present the first European-wide assessment of current and future flood risk to railway tracks for different global warming scenarios using an infrastructure-specific damage model. We find that the present risk, measured as expected annual damage, to railway networks in Europe is approx. (sic)581 million per year, with the highest risk relative to the length of the network in North Macedonia, Croatia, Norway, Portugal, and Germany. Based on an ensemble of climate projections for RCP8.5, we show that current risk to railway networks is projected to increase by 255% under a 1.5 degrees C, by 281% under a 2 degrees C, and by 310% under a 3 degrees C warming scenario. The largest increases in risk under a 3 degrees C scenario are projected for Slovakia, Austria, Slovenia, and Belgium. Our advances in the projection of flood risk to railway infrastructure are important given their criticality, and because losses to public infrastructure are usually not insured or even uninsurable in the private market. To cover the risk increase due to climate change, European member states would need to increase expenditure in transport by (sic)1.22 billion annually under a 3 degrees C warming scenario without further adaptation. Limiting global warming to the 1.5 degrees C goal of the Paris Agreement would result in avoided losses of (sic)317 million annually.
A digital filter is introduced which treats the problem of predictability versus time averaging in a continuous, seamless manner. This seamless filter (SF) is characterized by a unique smoothing rule that determines the strength of smoothing in dependence on lead time. The rule needs to be specified beforehand, either by expert knowledge or by user demand. As a result, skill curves are obtained that allow a predictability assessment across a whole range of time-scales, from daily to seasonal, in a uniform manner. The SF is applied to downscaled SEAS5 ensemble forecasts for two focus regions in or near the tropical belt, the river basins of the Karun in Iran and the Sao Francisco in Brazil. Both are characterized by strong seasonality and semi-aridity, so that predictability across various time-scales is in high demand. Among other things, it is found that from the start of the water year (autumn), areal precipitation is predictable with good skill for the Karun basin two and a half months ahead; for the Sao Francisco it is only one month, longer-term prediction skill is just above the critical level.
A digital filter is introduced which treats the problem of predictability versus time averaging in a continuous, seamless manner. This seamless filter (SF) is characterized by a unique smoothing rule that determines the strength of smoothing in dependence on lead time. The rule needs to be specified beforehand, either by expert knowledge or by user demand. As a result, skill curves are obtained that allow a predictability assessment across a whole range of time-scales, from daily to seasonal, in a uniform manner. The SF is applied to downscaled SEAS5 ensemble forecasts for two focus regions in or near the tropical belt, the river basins of the Karun in Iran and the Sao Francisco in Brazil. Both are characterized by strong seasonality and semi-aridity, so that predictability across various time-scales is in high demand. Among other things, it is found that from the start of the water year (autumn), areal precipitation is predictable with good skill for the Karun basin two and a half months ahead; for the Sao Francisco it is only one month, longer-term prediction skill is just above the critical level.
Continental rift systems form by propagation of isolated rift segments that interact, and eventually evolve into continuous zones of deformation. This process impacts many aspects of rifting including rift morphology at breakup, and eventual ocean-ridge segmentation. Yet, rift segment growth and interaction remain enigmatic. Here we present geological data from the poorly documented Ririba rift (South Ethiopia) that reveals how two major sectors of the East African rift, the Kenyan and Ethiopian rifts, interact. We show that the Ririba rift formed from the southward propagation of the Ethiopian rift during the Pliocene but this propagation was short-lived and aborted close to the Pliocene-Pleistocene boundary. Seismicity data support the abandonment of laterally offset, overlapping tips of the Ethiopian and Kenyan rifts. Integration with new numerical models indicates that rift abandonment resulted from progressive focusing of the tectonic and magmatic activity into an oblique, throughgoing rift zone of near pure extension directly connecting the rift sectors.
Accurate weather observations are the keystone to many quantitative applications, such as precipitation monitoring and nowcasting, hydrological modelling and forecasting, climate studies, as well as understanding precipitation-driven natural hazards (i.e. floods, landslides, debris flow). Weather radars have been an increasingly popular tool since the 1940s to provide high spatial and temporal resolution precipitation data at the mesoscale, bridging the gap between synoptic and point scale observations. Yet, many institutions still struggle to tap the potential of the large archives of reflectivity, as there is still much to understand about factors that contribute to measurement errors, one of which is calibration. Calibration represents a substantial source of uncertainty in quantitative precipitation estimation (QPE). A miscalibration of a few dBZ can easily deteriorate the accuracy of precipitation estimates by an order of magnitude. Instances where rain cells carrying torrential rains are misidentified by the radar as moderate rain could mean the difference between a timely warning and a devastating flood.
Since 2012, the Philippine Atmospheric, Geophysical, and Astronomical Services Administration (PAGASA) has been expanding the country’s ground radar network. We had a first look into the dataset from one of the longest running radars (the Subic radar) after devastating week-long torrential rains and thunderstorms in August 2012 caused by the annual southwestmonsoon and enhanced by the north-passing Typhoon Haikui. The analysis of the rainfall spatial distribution revealed the added value of radar-based QPE in comparison to interpolated rain gauge observations. However, when compared with local gauge measurements, severe miscalibration of the Subic radar was found. As a consequence, the radar-based QPE would have underestimated the rainfall amount by up to 60% if they had not been adjusted by rain gauge observations—a technique that is not only affected by other uncertainties, but which is also not feasible in other regions of the country with very sparse rain gauge coverage.
Relative calibration techniques, or the assessment of bias from the reflectivity of two radars, has been steadily gaining popularity. Previous studies have demonstrated that reflectivity observations from the Tropical Rainfall Measuring Mission (TRMM) and its successor, the Global Precipitation Measurement (GPM), are accurate enough to serve as a calibration reference for ground radars over low-to-mid-latitudes (± 35 deg for TRMM; ± 65 deg for GPM). Comparing spaceborne radars (SR) and ground radars (GR) requires cautious consideration of differences in measurement geometry and instrument specifications, as well as temporal coincidence. For this purpose, we implement a 3-D volume matching method developed by Schwaller and Morris (2011) and extended by Warren et al. (2018) to 5 years worth of observations from the Subic radar. In this method, only the volumetric intersections of the SR and GR beams are considered.
Calibration bias affects reflectivity observations homogeneously across the entire radar domain. Yet, other sources of systematic measurement errors are highly heterogeneous in space, and can either enhance or balance the bias introduced by miscalibration. In order to account for such heterogeneous errors, and thus isolate the calibration bias, we assign a quality index to each matching SR–GR volume, and thus compute the GR calibration bias as a qualityweighted average of reflectivity differences in any sample of matching SR–GR volumes. We exemplify the idea of quality-weighted averaging by using beam blockage fraction (BBF) as a quality variable. Quality-weighted averaging is able to increase the consistency of SR and GR observations by decreasing the standard deviation of the SR–GR differences, and thus increasing the precision of the bias estimates.
To extend this framework further, the SR–GR quality-weighted bias estimation is applied to the neighboring Tagaytay radar, but this time focusing on path-integrated attenuation (PIA) as the source of uncertainty. Tagaytay is a C-band radar operating at a lower wavelength and is therefore more affected by attenuation. Applying the same method used for the Subic radar, a time series of calibration bias is also established for the Tagaytay radar.
Tagaytay radar sits at a higher altitude than the Subic radar and is surrounded by a gentler terrain, so beam blockage is negligible, especially in the overlapping region. Conversely, Subic radar is largely affected by beam blockage in the overlapping region, but being an SBand radar, attenuation is considered negligible. This coincidentally independent uncertainty contributions of each radar in the region of overlap provides an ideal environment to experiment with the different scenarios of quality filtering when comparing reflectivities from the two ground radars. The standard deviation of the GR–GR differences already decreases if we consider either BBF or PIA to compute the quality index and thus the weights. However, combining them multiplicatively resulted in the largest decrease in standard deviation, suggesting that taking both factors into account increases the consistency between the matched samples.
The overlap between the two radars and the instances of the SR passing over the two radars at the same time allows for verification of the SR–GR quality-weighted bias estimation method. In this regard, the consistency between the two ground radars is analyzed before and after bias correction is applied. For cases when all three radars are coincident during a significant rainfall event, the correction of GR reflectivities with calibration bias estimates from SR overpasses dramatically improves the consistency between the two ground radars which have shown incoherent observations before correction. We also show that for cases where adequate SR coverage is unavailable, interpolating the calibration biases using a moving average can be used to correct the GR observations for any point in time to some extent. By using the interpolated biases to correct GR observations, we demonstrate that bias correction reduces the absolute value of the mean difference in most cases, and therefore improves the consistency between the two ground radars.
This thesis demonstrates that in general, taking into account systematic sources of uncertainty that are heterogeneous in space (e.g. BBF) and time (e.g. PIA) allows for a more consistent estimation of calibration bias, a homogeneous quantity. The bias still exhibits an unexpected variability in time, which hints that there are still other sources of errors that remain unexplored. Nevertheless, the increase in consistency between SR and GR as well as between the two ground radars, suggests that considering BBF and PIA in a weighted-averaging approach is a step in the right direction.
Despite the ample room for improvement, the approach that combines volume matching between radars (either SR–GR or GR–GR) and quality-weighted comparison is readily available for application or further scrutiny. As a step towards reproducibility and transparency in atmospheric science, the 3D matching procedure and the analysis workflows as well as sample data are made available in public repositories. Open-source software such as Python and wradlib are used for all radar data processing in this thesis. This approach towards open science provides both research institutions and weather services with a valuable tool that can be applied to radar calibration, from monitoring to a posteriori correction of archived data.
Study region: Tisza and Prut catchments, originating on the slopes of the Carpathian mountains. Study focus: The study reported here investigates (i) climate change impacts on flood risk in the region, and (ii) uncertainty related to hydrological modelling, downscaling techniques and climate projections. The climate projections used in the study were derived from five GCMs, downscaled either dynamically with RCMs or with the statistical downscaling model XDS. The resulting climate change scenarios were applied to drive the eco-hydrological model SWIM, which was calibrated and validated for the catchments in advance using observed climate and hydrological data. The changes in the 30-year flood hazards and 98 and 95 percentiles of discharge were evaluated for the far future period (2071-2100) in comparison with the reference period (1981-2010). New hydrological insights for the region: The majority of model outputs under RCP 4.5 show a small to strong increase of the 30-year flood level in the Tisza ranging from 4.5% to 62%, and moderate increase in the Prut ranging from 11% to 22%. The impact results under RCP 8.5 are more uncertain with changes in both directions due to high uncertainties in GCM-RCM climate projections, downscaling methods and the low density of available climate stations.
Spinning up large-scale coupled surface-subsurface numerical models can be a time and resource consuming task. If an uninformed initial condition is chosen, the spin-up can easily require 20 years of repeated simulations on high-performance computing machines. In this paper we compare the classical approach of starting from a fixed shallow depth to groundwater (here 3 m) with three more informed approaches for the definition of initial conditions in the spin up. In the first of these three approaches, we start from a known-steady state groundwater table, calculated with a 2-D groundwater model and the yearly net recharge, and combine it with an unsaturated zone that assumes hydrostatic conditions. In the second approach, we start from the same groundwater table combined with vertical profiles in the unsaturated zone with uniform vertical flow identical to the groundwater recharge. In the third approach we calculate a dynamic steady state from a simplified subsurface model combining a transient 2-D groundwater model with a limited number of 1-D transient unsaturated zone columns on top. Results for spinning-up a 3-D Parflow-CLM model using the different initial conditions show that large gains can be made by considering states in groundwater and the vadose zone that are consistent, i.e. where groundwater recharge and the vertical flux in the vadose zone agree. By this, the spin-up time was reduced from about 10 years to about 3 years of simulated time. In the light of seasonal fluctuations of net recharge, using the transient approach showed more stable results.
The Postmasburg Manganese Field (PMF), Northern Cape Province, South Africa, once represented one of the largest sources of manganese ore worldwide. Two belts of manganese ore deposits have been distinguished in the PMF, namely the Western Belt of ferruginous manganese ores and the Eastern Belt of siliceous manganese ores. Prevailing models of ore formation in these two belts invoke karstification of manganese-rich dolomites and residual accumulation of manganese wad which later underwent diagenetic and low-grade metamorphic processes. For the most part, the role of hydrothermal processes and metasomatic alteration towards ore formation has not been adequately discussed. Here we report an abundance of common and some rare Al-, Na-, K- and Ba-bearing minerals, particularly aegirine, albite, microcline, banalsite, serandite-pectolite, paragonite and natrolite in Mn ores of the PMF, indicative of hydrothermal influence. Enrichments in Na, K and/or Ba in the ores are generally on a percentage level for most samples analysed through bulk-rock techniques. The presence of As-rich tokyoite also suggests the presence of As and V in the hydrothermal fluid. The fluid was likely oxidized and alkaline in nature, akin to a mature basinal brine. Various replacement textures, particularly of Na- and K- rich minerals by Ba-bearing phases, suggest sequential deposition of gangue as well as ore-minerals from the hydrothermal fluid, with Ba phases being deposited at a later stage. The stratigraphic variability of the studied ores and their deviation from the strict classification of ferruginous and siliceous ores in the literature, suggests that a re-evaluation of genetic models is warranted. New Ar-Ar ages for K-feldspars suggest a late Neoproterozoic timing for hydrothermal activity. This corroborates previous geochronological evidence for regional hydrothermal activity that affected Mn ores at the PMF but also, possibly, the high-grade Mn ores of the Kalahari Manganese Field to the north. A revised, all-encompassing model for the development of the manganese deposits of the PMF is then proposed, whereby the source of metals is attributed to underlying carbonate rocks beyond the Reivilo Formation of the Campbellrand Subgroup. The main process by which metals are primarily accumulated is attributed to karstification of the dolomitic substrate. The overlying Asbestos Hills Subgroup banded iron formation (BIF) is suggested as a potential source of alkali metals, which also provides a mechanism for leaching of these BIFs to form high-grade residual iron ore deposits.
Large earthquakes initiate chains of surface processes that last much longer than the brief moments of strong shaking. Most moderate‐ and large‐magnitude earthquakes trigger landslides, ranging from small failures in the soil cover to massive, devastating rock avalanches. Some landslides dam rivers and impound lakes, which can collapse days to centuries later, and flood mountain valleys for hundreds of kilometers downstream. Landslide deposits on slopes can remobilize during heavy rainfall and evolve into debris flows. Cracks and fractures can form and widen on mountain crests and flanks, promoting increased frequency of landslides that lasts for decades. More gradual impacts involve the flushing of excess debris downstream by rivers, which can generate bank erosion and floodplain accretion as well as channel avulsions that affect flooding frequency, settlements, ecosystems, and infrastructure. Ultimately, earthquake sequences and their geomorphic consequences alter mountain landscapes over both human and geologic time scales. Two recent events have attracted intense research into earthquake‐induced landslides and their consequences: the magnitude M 7.6 Chi‐Chi, Taiwan earthquake of 1999, and the M 7.9 Wenchuan, China earthquake of 2008. Using data and insights from these and several other earthquakes, we analyze how such events initiate processes that change mountain landscapes, highlight research gaps, and suggest pathways toward a more complete understanding of the seismic effects on the Earth's surface.
Speleothems, secondary cave carbonates, are important tools for climate reconstruction, especially as they often can be very precisely dated with the UTh method. If the U-Th method fails, dating becomes difficult, and often results in abandonment of samples and study sites. Radiocarbon dating is the only other radiometric dating technique applicable to the last similar to 50 ka, but presents complexities related to temporal variability of the reservoir effect in speleothems. Thus, radiocarbon dating of speleothems is not straightforward, and there are currently no publicly available tools to define proper age-depth relationships with this method. Here, we present an improved version of a previously published radiocarbon based age-depth modelling approach (star, Lechleitner et al., 2016b), which is now made freely available. The software is easy to use and provides the possibility to obtain reliable age-depth relationships, without prior knowledge of reservoir effects and their variability. In addition, star is able to detect and handle growth stops and phases with different growth rates. We test star on artificially constructed data sets and illustrate steps to improve the model performance. Furthermore, we apply the new approach to published radiocarbon data of U-Th dated stalagmites. This offers the possibility to investigate the strengths and weaknesses of the new approach with respect to potentially significant long term trends in the radiocarbon reservoir effect, which might otherwise remain undetected. In summary, we have produced a valuable software, which easily enables to construct age-depth relationships on the basis of reservoir effect disturbed radiocarbon measurements.
During the earliest Triassic microbial mats flourished in the photic zones of marginal seas, generating widespread microbialites. It has been suggested that anoxic conditions in shallow marine environments, linked to the end-Permian mass extinction, limited mat-inhibiting metazoans allowing for this microbialite expansion. The presence of a diverse suite of proxies indicating oxygenated shallow sea-water conditions (metazoan fossils, biomarkers and redox proxies) from microbialite successions have, however, challenged the inference of anoxic conditions. Here, the distribution and faunal composition of Griesbachian microbialites from China, Iran, Turkey, Armenia, Slovenia and Hungary are investigated to determine the factors that allowed microbialite-forming microbial mats to flourish following the end-Permian crisis. The results presented here show that Neotethyan microbial buildups record a unique faunal association due to the presence of keratose sponges, while the Palaeotethyan buildups have a higher proportion of molluscs and the foraminifera Earlandia. The distribution of the faunal components within the microbial fabrics suggests that, except for the keratose sponges and some microconchids, most of the metazoans were transported into the microbial framework via wave currents. The presence of both microbialites and metazoan associations were limited to oxygenated settings, suggesting that a factor other than anoxia resulted in a relaxation of ecological constraints following the mass extinction event. It is inferred that the end-Permian mass extinction event decreased the diversity and abundance of metazoans to the point of significantly reducing competition, allowing photosynthesis-based microbial mats to flourish in shallow water settings and resulting in the formation of widespread microbialites.
Influence of the Main Border Faults on the 3D Hydraulic Field of the Central Upper Rhine Graben
(2019)
The Upper Rhine Graben (URG) is an active rift with a high geothermal potential. Despite being a well-studied area, the three-dimensional interaction of the main controlling factors of the thermal and hydraulic regime is still not fully understood. Therefore, we have used a data-based 3D structural model of the lithological configuration of the central URG for some conceptual numerical experiments of 3D coupled simulations of fluid and heat transport. To assess the influence of the main faults bordering the graben on the hydraulic and the deep thermal field, we carried out a sensitivity analysis on fault width and permeability. Depending on the assigned width and permeability of the main border faults, fluid velocity and temperatures are affected only in the direct proximity of the respective border faults. Hence, the hydraulic characteristics of these major faults do not significantly influence the graben-wide groundwater flow patterns. Instead, the different scenarios tested provide a consistent image of the main characteristics of fluid and heat transport as they have in common: (1) a topography-driven basin-wide fluid flow perpendicular to the rift axis from the graben shoulders to the rift center, (2) a N/NE-directed flow parallel to the rift axis in the center of the rift and, (3) a pronounced upflow of hot fluids along the rift central axis, where the streams from both sides of the rift merge. This upflow axis is predicted to occur predominantly in the center of the URG (northern and southern model area) and shifted towards the eastern boundary fault (central model area).
Southern Patagonia is a prime example of ongoing oceanic ridge collision and slab-window formation sustained over several million years. The impact of these phenomena on the thermal structure and exhumation of the crust have been mainly assessed with low-temperature thermochronology of bedrock samples. Here, we infer thermal histories from new and existing thermochronological data from the region of most recent ridge collision. In particular, we evaluate the potential far-reaching thermal effects of the evolving slab window, which have previously been considered responsible for patterns of late Miocene reheating associated with back-arc alkaline volcanism. Our model results define protracted cooling since similar to 15 Ma and stepwise exhumation since the late Miocene. The pattern of stepwise exhumation closely matches the onset of Patagonian glaciation at 7 Ma and the successive pulse of glacial incision coeval with neotectonic activity since 3-4 Ma that are also documented by independent geological and geomorphological evidence in the region. Importantly, our findings challenge the recently suggested lack of glacial erosion and incision since 5 Ma in this region. Furthermore, in contrast to previous modelling studies, we find that the available data do not evidence a previously proposed northward-propagating heating event associated with alkaline volcanism. We hypothesize that the anomalous alkaline volcanism in the Patagonian back-arc might be related to trench-orthogonal tears aligned with transform faults in the subducting plate. The substantial differences from the previous modelling procedure on some of the same samples is demonstrated to result from an important lack of convergence in model runs. (C) 2019 Elsevier B.V. All rights reserved.
Geophysical data acquisition in oceanic domains is challenging, implying measurements with low and/or nonhomogeneous spatial resolution. The evolution of satellite gravimetry and altimetry techniques allows testing 3-D density models of the lithosphere, taking advantage of the high spatial resolution and homogeneous coverage of satellites. However, it is not trivial to discretise the source of the gravity field at different depths. Here, we propose a new method for inferring tectonic boundaries at the crustal level. As a novelty, instead of modeling the gravity anomalies and assuming a flat Earth approximation, we model the vertical gravity gradients (VGG) in spherical coordinates, which are especially sensitive to density contrasts in the upper layers of the Earth. To validate the methodology, the complex oceanic domain of the Caribbean region is studied, which includes different crustal domains with a tectonic history since Late Jurassic time. After defining a lithospheric starting model constrained by up-to-date geophysical data sets, we tested several a-priory density distributions and selected the model with the minimum misfits with respect to the VGG calculated from the EIGEN-6C4 data set. Additionally, the density of the crystalline crust was inferred by inverting the VGG field. Our methodology enabled us not only to refine, confirm, and/or propose tectonic boundaries in the study area but also to identify a new anomalous buoyant body, located in the South Lesser Antilles subduction zone, and high-density bodies along the Greater, Lesser, and Leeward Antilles forearcs.
Core Ideas
3D MRI relaxation time maps reflect water mobility in root, rhizosphere, and soil.
3D NCT water content maps of the same plant complement relaxation time maps.
The relaxation time T1 decreases from soil to root, whereas water content increases.
Parameters together indicate modification of rhizosphere pore space by gel phase.
The zone of reduced T1 corresponds to the zone remaining dry after rewetting.
In situ investigations of the rhizosphere require high‐resolution imaging techniques, which allow a look into the optically opaque soil compartment. We present the novel combination of magnetic resonance imaging (MRI) and neutron computed tomography (NCT) to achieve synergistic information such as water mobility in terms of three‐dimensional (3D) relaxation time maps and total water content maps. Besides a stationary MRI scanner for relaxation time mapping, we used a transportable MRI system on site in the NCT facility to capture rhizosphere properties before desiccation and after subsequent rewetting. First, we addressed two questions using water‐filled test capillaries between 0.1 and 5 mm: which root diameters can still be detected by both methods, and to what extent are defined interfaces blurred by these imaging techniques? Going to real root system architecture, we demonstrated the sensitivity of the transportable MRI device by co‐registration with NCT and additional validation using X‐ray computed tomography. Under saturated conditions, we observed for the rhizosphere in situ a zone with shorter T1 relaxation time across a distance of about 1 mm that was not caused by reduced water content, as proven by successive NCT measurements. We conclude that the effective pore size in the pore network had changed, induced by a gel phase. After rewetting, NCT images showed a dry zone persisting while the MRI intensity inside the root increased considerably, indicating water uptake from the surrounding bulk soil through the still hydrophobic rhizosphere. Overall, combining NCT and MRI allows a more detailed analysis of the rhizosphere's functioning.
Meteoroids largely disintegrate during their entry into the atmosphere, contributing significantly to the input of cosmic material to Earth. Yet, their atmospheric entry is not well understood. Experimental studies on meteoroid material degradation in high-enthalpy facilities are scarce and when the material is recovered after testing, it rarely provides sufficient quantitative data for the validation of simulation tools. In this work, we investigate the thermochemical degradation mechanism of a meteorite in a high-enthalpy ground facility able to reproduce atmospheric entry conditions. A testing methodology involving measurement techniques previously used for the characterization of thermal protection systems for spacecraft is adapted for the investigation of ablation of alkali basalt (employed here as meteorite analog) and ordinary chondrite samples. Both materials are exposed to a cold-wall stagnation point heat flux of 1.2 MW m(-2). Numerous local pockets that formed on the surface of the samples by the emergence of gas bubbles reveal the frothing phenomenon characteristic of material degradation. Time-resolved optical emission spectroscopy data of ablated species allow us to identify the main radiating atoms and ions of potassium, calcium, magnesium, and iron. Surface temperature measurements provide maximum values of 2280 K for the basalt and 2360 K for the chondrite samples. We also develop a material response model by solving the heat conduction equation and accounting for evaporation and oxidation reaction processes in a 1D Cartesian domain. The simulation results are in good agreement with the data collected during the experiments, highlighting the importance of iron oxidation to the material degradation.
Humus forms are a distinctive morphological indicator of soil organic matter decomposition. The spatial distribution of humus forms depends on environmental factors such as topography, climate and vegetation. In montane and subalpine forests, environmental influences show a high spatial heterogeneity, which is reflected by a high spatial variability of humus forms. This study aims at examining spatial patterns of humus forms and their dependence on the spatial scale in a high mountain forest environment (Val di Sole/Val di Rabbi, Trentino, Italian Alps). On the basis of the distributions of environmental covariates across the study area, we described humus forms at the local scale (six sampling sites), slope scale (60 sampling sites) and landscape scale (30 additional sampling sites). The local variability of humus forms was analyzed with regard to the ground cover type. At the slope and landscape scale, spatial patterns of humus forms were modeled applying random forests and ordinary kriging of the model residuals. The results indicate that the occurrence of the humus form classes Mull, Mullmoder, Moder, Amphi and Eroded Moder generally depends on the topographical position. Local-scale patterns are mostly related to micro-topography (local accumulation and erosion sites) and ground cover, whereas slope-scale patterns are mainly connected with slope exposure and elevation. Patterns at the landscape scale show a rather irregular distribution, as spatial models at this scale do not account for local to slope-scale variations of humus forms. Moreover, models at the slope scale perform distinctly better than at the landscape scale. In conclusion, the results of this study highlight that landscape-scale predictions of humus forms should be accompanied by local- and slope-scale studies in order to enhance the general understanding of humus form patterns.
Land cover change is a dynamic phenomenon driven by synergetic biophysical and socioeconomic effects. It involves massive transitions from natural to less natural habitats and thereby threatens ecosystems and the services they provide. To retain intact ecosystems and reduce land cover change to a minimum of natural transition processes, a dense network of protected areas has been established across Europe. However, even protected areas and in particular the zones around protected areas have been shown to undergo land cover changes. The aim of our study was to compare land cover changes in protected areas, non-protected areas, and 1 km buffer zones around protected areas and analyse their relationship to climatic and socioeconomic factors across Europe between 2000 and 2012 based on earth observation data. We investigated land cover flows describing major change processes: urbanisation, afforestation, deforestation, intensification of agriculture, extensification of agriculture, and formation of water bodies. Based on boosted regression trees, we modelled correlations between land cover flows and climatic and socioeconomic factors. The results show that land cover changes were most frequent in 1 km buffer zones around protected areas (3.0% of all buffer areas affected). Overall, land cover changes within protected areas were less frequent than outside, although they still amounted to 18,800 km2 (1.5% of all protected areas) from 2000 to 2012. In some parts of Europe, urbanisation and intensification of agriculture still accounted for up to 25% of land cover changes within protected areas. Modelling revealed meaningful relationships between land cover changes and a combination of influencing factors. Demographic factors (accessibility to cities and population density) were most important for coarse-scale patterns of land cover changes, whereas fine-scale patterns were most related to longitude (representing the general east/west economic gradient) and latitude (representing the north/south climatic gradient).
Land cover change is a dynamic phenomenon driven by synergetic biophysical and socioeconomic effects. It involves massive transitions from natural to less natural habitats and thereby threatens ecosystems and the services they provide. To retain intact ecosystems and reduce land cover change to a minimum of natural transition processes, a dense network of protected areas has been established across Europe. However, even protected areas and in particular the zones around protected areas have been shown to undergo land cover changes. The aim of our study was to compare land cover changes in protected areas, non-protected areas, and 1 km buffer zones around protected areas and analyse their relationship to climatic and socioeconomic factors across Europe between 2000 and 2012 based on earth observation data. We investigated land cover flows describing major change processes: urbanisation, afforestation, deforestation, intensification of agriculture, extensification of agriculture, and formation of water bodies. Based on boosted regression trees, we modelled correlations between land cover flows and climatic and socioeconomic factors. The results show that land cover changes were most frequent in 1 km buffer zones around protected areas (3.0% of all buffer areas affected). Overall, land cover changes within protected areas were less frequent than outside, although they still amounted to 18,800 km2 (1.5% of all protected areas) from 2000 to 2012. In some parts of Europe, urbanisation and intensification of agriculture still accounted for up to 25% of land cover changes within protected areas. Modelling revealed meaningful relationships between land cover changes and a combination of influencing factors. Demographic factors (accessibility to cities and population density) were most important for coarse-scale patterns of land cover changes, whereas fine-scale patterns were most related to longitude (representing the general east/west economic gradient) and latitude (representing the north/south climatic gradient).
Most epiphytic bromeliads, especially those in the genus Tillandsia, lack functional roots and rely on the absorption of water and nutrients by large, multicellular trichomes on the epidermal surfaces of leaves and stems. Another important function of these structures is the spread of water over the epidermal surface by capillary action between trichome "wings" and epidermal surface. Although critical for the ultimate absorption by these plants, understanding of this function of trichomes is primarily based on light microscope observations. To better understand this phenomenon, the distribution of water was followed by its attenuation of cold neutrons following application of H2O to the cut end of Tillandsia usneoides shoots. Experiments confirmed the spread of added water on the external surfaces of this "atmospheric" epiphyte. In a morphologically and physiologically similar plant lacking epidermal trichomes, water added to the cut end of a shoot clearly moved via its internal xylem and not on its epidermis. Thus, in T. usneoides, water moves primarily by capillarity among the overlapping trichomes forming a dense indumentum on shoot surfaces, while internal vascular water movement is less likely. T. usneoides, occupying xeric microhabitats, benefits from reduction of water losses by low-shoot xylem hydraulic conductivities.
Proxy-based reconstructions and modeling of Holocene spatiotemporal precipitation patterns for China and Mongolia have hitherto yielded contradictory results indicating that the basic mechanisms behind the East Asian Summer Monsoon and its interaction with the westerly jet stream remain poorly understood. We present quantitative reconstructions of Holocene precipitation derived from 101 fossil pollen records and analyse them with the help of a minimal empirical model. We show that the westerly jet-stream axis shifted gradually southward and became less tilted since the middle Holocene. This was tracked by the summer monsoon rain band resulting in an early-Holocene precipitation maximum over most of western China, a mid-Holocene maximum in north-central and northeastern China, and a late-Holocene maximum in southeastern China. Our results suggest that a correct simulation of the orientation and position of the westerly jet stream is crucial to the reliable prediction of precipitation patterns in China and Mongolia.
Carbonate rocks, widely used for paleomagnetically quantifying the drift history of the Gondwana derived continental blocks of the Tibetan Plateau and evolution of the Paleo/Meso/Neo-Tethys Oceans, are prone to pervasive remagnetization. Identifying remagnetization is difficult because it is commonly undetectable through the classic paleomagnetic field tests. Here we apply comprehensive paleomagnetic, rock magnetic, and petrographic studies to upper Triassic limestones in the eastern Qiangtang block. Our results reveal that detrital/biogenic magnetite, which may carry the primary natural remanent magnetization (NRM), is rarely preserved in these rocks. In contrast, authigenic magnetite and hematite pseudomorphs after pyrite, and monoclinic pyrrhotite record three episodes of remagnetization. The earliest remagnetization was induced by oxidation of early diagenetic pyrite to magnetite, probably related to the collision between the northeastern Tibetan Plateau and the Qiangtang block after closure of the Paleo-Tethys Ocean in the Late Triassic. The second remagnetization, residing in hematite and minor goethite, which is the further subsurface oxidation product of pyrite/magnetite, is possibly related to the development of the localized Cenozoic basins soon after India-Asia collision in the Paleocene. The youngest remagnetization is a combination of thermoviscous and chemical remanent magnetization carried by authigenic magnetite and pyrrhotite, respectively. Our analyses suggest that a high supply of organic carbon during carbonate deposition, prevailing sulfate reducing conditions during early diagenesis, and widespread orogenic fluid migration related to crustal shortening during later diagenesis, have altered the primary remanence of the shallow-water Tethyan carbonate rocks of the Tibetan Plateau. We emphasize that all paleomagnetic results from these rocks must be carefully examined for remagnetization before being used for paleogeographic reconstructions. Future paleomagnetic investigations of the carbonate rocks in orogenic belts should be accompanied by thorough rock magnetic and petrographic studies to determine the origin of the NRM. (C) 2019 Elsevier B.V. All rights reserved.
Yukon’s Beaufort coast, Canada, is a highly dynamic landscape. Cultural sites, infrastructure, and travel routes used by the local population are particularly vulnerable to coastal erosion. To assess threats to these phenomena, rates of shoreline change for a 210 km length of the coast were analyzed and combined with socioeconomic and cultural information. Rates of shoreline change were derived from aerial and satellite imagery from the 1950s, 1970s, 1990s, and 2011. Using these data, conservative (S1) and dynamic (S2) shoreline projections were constructed to predict shoreline positions for the year 2100. The locations of cultural features in the archives of a Parks Canada database, the Yukon Archaeological Program, and as reported in other literature were combined with projected shoreline position changes. Between 2011 and 2100, approximately 850 ha (S1) and 2660 ha (S2) may erode, resulting in a loss of 45% (S1) to 61% (S2) of all cultural features by 2100. The last large, actively used camp area and two nearshore landing strips will likely be threatened by future coastal processes. Future coastal erosion and sedimentation processes are expected to increasingly threaten cultural sites and influence travelling and living along the Yukon coast.
Groundwater travel time distributions (TTDs) provide a robust description of the subsurface mixing behavior and hydrological response of a subsurface system. Lagrangian particle tracking is often used to derive the groundwater TTDs. The reliability of this approach is subjected to the uncertainty of external forcings, internal hydraulic properties, and the interplay between them. Here, we evaluate the uncertainty of catchment groundwater TTDs in an agricultural catchment using a 3-D groundwater model with an overall focus on revealing the relationship between external forcing, internal hydraulic properties, and TTD predictions. Eight recharge realizations are sampled from a high-resolution dataset of land surface fluxes and states. Calibration-constrained hydraulic conductivity fields (Ks fields) are stochastically generated using the null-space Monte Carlo (NSMC) method for each recharge realization. The random walk particle tracking (RWPT) method is used to track the pathways of particles and compute travel times. Moreover, an analytical model under the random sampling (RS) assumption is fit against the numerical solutions, serving as a reference for the mixing behavior of the model domain. The StorAge Selection (SAS) function is used to interpret the results in terms of quantifying the systematic preference for discharging young/old water. The simulation results reveal the primary effect of recharge on the predicted mean travel time (MTT). The different realizations of calibration-constrained Ks fields moderately magnify or attenuate the predicted MTTs. The analytical model does not properly replicate the numerical solution, and it underestimates the mean travel time. Simulated SAS functions indicate an overall preference for young water for all realizations. The spatial pattern of recharge controls the shape and breadth of simulated TTDs and SAS functions by changing the spatial distribution of particles' pathways. In conclusion, overlooking the spatial nonuniformity and uncertainty of input (forcing) will result in biased travel time predictions. We also highlight the worth of reliable observations in reducing predictive uncertainty and the good interpretability of SAS functions in terms of understanding catchment transport processes.
Cash vs. in-kind transfers
(2019)
Historically, India has relied on subsidizing staple food as a major instrument in improving food security. Recently, however, cash transfers have entered the debate as an alternative, as they are associated with lower market distortions, leakages and fiscal costs. This study contributes to this debate by analyzing India’s Targeted Public Distribution System (TPDS). Our main objective was to explain the under-purchase, or low take-up, from the TPDS, which is typically attributed to ‘leakage’, i.e. the diversion of food grains from eligible consumers. We provide an alternative solution based on self-targeting; while poorer households increase their consumption from the TPDS, wealthier households restrain from consuming subsidized commodities. Using a large household dataset, we estimated that such a voluntary opt-out system, based on income, would save a minimum of 6.5% of grains released through the TPDS. Besides these demand-driven aspects, our analysis indicates that poor regions perform better at lowering the diversion of grains and that large targeting errors exist among female-led households. Finally, we find substantial regional price differences that would benefit the poor and rural population under a uniform cash-transfer system that does not correct for regional price levels.
Armed conflicts trigger region-specific mechanisms that affect land use change. Deforestation is presented as one of the most common negative environmental impacts resulting from armed conflicts, with relevant consequences in terms of greenhouse gas emissions and loss of ecosystem services. However, the impact of armed conflict on forests is complex and may simultaneously lead to positive and negative environmental outcomes, i.e. forest regrowth and deforestation, in different regions even within a country. We investigate the impact that armed conflict exerted over forest dynamics at different spatial scales in Colombia and for the global tropics during the period 1992–2015. Through the analysis of its internally displaced population (departures) our results suggest that, albeit finding forest regrowth in some municipalities, the Colombian conflict predominantly exerted a negative impact on its forests. A further examination of georeferenced fighting locations in Colombia and across the globe shows that conflict areas were 8 and 4 times more likely to undergo deforestation, respectively, in the following years in relation to average deforestation rates. This study represents a municipality level, long-term spatial analysis of the diverging effects the Colombian conflict exerted over its forest dynamics over two distinct periods of increasing and decreasing conflict intensity. Moreover, it presents the first quantified estimate of conflict's negative impact on forest ecosystems across the globe. The relationship between armed conflict and land use change is of global relevance given the recent increase of armed conflicts across the world and the importance of a possible exacerbation of armed conflicts and migration as climate change impacts increase.
Soil landscape research is faced with wide-ranging questions of soil erosion, precision farming, and agricultural risk management. Digital Soil Morphometrics is a powerful tool to provide respective answers or recommendations but requires soil data from the pedon-to-field scale with high horizontal and vertical resolutions, including the subsoil. We present an efficient sampling and measurement method for easily obtainable soil driving cores with low-destructive preparation. Elemental contents and soil organic and mineral matter composition were measured rapidly and in large numbers using a multi-sensor approach, i.e., visible and near infrared (Vis-NIR), diffuse reflectance infrared Fourier transform (DRIFT), and X-ray fluorescence (XRF) spectroscopy. The suitability of the approach with respect to three-dimensional soil landscape models was tested using soils along a slope representing different stages of erosion and deposition in a hummocky landscape under arable land use (Calcaric Regosols, Calcic Luvisols, Luvic Stagnosols, Gleyic-Colluvic Regosols). The combination of soil core sampling, pedological description, and three spectroscopic techniques enabled rapid determination and interpretation of horizontal and vertical spatial distributions of soil organic carbon (SOC), soil organic and mineral matter composition, as well as CaCO3, Fe, and Mn contents. Depth profiles for SOC, CaCO3, and Fe contents were suitable indicators for site-specific degrees of erosion and matter transport processes at the pedon-to-field scale. Fe and Mn profiles helped identifying zones of reductive and oxic domains in subsoils (gleyzation). Further methodical developments should implement plant-availability of nutrients, characterization of Fe oxides, and calibration of the spectroscopic techniques to field-moist samples.
Mountains play a key role in the provision of nature’s contributions to people (NCP) worldwide that support societies’ quality of life. Simultaneously, mountains are threatened by multiple drivers of change. Due to the complex interlinkages between biodiversity, quality of life and drivers of change, research on NCP in mountains requires interdisciplinary approaches. In this study, we used the conceptual framework of the Intergovernmental Science-Policy Platform on Biodiversity and Ecosystem Services (IPBES) and the notion of NCP to determine to what extent previous research on ecosystem services in mountains has explored the different components of the IPBES conceptual framework. We conducted a systematic review of articles on ecosystem services in mountains published up to 2016 using the Web of Science and Scopus databases. Descriptive statistical and network analyses were conducted to explore the level of research on the components of the IPBES framework and their interactions. Our results show that research has gradually become more interdisciplinary by studying higher number of NCP, dimensions of quality of life, and indirect drivers of change. Yet, research focusing on biodiversity, regulating NCP and direct drivers has decreased over time. Furthermore, despite the fact that research on NCP in mountains becoming more policy-oriented over time, mainly in relation to payments for ecosystem services, institutional responses remained underexplored in the reviewed studies. Finally, we discuss the relevant knowledge gaps that should be addressed in future research in order to contribute to IPBES.
The North China Plain (NCP) is one of the most productive and intensive agricultural regions in China. High doses of mineral nitrogen (N) fertiliser, often combined with flood irrigation, are applied, resulting in N surplus, groundwater depletion and environmental pollution. The objectives of this thesis were to use the HERMES model to simulate the N cycle in winter wheat (Triticum aestivum L.)–summer maize (Zea mays L.) double crop rotations and show the performance of the HERMES model, of the new ammonia volatilisation sub-module and of the new nitrification inhibition tool in the NCP. Further objectives were to assess the models potential to save N and water on plot and county scale, as well as on short and long-term. Additionally, improved management strategies with the help of a model-based nitrogen fertiliser recommendation (NFR) and adapted irrigation, should be found.
Results showed that the HERMES model performed well under growing conditions of the NCP and was able to describe the relevant processes related to soil–plant interactions concerning N and water during a 2.5 year field experiment. No differences in grain yield between the real-time model-based NFR and the other treatments of the experiments on plot scale in Quzhou County could be found. Simulations with increasing amounts of irrigation resulted in significantly higher N leaching, higher N requirements of the NFR and reduced yields. Thus, conventional flood irrigation as currently practised by the farmers bears great uncertainties and exact irrigation amounts should be known for future simulation studies. In the best-practice scenario simulation on plot-scale, N input and N leaching, but also irrigation water could be reduced strongly within 2 years. Thus, the model-based NFR in combination with adapted irrigation had the highest potential to reduce nitrate leaching, compared to farmers practice and mineral N (Nmin)-reduced treatments. Also the calibrated and validated ammonia volatilisation sub-module of the HERMES model worked well under the climatic and soil conditions of northern China. Simple ammonia volatilisation approaches gave also satisfying results compared to process-oriented approaches. During the simulation with Ammonium sulphate Nitrate with nitrification inhibitor (ASNDMPP) ammonia volatilisation was higher than in the simulation without nitrification inhibitor, while the result for nitrate leaching was the opposite. Although nitrification worked well in the model, nitrification-born nitrous oxide emissions should be considered in future. Results of the simulated annual long-term (31 years) N losses in whole Quzhou County in Hebei Province were 296.8 kg N ha−1 under common farmers practice treatment and 101.7 kg N ha−1 under optimised treatment including NFR and automated irrigation (OPTai). Spatial differences in simulated N losses throughout Quzhou County, could only be found due to different N inputs. Simulations of an optimised treatment, could save on average more than 260 kg N ha−1a−1 from fertiliser input and 190 kg N ha−1a−1 from N losses and around 115.7 mm a−1 of water, compared to farmers practice. These long-term simulation results showed lower N and water saving potential, compared to short-term simulations and underline the necessity of long-term simulations to overcome the effect of high initial N stocks in soil.
Additionally, the OPTai worked best on clay loam soil except for a high simulated denitrification loss, while the simulations using farmers practice irrigation could not match the actual water needs resulting in yield decline, especially for winter wheat. Thus, a precise adaption of management to actual weather conditions and plant growth needs is necessary for future simulations. However, the optimised treatments did not seem to be able to maintain the soil organic matter pools, even with full crop residue input. Extra organic inputs seem to be required to maintain soil quality in the optimised treatments.
HERMES is a relatively simple model, with regard to data input requirements, to simulate the N cycle. It can offer interpretation of management options on plot, on county and regional scale for extension and research staff. Also in combination with other N and water saving methods the model promises to be a useful tool.
Arctic lowlands are characterized by large numbers of small waterbodies, which are known to affect surface energy budgets and the global carbon cycle. Statistical analysis of their size distributions has been hindered by the shortage of observations at sufficiently high spatial resolutions. This situation has now changed with the high-resolution (<5 m) circum-Arctic Permafrost Region Pond and Lake (PeRL) database recently becoming available. We have used this database to make the first consistent, high-resolution estimation of Arctic waterbody size distributions, with surface areas ranging from 0.0001 km(2) (100 m(2)) to 1 km(2). We found that the size distributions varied greatly across the thirty study regions investigated and that there was no single universal size distribution function (including power-law distribution functions) appropriate across all of the study regions. We did, however, find close relationships between the statistical moments (mean, variance, and skewness) of the waterbody size distributions from different study regions. Specifically, we found that the spatial variance increased linearly with mean waterbody size (R-2 = 0.97, p < 2.2e-16) and that the skewness decreased approximately hyperbolically. We have demonstrated that these relationships (1) hold across the 30 Arctic study regions covering a variety of (bio)climatic and permafrost zones, (2) hold over time in two of these study regions for which multi-decadal satellite imagery is available, and (3) can be reproduced by simulating rising water levels in a high-resolution digital elevation model. The consistent spatial and temporal relationships between the statistical moments of the waterbody size distributions underscore the dominance of topographic controls in lowland permafrost areas. These results provide motivation for further analyses of the factors involved in waterbody development and spatial distribution and for investigations into the possibility of using statistical moments to predict future hydrologic dynamics in the Arctic.
Fine particles or sediments are one of the important variables that should be considered for the proper management of water quality and aquatic ecosystems. In the present study, the effect of catchment characteristics on the performance of an already developed model for the estimation of fine sediments dynamics between the water column and sediment bed was tested, using 13 catchments distributed worldwide. The model was calibrated to determine two optimal model parameters. The first is the filtration parameter, which represents the filtration of fine sediments through pores of the stream bed during the recession period of a flood event. The second parameter is the bed erosion parameter that represents the active layer, directly related to the re-suspension of fine sediments during a flood event. A dependency of the filtration parameter with the catchment area was observed in catchments smaller than 100 km(2), whereas no particular relationship was observed for larger catchments (>100 km(2)). In contrast, the bed erosion parameter does not show a noticeable dependency with the area or other environmental characteristics. The model estimated the mass of fine sediments released from the sediment bed to the water column during flood events in the 13 catchments within 23% bias.
Through changes in policy and practice, the inherent intent of the ecosystem services (ES) concept is to safeguard ecosystems for human wellbeing. While impact is intrinsic to the concept, little is known about how and whether ES science leads to impact. Evidence of impact is needed. Given the lack of consensus on what constitutes impact, we differentiate between attributional impacts (transitional impacts on policy, practice, awareness or other drivers) and consequential impacts (real, on-the-ground impacts on biodiversity, ES, ecosystem functions and human wellbeing) impacts. We conduct rigorous statistical analyses on three extensive databases for evidence of attributional impact (the form most prevalently reported): the IPBES catalogue (n = 102), the Lautenbach systematic review (n = 504) and a 5-year in-depth survey of the OPERAs Exemplars (n = 13). To understand the drivers of impacts, we statistically analyse associations between study characteristics and impacts. Our findings show that there exists much confusion with regard to defining ES science impacts, and that evidence of attributional impact is scarce: only 25% of the IPBES assessments self-reported impact (7% with evidence); in our meta-analysis of Lautenbach’s systematic review, 33% of studies provided recommendations indicating intent of impacts. Systematic impact reporting was imposed by design on the OPERAs Exemplars: 100% reported impacts, suggesting the importance of formal impact reporting. The generalised linear models and correlations between study characteristics and attributional impact dimensions highlight four characteristics as minimum baseline for impact: study robustness, integration of policy instruments into study design, stakeholder involvement and type of stakeholders involved. Further in depth examination of the OPERAs Exemplars showed that study characteristics associated with impact on awareness and practice differ from those associated with impact on policy: to achieve impact along specific dimensions, bespoke study designs are recommended. These results inform targeted recommendations for ES science to break its impact glass ceiling.
Soil degradation by water is a serious environmental problem worldwide, with specific climatic factors being the major causes. We investigated the relationships between synoptic atmospheric patterns (i.e. weather types, WTs) and runoff, erosion and sediment yield throughout the Mediterranean basin by analyzing a large database of natural rainfall events at 68 research sites in 9 countries. Principal Component Analysis (PCA) was used to identify spatial relationships of the different WTs including three hydro-sedimentary variables: rainfall, runoff, and sediment yield (SY, used to refer to both soil erosion measured at plot scale and sediment yield registered at catchment scale). The results indicated 4 spatial classes of rainfall and runoff: (a) northern sites dependent on North (N) and North West (NW) flows; (b) eastern sites dependent on E and NE flows; (c) southern sites dependent on S and SE flows; and, finally, (d) western sites dependent on W and SW flows. Conversely, three spatial classes are identified for SY characterized by: (a) N and NE flows in northern sites (b) E flows in eastern sites, and (c) W and SW flows in western sites. Most of the rainfall, runoff and SY occurred during a small number of daily events, and just a few WTs accounted for large percentages of the total. Our results confirm that characterization by WT improves understanding of the general conditions under which runoff and SY occur, and provides useful information for understanding the spatial variability of runoff, and SY throughout the Mediterranean basin. The approach used here could be useful to aid of the design of regional water management and soil conservation measures.
Fire and grazing shape biodiversity in savannah landscapes. In land use management, knowing the effects of fire and grazing on biodiversity are important in order to ensure environmental sustainability. Beetles specifically are indicators of the biodiversity response to fire and grazing. A grazing exclusion and burning experiment in a split-plot design was used in order to investigate the interacting effects of fire and wildlife grazing on biomass, diversity, and species composition of darkling beetles (Coleoptera, Tenebrionidae) over time after fire. Darkling beetle species richness and diversity were responding in a three-way-interaction to fire, grazing, and time after fire, whereby biomass of darkling beetles remained unaffected and species compositional changes were attributed to seasonal changes of time only. Fire on ungrazed plots had a negative effect on species diversity and richness 2 weeks and 6 months post fire, whereas fire on grazed plots had no impact on species diversity and richness. Grazing only lowered species diversity and richness 6 months after fire treatments. Results suggest that grazing overrides the effects of fire and that the similar effects caused by fire and grazing are due to niche and assemblage simplification of the habitat.