Pore water concentrations and solid phase concentrations of sediment cores from the Ligurian continental slope

In October 1979, a period of heavy rainfall along the French Riviera was followed by the collapse of the Ligurian continental slope adjacent to the airport of Nice, France. A body of slope sediments, which was shortly beforehand affected by construction work south of the airport, was mobilized and traveled hundreds of kilometers downslope into the Var submarine canyon and, eventually, into the deep Ligurian basin. As a direct consequence, the construction was destroyed, seafloor cables were torn, and a small tsunami hit Antibes shortly after the failure. Hypotheses regarding the trigger mechanism include (i) vertical loading by construction of an embankment south of the airport, (ii) failure of a layer of sensitive clay within the slope sequence, and (iii) excess pore fluid pressures from charged aquifers in the underground. Over the previous decades, both the sensitive clay layers and the permeable sand and gravel layers were sampled to detect freshened waters. In 2007, the landslide scar and adjacent slopes were revisited for high-resolution seafloor mapping and systematic sampling. Results from half a dozen gravity and push cores in the shallow slope area reveal a limited zone of freshening (i.e. groundwater influence). A 100-250 m wide zone of the margin shows pore water salinities of 5-50% SW concentration and depletion in Cl, SO4, but Cr enrichment, while cores east or west of the landslide scar show regular SW profiles. Most interestingly, the three cores inside the landslide scar hint towards a complex hydrological system with at least two sources for groundwater. The aquifer system also showed strong freshening after a period of several months without significant precipitation. This freshening implies that charged coarse-grained layers represent a permanent threat to the slope's stability, not just after periods of major rainfall such as in October 1979.

137Cs and 90Sr depositions, and 137Cs/90Sr ratios of samples from the Kraton-3 underground nuclear explosion site, Russia

Strontium-90 activity concentrations in surface soils and areal deposition densities have been studied at a site contaminated by an accidental release to atmosphere from the underground nuclear explosion 'Kraton-3' conducted near the Polar Circle (65.9°N, 112.3°E) within the territory of the former USSR in 1978. In 2001-2002, the ground surface contamination at 14 plots studied ranged from 20 to 15000 kBq/m**2, which significantly exceeds the value of 0.44 kBq/m**2 deduced for three background plots. The zone with substantial radiostrontium contamination extends, at least, 2.5 km in a north-easterly direction from the borehole. The average 137Cs/90Sr ratio in the ground contamination originated from the 'Kraton-3' fallout was estimated to be 0.55, which is significantly different from the ratio of 2.05 evaluated for background plots contaminated mostly from global fallout. Although vertical migration of 90Sr in all undisturbed soil profiles studied is more rapid than that for 137Cs, the depth of percolation of both radionuclides into the ground is mostly limited to the top 10-20 cm, which may be explained, primarily, by permafrost conditions. The horizontal migration rate of radiostrontium in the aqueous phase exceeds the radiocaesium migration rate by many times. This phenomenon seems to be a reason for the significant enrichment of the soil surface layers by radiostrontium at some sites, with variations occurring in accordance with small-scale irregularities of landscape.

(Table 1) Isotopic and geochemical parameters of gas phase from underground fluids of the Baikal rift and its framing

