Zinc deposits and related mineralization of the Burketown mineral field, including the world-class Century Deposit, Northern Australia: Fluid inclusion and stable isotope evidence for basin fluid sources

The stratiform Century Zn-Pb deposit and the discordant Zn-Pb lode deposits of the Burketown mineral field, northern Australia, host ore and gangue minerals with primary fluid inclusions that have not been affected by the Isan orogeny, thus providing a unique opportunity to investigate the nature of the ore-forming brines. All of the deposits are hosted in shales and siltstones belonging to the Isa superbasin and comprise sphalerite, pyrite, carbonate, quartz, galena, minor chalcopyrite, and minor illite. According to Pb model ages, the main ore stage of mineralization at Century formed at I575 Ma, some 20 m.y. after deposition of the host shale sequence. Microthermometry on undeformed, primary fluid inclusions hosted in porous sphalerite shows that the Zn at Century was transported to the deposit by a homogeneous, Ca2+- and Na+-bearing brine with a salinity of 21.6 wt percent NaCl equiv. delta D-fluid of the fluid inclusion water ranges from -89 to -83 per mil, consistent with a basinal brine that evolved from meteoric water. Fluid inclusion homogenization temperatures range between 74 degrees and 125 degrees C, which are lower than the 120 degrees to 160 degrees C range calculated from vitrinite reflectance and illite crystallinity data from the deposit. This discrepancy indicates that mineralization likely formed at 50 to 85 Mpa, corresponding to a depth of 1,900 to 3,100 m. Transgressive galena-sphalerite veins that cut stratiform mineralization at Century and breccia-filled quartz-dolomite-sphalerite-galena veins in the discordant Zn-Pb lodes have Pb model ages between 1575 and 1485 Ma. Raman spectroscopy and microthermometry reveal that the primary fluid inclusions in these veins contain Ca2+, Na+. but they have lower salinities between 23 and 10 wt percent NaCl equiv and higher delta D-fluid values ranging from -89 to -61 per mil than fluid inclusions in porous sphalerite from Century. Fluid inclusion water from sphalerite in one of the lode deposits has delta O-18(fluid) values of 1.6 and 2.4 per mil, indistinguishable from delta O-18(fluid) values between -0.3 to +7.4 per mil calculated from the isotopic composition of coexisting quartz, dolomite, and illite. The trend toward lower salinities and higher delta D-fluid values relative to the earlier mineralizing fluids is attributed to mixing between the fluid that formed Century and a seawater-derived fluid from a different source. Based on seismic data from the Lawn Hill platform and paragenetic and geochemical results from the Leichhardt River fault trough to the south, diagenetic aquifers in the Underlying Calvert superbasin appear to have been the most likely sources for the fluids that formed Century and the discordant lode deposits. Paragenetically late sphalerite and calcite cut sphalerite, quartz, and dolomite in the lode deposits and contain Na+-dominated fluid inclusions with much lower salinities than their older counterparts. The isotopic composition of calcite also indicates delta O-18(fluid) from 3.3 to 10.7 per mil, which is larger than the range obtained from synmineralization minerals, supporting the idea that a unique fluid source was involved. The absolute timing of this event is unclear, but a plethora of Pb model, K-Ar, and Ar-40/Ar-39 ages between 1440 and 1300 Ma indicate that a significant volume of fluid was mobilized at this time. The deposition of the Roper superbasin from ca. 1492 +/- 4 Ma suggests that these late veins formed from fluids that may have been derived from aquifers in overlying sediments of the Roper superbasin. Clear, buck, and drusy quartz in veins unrelated to any form of Pb-Zn mineralization record the last major fluid event in the Burketown mineral field and form distinct outcrops and ridges in the district. Fluid inclusions in these veins indicate formation from a low-salinity, 300 degrees +/- 80 degrees C fluid. Temperatures approaching 300 degrees C recorded in organic matter adjacent to faults and at sequence boundaries correspond to K-Ar ages spanning 1300 to 1100 Ma, which coincides with regional hydrothermal activity in the northern Lawn Hill platform and the emplacement of the Lakeview Dolerite at the time of assemblage of the Rodinia supercontinent.

