Refine
Year of publication
- 2019 (80) (remove)
Document Type
- Article (55)
- Doctoral Thesis (8)
- Postprint (7)
- Other (6)
- Part of a Book (1)
- Habilitation Thesis (1)
- Report (1)
- Review (1)
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 (80) (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.