Results of studying isotopic composition of helium in underground fluids of the Baikal-Mongolian region during the last quarter of XX century are summarized. Determinations of 3He/4He ratio in 139 samples of gas phase from fluids, collected at 104 points of the Baikal rift zone and adjacent structures are given. 3He/4He values lie within the range from 1x10**-8 (typical for crustal radiogenic helium) to 1.1x10**-5 (close to typical MORB reservoir). Repeated sampling in some points during more than 20 years showed stability of helium isotopic composition in time in each of them at any level of 3He/4He values. There is no systematic differences of 3He/4He in samples from surface water sources and deeper intervals of boreholes in the same areas. Universal relationship between isotopic composition of helium and general composition of gas phase is absent either, but the minimum 3He/4He values occurred in methane gas of hydrocarbon deposits, whereas in nitrogen and carbon dioxide gases of helium composition varied (in the latter maximum 3He/4He values have been measured). According to N2/Ar_atm ratio nitrogen gases are atmospheric. In carbonic gas fN2/fNe ratio indicates presence of excessive (non-atmogenic) nitrogen, but the attitude CO2/3He differs from one in MORB. Comparison of helium isotopic composition with its concentration and composition of the main components of gas phase from fluids shows that it is formed under influence of fractionation of components with different solubility in the gas-water system and generation/consumption of reactive gases in the crust. Structural and tectonic elements of the region differ from the spectrum of 3He/4He values. At the pre-Riphean Siberian Platform the mean 3He/4He = (3.6+/-0.9)x10**- 8 is very close to radiogenic one. In the Paleozoic crust of Khangay 3He/4He = (16.3+/-4.6)x10**-8, and the most probable estimate is (12.3+/-2.9)x10**-8. In structures of the eastern flank of the Baikal rift zone (Khentei, Dauria) affected by the Mz-Kz activization 3He/4He values range from 4.4x10**-8 to 2.14x10**-6 (average 0.94x10**-6). Distribution of 3He/4He values across the strike of the Baikal rift zone indicates advective heat transfer from the mantle not only in the rift zone, but also much further to the east. In fluids of the Baikal rift zone range of 3He/4He values is the widest: from 4x10**-8 to 1.1x10**-5. Their variations along the strike of the rift zone are clearly patterned, namely, decrease of 3He/4He values in both directions from the Tunka depression. Accompanied by decrease in density of conductive heat flow and in size of rift basins, this trend indicates decrease in intensity of advective heat transfer from the mantle to peripheral segments of the rift zone. Comparing this trend with data on other continental rift zones and mid-ocean ridges leads to the conclusion about fundamental differences in mechanisms of interaction between the crust and the mantle in these environments.

Full scale test for explosion water barriers in small cross-section galleries

Underground coal mines explosions generally arise from the inflammation of a methane/air mixture. This explosion can also generate a subsequent coal dust explosion. Traditionally such explosions have being fought eliminating one or several of the factors needed by the explosion to take place. Although several preventive measures are taken to prevent explosions, other measures should be considered to reduce the effects or even to extinguish the flame front. Unlike other protection methods that remove one or two of the explosion triangle elements, namely; the ignition source, the oxidizing agent and the fuel, explosion barriers removes all of them: reduces the quantity of coal in suspension, cools the flame front and the steam generated by vaporization removes the oxygen present in the flame. Passive water barriers are autonomous protection systems against explosions that reduce to a satisfactory safety level the effects of methane and/or flammable dust explosions. The barriers are activated by the pressure wave provoked in the explosion destroying the barrier troughs and producing a uniform dispersion of the extinguishing agent throughout the gallery section in quantity enough to extinguish the explosion flame. Full scale tests have been carried out in Polish Barbara experimental mine at GIG Central Mining Institute in order to determine the requirements and the optimal installation conditions of these devices for small sections galleries which are very frequent in the Spanish coal mines. Full scale tests results have been analyzed to understand the explosion timing and development, in order to assess on the use of water barriers in the typical small crosssection Spanish galleries. Several arrangements of water barriers have been designed and tested to verify the effectiveness of the explosion suppression in each case. The results obtained demonstrate the efficiency of the water barriers in stopping the flame front even with smaller amounts of water than those established by the European standard. According to the tests realized, water barriers activation times are between 0.52 s and 0.78 s and the flame propagation speed are between 75 m/s and 80 m/s. The maximum pressures (Pmax) obtained in the full scale tests have varied between 0.2 bar and 1.8 bar. Passive barriers protect effectively against the spread of the flame but cannot be used as a safeguard of the gallery between the ignition source and the first row of water troughs or bags, or even after them, as the pressure could remain high after them even if the flame front has been extinguished.

Implementation of the Heated Pulsed Theory using actively heated fiber optics: measurements of soil volumetric water content and soil volumetric heat capacity

