Coupling hydrogeological with surface runoff model in a Poltva case study in Western Ukraine

This paper presents the hydrological coupling of the software framework OpenGeoSys (OGS) with the EPA Storm Water Management Model (SWMM). Conceptual models include the Saint Venant equation for river flow, the 2D Darcy equations for confined and unconfined groundwater flow, a two-way hydrological coupling flux in a compartment coupling approach (conductance concept), and Lagrangian particles for solute transport in the river course. A SWMM river–OGS aquifer inter-compartment coupling flux is examined for discharging groundwater in a systematic parameter sensitivity analysis. The parameter study involves a small perturbation (first-order) sensitivity analysis and is performed for a synthetic test example base-by-base through a comprehensive range of aquifer parametrizations. Through parametrization, the test cases enables to determine the leakance parameter for simulating streambed clogging and non-ocillatory river-aquifer water exchange rates with the sequential (partitioned) coupling scheme. The implementation is further tested with a hypothetical but realistic 1D river–2D aquifer model of the Poltva catchment, where discharging groundwater in the upland area affects the river–aquifer coupling fluxes downstream in the river course (propagating feedbacks). Groundwater contribution in the moving river water is numerically determined with Lagrangian particles. A numerical experiment demonstrates that the integrated river–aquifer model is a serviceable and realistic constituent in a complete compartment model of the Poltva catchment.     

Uncertainties of geochemical modeling during CO<Subscript>2</Subscript> sequestration applying batch equilibrium calculations

One of the uncertainties in the field of carbon dioxide capture and storage (CCS) is caused by the parameterization of geochemical models. The application of geochemical models contributes significantly to calculate the fate of the CO₂ after its injection. The choice of the thermodynamic database used, the selection of the secondary mineral assemblage as well as the option to calculate pressure dependent equilibrium constants influence the CO₂ trapping potential and trapping mechanism. Scenario analyses were conducted applying a geochemical batch equilibrium model for a virtual CO₂ injection into a saline Keuper aquifer. The amount of CO₂ which could be trapped in the formation water and in the form of carbonates was calculated using the model code PHREEQC. Thereby, four thermodynamic datasets were used to calculate the thermodynamic equilibria. Furthermore, the equilibrium constants were re-calculated with the code SUPCRT92, which also applied a pressure correction to the equilibrium constants. Varying the thermodynamic database caused a range of 61% in the amount of trapped CO₂ calculated. Simultaneously, the assemblage of secondary minerals was varied, and the potential secondary minerals dawsonite and K-mica were included in several scenarios. The selection of the secondary mineral assemblage caused a range of 74% in the calculated amount of trapped CO₂. Correcting the equilibrium constants with respect to a pressure of 125 bars had an influence of 11% on the amount of trapped CO₂. This illustrates the need for incorporating sensitivity analyses into reaction pathway modeling.     

Evaluation of the water resource reproducible ability on Tarim River Basin in south of Xinjiang,northwest China

With the over-exploitation of water resources, water pollution and poor management of water infrastructures are exacerbated. Ecosystem degradation is apparent at the basin level. The Tarim River Basin in northwest China has seen intensive confrontation between environmental protection and economic development over the past five decades. Ambitious agricultural development and land reclamation projects implemented by the Chinese government in the early 1960s led to several influences. For example, the construction of dams like the Daxihaizi Reservoir disrupted the stream-flow to the lower reaches. Water resource reproducible ability (WRRA) refers to the ability of water resources to be continually added by the natural water cycle. It includes the supplementation of water quantity and the self-purification of water quality in the natural cycle of water resources. This study discusses the WRRA index and introduces the computational method for calculating the WRRA index for the Tarim River Basin. The following conclusions are observed: (1) from 1956 to 2005, the indices of WRRA in the Hotan River Basin, Yarkand and Aksu River Basin are 0.26, 0.55, 0.58, respectively, which are between 0 and 1. The results indicate that the hydrological cycle in these three sub-basins is in a reproducible state. (2) The WRRA index in the Kaidu-Kong River Basin is 1.23 > 1, which indicates that floods may occur in the Kaidu-Kong River Basin. (3) The index of WRRA in the main stream is 0, which indicates that the WRRA is very weak at this location, and zero-flow may occur. Calculating the WRRA of a basin can provide a basis for corresponding basin water resources management.     

Challenges of an integrated water resource management for the Distrito Federal, Western Central Brazil: climate, land-use and water resources

Land-use/cover change (LUCC) and climate change are major controlling factors for water resources in the Distrito Federal in Western Central Brazil. Dynamic LUCC in the region has severe impacts on water resources, while climate changes during the last three decades is thought to have only moderate effects. LUCC affects water quantity mostly during base flow conditions. River basins with substantial expansion of agriculture since the end 1970s show a dramatic decrease of base flow discharge by 40?C70%, presumably due to irrigation. In contrast, the effects of urbanization on runoff are less distinct, since factors controlling runoff generation might be more variable. For water quality, we found urban areas to have a strong influence on the parameters CSB, NH₄ ⁺, and suspended solids. In addition, we assume emerging pollutants, e.g. organic (micro)pollutants, might play a major role in the future. The project IWAS-áGUA DF focuses on creating the scientific base to face these problems in frame of an IWRM concept for the region. Results of our study will be a contribution to an IRWM concept for the Distrito Federal and will help to maintain high standards in water supply for the region.     

Electrical resistivity response due to variation in embankment shape and reservoir levels