Stochastic modelling and simulations for solute transport in porous media

A stochastic model for solute transport in aquifers is studied based on the concepts of stochastic velocity and stochastic diffusivity. By applying finite difference techniques to the spatial variables of the stochastic governing equation, a system of stiff stochastic ordinary differential equations is obtained. Both the semi-implicit Euler method and the balanced implicit method are used for solving this stochastic system. Based on the Karhunen-Loeve expansion, stochastic processes in time and space are calculated by means of a spatial correlation matrix. Four types of spatial correlation matrices are presented based on the hydraulic properties of physical parameters. Simulations with two types of correlation matrices are presented.

Indagine modellistica sulla tenuta idraulica delle difese arginali deputate alla protezione idraulica del territorio

La prima parte di questo lavoro di tesi tratta dell’interazione tra un bacino di laminazione e il sottostante acquifero: è in fase di progetto, infatti, la costruzione di una cassa di espansione sul torrente Baganza, a monte della città di Parma. L’obiettivo di tale intervento è di ridurre il rischio di esondazione immagazzinando temporaneamente, in un serbatoio artificiale, la parte più pericolosa del volume di piena che verrebbe rilasciata successivamente con portate che possono essere agevolmente contenute nel tratto cittadino del torrente. L’acquifero è stato preliminarmente indagato e monitorato permettendone la caratterizzazione litostratigrafica. La stratigrafia si può riassumere in una sequenza di strati ghiaioso-sabbiosi con successione di lenti d’argilla più o meno spesse e continue, distinguendo due acquiferi differenti (uno freatico ed uno confinato). Nel presente studio si fa riferimento al solo acquifero superficiale che è stato modellato numericamente, alle differenze finite, per mezzo del software MODFLOW_2005. L'obiettivo del presente lavoro è di rappresentare il sistema acquifero nelle condizioni attuali (in assenza di alcuna opera) e di progetto. La calibrazione è stata condotta in condizioni stazionarie utilizzando i livelli piezometrici raccolti nei punti d’osservazione durante la primavera del 2013. I valori di conducibilità idraulica sono stati stimati per mezzo di un approccio geostatistico Bayesiano. Il codice utilizzato per la stima è il bgaPEST, un software gratuito per la soluzione di problemi inversi fortemente parametrizzati, sviluppato sulla base dei protocolli del software PEST. La metodologia inversa stima il campo di conducibilità idraulica combinando osservazioni sullo stato del sistema (livelli piezometrici nel caso in esame) e informazioni a-priori sulla struttura dei parametri incogniti. La procedura inversa richiede il calcolo della sensitività di ciascuna osservazione a ciascuno dei parametri stimati; questa è stata valutata in maniera efficiente facendo ricorso ad una formulazione agli stati aggiunti del codice in avanti MODFLOW_2005_Adjoint. I risultati della metodologia sono coerenti con la natura alluvionale dell'acquifero indagato e con le informazioni raccolte nei punti di osservazione. Il modello calibrato può quindi essere utilizzato come supporto alla progettazione e gestione dell’opera di laminazione. La seconda parte di questa tesi tratta l'analisi delle sollecitazioni indotte dai percorsi di flusso preferenziali causati da fenomeni di piping all’interno dei rilevati arginali. Tali percorsi preferenziali possono essere dovuti alla presenza di gallerie scavate da animali selvatici. Questo studio è stato ispirato dal crollo del rilevato arginale del Fiume Secchia (Modena), che si è verificato in gennaio 2014 a seguito di un evento alluvionale, durante il quale il livello dell'acqua non ha mai raggiunto la sommità arginale. La commissione scientifica, la cui relazione finale fornisce i dati utilizzati per questo studio, ha attribuito, con molta probabilità, il crollo del rilevato alla presenza di tane di animali. Con lo scopo di analizzare il comportamento del rilevato in condizioni integre e in condizioni modificate dall'esistenza di un tunnel che attraversa il manufatto arginale, è stato realizzato un modello numerico 3D dell’argine mediante i noti software Femwater e Feflow. I modelli descrivono le infiltrazioni all'interno del rilevato considerando il terreno in entrambe le porzioni sature ed insature, adottando la tecnica agli elementi finiti. La tana è stata rappresentata da elementi con elevata permeabilità e porosità, i cui valori sono stati modificati al fine di valutare le diverse influenze sui flussi e sui contenuti idrici. Per valutare se le situazioni analizzate presentino o meno il verificarsi del fenomeno di erosione, sono stati calcolati i valori del fattore di sicurezza. Questo è stato valutato in differenti modi, tra cui quello recentemente proposto da Richards e Reddy (2014), che si riferisce al criterio di energia cinetica critica. In ultima analisi è stato utilizzato il modello di Bonelli (2007) per calcolare il tempo di erosione ed il tempo rimanente al collasso del rilevato.

