Coupled Process Models as a Tool for Analysing Hydrothermal Systems

Hydrothermal systems are characterised by complex interactions between heat transfer, fluid flow, deformation, species transport and chemical reactions. Numerical models can provide quantitatively constrained information in regions where acquisition of new data is difficult or expensive thus providing a means for reducing risks, costs, and effort during targeting, production, and management of resources linked to hydrothermal systems. Here we show how numerical simulations of hydrothermal processes can be used to better understand coupled reactive transport in modern geothermal systems and in ancient hydrothermal ore deposits. We give examples based on the Enhanced Geothermal System at Soultz-sous-Forêts in France, hydrothermal mineralisation at Mount Isa in Australia, and the geothermal resource at Hamburg-Allermöhe in Germany.     

3D finite volume groundwater and heat transport modeling with non-orthogonal grids,using a coordinate transformation method

Many popular groundwater modeling codes are based on the finite differences or finite volume method for orthogonal grids. In cases of complex subsurface geometries this type of grid either leads to coarse geometric representations or to extremely fine meshes. We use a coordinate transformation method (CTM) to circumvent this shortcoming. In computational fluid dynamics (CFD), this method has been applied successfully to the general Navier–Stokes equation. The method is based on tensor analysis and performs a transformation of a curvilinear into a rectangular unit grid, on which a modified formulation of the differential equations is applied. Therefore, it is not necessary to reformulate the code in total. We applied the CTM to an existing three-dimensional code (SHEMAT), a simulator for heat conduction and advection in porous media. The finite volume discretization scheme for the non-orthogonal, structured, hexahedral grid leads to a 19-point stencil and a correspondingly banded system matrix. The implementation is straightforward and it is possible to use some existing routines without modification. The accuracy of the modified code is demonstrated for single phase flow on a two-dimensional analytical solution for flow and heat transport. Additionally, a simple case of potential flow is shown for a two-dimensional grid which is increasingly deformed. The result reveals that the corresponding error increases only slightly. Finally, a thermal free-convection benchmark is discussed. The result shows, that the solution obtained with the new code is in good agreement with the ones obtained by other codes.     

A weighing lysimeter in a regenerating eucalypt forest: Design,construction, and performance

A 36 tonne monolithic weighing lysimeter (3.7 m diameter and 1.5 m deep) was installed in Kioloa State Forest near Batemans Bay, New South Wales, Australia, to provide a continuous record of water use from a regenerating natural eucaplyt community, five years old at installation. The resolution is equivalent to 0.05 mm of evaporation and provides hourly information for diurnal variation of forest evaporation. The sensitivity of 14 ppm compares favourably with those reported for other lysimeters, similar in magnitude, in the U.S.A. Design criteria for selecting, isolating and suspending a sample community ‘in situ’ with an electro-mechanical balance are outlined. Procedures adopted during installation are described with emphasis placed on precautions taken to minimize disturbance and to avert possible changes in character of the isolated natural community. Lysimeter performance was assessed both in terms of its accuracy and its operation as a remote facility with infrequent attention.     

Validating quantitative precipitation forecast for the Flood Meteorological Office,Patna region during 2011–2014

In order to issue an accurate warning for flood, a better or appropriate quantitative forecasting of precipitation is required. In view of this, the present study intends to validate the quantitative precipitation forecast (QPF) issued during southwest monsoon season for six river catchments (basin) under the flood meteorological office, Patna region. The forecast is analysed statistically by computing various skill scores of six different precipitation ranges during the years 2011–2014. The analysis of QPF validation indicates that the multi-model ensemble (MME) based forecasting is more reliable in the precipitation ranges of 1–10 and 11–25 mm. However, the reliability decreases for higher ranges of rainfall and also for the lowest range, i.e., below 1 mm. In order to testify synoptic analogue method based MME forecasting for QPF during an extreme weather event, a case study of tropical cyclone Phailin is performed. It is realized that in case of extreme events like cyclonic storms, the MME forecasting is qualitatively useful for issue of warning for the occurrence of floods, though it may not be reliable for the QPF. However, QPF may be improved using satellite and radar products.     

Screening of Microalgae for Feasible Mass Production in Industrial Hypersaline Wastewater Using Disposable Bioreactors

Ten algae species were analyzed by comparing their growth in specific hypersaline industrial wastewater. It was a by‐product of fertilizer production which was released by K + S Aktiengesellschaft, Germany. Due to processing, brine water contains a high amount of salts ( 200 g L^?1). A successful algal biotechnology mainly depends on choosing and screening the adequate algae for a specific application along with the design of optimal culture conditions with comparable photo bioreactor technologies. Therefore, a high throughput screening technology was developed. In comparison to glass flasks or flat panel reactors this system was eligible for screening applications because of disposable characteristics and the equability of each culture tube. Dunaliella salina, Tetraselmis tetrathele, and Nannochloropsis salina grew in the presence of hypersaline wastewater where T. tetrathele grew best to a wastewater concentration of 75% by salt shock experiments. D. salina tolerates a wastewater level up to 80% by gradual increase. Intracellular ion contents of lyophilized algae samples were measured. They feature special transporter to either exclude ions, i.e., sodium from the cell, or to include ions like potassium and magnesium in order to secure functionality of sensitive enzymes. Under saline stress conditions these transport systems as well as metabolic pathways leading to the production of compatible osmolytes could be induced. Stress tolerance mechanisms developed in initially unstressed culture either by stepwise adaptation or by shock exposure to harsh salt condition. For this reason a feasible mass production in industrial hypersaline wastewater was possible.     