Existe una creciente necesidad de hacer el mejor uso del agua para regadío. Una alternativa eficiente consiste en la monitorización del contenido volumétrico de agua (θ), utilizando sensores de humedad. A pesar de existir una gran diversidad de sensores y tecnologías disponibles, actualmente ninguna de ellas permite obtener medidas distribuidas en perfiles verticales de un metro y en escalas laterales de 0.1-1,000 m. En este sentido, es necesario buscar tecnologías alternativas que sirvan de puente entre las medidas puntuales y las escalas intermedias. Esta tesis doctoral se basa en el uso de Fibra Óptica (FO) con sistema de medida de temperatura distribuida (DTS), una tecnología alternativa de reciente creación que ha levantado gran expectación en las últimas dos décadas. Específicamente utilizamos el método de fibra calentada, en inglés Actively Heated Fiber Optic (AHFO), en la cual los cables de Fibra Óptica se utilizan como sondas de calor mediante la aplicación de corriente eléctrica a través de la camisa de acero inoxidable, o de un conductor eléctrico simétricamente posicionado, envuelto, alrededor del haz de fibra óptica. El uso de fibra calentada se basa en la utilización de la teoría de los pulsos de calor, en inglés Heated Pulsed Theory (HPP), por la cual el conductor se aproxima a una fuente de calor lineal e infinitesimal que introduce calor en el suelo. Mediante el análisis del tiempo de ocurrencia y magnitud de la respuesta térmica ante un pulso de calor, es posible estimar algunas propiedades específicas del suelo, tales como el contenido de humedad, calor específico (C) y conductividad térmica. Estos parámetros pueden ser estimados utilizando un sensor de temperatura adyacente a la sonda de calor [método simple, en inglés single heated pulsed probes (SHPP)], ó a una distancia radial r [método doble, en inglés dual heated pulsed probes (DHPP)]. Esta tesis doctoral pretende probar la idoneidad de los sistemas de fibra óptica calentada para la aplicación de la teoría clásica de sondas calentadas. Para ello, se desarrollarán dos sistemas FO-DTS. El primero se sitúa en un campo agrícola de La Nava de Arévalo (Ávila, España), en el cual se aplica la teoría SHPP para estimar θ. El segundo sistema se desarrolla en laboratorio y emplea la teoría DHPP para medir tanto θ como C. La teoría SHPP puede ser implementada con fibra óptica calentada para obtener medidas distribuidas de θ, mediante la utilización de sistemas FO-DTS y el uso de curvas de calibración específicas para cada suelo. Sin embargo, la mayoría de aplicaciones AHFO se han desarrollado exclusivamente en laboratorio utilizando medios porosos homogéneos. En esta tesis se utiliza el programa Hydrus 2D/3D para definir tales curvas de calibración. El modelo propuesto es validado en un segmento de cable enterrado en una instalación de fibra óptica y es capaz de predecir la respuesta térmica del suelo en puntos concretos de la instalación una vez que las propiedades físicas y térmicas de éste son definidas. La exactitud de la metodología para predecir θ frente a medidas puntuales tomadas con sensores de humedad comerciales fue de 0.001 a 0.022 m3 m-3 La implementación de la teoría DHPP con AHFO para medir C y θ suponen una oportunidad sin precedentes para aplicaciones medioambientales. En esta tesis se emplean diferentes combinaciones de cables y fuentes emisoras de calor, que se colocan en paralelo y utilizan un rango variado de espaciamientos, todo ello en el laboratorio. La amplitud de la señal y el tiempo de llegada se han observado como funciones del calor específico del suelo. Medidas de C, utilizando esta metodología y ante un rango variado de contenidos de humedad, sugirieron la idoneidad del método, aunque también se observaron importantes errores en contenidos bajos de humedad de hasta un 22%. La mejora del método requerirá otros modelos más precisos que tengan en cuenta el diámetro del cable, así como la posible influencia térmica del mismo. ABSTRACT There is an increasing need to make the most efficient use of water for irrigation. A good approach to make irrigation as efficient as possible is to monitor soil water content (θ) using soil moisture sensors. Although, there is a broad range of different sensors and technologies, currently, none of them can practically and accurately provide vertical and lateral moisture profiles spanning 0-1 m depth and 0.1-1,000 m lateral scales. In this regard, further research to fulfill the intermediate scale and to bridge single-point measurement with the broaden scales is still needed. This dissertation is based on the use of Fiber Optics with Distributed Temperature Sensing (FO-DTS), a novel approach which has been receiving growing interest in the last two decades. Specifically, we employ the so called Actively Heated Fiber Optic (AHFO) method, in which FO cables are employed as heat probe conductors by applying electricity to the stainless steel armoring jacket or an added conductor symmetrically positioned (wrapped) about the FO cable. AHFO is based on the classic Heated Pulsed Theory (HPP) which usually employs a heat probe conductor that approximates to an infinite line heat source which injects heat into the soil. Observation of the timing and magnitude of the thermal response to the energy input provide enough information to derive certain specific soil thermal characteristics such as the soil heat capacity, soil thermal conductivity or soil water content. These parameters can be estimated by capturing the soil thermal response (using a thermal sensor) adjacent to the heat source (the heating and the thermal sources are mounted together in the so called single heated pulsed probe (SHPP)), or separated at a certain distance, r (dual heated pulsed method (DHPP) This dissertation aims to test the feasibility of heated fiber optics to implement the HPP theory. Specifically, we focus on measuring soil water content (θ) and soil heat capacity (C) by employing two types of FO-DTS systems. The first one is located in an agricultural field in La Nava de Arévalo (Ávila, Spain) and employ the SHPP theory to estimate θ. The second one is developed in the laboratory using the procedures described in the DHPP theory, and focuses on estimating both C and θ. The SHPP theory can be implemented with actively heated fiber optics (AHFO) to obtain distributed measurements of soil water content (θ) by using reported soil thermal responses in Distributed Temperature Sensing (DTS) and with a soil-specific calibration relationship. However, most reported AHFO applications have been calibrated under laboratory homogeneous soil conditions, while inexpensive efficient calibration procedures useful in heterogeneous soils are lacking. In this PhD thesis, we employ the Hydrus 2D/3D code to define these soil-specific calibration curves. The model is then validated at a selected FO transect of the DTS installation. The model was able to predict the soil thermal response at specific locations of the fiber optic cable once the surrounding soil hydraulic and thermal properties were known. Results using electromagnetic moisture sensors at the same specific locations demonstrate the feasibility of the model to detect θ within an accuracy of 0.001 to 0.022 m3 m-3. Implementation of the Dual Heated Pulsed Probe (DPHP) theory for measurement of volumetric heat capacity (C) and water content (θ) with Distributed Temperature Sensing (DTS) heated fiber optic (FO) systems presents an unprecedented opportunity for environmental monitoring. We test the method using different combinations of FO cables and heat sources at a range of spacings in a laboratory setting. The amplitude and phase-shift in the heat signal with distance was found to be a function of the soil volumetric heat capacity (referred, here, to as Cs). Estimations of Cs at a range of θ suggest feasibility via responsiveness to the changes in θ (we observed a linear relationship in all FO combinations), though observed bias with decreasing soil water contents (up to 22%) was also reported. Optimization will require further models to account for the finite radius and thermal influence of the FO cables, employed here as “needle probes”. Also, consideration of the range of soil conditions and cable spacing and jacket configurations, suggested here to be valuable subjects of further study and development.

When Wells Run Dry: Water and Tourism Along the Western Coast of Nicaragua

This dissertation uses a political ecology approach to examine the relationship between tourism development and groundwater in southwest Nicaragua. Tourism in Nicaragua is a booming industry bolstered by ‘unspoiled’ natural beauty, low crime rates, and government incentives. This growth has led to increased infrastructure, revenue, and employment opportunities for many local communities along the Pacific coast. Not surprisingly, it has also brought concomitant issues of deeper poverty, widening gaps between rich and poor, and competition over natural resources. Adequate provisions of freshwater are necessary to sustain the production and reproduction of tourism; however, it remains uncertain if groundwater supplies can keep pace with demand. The objective of this research is to assess water supply availability amidst tourism development in the Playa Gigante area. It addresses the questions: 1) are local groundwater supplies sufficient to sustain the demand for freshwater imposed by increased tourism development? and 2) is there a power relationship between tourism development and control over local freshwater that would prove inequitable to local populations? Integrating the findings of groundwater monitoring, geological mapping, and ethnographic and survey research from a representative stretch of Pacific coastline, this dissertation shows that diminishing recharge and increased groundwater consumption is creating conflict between stakeholders with various levels of knowledge, power, and access. Although national laws are structured to protect the environment and ensure equitable access to groundwater, the current scramble to secure water has powerful implications on social relations and power structures associated with tourism development. This dissertation concludes that marginalization due to environmental degradation is attributable to the nexus of a political promotion of tourism, poorly enforced state water policies, insufficient water research, and climate change. Greater technical attention to hydrological dynamics and collaboration amongst stakeholders are necessary for equitable access to groundwater, environmental sustainability, and profitability of tourism.