Theory and application of inverse problems in groundwater: numerical, laboratory and field studies.

Aquifers are a vital water resource whose quality characteristics must be safeguarded or, if damaged, restored. The extent and complexity of aquifer contamination is related to characteristics of the porous medium, the influence of boundary conditions, and the biological, chemical and physical processes. After the nineties, the efforts of the scientists have been increased exponentially in order to find an efficient way for estimating the hydraulic parameters of the aquifers, and thus, recover the contaminant source position and its release history. To simplify and understand the influence of these various factors on aquifer phenomena, it is common for researchers to use numerical and controlled experiments. This work presents some of these methods, applying and comparing them on data collected during laboratory, field and numerical tests. The work is structured in four parts which present the results and the conclusions of the specific objectives.

Green technologies in the automotive industry

In most of the discussions about environmental issues and policies, transportation is highlighted as one of the main sources of pollutant emissions and energy consumption. The attention given to the automotive industry is understandable in this context due to its size, expansion, presence in our daily lives, and of course its environmental impact. If we scrutinize the “greenness” of car manufacturers we will find issues of concern from the raw material use, production processes, use, and end-of-life of vehicles. The main issues for production are high consumption of energy, raw materials, water and the waste stream, which contains the four substances of concern (cadmium, lead, hexavalent chromium, mercury). In respect of carbon emissions and energy use the use of cars is the main phase of its life-cycle due to the combination of internal combustion engines with fossil fuels. The most recent pressure is aimed at the end-of-life vehicles (ELV). In addition to the pollution from vehicle use, traffic jams and car accidents continue to be part of the downside of a car culture. Landfills sites are becoming scarce and the contamination of soil and aquifers completes the picture.

Adsorption and desorption of phosphate on limestone in experiments simulating seawater intrusion

This study investigates the potential release of from carbonate aquifers exposed to seawater intrusion. Adsorption and desorption of in the presence of deionized water (DIW) and seawater were conducted on a large block of Pleistocene age limestone to simulate the effects of seawater intrusion into a coastal carbonate aquifer at the laboratory scale. The limestone showed strong adsorption of in DIW, while adsorption was significantly less in the presence of seawater. Dissolution of CaCO3 was found to prevent adsorption at salinities less than 30 psu. Adsorption of was limited at higher salinities (30–33 psu), due to competition with ions for adsorption sites. At a salinity3 precipitated. Concentrations of between 2 and 5 μmol/L were released by desorption when the limestone was exposed to seawater. The results of this study suggest that as seawater intrudes into an originally freshwater coastal aquifer, adsorbed may be released into the groundwater. Consequently, adsorbed is expected to be released from coastal carbonate aquifers world-wide as sea level continues to rise exposing more of the freshwater aquifer to seawater.

Geochemical determination of the fate and transport of injected fresh wastewater to a deep saline aquifer