Geomorphic signatures of ore deposits — A case study from sukinda chromite and nickel complex,Orissa

The present scientific pursuit attempts to quantify various geomorphic parameters associated with the Chromite and Nickel ores of Sukinda valley and to establish geomorphic features like land forms, slope, relative relief, drainage pattern, density and frequency, lineaments and erosion surfaces. The present morphologic features of the valley were studied by taking inputs from relevant False Colour Composite, Toposheets and by undertaking ground truth checks. It is observed that both these ores are associated with parallel drainage pattern, which is prevalent in the valley, and are characterised by low drainage density and frequency. However Nickel ore exhibits a higher range in drainage density i.e. from 0.4 to 3.0 km./sq. km. in comparison to Chromite. As regards to erosion surface Cr is found to be associated with 140m and 200m end surfaces while Ni is confined to only 140 m. Though the range of slope varies for both the ore types from 0 to 5 degrees, Cr is reported to be associated mainly with foothills and Ni with valley flats.     

Joint three-dimensional inversion of coupled groundwater flow and heat transfer based on automatic differentiation: sensitivity calculation, verification, and synthetic examples

Inverse methods are useful tools not only for deriving estimates of unknown parameters of the subsurface, but also for appraisal of the thus obtained models. While not being neither the most general nor the most efficient methods, Bayesian inversion based on the calculation of the Jacobian of a given forward model can be used to evaluate many quantities useful in this process. The calculation of the Jacobian, however, is computationally expensive and, if done by divided differences, prone to truncation error. Here, automatic differentiation can be used to produce derivative code by source transformation of an existing forward model. We describe this process for a coupled fluid flow and heat transport finite difference code, which is used in a Bayesian inverse scheme to estimate thermal and hydraulic properties and boundary conditions form measured hydraulic potentials and temperatures. The resulting derivative code was validated by comparison to simple analytical solutions and divided differences. Synthetic examples from different flow regimes demonstrate the use of the inverse scheme, and its behaviour in different configurations.     

Characterising thermal water circulation in fractured bedrock using a multidisciplinary approach: a case study of St. Gorman’s Well,Ireland

A hydrogeological conceptual model of the source, circulation pathways and temporal variation of a low-enthalpy thermal spring in a fractured limestone setting is derived from a multidisciplinary approach. St. Gorman’s Well is a thermal spring in east-central Ireland with a complex and variable temperature profile (maximum of 21.8 °C). Geophysical data from a three-dimensional(3D)audio-magnetotelluric(AMT) survey are combined with time-lapse hydrogeological data and information from a previously published hydrochemical analysis to investigate the operation of this intriguing hydrothermal system. Hydrochemical analysis and time-lapse measurements suggest that the thermal waters flow within the fractured limestones of the Carboniferous Dublin Basin at all times but display variability in discharge and temperature. The 3D electrical resistivity model of the subsurface revealed two prominent structures: (1) a NW-aligned faulted contact between two limestone lithologies; and (2) a dissolutionally enhanced, N-aligned, fault of probable Cenozoic age. The intersection of these two structures, which has allowed for karstification of the limestone bedrock, has created conduits facilitating the operation of relatively deep hydrothermal circulation (likely estimated depths between 240 and 1,000 m) within the limestone succession of the Dublin Basin. The results of this study support a hypothesis that the maximum temperature and simultaneous increased discharge observed at St. Gorman’s Well each winter is the result of rapid infiltration, heating and recirculation of meteoric waters within a structurally controlled hydrothermal circulation system.

     

Freshwater flow to the San Francisco Bay‐Delta estuary over nine decades (Part 2): Change attribution

This paper explains observed trends in freshwater flow to the San Francisco Bay‐Delta estuary as reported in a companion paper (Hutton, Rath, & Roy, 2017 ). We employ a historical hydrologic record spanning nine decades and define a set of idealized flow scenarios to identify drivers of change in delta outflow and consequent salinity regime. Flow changes are measured against a baseline scenario representing 1920‐level land use and water management conditions. Additional scenarios are defined to represent the system absent state and federal water project reservoir and export operations, absent key non‐project reservoir operations, and absent historically‐observed sea level rise. These scenarios, in conjunction with the principle of superposition, are used to ascribe outflow and salinity trends to different anthropogenic and natural causes. We find that project and non‐project water management are attributed similar responsibility for decreasing outflow trends in April and May and consequent increasing spring salinity trends. In contrast, we find that increasing July and August outflow trends (and lagged decreasing salinity trends) are attributed to flow contributions from project water management; these contributions more than fully attenuate impacts associated with non‐project water management.     

Distribution of geochemical fractions of phosphorus and its ecological risk in sediment cores of a largest brackish water lake,South Asia

The current study focuses on the vertical profile of different geochemical fractions of phosphorus-loosely bound(Lo–P),aluminium bound(Al–P),iron bound(Fe–P),calcium bound(Ca–P),and organic bound phosphorus(O–P)along with ecological risk assessment of sediment cores from Chilika Lake,eastcoast of India.The percentage contribution of the different fractions to the sedimentary phosphorus in the sediment column of the whole lake are on the order:O–P(33.2%)>Ca–P(20.3%)>Fe–P(18%)>Al–P(6.7%)>Lo–P(0.35%).The Phosphorus Pollution Index(PPI)revealed the contamination of lake sediment with phosphorus.The principal component and cluster analyses highlighted the anthropogenic contribution of phosphorus.The negative loading of Ca–P with Ca points towards its origin from marine shells.The discriminate analysis showed that the variables like Ca–P,bio-available phosphorus(BAP),and pH were able to effectively discriminate the sectors in a significant manner.