The safety of the center-core type of fill dam structure was assessed by examining the effects of the distortion of electrical response verified in terms of two-dimensional (2D) apparent resistivity and its inverted sections from three-dimensional (3D) modeling for the embankment. The distortion effect is due to 2D interpretation of the 3D structure of the embankment. From the analysis, it was found that water level was correctly described by the resistivity section around the middle part rather than each side at the end of the embankment. This is due to the 3D terrain effect when the material of the embankment is assumed to be horizontally uniform. In addition, when the slope of the outer rock-fill section is set as uniform, the resistivity section is more similar for sharper center-core slopes. On the other hand, when the rock-fill slope is steep, the resistivity section shows the water level at a lower position than the real one, and the 3D distortion effect at the end of the embankment is enhanced.     

Heavy metal pollution in the surrounding ore deposits and mining activity: a case study from Koyulhisar (Sivas-Turkey)

The oxidation of sulfide minerals generates acidic waters containing high levels of SO₄ and Fe. The study area has active Pb?CZn?CCu mining. It is thought that the surface/subsurface/underground sulfide minerals in the region generally contribute to the acidification of groundwater. Low pH values are also responsible for dissolved metals (Al, Fe, Mn, SO₄, Pb, Zn) in the groundwater and river. Furthermore, current mining wastes have affected concentrations of trace elements in the water. High Fe and Mn concentrations are generally found in the spring which has acidic and low Eh values, while Al, Fe and Mn concentrations in the acidic waters show notable increases with the maximum values reaching 8,829, 19,084 and 1,708?ppb, respectively. These values exceed the Turkish drinking water standard of 200, 200 and 50?ppb, respectively.     

Practical evaluation of carbon sources of forest soils in Slovenia from stable and radio-carbon isotope measurements

Carbon sources were estimated by measuring carbon isotope ratios (??¹³C and ??¹⁴C) with accelerator mass spectrometry (AMS) in forest soils of different lithology. Six locations were selected in temperate deciduous and coniferous stands in Slovenia (?irovski vrh, Idrija, Ko?evski Rog, Pohorje, Gori?nica, and Rakitna), where carbonate rocks consisting of limestone and dolomite are abundant as underlying bedrock. Carbon isotope fractionation would not have occurred in two carbonaceous soils, since the values of both ??¹³C and ??¹⁴C changed consistently in these soils after thermal (550°C, 2?h) or chemical (1?M HCl, 24?h) treatments. Organic components were found to be predominant carbon sources (70?C100%) in the uppermost portions (0?C2?cm in depth). In deeper portions at a depth of about 30?C35?cm, soil carbon may be derived completely from underlying carbonate minerals in Idria, western part of Slovenia. The Combination of heat and chemical treatments with AMS provides practical information on soil carbon sources in carbonaceous soils.     

Numerical analysis of thermal effects during carbon dioxide injection with enhanced gas recovery: a theoretical case study for the Altmark gas field

Prediction about reservoir temperature change during carbon dioxide injection requires consideration of all, often subtle, thermal effects. In particular, Joule?CThomson cooling (JTC) and the viscous heat dissipation (VHD) effect are factors that cause flowing fluid temperature to differ from the static formation temperature. In this work, warm-back behavior (thermal recovery after injection completed), as well as JTC and VHD effects, at a multi-layered depleted gas reservoir are demonstrated numerically. OpenGeoSys (OGS) is able to solve coupled partial differential equations for pressure, temperature and mole-fraction of each component of the mixture with a combination of monolithic and staggered approaches. The Galerkin finite element approach is adapted for space discretization of governing equations, whereas for temporal discretization, a generalized implicit single-step scheme is used. For numerical modeling of warm-back behavior, we chose a simplified test case of carbon dioxide injection. This test case is numerically solved by using OGS and FeFlow simulators independently. OGS differs from FeFlow in the capability of representing multi-componential effects on warm-back behavior. We verify both code results by showing the close comparison of shut-in temperature profiles along the injection well. As the JTC cooling rate is inversely proportional to the volumetric heat capacity of the solid matrix, the injection layers are cooled faster as compared to the non-injection layers. The shut-in temperature profiles are showing a significant change in reservoir temperature; hence it is important to account for thermal effects in injection monitoring.     

Europe’s longest-operating on-shore CO2 storage site at Ketzin, Germany: a progress report after three years of injection

The Ketzin pilot site, led by the GFZ German Research Centre for Geosciences, is Europe??s longest-operating on-shore CO₂ storage site with the aim of increasing the understanding of geological storage of CO₂ in saline aquifers. Located near Berlin, the Ketzin pilot site is an in situ laboratory for CO₂ storage in an anticlinal structure in the Northeast German Basin. Starting research within the framework of the EU project CO₂SINK in 2004, Ketzin is Germany??s first CO₂ storage site and fully in use since the injection began in June 2008. After 39?months of operation, about 53,000 tonnes of CO₂ have been stored in 630?C650?m deep sandstone units of the Upper Triassic Stuttgart Formation. An extensive monitoring program integrates geological, geophysical and geochemical investigations at Ketzin for a comprehensive characterization of the reservoir and the CO₂ migration at various scales. Integrating a unique field and laboratory data set, both static geological modeling and dynamic simulations are regularly updated. The Ketzin project successfully demonstrates CO₂ storage in a saline aquifer on a research scale. The results of monitoring and modeling can be summarized as follows: (1) Since the start of the CO₂ injection in June 2008, the operation has been running reliably and safely. (2) Downhole pressure data prove correlation between the injection rate and the reservoir pressure and indicates the presence of an overall dynamic equilibrium within the reservoir. (3) The extensive geochemical and geophysical monitoring program is capable of detecting CO₂ on different scales and gives no indication for any leakage. (4) Numerical simulations (history matching) are in good agreement with the monitoring results.