Deep well injection into non-potable saline aquifers of treated domestic wastewater has been used in Florida for decades as a safe and effective alternative to ocean outfall disposal. The objectives of this study were to determine the fate and transport of injected wastewater at two deep well injection sites in Miami Dade County, Florida, USA. Detection of ammonium in the Middle Confining units of the Floridan aquifer above the injection zone at both sites has been interpreted as evidence of upward migration of injected wastewater, posing a risk to underground sources of drinking water. Historical water quality data, including ammonia, chloride, temperature, and pH from existing monitoring wells at both sites from 1983 to 2008, major ions collected monthly from 2006 and 2008, and a synoptic sampling event for stable isotopes, tritium, and dissolved gases in 2008, were used to determine the source of ammonium in groundwater and possible migration pathways. Geochemical modeling was used to determine possible effects of injected wastewater on native water and aquifer matrix geochemistry. Injected wastewater was determined to be the source of elevated ammonium concentrations above ambient water levels, based on the results of major ion concentrations, tritium, dissolved noble gases and 15N isotopes analyses. Various possible fluid migration pathways were identified at the sites. Data for the south site suggest buoyancy-driven vertical pathways to overlying aquifers bypassing the confining units, with little mixing of injected wastewater with native water as it migrated upward. Once it is introduced into an aquifer, the injectate appeared to migrate advectively with the regional groundwater flow. Geochemical modeling indicated that CO 2-enriched injected wastewater allowed for carbonate dissolution along the vertical pathways, enhancing permeability along these flowpaths. At the north site, diffusive upward flow through the confining units or offsite vertical pathways were determined to be possible, however no evidence was detected for any on-site confining unit bypass pathway. No evidence was observed at either site of injected wastewater migration to the Upper Floridan aquifer, which is used as a municipal water supply and for aquifer storage and recovery.

Determination of vertical and horizontal pathways of injected fresh wastewater into a deep saline aquifer (Florida, USA) using natural chemical tracers

Two deep-well injection sites in south Florida, USA, inject an average of 430 million liters per day (MLD) of treated domestic fresh wastewater into a deep saline aquifer 900 m below land surface. Elevated levels of NH3 (highest concentration 939 µmol) in the overlying aquifer above ambient concentrations (concentration less than 30 µmol) were evidence of the upward migration of injected fluids. Three pathways were distinguished based on ammonium, chloride and bromide ratios, and temperature. At the South District Wastewater Treatment Plant, the tracer ratios showed that the injectate remained chemically distinct as it migrated upwards through rapid vertical pathways via density-driven buoyancy. The warmer injectate (mean 28°C) retained the temperature signal as it vertically migrated upwards; however, the temperature signal did not persist as the injectate moved horizontally into the overlying aquifers. Once introduced, the injectate moved slowly horizontally through the aquifer and mixed with ambient water. At the North District Wastewater Treatment Plant, data provide strong evidence of a one-time pulse of injectate into the overlying aquifers due to improper well construction. No evidence of rapid vertical pathways was observed at the North District Wastewater Treatment Plant.

The Influence of Hydrologic Restoration on Groundwater-Surface Water Interactions in a Karst Wetland, the Everglades (FL, USA)

Efforts to rehydrate and restore surface water flow in karst wetlands can have unintended consequences, as these highly conductive and heterogeneous aquifers create a close connection between groundwater and surface water. Recently, hydrologic restoration efforts in the karstic Taylor Slough portion of the Everglades has changed from point source delivery of canal water (direct restoration), to the use of a series of surface water recharge retention basins (diffuse restoration). To determine the influence of restoration on groundwater-surface water interactions in the Taylor Slough headwaters, a water budget was constructed for 1997–2011 using 70 hydro-meteorological stations. With diffuse restoration, groundwater seepage from the Everglades toward the urban boundary increased, while the downstream delivery of surface water to the main portion of the slough declined. The combined influence of diffuse restoration and climate led to increased intra-annual variability in the volume of groundwater and surface water in storage but supported a more seasonally hydrated wetland compared to the earlier direct tactics. The data further indicated that hydrologic engineering in karst wetland landscapes enhances groundwater-surface water interactions, even those designed for restoration purposes.

Evolução espaço-temporal da ocupação urbana sobre áreas naturais no município de Natal-RN