Effect of water vapour condensation on the radon content in subsurface air in a hypogeal inactive-volcanic environment in Galdar cave, Spain

Fluctuations of trace gas activity as a response to variations in weather and microclimate conditions were monitored over a year in a shallow volcanic cave (Painted Cave, Galdar, Canary Islands, Spain). 222Rn concentration was used due to its greater sensitivity to hygrothermal variations than CO2 concentration. Radon concentration in the cave increases as effective vapour condensation within the porous system of the rock surfaces inside the cave increases due to humidity levels of more than 70%. Condensed water content in pores was assessed and linked to a reduction in the direct passage of trace gases. Fluctuations in radon activity as a response to variations in weather and microclimate conditions were statistically identified by clustering entropy changes on the radon signal and parameterised to predict radon concentration anomalies. This raises important implications for other research fields, including the surveillance of shallow volcanic and seismic activity, preventive conservation of cultural heritage in indoor spaces, indoor air quality control and studies to improve understanding of the role of subterranean terrestrial ecosystems as reservoirs and/or temporary sources of trace gases.

The Modelling of Photovoltaic, Solar Thermal, and Photovoltaic/Thermal Domestic Hot Water Systems

Solar heating of potable water has traditionally been accomplished through the use of solar thermal (ST) collectors. With the recent increases in availability and lower cost of photovoltaic (PV) panels, the potential of coupling PV solar arrays to electrically heated domestic hot water (DHW) tanks has been considered. Additionally, innovations in the SDHW industry have led to the creation of photovoltaic/thermal (PV/T) collectors, which heat water using both electrical and thermal energy. The current work compared the performance and cost-effectiveness of a traditional solar thermal (ST) DHW system to PV-solar-electric DHW systems and a PV/T DHW system. To accomplish this, a detailed TRNSYS model of the solar hot water systems was created and annual simulations were performed for 250 L/day and 325 L/day loads in Toronto, Vancouver, Montreal, Halifax, and Calgary. It was shown that when considering thermal performance, PV-DHW systems were not competitive when compared to ST-DHW and PVT-DHW systems. As an example, for Toronto the simulated annual solar fractions of PV-DHW systems were approximately 30%, while the ST-DHW and PVT-DHW systems achieved 65% and 71% respectively. With current manufacturing and system costs, the PV-DHW system was the most cost-effective system for domestic purposes. The capital cost of the PV-DHW systems were approximately $1,923-$2,178 depending on the system configuration, and the ST-DHW and PVT system were estimated to have a capital cost of $2,288 and $2,373 respectively. Although the capital cost of the PVT-DHW system was higher than the other systems, a Present Worth analysis for a 20-year period showed that for a 250 L/day load in Toronto the Present Worth of the PV/T system was approximately $4,597, with PV-DHW systems costing approximately $7,683-$7,816 and the ST-DHW system costing $5,238.