In this work the landscape morphodynamics was used to check the strength and importance of the changes carried out by man on the environment over time, in Natal-RN municipality. The occupation of partially preserved natural areas was analyzed, but environmentally fragile, such as riparian forests, vegetation on the banks of waterways, which play regulatory role of the water flow, and the dunes, which guarantee the rapid recharge of aquifers. The impacts of urban sprawl in Natal Southern and West zones Were identified and characterized, through a detailed mapping in the period between 1969 and 2013 the main Permanent Preservation Areas - PPA (banks of rivers and lagoons, and dunes remaining) and their temporal changes. For this were used aerial photographs and satellite imagery, altimetry data, and pre-existing information, which allowed the creation of a spatial database, and evolution of maps of impervious areas, evolution of the use and occupation and Digital Terrain Model (DTM) from contour lines with contour interval of 1 meter. Based on this study presents a diagnosis of the environmental situation and the state of conservation of natural areas, over the last 44 years, compared to human pressures. In general, it was found that the urban settlement has advanced about 60% of studied natural areas. This advance was growing by the year 2006, when there was a slowdown in the process, except for the Environmental Protection Zone (EPZ) 03, where the river Pitimbú and your PPA, which experienced a more significant loss area. The urban occupation affected the natural drainage and contributed to the contamination of groundwater Natal, due to increased sealed area, the release of liquid and solid waste, as well as the removal of riparian vegetation. Changed irreversibly the natural landscape, and reduced the quality and quantity of water resources necessary for the population. Thus, it is necessary to stimulate the adoption of use and protection of PPA planning measures, to the preservation of the San Valley Region inserted into the EPZ 01, and integrate more remaining dunes, in good condition, this EPZ, due to the importance of those remaining on the environment and the maintenance of quality of life. It is suggested, also, protection of catchment areas, such as PPA ponds and Pitimbú River. Finally, it is expected that this study can assist the managers in making decisions in urban and environmental planning of the municipality

Análise computacional de imagens de lâminas delgadas na predição dos parâmetros hidráulicos porosidade e condutividade hidráulica do aquífero Barreiras-RN

The understanding of the occurrence and flow of groundwater in the subsurface is of fundamental importance in the exploitation of water, just like knowledge of all associated hydrogeological context. These factors are primarily controlled by geometry of a certain pore system, given the nature of sedimentary aquifers. Thus, the microstructural characterization, as the interconnectivity of the system, it is essential to know the macro properties porosity and permeability of reservoir rock, in which can be done on a statistical characterization by twodimensional analysis. The latter is being held on a computing platform, using image thin sections of reservoir rock, allowing the prediction of the properties effective porosity and hydraulic conductivity. For Barreiras Aquifer to obtain such parameters derived primarily from the interpretation of tests of aquifers, a practice that usually involves a fairly complex logistics in terms of equipment and personnel required in addition to high cost of operation. Thus, the analysis and digital image processing is presented as an alternative tool for the characterization of hydraulic parameters, showing up as a practical and inexpensive method. This methodology is based on a flowchart work involving sampling, preparation of thin sections and their respective images, segmentation and geometric characterization, three-dimensional reconstruction and flow simulation. In this research, computational image analysis of thin sections of rocks has shown that aquifer storage coefficients ranging from 0,035 to 0,12 with an average of 0,076, while its hydrogeological substrate (associated with the top of the carbonate sequence outcropping not region) presents effective porosities of the order of 2%. For the transport regime, it is evidenced that the methodology presents results below of those found in the bibliographic data relating to hydraulic conductivity, mean values of 1,04 x10-6 m/s, with fluctuations between 2,94 x10-6 m/s and 3,61x10-8 m/s, probably due to the larger scale study and the heterogeneity of the medium studied.

Table 1. pH, Temperature, conductivity and hydrochemistry of samples

Carbon dioxide deep geological storage, especially in deep saline aquifers, is one of the preferred technological options to mitigate the effects of greenhouse gases emissions. Thus, in the last decade, studies characterising the behaviour of potential CO2 deep geological storage sites along with thorough safety assessments have been considered essential in order to minimise the risks associated with these sites. The study of natural analogues represents the best source of reliable information about the expected hydrogeochemical processes involved in the CO2 storage in such deep saline aquifers. In this work, a comprehensive study of the hydrogeochemical features and processes taking place at the natural analogue of the Alicún de las Torres thermal system (Betic Cordillera) has been conducted. Thus, the main water/CO2/rock interaction processes occurring at the thermal system have been identified, quantified and modelled, and a principle conclusion is that the hydrogeochemical evolution of the thermal system is controlled by a global dedolomitization process triggered by gypsum dissolution. This geochemical process generates a different geochemical environment to that which would result from the exclusive dissolution of carbonates from the deep aquifer, which is generally considered as the direct result of CO2 injection in a deep carbonate aquifer. Therefore, discounting of the dedolomitization process in any CO2 deep geological storage may lead to erroneous conclusions. This process will also influence the porosity evolution of the CO2 storage formation, which is a very relevant parameter when evaluating a reservoir for CO2 storage. The geothermometric calculation performed in this work leads to estimate that the thermal water reservoir is located between 650 and 800 m depth, which is very close to the minimum required to inject CO2 in a deep geological storage. It is clear that the proper characterisation of the features and hydrogeochemical processes taking place at a natural system analogous to a man-made deep geological storage will provide useful conceptual, semi-quantitative and even quantitative information about the processes and consequences that may occur at the artificial storage system.

Hydrochemical and microbiological properties of groundwater from pristine aquifer sampled at the Lower Rhine Embayment

The hydrochemistry and the microbial diversity of a pristine aquifer system near Garzweiler, Germany next to the open-pit lignite mine Garzweiler 1, were characterized. Hydrogeochemical and isotopic data indicate a recent activity of sulfate-reducing bacteria in the Tertiary marine sands. The community structure in the aquifer was studied by fluorescence in situ hybridization (FISH). Up to 7.3 x 10**5 cells/ml were detected by DAPIstaining. Bacteria (identified by the probe EUB338) were dominant, representing 51.9% of the total cell number (DAPI). Another 25.7% of total cell were affiliated with the domain Archaea as identified by the probe ARCH915. Within the domain Bacteria, the beta-Proteobacteria were most abundant (21.0% of total cell counts). Using genusspecific probes for sulfate-reducing bacteria (SRB), 2.5% of the total cells were identified as members of the genus Desulfotomaculum. This reflects the predominant role these microorganisms have been found to play in sulfatereducing zones of aquifers at other sites. Previously, all SRB cultured from this site were from the spore-forming genera Desulfotomaculum and Desulfosporosinus. Samples were taken after pumping for >= 40 min and after parameters such as temperature, pH, redox potential, oxygen and conductivity of the groundwater had remained stable for >= 15 min due to recharge of aquifer water. Hybridization and microscopy counts of hybridized and 4',6'-diamidino-2-phenylindole (DAPI)- stained cells were performed as described in Snaidr et al., (1997, http://aem.asm.org/content/63/7/2884.full.pdf). Means were calculated from 10 to 20 randomly chosen fields on each filter section, corresponding to 800-1000 DAPI stained cells. Counting results were always corrected by subtracting signals observed with the probe NON338. Formamide concentrations and oligonucleotide probes used please see further details.

Hydrochemical properties of groundwater samples in Semarang area, Java Island, Indonesia (1992, 1993, 2003, 2006, and 2007)

In Semarang City, groundwater has been exploited as a natural resource since 1841. The groundwater exploited in deep wells is concentrated in confined aquifers. The previous hydrogeological model was developed in one unit of aquifer and refined then by using several hydrostratigraphical units following a regional hydrogeological map without any further analysis. At present, there is a lack of precise hydrogeological model which integrates geological and hydrogeological data, in particular for multiple aquifers in Semarang. Thus, the aim of this paper is to develop a hydrogeological model for the multiple aquifers in Semarang using an integrated data approach. Groundwater samples in the confined aquifers have been analyzed to define the water type and its lateral distribution. Two hydrogeological cross sections were then created based on several borelog data to define a hydrostratigraphical unit (HSU). The HSU result indicates the hydrogeological model of Semarang consists of two aquifers, three aquitards, and one aquiclude. Aquifer 1 is unconfined, while Aquifer 2 is confined. Aquifer 2 is classified into three groups (2a, 2b, and 2c) based on analyses of major ion content and hydrostratigraphical cross sections.