Dodecyltrimethylammonium chloride adsorption at the silica/water interface studied by Sum Frequency Generation

La génération des fréquences somme (SFG), une technique spectroscopique spécifique aux interfaces, a été utilisée pour caractériser les changements de la structure macromoléculaire du surfactant cationique chlorure de dodécyltriméthylammonium (DTAC) à l’interface silice/eau dans une plage de pH variant entre 3 et 11. Les conditions expérimentales ont été choisies pour imiter les conditions les plus communes trouvées pendant les opérations de récupération assistée du pétrole. Particulièrement, la silice a été étudiée, car elle est un des composantes des surfaces minérales des réservoirs de grès, et l’adsorption du surfactant a été étudiée avec une force ionique pertinente pour les fluides de la fracturation hydraulique. Les spectres SFG ont présenté des pics détectables avec une amplitude croissante dans la région des étirements des groupes méthylène et méthyle lorsque le pH est diminué jusqu’à 3 ou augmenté jusqu’à 11, ce qui suggère des changements de la structure des agrégats de surfactant à l’interface silice/eau à une concentration de DTAC au-delà de la concentration micellaire critique. De plus, des changements dans l’intensité SFG ont été observés pour le spectre de l’eau quand la concentration de DTAC augmente de 0,2 à 50 mM dans les conditions acide, neutre et alcaline. À pH 3, près du point de charge zéro de la surface de silice, l’excès de charge positive en raison de l’adsorption du surfactant cationique crée un champ électrostatique qui oriente les molécules d’eau à l’interface. À pH 7 et 11, ce qui sont des valeurs au-dessus du point de charge zéro de la surface de silice, le champ électrostatique négatif à l’interface silice/eau diminue par un ordre de grandeur avec l’adsorption du surfactant comme résultat de la compensation de la charge négative à la surface par la charge positive du DTAC. Les résultats SFG ont été corrélés avec des mesures de l’angle de contact et de la tension interfaciale à pH 3, 7 et 11.

Alteration of natural ³⁷Ar activity concentration in the subsurface by gas transport and water infiltration

High ³⁷Ar activity concentration in soil gas is proposed as a key evidence for the detection of underground nuclear explosion by the Comprehensive Nuclear Test-Ban Treaty. However, such a detection is challenged by the natural background of ³⁷Ar in the subsurface, mainly due to Ca activation by cosmic rays. A better understanding and improved capability to predict ³⁷Ar activity concentration in the subsurface and its spatial and temporal variability is thus required. A numerical model integrating ³⁷Ar production and transport in the subsurface is developed, including variable soil water content and water infiltration at the surface. A parameterized equation for ³⁷Ar production in the first 15 m below the surface is studied, taking into account the major production reactions and the moderation effect of soil water content. Using sensitivity analysis and uncertainty quantification, a realistic and comprehensive probability distribution of natural ³⁷Ar activity concentrations in soil gas is proposed, including the effects of water infiltration. Site location and soil composition are identified as the parameters allowing for a most effective reduction of the possible range of ³⁷Ar activity concentrations. The influence of soil water content on ³⁷Ar production is shown to be negligible to first order, while ³⁷Ar activity concentration in soil gas and its temporal variability appear to be strongly influenced by transient water infiltration events. These results will be used as a basis for practical CTBTO concepts of operation during an OSI.

Dodecyltrimethylammonium chloride adsorption at the silica/water interface studied by Sum Frequency Generation

La génération des fréquences somme (SFG), une technique spectroscopique spécifique aux interfaces, a été utilisée pour caractériser les changements de la structure macromoléculaire du surfactant cationique chlorure de dodécyltriméthylammonium (DTAC) à l’interface silice/eau dans une plage de pH variant entre 3 et 11. Les conditions expérimentales ont été choisies pour imiter les conditions les plus communes trouvées pendant les opérations de récupération assistée du pétrole. Particulièrement, la silice a été étudiée, car elle est un des composantes des surfaces minérales des réservoirs de grès, et l’adsorption du surfactant a été étudiée avec une force ionique pertinente pour les fluides de la fracturation hydraulique. Les spectres SFG ont présenté des pics détectables avec une amplitude croissante dans la région des étirements des groupes méthylène et méthyle lorsque le pH est diminué jusqu’à 3 ou augmenté jusqu’à 11, ce qui suggère des changements de la structure des agrégats de surfactant à l’interface silice/eau à une concentration de DTAC au-delà de la concentration micellaire critique. De plus, des changements dans l’intensité SFG ont été observés pour le spectre de l’eau quand la concentration de DTAC augmente de 0,2 à 50 mM dans les conditions acide, neutre et alcaline. À pH 3, près du point de charge zéro de la surface de silice, l’excès de charge positive en raison de l’adsorption du surfactant cationique crée un champ électrostatique qui oriente les molécules d’eau à l’interface. À pH 7 et 11, ce qui sont des valeurs au-dessus du point de charge zéro de la surface de silice, le champ électrostatique négatif à l’interface silice/eau diminue par un ordre de grandeur avec l’adsorption du surfactant comme résultat de la compensation de la charge négative à la surface par la charge positive du DTAC. Les résultats SFG ont été corrélés avec des mesures de l’angle de contact et de la tension interfaciale à pH 3, 7 et 11.

More about leaking underground storage tanks : a background booklet for the chemical advisory /

Mode of access: Internet.

Evaluation of underground injection of industrial waste in Illinois /

Published cooperatively by the Department of Energy and Natural Resources, Illinois State Geological Survey, Illinois State Water Survey, and Hazardous Waste Research and Information Center.

Evaluation of current underground injection of industrial waste in Illinois : final draft report /

"HWRIC RR 008." -- Cover.

Replenishment of ground water supplies by artificial means /

"Based on Technical bulletin no. 578, Spreading water for storage underground."