On-site wastewater treatment using subsurface flow constructed wetlands in Ireland.

PubMed

Gill, Laurence W; O'Luanaigh, Niall; Johnston, Paul M

2011-01-01

The results from an Irish EPA-funded project on the effectiveness of using constructed wetlands for treating wastewater from single households is presented, which has contributed to the design guidelines included in the new EPA Code of Practice. Three subsurface flow gravel-filled wetlands were constructed on separate sites--one to provide secondary treatment and the other two to provide tertiary treatment stages for the domestic effluent. A comprehensive analysis over three years was then conducted to provide a robust characterization of the internal dynamics of the systems, particularly with respect to N and P removal as well as evaluating the temporal water balance across the different seasons. The removal of Total N was only 29% and 30% in the secondary and tertiary treatment wetlands, respectively; particularly disappointing for the tertiary treatment process, which was receiving nitrified effluent. Studies on the (15)N stable isotope confirmed that 35% of the ammonium from the septic tank was passing straight through the process without taking part in any biogeochemical processes. However, influent N in the wetlands was shown to be biologically assimilated into organic nitrogen and then released again as soluble ammonium--so-called nitrogen "spiraling." Removal of Total P in the wetlands averaged from 28% to 45% with higher P removals measured during summer periods, although the effluent concentrations were still found to be high (> 5 mg/l on average). The phosphorus in the plant material was also analysed revealing that the annual above-ground stem matter only accounted for 1.3% to 8.4% of the annual total P-load in the wetlands. Finally, the water balance analyses showed that the mean flow discharging from both the secondary and tertiary treatment wetlands was slightly greater than the mean flow to the reed bed over the trial period, with rainfall acting to increase flows by 13% and 5%, respectively, on average in winter while just about balancing

Treatment of landfill leachate using an aerated, horizontal subsurface-flow constructed wetland.

PubMed

Nivala, J; Hoos, M B; Cross, C; Wallace, S; Parkin, G

2007-07-15

A pilot-scale subsurface-flow constructed wetland was installed at the Jones County Municipal Landfill, near Anamosa, Iowa, in August 1999 to demonstrate the use of constructed wetlands as a viable low-cost treatment option for leachate generated at small landfills. The system was equipped with a patented wetland aeration process to aid in removal of organic matter and ammonia nitrogen. The high iron content of the leachate caused the aeration system to cease 2 years into operation. Upon the installation of a pretreatment chamber for iron removal and a new aeration system, treatment efficiencies dramatically improved. Seasonal performance with and without aeration is reported for 5-day biochemical oxygen demand (BOD(5)), chemical oxygen demand (COD), ammonia nitrogen (NH(4)-N), and nitrate nitrogen (NO(3)-N). Since winter air temperatures in Iowa can be very cold, a layer of mulch insulation was installed on top of the wetland bed to keep the system from freezing. When the insulation layer was properly maintained (either through sufficient litterfall or replenishing the mulch layer), the wetland sustained air temperatures of as low as -26 degrees C without freezing problems.

Design, construction and performance of a horizontal subsurface flow wetland system in Australia.

PubMed

Bolton, Lise M W; Bolton, Keith G E

2013-01-01

Malabugilmah is a remote Aboriginal community located in Clarence Valley, Northern NSW, Australia. In 2006, seven horizontal subsurface flow wetland clusters consisting of 3 m × 2 m wetland cells in series were designed and constructed to treat septic tank effluent to a secondary level (Total Suspended Solids (TSS) < 30 mg/L and Biochemical Oxygen Demand (BOD5) <20 mg/L) and achieve >50% Total Nitrogen (TN) reduction, no net Total Phosphorus (TP) export and ≥99.9% Faecal Coliform (FC) reduction. The wetland cell configuration allowed the wetlands to be located on steeper terrain, enabling effluent to be treated to a secondary level without the use of pumps. In addition to the water quality targets, the wetlands were designed and constructed to satisfy environmental, economic and social needs of the community. The wetland systems were planted with a local Australian wetland tree species which has become well established. Two wetland clusters have been monitored over the last 4 years. The wetlands have demonstrated to be robust over time, providing a high level of secondary treatment over an extended period.

Design and monitoring of horizontal subsurface-flow constructed wetlands for treating nursery leachates.

PubMed

Narváez, Lola; Cunill, Conrad; Cáceres, Rafaela; Marfà, Oriol

2011-06-01

Nursery leachates usually contain high concentrations of nitrates, phosphorus and potassium, so discharging them into the environment often causes pollution. Single-stage or two-stage horizontal subsurface flow constructed wetlands (HSSCW) filled with different substrates were designed to evaluate the effect and evolution over time of the removal of nitrogen and other nutrients contained in nursery leachates. The addition of sodium acetate to achieve a C:NO(3)(-)-N ratio of 3:1 was sufficient to reach complete denitrification in all HSSCW. The removal rate of nitrate was high throughout the operation period (over 98%). Nevertheless, the removal rate of ammonium decreased about halfway through the operation. Removal of the COD was enhanced by the use of two-stage HSSCW. In general, the substrates and the number of stages of the wetlands did not affect the removal of nitrogen, total phosphorus and potassium. Copyright © 2011 Elsevier Ltd. All rights reserved.

Design configurations affecting flow pattern and solids accumulation in horizontal free water and subsurface flow constructed wetlands.

PubMed

Pedescoll, A; Sidrach-Cardona, R; Sánchez, J C; Carretero, J; Garfi, M; Bécares, E

2013-03-01

The aim of this study was to evaluate the effect of different horizontal constructed wetland (CW) design parameters on solids distribution, loss of hydraulic conductivity over time and hydraulic behaviour, in order to assess clogging processes in wetlands. For this purpose, an experimental plant with eight CWs was built at mesocosm scale. Each CW presented a different design characteristic, and the most common CW configurations were all represented: free water surface flow (FWS) with different effluent pipe locations, FWS with floating macrophytes and subsurface flow (SSF), and the presence of plants and specific species (Typha angustifolia and Phragmites australis) was also considered. The loss of the hydraulic conductivity of gravel was greatly influenced by the presence of plants and organic load (representing a loss of 20% and c.a. 10% in planted wetlands and an overloaded system, respectively). Cattail seems to have a greater effect on the development of clogging since its below-ground biomass weighed twice as much as that of common reed. Hydraulic behaviour was greatly influenced by the presence of a gravel matrix and the outlet pipe position. In strict SSF CW, the water was forced to cross the gravel and tended to flow diagonally from the top inlet to the bottom outlet (where the inlet and outlet pipes were located). However, when FWS was considered, water preferentially flowed above the gravel, thus losing half the effective volume of the system. Only the presence of plants seemed to help the water flow partially within the gravel matrix. Copyright © 2012 Elsevier Ltd. All rights reserved.

Limestone and Zeolite as Alternative Media in Horizontal Subsurface Flow Constructed Wetlands: Laboratory-Scale Studies

NASA Astrophysics Data System (ADS)

Lizama, K.; Jaque, I.; Ayala, J.

2016-12-01

Arsenic is well known for its chronic toxicity. Millions of people around the world are currently at risk, drinking water with As concentrations above 10 ppb, the WHO drinking water guideline. Although different treatment options exist, they are often limited by elevated costs and maintenance requirements. Constructed wetlands are a natural water treatment system, capable to remove metals and metalloids -including As- via different physical, chemical and biological processes. The use of alternative supporting media to enhance As removal in subsurface flow wetlands has been recommended, but not sufficiently studied. Limestone and zeolite have been identified as effective supporting media in subsurface flow wetlands aiming As removal. However, there are still key aspects to be addressed, such as the implications of using these media, the speciation in the solid phase, the role of vegetation, etc. This study investigated the performance of limestone and zeolite in three types of experiments: batch, column and as main supporting media in a bench scale horizontal subsurface flow wetland system. Synthetic water resembling a contaminated river in Chile (As concentration=3 mg/L, Fe concentration= 100 mg/L, pH=2) was used in all experiments. In the batch experiments, the As concentration, the mass of media and the contact time were varied. The column system consisted of three limestone columns and three zeolite columns, operated under a hydraulic loading of 20 mm/d. The wetland system consisted of twelve PVC cells: six filled with zeolite and six with limestone. Phragmites australis were planted in three cells of each media type, as control cells. From the batch experiments, maximum As sorption capacities as indicated by Langmuir model were 1.3 mg/g for limestone and 0.17 mg/g for zeolite, at 18 h contact time and 6.3 g/L medium concentration. EDS and XPS analyses revealed that As and Fe were retained in zeolite at the end of the batch experiments. Zeolite and limestone

Optimal conditions for chlorothalonil and dissolved organic carbon in horizontal subsurface flow constructed wetlands.

PubMed

Rìos-Montes, Karina A; Casas-Zapata, Juan C; Briones-Gallardo, Roberto; Peñuela, Gustavo

2017-04-03

The most efficient system of horizontal subsurface flow constructed wetlands (HSSFCW) for removing dissolved organic carbon (DOC) in the presence of chlorothalonil pesticide (CLT) present in synthetic domestic wastewater was determined using the macrophyte Phragmites australis. Two concentrations of CLT (85 and 385 μg L -1 ) and one concentration of glucose (20 mg L -1 ) were evaluated in four pilot scale horizontal surface flow constructed wetlands coupled with two sizes of silica gravel, igneous gravel, fine chalky gravel (3.18-6.35 mm), coarse gravel (12.70-25.40 mm) and two water surface heights (20 and 40 cm). For a month, wetlands were acclimated with domestic wastewater. Some groups of bacteria were also identified in the biofilm attached to the gravel. In each treatment periodic samplings were conducted in the influent and effluent. Chlorothalonil was quantified by gas chromatography (GC-ECD m), DOC by an organic carbon analyzer and bacterial groups using conventional microbiology in accordance with Standard Methods. The largest removals of DOC (85.82%-85.31%) were found when using fine gravel (3.18-6.35 mm) and the lower layer of water (20 cm). The bacterial groups quantified in the biofilm were total heterotrophic, revivable heterotrophic, Pseudomonas and total coliforms. The results of this study indicate that fine grain gravel (3.18-6.35 mm) and both water levels (20 to 40 cm) can be used in the removal of organic matter and for the treatment of agricultural effluents contaminated with organo-chloride pesticides like CLT in HSSFCW.

Reconstruction of a constructed wetland with horizontal subsurface flow after 18 years of operation.

PubMed

Hudcová, Tereza; Vymazal, Jan; Dunajský, Michal Kriška

2013-01-01

The constructed wetland (CW) for 326 PE with horizontal subsurface flow at Kotenčice, Central Bohemia, Czech Republic, was built in 1994. Despite the relatively high efficiency of the CW, the filtration beds suffered from clogging, and therefore it was decided in 2011 to rebuild the whole system. The new treatment system was built as an experimental system consisting of four different combinations of horizontal and vertical beds. The major aim of the design was to determine the best hybrid combination which then could be used in the future for refurbishment of older horizontal flow CWs or for the new systems. The mechanical pretreatment consists of mechanical bar screens, a new Imhoff tank, and the original settling tank which has been converted into the accumulation tank from where the wastewater is pumped into the wetlands. The filters are planted with Phragmites australis, Phalaris arundinacea, Iris pseudacorus, Iris sibirica, Glyceria maxima and Lythrum salicaria in order to evaluate and compare various plant species' effect on the treatment process. The new technology includes a tertiary treatment which consists of a greenhouse with a photo-reactor for the cultivation of algae and hydroponic systems (residual nutrients removal), sludge reed-beds and a composting field.

Nitrogen and COD Removal from Septic Tank Wastewater in Subsurface Flow Constructed Wetlands: Plants Effects.

PubMed

Collison, R S; Grismer, M E

2015-11-01

We evaluated subsurface flow (SSF) constructed wetland treatment performance with respect to organics (COD) and nitrogen (ammonium and nitrate) removal from domestic (septic tank) wastewater as affected by the presence of plants, substrate "rock" cation exchange capacity (CEC), laboratory versus field conditions and use of synthetic as compared to actual domestic wastewater. This article considers the effects of plants on constructed wetland treatment in the field. Each constructed wetland system was comprised of two beds (2.6 m long by 0.28 m wide and deep filled with ~18 mm crushed lava rock) separated by an aeration tank connected in series. The lava rock had a porosity of ~47% and a CEC of 4 meq/100 gm. One pair of constructed wetland systems was planted with cattails in May 2008, while an adjacent pair of systems remained un-planted. Collected septic tank or synthesized wastewater was allowed to gravity feed each constructed wetland system and effluent samples were regularly collected and tested for COD and nitrogen species during four time periods spanning November 2008 through June 2009. These effluent concentrations were tested for statistical differences at the 95% level for individual time periods as well as the overall 6-month period. Organics removal from domestic wastewater was 78.8% and 76.1% in the planted and un-planted constructed wetland systems, respectively, while ammonium removal was 94.5% and 90.2%, respectively. Similarly, organics removal from the synthetic wastewater of equivalent strength was 88.8% and 90.1% for planted and un-planted constructed wetland systems, respectively, while ammonium removal was 96.9% and 97.3%, respectively.

Removal of total suspended solids from wastewater in constructed horizontal flow subsurface wetlands.

PubMed

Manios, T; Stentiford, E I; Millner, P

2003-06-01

Subsurface horizontal flow experimental wetlands (reed beds), were designed and built based on a combination of two design methodologies, that of the WRc and Severn Trent Water plc (1996) and that of the USA, EPA (1988). Four different growing media were used with a combination of top soil, gravel, river sand, and mature sewage sludge compost, to determine the best substrate for total suspended solids (TSS) removal. Eight units were constructed, two for each growing media. One bed for each pair was planted with Typha latifolia plants commonly known as cattails. Primary treated domestic wastewater, was continuously fed to the beds for more than six months. All eight beds performed very well. The best performance was achieved by the gravel reed beds with an almost constant removal rate above 95% and an average effluent concentration of less than 10 mg/L. Soil based beds containing top soil and sand, managed to reach values of removal around 90%. The wetlands containing compost in their substrate, produced an effluent with average concentration of less than 30 mg/L and a percentage removal between 80% and 90%. As expected, there was no significant difference in the performance of planted and unplanted wetlands.

Vertical Subsurface Flow (VSSF) constructed wetland for domestic wastewater treatment

NASA Astrophysics Data System (ADS)

Perdana, M. C.; Sutanto, H. B.; Prihatmo, G.

2018-04-01

Vertical Subsurface Flow Constructed Wetland (VSSF) is appraised to become an alternative solution for treating domestic wastewater effectively and efficiently. The system which imitates the natural wetland concept is able to reduce organic material and nutrients in wastewater; therefore, it will be more feasible to be discharged to the environment. This study aimed to compare which species is more recommended to be applied for reducing organic material and nutrients in domestic wastewater. This experimental study applied four treatments, i.e 1) control (unplanted), 2) single species Iris pseudacorus, 3) single species Echinodorus palaefolius, and 4) combination (Iris pseudacorus and Echinodorus palaefolius) with three days of retention time. The application of those plants aims for holding the role in increasing wastewater quality and adding aesthetic impression at once. The plants were planted on VSSF media, in relatively same of weight and size to compare their effectiveness in decreasing organic and inorganic load. The parameters measured pervade TDS, pH, BOD5, COD, Nitrate, and Phosphate. The plants’ condition was also observed during and after the system worked. The result showed that the best average value of effectiveness for each of parameters: COD by combination treatment (50.76%), BOD5 by single I. pseudacorus (30.15%), Nitrate by single E. palaefolius (58.06%), Phosphate by single E. palaefolius (99.5%), and TDS by E.palaefolius (3.25%). The result showed that there was a significant difference of Nitrate and Phosphate reduction between control and three other treatments, while pH parameter showed non-significant change among them. In term of performance, I.pseudacorus seemed showed a preferable achievement.

Phytoremediation of selenium by two helophyte species in subsurface flow constructed wetland.

PubMed

Shardendu; Salhani, N; Boulyga, S F; Stengel, E

2003-03-01

The phytoremediation of selenium by two different wetland species was investigated. Selenium (20.4 microg/l) was supplied continuously to subsurface flow constructed wetlands, one vegetated with Typha latifolia L. and the other with Phragmites australis (Cav.) Trin. ex Steud. The beds of both species had same hydraulic loading rate (0.079 m(3)/m(2)/d) and water retention time (24 h). However, the mass loading rate was 1.27 mg Se/m(2)/d for Phragmites and 1.35 mg Se/m(2)/d for Typha. In the Typha bed Se migrated faster than in the Phragmites bed. After 25 d of Se supplementation in the Typha bed about 54% of the Se inlet concentration remained in the outlet water. In the Phragmites bed Se was removed completely from the water after passing through 3/4 of the bed length. After 65 d of Se supplementation the highest amount of Se (2.8 microg/g dry matter) was determined in the organic material of the Typha bed. Roots and rhizomes accumulated 2.2 and 1.8 microg/g dry matter respectively. Phragmites accumulated Se in the leaves and stems, but not in the rhizomes. The accumulation in the leaves (1.8 microg Se/g dry matter) was three times higher than in the stems (0.6 microg Se/g dry matter). Copyright 2002 Elsevier Science Ltd.

Performance assessment and microbial diversity of two pilot scale multi-stage sub-surface flow constructed wetland systems.

PubMed

Babatunde, A O; Miranda-CasoLuengo, Raul; Imtiaz, Mehreen; Zhao, Y Q; Meijer, Wim G

2016-08-01

This study assessed the performance and diversity of microbial communities in multi-stage sub-surface flow constructed wetland systems (CWs). Our aim was to assess the impact of configuration on treatment performance and microbial diversity in the systems. Results indicate that at loading rates up to 100gBOD5/(m(2)·day), similar treatment performances can be achieved using either a 3 or 4 stage configuration. In the case of phosphorus (P), the impact of configuration was less obvious and a minimum of 80% P removal can be expected for loadings up to 10gP/(m(2)·day) based on the performance results obtained within the first 16months of operation. Microbial analysis showed an increased bacterial diversity in stage four compared to the first stage. These results indicate that the design and configuration of multi-stage constructed wetland systems may have an impact on the treatment performance and the composition of the microbial community in the systems, and such knowledge can be used to improve their design and performance. Copyright © 2016. Published by Elsevier B.V.

The Cartridge Theory: a description of the functioning of horizontal subsurface flow constructed wetlands for wastewater treatment, based on modelling results.

PubMed

Samsó, Roger; García, Joan

2014-03-01

Despite the fact that horizontal subsurface flow constructed wetlands have been in operation for several decades now, there is still no clear understanding of some of their most basic internal functioning patterns. To fill this knowledge gap, on this paper we present what we call "The Cartridge Theory". This theory was derived from simulation results obtained with the BIO_PORE model and explains the functioning of urban wastewater treatment wetlands based on the interaction between bacterial communities and the accumulated solids leading to clogging. In this paper we start by discussing some changes applied to the biokinetic model implemented in BIO_PORE (CWM1) so that the growth of bacterial communities is consistent with a well-known population dynamics models. This discussion, combined with simulation results for a pilot wetland system, led to the introduction of "The Cartridge Theory", which states that the granular media of horizontal subsurface flow wetlands can be assimilated to a generic cartridge which is progressively consumed (clogged) with inert solids from inlet to outlet. Simulations also revealed that bacterial communities are poorly distributed within the system and that their location is not static but changes over time, moving towards the outlet as a consequence of the progressive clogging of the granular media. According to these findings, the life-span of constructed wetlands corresponds to the time when bacterial communities are pushed as much towards the outlet that their biomass is not anymore sufficient to remove the desirable proportion of the influent pollutants. Copyright © 2013 Elsevier B.V. All rights reserved.

Advantages of using subsurface flow constructed wetlands for wastewater treatment in space applications: Ground-based mars base prototype

NASA Astrophysics Data System (ADS)

Nelson, M.; Alling, A.; Dempster, W. F.; van Thillo, M.; Allen, John

Research and design of subsurface flow wetland wastewater treatment systems for a ground-based experimental prototype Mars Base facility has been carried out, using a subsurface flow approach. These systems have distinct advantages in planetary exploration scenarios: they are odorless, relatively low-labor and low-energy, assist in purification of water and recycling of atmospheric CO2, and will support some food crops. An area of 6-8 m2 may be sufficient for integration of wetland wastewater treatment with a prototype Mars Base supporting 4-5 people. Discharge water from the wetland system will be used as irrigation water for the agricultural crop area, thus ensuring complete recycling and utilization of nutrients. Since the primary requirements for wetland treatment systems are warm temperatures and lighting, such bioregenerative systems may be integrated into early Mars base habitats, since waste heat from the lights may be used for temperature maintenance in the human living environment. "Wastewater gardens ™" can be modified for space habitats to lower space and mass requirements. Many of its construction requirements can eventually be met with use of in-situ materials, such as gravel from the Mars surface. Because the technology requires little machinery and no chemicals, and relies more on natural ecological mechanisms (microbial and plant metabolism), maintenance requirements are minimized, and systems can be expected to have long operating lifetimes. Research needs include suitability of Martian soil and gravel for wetland systems, system sealing and liner options in a Mars Base, and wetland water quality efficiency under varying temperature and light regimes.

Bacterial transformation and biodegradation processes simulation in horizontal subsurface flow constructed wetlands using CWM1-RETRASO.

PubMed

Llorens, Esther; Saaltink, Maarten W; Poch, Manel; García, Joan

2011-01-01

The performance and reliability of the CWM1-RETRASO model for simulating processes in horizontal subsurface flow constructed wetlands (HSSF CWs) and the relative contribution of different microbial reactions to organic matter (COD) removal in a HSSF CW treating urban wastewater were evaluated. Various different approaches with diverse influent configurations were simulated. According to the simulations, anaerobic processes were more widespread in the simulated wetland and contributed to a higher COD removal rate [72-79%] than anoxic [0-1%] and aerobic reactions [20-27%] did. In all the cases tested, the reaction that most contributed to COD removal was methanogenesis [58-73%]. All results provided by the model were in consonance with literature and experimental field observations, suggesting a good performance and reliability of CWM1-RETRASO. According to the good simulation predictions, CWM1-RETRASO is the first mechanistic model able to successfully simulate the processes described by the CWM1 model in HSSF CWs. Copyright © 2010 Elsevier Ltd. All rights reserved.

Determination and removal of antibiotics in secondary effluent using a horizontal subsurface flow constructed wetland.

PubMed

Zhang, Chunhui; Ning, Ke; Zhang, Wenwen; Guo, Yuanjie; Chen, Jun; Liang, Chen

2013-04-01

Increased attention is currently being directed towards the potential negative effects of antibiotics and other PPCPs discharged into the aquatic environment via municipal WWTP secondary effluents. A number of analytical methods, such as high performance liquid chromatography technologies, including a high performance liquid chromatography-fluorescence method (HPLC-FLD), high performance liquid chromatography-UV detection method (HPLC-UV) and high performance liquid chromatography-mass spectrometry method (HPLC-MS), have been suggested as determination technologies for antibiotic residues in water. In this study, we implement a HPLC-MS/MS combined method to detect and analyze antibiotics in WWTP secondary effluent and apply a horizontal subsurface flow constructed wetland (CW) as an advanced wastewater treatment for removing antibiotics in the WWTP secondary effluent. The results show that there were 2 macrolides, 2 quinolones and 5 sulfas in WWTP secondary effluent among all the 22 antibiotics considered. After the CW advanced treatment, the concentration removal efficiencies and removal loads of 9 antibiotics were 53-100% and 0.004-0.7307 μg m(-2) per day, respectively.

A field study of colloid transport in surface and subsurface flows

NASA Astrophysics Data System (ADS)

Zhang, Wei; Tang, Xiang-Yu; Xian, Qing-Song; Weisbrod, Noam; Yang, Jae E.; Wang, Hong-Lan

2016-11-01

Colloids have been recognized to enhance the migration of strongly-sorbing contaminants. However, few field investigations have examined combined colloid transport via surface runoff and subsurface flows. In a headwater catchment of the upper Yangtze River, a 6 m (L) by 4 m (W) sloping (6°) farmland plot was built by cement walls to form no-flow side boundaries. The plot was monitored in the summer of 2014 for the release and transport of natural colloids via surface runoff and subsurface flows (i.e., the interflow from the soil-mudrock interface and fracture flow from the mudrock-sandstone interface) in response to rain events. The water sources of the subsurface flows were apportioned to individual rain events using a two end-member model (i.e., mobile pre-event soil water extracted by a suction-cup sampler vs. rainwater (event water)) based on δ18O measurements. For rain events with high preceding soil moisture, mobile pre-event soil water was the main contributor (generally >60%) to the fracture flow. The colloid concentration in the surface runoff was 1-2 orders of magnitude higher than that in the subsurface flows. The lowest colloid concentration was found in the subsurface interflow, which was probably the result of pore-scale colloid straining mechanisms. The rainfall intensity and its temporal variation govern the dynamics of the colloid concentrations in both surface runoff and subsurface flows. The duration of the antecedent dry period affected not only the relative contributions of the rainwater and the mobile pre-event soil water to the subsurface flows but also the peak colloid concentration, particularly in the fracture flow. The <10 μm fine colloid size fraction accounted for more than 80% of the total suspended particles in the surface runoff, while the colloid size distributions of both the interflow and the fracture flow shifted towards larger diameters. These results highlight the need to avoid the application of strongly

Advantages of using subsurface flow constructed wetlands for wastewater treatment in space applications: ground-based Mars Base prototype.

PubMed

Nelson, M; Alling, A; Dempster, W F; van Thillo, M; Allen, John

2003-01-01

Research and design of subsurface flow wetland wastewater treatment systems for a ground-based experimental prototype Mars Base facility has been carried out, using a subsurface flow approach. These systems have distinct advantages in planetary exploration scenarios: they are odorless, relatively low-labor and low-energy, assist in purification of water and recycling of atmospheric CO2, and will support some food crops. An area of 6-8 m2 may be sufficient for integration of wetland wastewater treatment with a prototype Mars Base supporting 4-5 people. Discharge water from the wetland system will be used as irrigation water for the agricultural crop area, thus ensuring complete recycling and utilization of nutrients. Since the primary requirements for wetland treatment systems are warm temperatures and lighting, such bioregenerative systems may be integrated into early Mars base habitats, since waste heat from the lights may be used for temperature maintenance in the human living environment. "Wastewater gardens (TM)" can be modified for space habitats to lower space and mass requirements. Many of its construction requirements can eventually be met with use of in-situ materials, such as gravel from the Mars surface. Because the technology requires little machinery and no chemicals, and relies more on natural ecological mechanisms (microbial and plant metabolism), maintenance requirements are minimized, and systems can be expected to have long operating lifetimes. Research needs include suitability of Martian soil and gravel for wetland systems, system sealing and liner options in a Mars Base, and wetland water quality efficiency under varying temperature and light regimes. c2003 COSPAR. Published by Elsevier Science Ltd. All rights reserved.

Understanding Subsurface Flow Mechanisms by Studying Recession Flow Curves

NASA Astrophysics Data System (ADS)

patnaik, S.; Biswal, B.; D, N.

2013-12-01

The recession flows offer valuable information on the subsurface systems of the drainage which cannot be observed due to technological limitations. Many analytical frameworks have been proposed in the past to analyze recession flow curves assess. Among them the most widely used one is Brutsaert-Neiber method of expressing negative time derivative of Q (discharge at the basin outlet at time t), -dQ/dt, as a function of Q itself, which eliminates the need of finding a reference time. Typically, basins across geographical regions display a power law relationship of the type: -dQ/dt = kQ^α. For a particular basin, the exponent α remains fairly constant recession events while the coefficient k varies greatly from one recession event to another, indicating the dynamic nature -dQ/dt-Q relationship. Recent observations show that subsurface storage in a basin mainly controls the dynamic parameter k. As subsurface water takes long time to fully drain, k of a recession event can also be influenced by the storage that occurred during the past rainfall events. We indirectly analyze the effect of past storage on recession flow by considering past streamflow as a proxy of past storage. A stronger relationship implies that the basin is able to store water for longer duration, and vice versa. In this study, we used streamflow data from 388 USGS basins that are relatively unaffected by human activities to find out the factors that affect the relationship between the power law correlation (R^2_PN) between past discharge and k, where the subscript N is the number of days of past streamflow observations considered for the recession event. For most of the basins R^2_PN decreases with N. We then selected 18 physical and climatological parameters for each study basin and investigated how they influence the value of R^2_PN for each N. We followed multiple linear regression method and found that R^2_PN is strongly influenced by the selected parameters (R^2 = 0.58) for N =30 days. We also

Gaseous fluxes from subsurface flow constructed wetlands for wastewater treatment.

PubMed

Mander, Ulo; Lõhmus, Krista; Teiter, Sille; Nurk, Kaspar; Mauring, Tõnu; Augustin, Jürgen

2005-01-01

We measured nitrous oxide (N2O), dinitrogen (N2), and methane (CH4) fluxes in two constructed wetlands (CW) in Estonia using the closed chamber method and the He-O method in the period from October 2000 to March 2003. Emission rates of N2O-N, N2-N and CH4-C from both CWs varied significantly on a both spatial and temporal scale, ranging from 1 to 2,600, 170 to 130,000, and -1.7 to 87,200 microg m(-2) h(-1) respectively. The average flux of N2O from the microsites in the Kodijärve horizontal subsurface flow (HSSF) CW and Kõo hybrid CW ranged from 27 to 370 and from 72 to 500 microg N2O-N m(-2) h(-1), respectively, whereas the average dinitrogen flux from the microsites in the HSSF CW in Kodijärve was 2-3 magnitudes higher than the N2O flux, ranging from 19,500 to 33,300 microg N2-N m(-2) h(-1). The average methane emissions from the microsites in the Kodijärve HSSF CW and the Kõo hybrid CW ranged from 31 to 12,100 and from 950 to 5,750 microg CH4-C m(-2) h(-1), respectively. The highest emission values for all three gases were observed in the warm period. There was a significant relationship between emission rates and water table depth: CH4 and N2 emission increased and N2O emission decreased when the water table did rise. Although the emission of N2O and CH4 from CWs was found to be relatively high, their global warming potential (GWP) in the time horizon of 100 years is not significant, ranging from 4.5 to 16.3 tonnes of CO2 equivalents per ha per year in Kodijärve and from 12.1 to 17.3 t CO2 equivalents ha(-1) yr(-1) in Kõo.

Influence of chlorothalonil on the removal of organic matter in horizontal subsurface flow constructed wetlands.

PubMed

Casas-Zapata, Juan C; Ríos, Karina; Florville-Alejandre, Tomás R; Morató, Jordi; Peñuela, Gustavo

2013-01-01

This study investigates the effects of chlorothalonil (CLT) on chemical oxygen demand (COD) and dissolved organic carbon (DOC) in pilot-scale horizontal subsurface flow constructed wetlands (HSSFCW) planted with Phragmites australis. Physicochemical parameters of influent and effluent water samples, microbial population counting methods and statistical analysis were used to evaluate the influence of CLT on organic matter removal efficiency. The experiments were conducted on four planted replicate wetlands (HSSFCW-Pa) and one unplanted control wetland (HSSFCW-NPa). The wetlands exhibited high average organic matter removal efficiencies (HSSFCW-Pa: 80.6% DOC, 98.0% COD; HSSFCW-NPa: 93.2% DOC, 98.4% COD). The addition of CLT did not influence organic removal parameters. In all cases CLT concentrations in the effluent occurred in concentrations lower than the detection limit of the analytical method. Microbial population counts from HSSFCW-Pa showed significant correlations among different microbial groups and with different physicochemical variables. The apparent independence of organic matter removal and CLT inputs, along with the CLT depletion observed in effluent samples demonstrated that HSSFCW are a viable technology for the treatment of agricultural effluents contaminated with organo-chloride pesticides like CLT.

Multi-stage hybrid subsurface flow constructed wetlands for treating piggery and dairy wastewater in cold climate.

PubMed

Zhang, Xiaomeng; Inoue, Takashi; Kato, Kunihiko; Izumoto, Hayato; Harada, June; Wu, Da; Sakuragi, Hiroaki; Ietsugu, Hidehiro; Sugawara, Yasuhide

2017-01-01

This study followed three field-scale hybrid subsurface flow constructed wetland (CW) systems constructed in Hokkaido, northern Japan: piggery O (2009), dairy G (2011), and dairy S (2006). Treatment performance was monitored from the outset of operation for each CW. The ranges of overall purification efficiency for these systems were 70-86%, 40-85%, 71-90%, 91-96%, 94-98%, 84-97%, and 70-97% for total N (TN), NH 4 -N, total P, chemical oxygen demand (COD), biochemical oxygen demand, suspended solid, and total Coliform, respectively. The hybrid system's removal rates were highest when influent loads were high. COD removal rates were 46.4 ± 49.2, 94.1 ± 36.6, and 25.1 ± 15.5 g COD m -2 d -1 in piggery O, dairy G, and dairy S, with average influent loads of 50.5 ± 51.5, 98.9 ± 37.1, and 26.9 ± 16.0 g COD m -2 d -1 , respectively. The systems had overall COD removal efficiencies of around 90%. TN removal efficiencies were 62 ± 19%, 82 ± 9%, and 82 ± 15% in piggery O, dairy G, and dairy S, respectively. NH 4 -N removal efficiency was adversely affected by the COD/TN ratio. Results from this study prove that these treatment systems have sustained and positive pollutant removal efficiencies, which were achieved even under extremely cold climate conditions and many years after initial construction.

Spatial Variation of Phosphorous Retention Capacity in Subsurface Flow Constructed Wetlands: Effect of Wetland Type and Inflow Loading.

PubMed

Yu, Guangwei; Tan, Meijuan; Chong, Yunxiao; Long, Xinxian

2015-01-01

For verification of spatial distribution of phosphorous retention capacity in constructed wetlands systems(CWs), two horizontal subsurface flow(HSSF) CWs and two vertical subsurface flow(VSSF) CWs, using sand as substrate and Typha latifolia as wetland plants, were constructed and put into use for synthetic wastewater treatment. Five months later, significant spatial variations of TP and inorganic phosphorus(Ca-P, Fe-P and Al-P) were observed, which were found to be greatly affected by CWs type and hydraulic loading. The results revealed that though spatial distribution of Fe-P and Al-P displayed a similar order of substrate content as "rhizosphere" > "near-rhizosphere" > "non-rhizosphere" and "inflow section" > "outflow section" regardless of types and loading, the distribution of Ca-P was positively correlated to that of Fe-P and Al-P in HSSF CWs, while negative correlation was shown in VSSF CWs. As a result, TP spatial distribution in HSSF CWs demonstrated a greater dissimilarity than that in VSSF CWs. For HSSF CWs with low hydraulic loading, the lowest TP content was found in non-rhizosphere substrate of outflow section, while the highest one was discovered in rhizonsphere substrate of inflow section. The values in 6 parts of areas ranged from 0.138 g·kg-1 to 2.710 g·kg-1, which also were from -33.5% to 1209% compared to the control value. On contrast, spatial difference of TP content in substrates of VSSF CWs was insignificant, with a variation ranging from 0.776 g·kg-1 to 1.080 g·kg-1, that was 275% to 421% higher than the control value. In addition, when hydraulic loading was increased, TP content in VSSF CWs sharply decreased, ranging from 0.210 g·kg-1 to 0.634 g·kg-1. Meanwhile, dissimilarity of TP spatial distribution in HSSF CWs was reduced, with TP content ranging from 0.258 g·kg-1 to 2.237 g·kg-1. The results suggested that P spatial distribution should be taken into account for CWs design and operation.

Spatial Variation of Phosphorous Retention Capacity in Subsurface Flow Constructed Wetlands: Effect of Wetland Type and Inflow Loading

PubMed Central

Yu, Guangwei; Tan, Meijuan; Chong, Yunxiao; Long, Xinxian

2015-01-01

For verification of spatial distribution of phosphorous retention capacity in constructed wetlands systems(CWs), two horizontal subsurface flow(HSSF) CWs and two vertical subsurface flow(VSSF) CWs, using sand as substrate and Typha latifolia as wetland plants, were constructed and put into use for synthetic wastewater treatment. Five months later, significant spatial variations of TP and inorganic phosphorus(Ca-P, Fe-P and Al-P) were observed, which were found to be greatly affected by CWs type and hydraulic loading. The results revealed that though spatial distribution of Fe-P and Al-P displayed a similar order of substrate content as "rhizosphere" > "near-rhizosphere" > "non-rhizosphere" and "inflow section" > "outflow section" regardless of types and loading, the distribution of Ca-P was positively correlated to that of Fe-P and Al-P in HSSF CWs, while negative correlation was shown in VSSF CWs. As a result, TP spatial distribution in HSSF CWs demonstrated a greater dissimilarity than that in VSSF CWs. For HSSF CWs with low hydraulic loading, the lowest TP content was found in non-rhizosphere substrate of outflow section, while the highest one was discovered in rhizonsphere substrate of inflow section. The values in 6 parts of areas ranged from 0.138 g·kg-1 to 2.710 g·kg-1, which also were from -33.5% to 1209% compared to the control value. On contrast, spatial difference of TP content in substrates of VSSF CWs was insignificant, with a variation ranging from 0.776 g·kg-1 to 1.080 g·kg-1, that was 275% to 421% higher than the control value. In addition, when hydraulic loading was increased, TP content in VSSF CWs sharply decreased, ranging from 0.210 g·kg-1 to 0.634 g·kg-1. Meanwhile, dissimilarity of TP spatial distribution in HSSF CWs was reduced, with TP content ranging from 0.258 g·kg-1 to 2.237 g·kg-1. The results suggested that P spatial distribution should be taken into account for CWs design and operation. PMID:26218872

Integrated Coupling of Surface and Subsurface Flow with HYDRUS-2D

NASA Astrophysics Data System (ADS)

Hartmann, Anne; Šimůnek, Jirka; Wöhling, Thomas; Schütze, Niels

2016-04-01

Describing interactions between surface and subsurface flow processes is important to adequately define water flow in natural systems. Since overland flow generation is highly influenced by rainfall and infiltration, both highly spatially heterogeneous processes, overland flow is unsteady and varies spatially. The prediction of overland flow needs to include an appropriate description of the interactions between the surface and subsurface flow. Coupling surface and subsurface water flow is a challenging task. Different approaches have been developed during the last few years, each having its own advantages and disadvantages. A new approach by Weill et al. (2009) to couple overland flow and subsurface flow based on a generalized Richards equation was implemented into the well-known subsurface flow model HYDRUS-2D (Šimůnek et al., 2011). This approach utilizes the one-dimensional diffusion wave equation to model overland flow. The diffusion wave model is integrated in HYDRUS-2D by replacing the terms of the Richards equation in a pre-defined runoff layer by terms defining the diffusion wave equation. Using this approach, pressure and flux continuity along the interface between both flow domains is provided. This direct coupling approach provides a strong coupling of both systems based on the definition of a single global system matrix to numerically solve the coupled flow problem. The advantage of the direct coupling approach, compared to the loosely coupled approach, is supposed to be a higher robustness, when many convergence problems can be avoided (Takizawa et al., 2014). The HYDRUS-2D implementation was verified using a) different test cases, including a direct comparison with the results of Weill et al. (2009), b) an analytical solution of the kinematic wave equation, and c) the results of a benchmark test of Maxwell et al. (2014), that included several known coupled surface subsurface flow models. Additionally, a sensitivity analysis evaluating the effects

Nitrogen removal in wood chip combined substrate baffled subsurface-flow constructed wetlands: impact of matrix arrangement and intermittent aeration.

PubMed

Li, Huai; Chi, Zifang; Yan, Baixing; Cheng, Long; Li, Jianzheng

2017-02-01

In this study, two lab-scale baffled subsurface-flow constructed wetlands (BSFCWs), including gravel-wood chips-slag and gravel-slag-wood chips, were operated at different intermittent aeration to evaluate the effect of artificial aeration and slow-released carbon source on the treatment efficiency of high-strength nitrogen wastewater. Results indicated that gravel-slag-wood chips extended aerobic/anaerobic alternating environment to gravel and slag zones and maintained anaerobic condition in the subsequent wood chip section. The order of gravel-slag-wood chip was more beneficial to pollutant removal. Sufficient carbon source supply resulted from wood-chip-framework substrate simultaneously obtained high removals of COD (97%), NH 4 + -N (95%), and TN (94%) in BSFCWs at 2 h aeration per day. The results suggest that intermittent aeration combined with wood chips could achieve high nitrogen removal in BSFCWs.

A Generalized Subsurface Flow Parameterization Considering Subgrid Spatial Variability of Recharge and Topography

DOE Office of Scientific and Technical Information (OSTI.GOV)

Huang, Maoyi; Liang, Xu; Leung, Lai R.

2008-12-05

Subsurface flow is an important hydrologic process and a key component of the water budget, especially in humid regions. In this study, a new subsurface flow formulation is developed that incorporates spatial variability of both topography and recharge. It is shown through theoretical derivation and case studies that the power law and exponential subsurface flow parameterizations and the parameterization proposed by Woods et al.[1997] are all special cases of the new formulation. The subsurface flows calculated using the new formulation compare well with values derived from observations at the Tulpehocken Creek and Walnut Creek watersheds. Sensitivity studies show that whenmore » the spatial variability of topography or recharge, or both is increased, the subsurface flows increase at the two aforementioned sites and the Maimai hillslope. This is likely due to enhancement of interactions between the groundwater table and the land surface that reduce the flow path. An important conclusion of this study is that the spatial variability of recharge alone, and/or in combination with the spatial variability of topography can substantially alter the behaviors of subsurface flows. This suggests that in macroscale hydrologic models or land surface models, subgrid variations of recharge and topography can make significant contributions to the grid mean subsurface flow and must be accounted for in regions with large surface heterogeneity. This is particularly true for regions with humid climate and relatively shallow groundwater table where the combined impacts of spatial variability of recharge and topography are shown to be more important. For regions with arid climate and relatively deep groundwater table, simpler formulations, especially the power law, for subsurface flow can work well, and the impacts of subgrid variations of recharge and topography may be ignored.« less

Researchers Mine Information from Next-Generation Subsurface Flow Simulations

DOE PAGES

Gedenk, Eric D.

2015-12-01

A research team based at Virginia Tech University leveraged computing resources at the US Department of Energy's (DOE's) Oak Ridge National Laboratory to explore subsurface multiphase flow phenomena that can't be experimentally observed. Using the Cray XK7 Titan supercomputer at the Oak Ridge Leadership Computing Facility, the team took Micro-CT images of subsurface geologic systems and created two-phase flow simulations. The team's model development has implications for computational research pertaining to carbon sequestration, oil recovery, and contaminant transport.

Effect of water depth on the removal of organic matter in horizontal subsurface flow constructed wetlands.

PubMed

Aguirre, Paula; Ojeda, Esther; García, Joan; Barragán, Jesús; Mujeriego, Rafael

2005-01-01

The objective of this article is to evaluate the effect of water depth on organic matter removal efficiency in horizontal subsurface flow constructed wetlands (SSFs). Experiments were carried out in a pilot plant comprising eight parallel SSF of almost equal surface area (54-56 m2 each) and treating urban wastewater. Each SSF differs from the others in the aspect ratio or the size of the granular medium or the water depth. During a period of two years, the shallow SSFs (0.27 m water depth) removed more chemical oxygen demand (COD) (72-81%), biochemical oxygen demand (BOD)5 (72-85%), ammonia (35-56%), and dissolved reactive phosphorus (DRP) (8-23%) than deep SSFs (0.5 m water depth) (59-64% for COD; 51-57% for BOD5; 18-29% for ammonia; and 0-7% for DRP). Experiments carried out during the summer indicated that sulphate reduction accounted for a clearly higher organic matter removal in the deep SSFs than in the shallow ones. Denitrification seemed to be a significant mechanism for organic matter removal to occur in shallow SSFs. The results suggest that the relative contribution of different metabolic pathways varies with depth.

Treatment of domestic wastewater by subsurface flow constructed wetlands filled with gravel and tire chip media.

PubMed

Richter, A Y; Weaver, R W

2003-12-01

Subsurface flow constructed wetlands (SFCWs) are becoming increasingly common in on-site treatment of wastewater. Gravel is the most popular form of wetland fill medium, but tire chips provide more porosity, are less dense, and less expensive. This study determines the treatment efficiency of SFCWs filled with gravel or tire chip media to treat domestic wastewater. The influent and effluent of six SFCWs filled with tire chip medium and six SFCWs filled with gravel were monitored for 5 to 16 consecutive months. Parameters measured included pH, biochemical oxygen demand (BOD5), total and volatile suspended solids, NH4, P, and fecal and total coliforms. The only clear difference between medium types in wetland performance was for P. Soluble P in the effluent averaged 1.6 +/- 1.0 mg l(-1) in the tire chip-filled wetlands and 4.8 +/- 3.2 mg l(-1) in the gravel-filled wetlands. Most likely, Fe from exposed wires in shredded steel-belted tires complexed with P to create an insoluble compound. Tire chips may be a better fill medium for SFCWs than gravel because of higher porosity, lower cost, and greater reduction of P in effluent.

Simulation and validation of concentrated subsurface lateral flow paths in an agricultural landscape

NASA Astrophysics Data System (ADS)

Zhu, Q.; Lin, H. S.

2009-08-01

The importance of soil water flow paths to the transport of nutrients and contaminants has long been recognized. However, effective means of detecting concentrated subsurface flow paths in a large landscape are still lacking. The flow direction and accumulation algorithm based on single-direction flow algorithm (D8) in GIS hydrologic modeling is a cost-effective way to simulate potential concentrated flow paths over a large area once relevant data are collected. This study tested the D8 algorithm for simulating concentrated lateral flow paths at three interfaces in soil profiles in a 19.5-ha agricultural landscape in central Pennsylvania, USA. These interfaces were (1) the interface between surface plowed layers of Ap1 and Ap2 horizons, (2) the interface with subsoil water-restricting clay layer where clay content increased to over 40%, and (3) the soil-bedrock interface. The simulated flow paths were validated through soil hydrologic monitoring, geophysical surveys, and observable soil morphological features. The results confirmed that concentrated subsurface lateral flow occurred at the interfaces with the clay layer and the underlying bedrock. At these two interfaces, the soils on the simulated flow paths were closer to saturation and showed more temporally unstable moisture dynamics than those off the simulated flow paths. Apparent electrical conductivity in the soil on the simulated flow paths was elevated and temporally unstable as compared to those outside the simulated paths. The soil cores collected from the simulated flow paths showed significantly higher Mn content at these interfaces than those away from the simulated paths. These results suggest that (1) the D8 algorithm is useful in simulating possible concentrated subsurface lateral flow paths if used with appropriate threshold value of contributing area and sufficiently detailed digital elevation model (DEM); (2) repeated electromagnetic surveys can reflect the temporal change of soil water storage

Need to improve SWMM's subsurface flow routing algorithm for green infrastructure modeling

EPA Science Inventory

SWMM can simulate various subsurface flows, including groundwater (GW) release from a subcatchment to a node, percolation out of storage units and low impact development (LID) controls, and rainfall derived inflow and infiltration (RDII) at a node. Originally, the subsurface flow...

Nitrogen and COD removal from domestic and synthetic wastewater in subsurface-flow constructed wetlands.

PubMed

Collison, R S; Grismer, M E

2013-09-01

Comparisons of the performance of constructed-wetland systems (CWs) for treating domestic wastewater in the laboratory and field may use pathogen-free synthetic wastewater to avoid regulatory health concerns. However, little to no data are available describing the relative treatment efficiencies of CWs to both actual and synthetic domestic wastewaters so as to enable such comparison. To fill this gap, treatment performances with respect to organics (chemical organic demand; COD) and nitrogen (ammonium and nitrate) removal from domestic (septic tank) and a similar-strength synthetic wastewater under planted and non-planted subsurface-flow CWs are determined. One pair of CWs was planted with cattails in May 2008, whereas the adjacent system was non-planted. Collected septic tank or synthesized wastewater was allowed to gravity feed each CWs, and effluent samples were collected and tested for COD and nitrogen species regularly during four different periods over six months. Overall, statistically significant greater removal of COD (-12%) and nitrogen (-5%) occurred from the synthetic as compared with the domestic wastewater from the planted and non-planted CWs. Effluent BOD5/COD ratios from the synthetic wastewater CWs averaged nearly twice that from the domestic wastewater CWs (0.17 vs 0.10), reflecting greater concentrations of readily degraded compounds. That removal fractions were consistent across the mid-range loading rates to the CWs suggests that the synthetic wastewater can be used in testing laboratory CWs with reasonable success in application of their results to the field.

Fractal topography and subsurface water flows from fluvial bedforms to the continental shield

USGS Publications Warehouse

Worman, A.; Packman, A.I.; Marklund, L.; Harvey, J.W.; Stone, S.H.

2007-01-01

Surface-subsurface flow interactions are critical to a wide range of geochemical and ecological processes and to the fate of contaminants in freshwater environments. Fractal scaling relationships have been found in distributions of both land surface topography and solute efflux from watersheds, but the linkage between those observations has not been realized. We show that the fractal nature of the land surface in fluvial and glacial systems produces fractal distributions of recharge, discharge, and associated subsurface flow patterns. Interfacial flux tends to be dominated by small-scale features while the flux through deeper subsurface flow paths tends to be controlled by larger-scale features. This scaling behavior holds at all scales, from small fluvial bedforms (tens of centimeters) to the continental landscape (hundreds of kilometers). The fractal nature of surface-subsurface water fluxes yields a single scale-independent distribution of subsurface water residence times for both near-surface fluvial systems and deeper hydrogeological flows. Copyright 2007 by the American Geophysical Union.

Dynamic coupling of subsurface and seepage flows solved within a regularized partition formulation

NASA Astrophysics Data System (ADS)

Marçais, J.; de Dreuzy, J.-R.; Erhel, J.

2017-11-01

Hillslope response to precipitations is characterized by sharp transitions from purely subsurface flow dynamics to simultaneous surface and subsurface flows. Locally, the transition between these two regimes is triggered by soil saturation. Here we develop an integrative approach to simultaneously solve the subsurface flow, locate the potential fully saturated areas and deduce the generated saturation excess overland flow. This approach combines the different dynamics and transitions in a single partition formulation using discontinuous functions. We propose to regularize the system of partial differential equations and to use classic spatial and temporal discretization schemes. We illustrate our methodology on the 1D hillslope storage Boussinesq equations (Troch et al., 2003). We first validate the numerical scheme on previous numerical experiments without saturation excess overland flow. Then we apply our model to a test case with dynamic transitions from purely subsurface flow dynamics to simultaneous surface and subsurface flows. Our results show that discretization respects mass balance both locally and globally, converges when the mesh or time step are refined. Moreover the regularization parameter can be taken small enough to ensure accuracy without suffering of numerical artefacts. Applied to some hundreds of realistic hillslope cases taken from Western side of France (Brittany), the developed method appears to be robust and efficient.

Subsurface flow and vegetation patterns in tidal environments

NASA Astrophysics Data System (ADS)

Ursino, Nadia; Silvestri, Sonia; Marani, Marco

2004-05-01

Tidal environments are characterized by a complex interplay of hydrological, geomorphic, and biological processes, and their understanding and modeling thus require the explicit description of both their biotic and abiotic components. In particular, the presence and spatial distribution of salt marsh vegetation (a key factor in the stabilization of the surface soil) have been suggested to be related to topographic factors and to soil moisture patterns, but a general, process-based comprehension of this relationship has not yet been achieved. The present paper describes a finite element model of saturated-unsaturated subsurface flow in a schematic salt marsh, driven by tidal fluctuations and evapotranspiration. The conditions leading to the establishment of preferentially aerated subsurface zones are studied, and inferences regarding the development and spatial distribution of salt marsh vegetation are drawn, with important implications for the overall ecogeomorphological dynamics of tidal environments. Our results show that subsurface water flow in the marsh induces complex water table dynamics, even when the tidal forcing has a simple sinusoidal form. The definition of a space-dependent aeration time is then proposed to characterize root aeration. The model shows that salt marsh subsurface flow depends on the distance from the nearest creek or channel and that the subsurface water movement near tidal creeks is both vertical and horizontal, while farther from creeks, it is primarily vertical. Moreover, the study shows that if the soil saturated conductivity is relatively low (10-6 m s-1, values quite common in salt marsh areas), a persistently unsaturated zone is present below the soil surface even after the tide has flooded the marsh; this provides evidence of the presence of an aerated layer allowing a prolonged presence of oxygen for aerobic root respiration. The results further show that plant transpiration increases the extent and persistence of the aerated

Characteristics of Nitrogen Loss through Surface-Subsurface Flow on Red Soil Slopes of Southeast China

NASA Astrophysics Data System (ADS)

Zheng, Haijin; Liu, Zhao; Zuo, Jichao; Wang, Lingyun; Nie, Xiaofei

2017-12-01

Soil nitrogen (N) loss related to surface flow and subsurface flow (including interflow and groundwater flow) from slope lands is a global issue. A lysimetric experiment with three types of land cover (grass cover, GC; litter cover, LC; and bare land, BL) were carried out on a red soil slope land in southeast China. Total Nitrogen (TN) loss through surface flow, interflow and groundwater flow was observed under 28 natural precipitation events from 2015 to 2016. TN concentrations from subsurface flow on BL and LC plots were, on average, 2.7-8.2 and 1.5-4.4 times greater than TN concentrations from surface flow, respectively; the average concentration of TN from subsurface flow on GC was about 36-56% of that recorded from surface flow. Surface flow, interflow and groundwater flow contributed 0-15, 2-9 and 76-96%, respectively, of loss load of TN. Compared with BL, GC and LC intercepted 83-86% of TN loss through surface runoff; GC intercepted 95% of TN loss through subsurface flow while TN loss through subsurface flow on LC is 2.3 times larger than that on BL. In conclusion, subsurface flow especially groundwater flow is the dominant hydrological rout for N loss that is usually underestimated. Grass cover has the high retention of N runoff loss while litter mulch will increase N leaching loss. These findings provide scientific support to control N runoff loss from the red soil slope lands by using suitable vegetation cover and mulching techniques.

Effects of rainfall patterns and land cover on the subsurface flow generation of sloping Ferralsols in southern China

PubMed Central

Yang, Jie; Tang, Chongjun; Chen, Lihua; Liu, Yaojun; Wang, Lingyun

2017-01-01

Rainfall patterns and land cover are two important factors that affect the runoff generation process. To determine the surface and subsurface flows associated with different rainfall patterns on sloping Ferralsols under different land cover types, observational data related to surface and subsurface flows from 5 m × 15 m plots were collected from 2010 to 2012. The experiment was conducted to assess three land cover types (grass, litter cover and bare land) in the Jiangxi Provincial Soil and Water Conservation Ecological Park. During the study period, 114 natural rainfall events produced subsurface flow and were divided into four groups using k-means clustering according to rainfall duration, rainfall depth and maximum 30-min rainfall intensity. The results showed that the total runoff and surface flow values were highest for bare land under all four rainfall patterns and lowest for the covered plots. However, covered plots generated higher subsurface flow values than bare land. Moreover, the surface and subsurface flows associated with the three land cover types differed significantly under different rainfall patterns. Rainfall patterns with low intensities and long durations created more subsurface flow in the grass and litter cover types, whereas rainfall patterns with high intensities and short durations resulted in greater surface flow over bare land. Rainfall pattern I had the highest surface and subsurface flow values for the grass cover and litter cover types. The highest surface flow value and lowest subsurface flow value for bare land occurred under rainfall pattern IV. Rainfall pattern II generated the highest subsurface flow value for bare land. Therefore, grass or litter cover are able to convert more surface flow into subsurface flow under different rainfall patterns. The rainfall patterns studied had greater effects on subsurface flow than on total runoff and surface flow for covered surfaces, as well as a greater effect on surface flows associated

Let's Go Off the Grid: Subsurface Flow Modeling With Analytic Elements

NASA Astrophysics Data System (ADS)

Bakker, M.

2017-12-01

Subsurface flow modeling with analytic elements has the major advantage that no grid or time stepping are needed. Analytic element formulations exist for steady state and transient flow in layered aquifers and unsaturated flow in the vadose zone. Analytic element models are vector-based and consist of points, lines and curves that represent specific features in the subsurface. Recent advances allow for the simulation of partially penetrating wells and multi-aquifer wells, including skin effect and wellbore storage, horizontal wells of poly-line shape including skin effect, sharp changes in subsurface properties, and surface water features with leaky beds. Input files for analytic element models are simple, short and readable, and can easily be generated from, for example, GIS databases. Future plans include the incorporation of analytic element in parts of grid-based models where additional detail is needed. This presentation will give an overview of advanced flow features that can be modeled, many of which are implemented in free and open-source software.

Pathways of nitrobenzene degradation in horizontal subsurface flow constructed wetlands: Effect of intermittent aeration and glucose addition.

PubMed

Kirui, Wesley K; Wu, Shubiao; Kizito, Simon; Carvalho, Pedro N; Dong, Renjie

2016-01-15

Intermittent aeration and addition of glucose were applied to horizontal subsurface flow constructed wetlands in order to investigate the effect on pathways of nitrobenzene (NB) degradation and interactions with microbial nitrogen and sulphur transformations. The experiment was carried out in three phases A, B and C consisting of different NB loading and glucose dosing. For each phase, the effect of aeration was assessed by intermittently aerating one wetland and leaving one unaerated. Regardless of whether or not the wetland was aerated, at an influent NB concentration of 140 mg/L, both wetlands significantly reduced NB to less than 2 mg/L, a reduction efficiency of 98%. However, once the influent NB concentration was increased to 280 mg/L, the aerated wetland had a higher removal performance 82% compared to that of the unaerated wetland 71%. Addition of glucose further intensified the NB removal to 95% in the aerated wetlands and 92% in the unaerated. Aeration of wetlands enhanced NB degradation, but also resulted in higher NB volatilization of 6 mg m(-2) d(-1). The detected high concentration of sulphide 20-60 mg/L in the unaerated wetland gave a strong indication that NB may act as an electron donor to sulphate-reducing bacteria, but this should be further investigated. Aeration positively improved NB removal in constructed wetlands, but resulted in higher NB volatilization. Glucose addition induced co-metabolism to enhance NB degradation. Copyright © 2015 Elsevier Ltd. All rights reserved.

Performance assessment of pilot horizontal sub-surface flow constructed wetlands for removal of diesel from wastewater by Scirpus grossus.

PubMed

Al-Baldawi, Israa Abdulwahab; Sheikh Abdullah, Siti Rozaimah; Anuar, Nurina; Suja, Fatihah; Idris, Mushrifah

2013-01-01

One of the appropriate development technology options for the treatment of wastewater contaminated with diesel is constructed wetlands (CWs). Throughout 72 days of exposure, sampling was carried out for monitoring of physical parameters, plant growth and the efficiency of total petroleum hydrocarbon (TPH) removal, as an indication for diesel contamination, to assess the pilot-scale performance. Four pilot CWs with a horizontal sub-surface flow system were applied using the bulrush of Scirpus grossus. The CWs were loaded with different diesel concentrations of 0, 0.1, 0.2 and 0.25% (Vdiesel/Vwater). The TPH removal efficiencies were 82, 71, and 67% at the end of 72 days for diesel concentrations of 0.1, 0.2, and 0.25% respectively. In addition, the high removal efficiency of total suspended solids and chemical oxygen demand (COD) were 100 and 75.4% respectively, for a diesel concentration of 0.1%. It was concluded that S. grossus is a potential plant that can be used in a well-operated CW for restoring 0.1% diesel-contaminated water.

Characterization of enterococci populations collected from a subsurface flow constructed wetland.

PubMed

Graves, A K; Weaver, R W

2010-04-01

The aim of this study was to identify and characterize the population of Enterococcus sp. in domestic wastewater as it flows through a constructed wetland. Four hundred and eighty-four Enterococcus isolates were collected from the inlet, various sites within and from the outlet of a plastic lined constructed wetland in College Station, TX. The wetland treated septic tank effluent that passed sequentially through two 1.89 m(3) septic tanks and a 1.89 m(3) pump tank allowing 48 l doses at a 24 l min(-1) rate. The Enterococcus isolates were identified to species using the commercial Biolog system. The 484 Enterococcus isolates were comprised of ten different species, including Enterococcus faecalis (30.6%), Enterococcus pseudoavium (24.0%), Enterococcus casseliflavus (12.8%), Enterococcus faecium (11.2%), Enterococcus mundtii (7.9%), Enterococcus gallinarum (6.2%), Enterococcus dispar (3.7%), Enterococcus hirae (2.1%), Enterococcus durans and Enterococcus flavescens both 0.8%. Of the 88 isolates collected from the inlet, only 9.1% of the isolates were identified as Ent. faecalis and Ent. pseudoavium (36.4%) was identified as the predominant species. Whereas of the 74 isolates collected from the outlet, the predominant species were identified as Ent. faecalis (29.7%). Species identification varied among sites within the wetland, but often Ent. faecalis was the predominant species. Our data suggest that while Ent. faecalis is the predominant species of Enterococcus found in domestic wastewater, the populations may shift during treatment as the wastewater flows through the constructed wetland. We found that shifts in Enterococcus species composition occurred during domestic wastewater treatment. This has implications for the identification of faecal pollution based on the presence of specific bacterial types associated with domestic wastewater.

Imaging lateral groundwater flow in the shallow subsurface using stochastic temperature fields

NASA Astrophysics Data System (ADS)

Fairley, Jerry P.; Nicholson, Kirsten N.

2006-04-01

Although temperature has often been used as an indication of vertical groundwater movement, its usefulness for identifying horizontal fluid flow has been limited by the difficulty of obtaining sufficient data to draw defensible conclusions. Here we use stochastic simulation to develop a high-resolution image of fluid temperatures in the shallow subsurface at Borax Lake, Oregon. The temperature field inferred from the geostatistical simulations clearly shows geothermal fluids discharging from a group of fault-controlled hydrothermal springs, moving laterally through the subsurface, and mixing with shallow subsurface flow originating from nearby Borax Lake. This interpretation of the data is supported by independent geochemical and isotopic evidence, which show a simple mixing trend between Borax Lake water and discharge from the thermal springs. It is generally agreed that stochastic simulation can be a useful tool for extracting information from complex and/or noisy data and, although not appropriate in all situations, geostatistical analysis may provide good definition of flow paths in the shallow subsurface. Although stochastic imaging techniques are well known in problems involving transport of species, e.g. delineation of contaminant plumes from soil gas survey data, we are unaware of previous applications to the transport of thermal energy for the purpose of inferring shallow groundwater flow.

Evaluation of the hydrological flow paths in a gravel bed filter modeling a horizontal subsurface flow wetland by using a multi-tracer experiment.

PubMed

Birkigt, Jan; Stumpp, Christine; Małoszewski, Piotr; Nijenhuis, Ivonne

2018-04-15

In recent years, constructed wetland systems have become into focus as means of cost-efficient organic contaminant management. Wetland systems provide a highly reactive environment in which several removal pathways of organic chemicals may be present at the same time; however, specific elimination processes and hydraulic conditions are usually separately investigated and thus not fully understood. The flow system in a three dimensional pilot-scale horizontal subsurface constructed wetland was investigated applying a multi-tracer test combined with a mathematical model to evaluate the flow and transport processes. The results indicate the existence of a multiple flow system with two distinct flow paths through the gravel bed and a preferential flow at the bottom transporting 68% of tracer mass resulting from the inflow design of the model wetland system. There the removal of main contaminant chlorobenzene was up to 52% based on different calculation approaches. Determined retention times in the range of 22d to 32.5d the wetland has a heterogeneous flow pattern. Differences between simulated and measured tracer concentrations in the upper sediment indicate diffusion dominated processes due to stagnant water zones. The tracer study combining experimental evaluation with mathematical modeling demonstrated the complexity of flow and transport processes in the constructed wetlands which need to be taken into account during interpretation of the determining attenuation processes. Copyright © 2017 Elsevier B.V. All rights reserved.

Dissolved phosphorus retention of light-weight expanded shale and masonry sand used in subsurface flow treatment wetlands.

PubMed

Forbes, Margaret G; Dickson, Kenneth R; Golden, Teresa D; Hudak, Paul; Doyle, Robert D

2004-02-01

Using surface flow constructed wetlands for long-term phosphorus (P) retention presents a challenge due to the fact that P is stored primarily in the sediments. Subsurface flow wetlands have the potential to greatly increase P retention; however, the substrate needs to have both high hydraulic conductivity and high P sorption capacity. The objective of our study was to assess the P retention capacity of two substrates, masonry sand and lightweight expanded shale. We used sorption/desorption isotherms, flow-through column experiments, and pilot-scale wetlands to quantify P retained from treated municipal wastewater. Langmuir sorption isotherms predicted that the expanded shale has a maximum sorption capacity of 971 mg/kg and the masonry sand 58.8 mg/kg. In column desorption and column flow-through experiments, the masonry sand desorbed P when exposed to dilute P solutions. The expanded shale, however, had very little desorption and phosphorus did not break through the columns during our experiment. In pilot cells, masonry sand retained (mean +/- standard deviation) 45 +/- 62 g P/m2/yr and expanded shale retained 164 +/- 110 g P/m2/yr. We conclude that only the expanded shale would be a suitable substrate for retaining P in a subsurface flow wetland.

Investigation of lab-scale horizontal subsurface flow constructed wetlands treating industrial cork boiling wastewater.

PubMed

Gomes, Arlindo C; Silva, Lúcia; Albuquerque, António; Simões, Rogério; Stefanakis, Alexandros I

2018-09-01

The feasibility and treatment efficiency of horizontal subsurface flow constructed wetlands (HSFCW) was assessed for the first time for cork boiling wastewater (CBW) through laboratory experiments. CBW is known for its high content of phenolic compounds, complex composition of biorecalcitrant and toxic nature. Two lab-scale units, one planted with Phragmites australis (CWP) and one unplanted (CWC), were used to evaluate the removals of COD, BOD, total phenolic compounds (TPh) and decolourization over a 2.5-years monitoring period under Mediterranean climatic conditions. Seven organic and hydraulic loading rates ranging from 2.6 to 11.5 g COD/m 2 /d and 5.7-9.1 L/m 2 /d were tested under average hydraulic retention time (HRT) of 5 ± 1 days required due to the CWB limited biodegradability (i.e., BOD 5 /COD of 0.19). Average removals of the CWP exceeded those of the CWC and reached 74.6%, 91.7% and 69.1% for COD, BOD 5 and TPh, respectively, with respective mass removals rates up to 7.0, 1.7 and 0.5 (in g/m 2 /d). Decolourization was limited to 35%, since it mainly depends on physical processes rather than biodegradation. CBW concentration of nine phenolic compounds ranged from 1.2 to 38.4 mg/L (for the syringic and ellagic acids, respectively) in the raw CBW, with respective removals in the CWP unit ranging from 41.8 to 76.3%, higher than those in the control unit. Despite CBW high concentration of TPhs (average of 116.3 mg/L), the HSFCW reached organic load removals higher than those of conventional biological treatment methods. Copyright © 2018 Elsevier Ltd. All rights reserved.

Phytoextraction, phytotransformation and rhizodegradation of ibuprofen associated with Typha angustifolia in a horizontal subsurface flow constructed wetland.

PubMed

Li, Yifei; Zhang, Jiefeng; Zhu, Guibing; Liu, Yu; Wu, Bing; Ng, Wun Jern; Appan, Adhityan; Tan, Soon Keat

2016-10-01

Widespread occurrence of trace pharmaceutical residues in aquatic environments is of great concerns due to the potential chronic toxicity of certain pharmaceuticals including ibuprofen on aquatic organisms even at environmental levels. In this study, the phytoextraction, phytotransformation and rhizodegradation of ibuprofen associated with Typha angustifolia were investigated in a horizontal subsurface flow constructed wetland system. The experimental wetland system consisted of a planted bed with Typha angustifolia and an unplanted bed (control) to treat ibuprofen-loaded wastewater (∼107.2 μg L(-1)). Over a period of 342 days, ibuprofen was accumulated in leaf sheath and lamina tissues at a mean concentration of 160.7 ng g(-1), indicating the occurrence of the phytoextraction of ibuprofen. Root-uptake ibuprofen was partially transformed to ibuprofen carboxylic acid, 2-hydroxy ibuprofen and 1-hydroxy ibuprofen which were found to be 1374.9, 235.6 and 301.5 ng g(-1) in the sheath, respectively, while they were 1051.1, 693.6 and 178.7 ng g(-1) in the lamina. The findings from pyrosequencing analysis of the rhizosphere bacteria suggest that the Dechloromonas sp., the Clostridium sp. (e.g. Clostridium saccharobutylicum), the order Sphingobacteriales, and the Cytophaga sp. in the order Cytophagales were most probably responsible for the rhizodegradation of ibuprofen. Copyright © 2016 Elsevier Ltd. All rights reserved.

[Community Characteristics of ANAMMOX Bacteria in Subsurface Flow Constructed Wetland (SSFCW) for Processing of Aquaculture Waster Water].

PubMed

Zeng, Xian-lei; Liu, Xing-guo; Wu, Zong-fan; Shi, Xu; Lu, Shi-min

2016-02-15

Anaerobic ammonium oxidation (ANAMMOX) is one of the important functions in waste water treatment by subsurface flow constructed wetland (SSFCW), however, there are few studies on ANAMMOX in SSFCW environment at present. The community characteristics of ANAMMOX in the SSFCW of processing aquaculture waste water were explored in this study. In order to analyze the structure, diversity and abundance of ANAMMOX bacteria, several 16S rRNA clone libraries were constructed and real-time PCR targeting specific 16S rRNA genes together with diversity analysis was adopted. The obtained results showed that the SSFCW identified a total of three unknown clusters and two known clusters including Candidatus brocadia and Candidatus kuenenia. The dominant cluster was Candidatus brocadia. The highest diversity levels of ANAMMOX bacteria occurred in autumn (H', 1.21), while the lowest in spring (H', 0.64). The abundance of ANAMMOX bacteria in SSFCW environment ranged from 8.0 x 10(4) to 9.4 x 10(6) copies x g(-1) of fresh weight and the copy number of total bacterial 16S rRNA genes ranged from 7.3 x 10(9) to 9.1 x 10(10) copies x g(-1) of fresh weight during culture cycle. There were significant differences in the ANAMMOX bacteria abundances of different stratum and seasons in SSFCW environment, but the differences in total bacterial abundances were not obvious. In addition, the differences in ANAMMOX bacteria abundances in different stratum and seasons in SSFCW environment were irregular in different culture cycle. According to the distribution characteristics of ANAMMOX bacteria in the wetland, the denitrification effect of SSFCW could be improved by changing the supplying manners of aquaculture wastewater and adjusting the structure of wetland. The research results will provide reference for further optimizing the SSFCW and improving the efficiency of purification.

Controlling Subsurface Fractures and Fluid Flow: A Basic Research Agenda

DOE Office of Scientific and Technical Information (OSTI.GOV)

Pyrak-Nolte, Laura J; DePaolo, Donald J.; Pietraß, Tanja

2015-05-22

From beneath the surface of the earth, we currently obtain about 80-percent of the energy our nation consumes each year. In the future we have the potential to generate billions of watts of electrical power from clean, green, geothermal energy sources. Our planet’s subsurface can also serve as a reservoir for storing energy produced from intermittent sources such as wind and solar, and it could provide safe, long-term storage of excess carbon dioxide, energy waste products and other hazardous materials. However, it is impossible to underestimate the complexities of the subsurface world. These complexities challenge our ability to acquire themore » scientific knowledge needed for the efficient and safe exploitation of its resources. To more effectively harness subsurface resources while mitigating the impacts of developing and using these resources, the U.S. Department of Energy established SubTER – the Subsurface Technology and Engineering RD&D Crosscut team. This DOE multi-office team engaged scientists and engineers from the national laboratories to assess and make recommendations for improving energy-related subsurface engineering. The SubTER team produced a plan with the overall objective of “adaptive control of subsurface fractures and fluid flow.”This plan revolved around four core technological pillars—Intelligent Wellbore Systems that sustain the integrity of the wellbore environment; Subsurface Stress and Induced Seismicity programs that guide and optimize sustainable energy strategies while reducing the risks associated with subsurface injections; Permeability Manipulation studies that improve methods of enhancing, impeding and eliminating fluid flow; and New Subsurface Signals that transform our ability to see into and characterize subsurface systems. The SubTER team developed an extensive R&D plan for advancing technologies within these four core pillars and also identified several areas where new technologies would require additional basic

Nitrogen removal and its relationship with the nitrogen-cycle genes and microorganisms in the horizontal subsurface flow constructed wetlands with different design parameters.

PubMed

Chen, Jun; Ying, Guang-Guo; Liu, You-Sheng; Wei, Xiao-Dong; Liu, Shuang-Shuang; He, Liang-Ying; Yang, Yong-Qiang; Chen, Fan-Rong

2017-07-03

This study aims to investigate nitrogen removal and its relationship with the nitrogen-cycle genes and microorganisms in the horizontal subsurface flow constructed wetlands (CWs) with different design parameters. Twelve mesocosm-scale CWs with four substrates and three hydraulic loading rates were set up in the outdoor. The result showed the CWs with zeolite as substrate and HLR of 20 cm/d were selected as the best choice for the TN and NH 3 -N removal. It was found that the single-stage mesocosm-scale CWs were incapable to achieve high removals of TN and NH 3 -N due to inefficient nitrification process in the systems. This was demonstrated by the lower abundance of the nitrification genes (AOA and AOB) than the denitrification genes (nirK and nirS), and the less diverse nitrification microorganisms than the denitrification microorganisms in the CWs. The results also show that microorganism community structure including nitrogen-cycle microorganisms in the constructed wetland systems was affected by the design parameters especially the substrate type. These findings show that nitrification is a limiting factor for the nitrogen removal by CWs.

Mechanistic assessment of hillslope transpiration controls of diel subsurface flow: a steady-state irrigation approach

Treesearch

H.R. Barnard; C.B. Graham; W.J. van Verseveld; J.R. Brooks; B.J. Bond; J.J. McDonnell

2010-01-01

Mechanistic assessment of how transpiration influences subsurface flow is necessary to advance understanding of catchment hydrology. We conducted a 24-day, steady-state irrigation experiment to quantify the relationships among soil moisture, transpiration and hillslope subsurface flow. Our objectives were to: (1) examine the time lag between maximum transpiration and...

A one-dimensional model of subsurface hillslope flow

Treesearch

Jason C. Fisher

1997-01-01

Abstract - A one-dimensional, finite difference model of saturated subsurface flow within a hillslope was developed. The model uses rainfall, elevation data, a hydraulic conductivity, and a storage coefficient to predict the saturated thickness in time and space. The model was tested against piezometric data collected in a swale located in the headwaters of the North...

Discharge-nitrate data clustering for characterizing surface-subsurface flow interaction and calibration of a hydrologic model

NASA Astrophysics Data System (ADS)

Shrestha, R. R.; Rode, M.

2008-12-01

Concentration of reactive chemicals has different chemical signatures in baseflow and surface runoff. Previous studies on nitrate export from a catchment indicate that the transport processes are driven by subsurface flow. Therefore nitrate signature can be used for understanding the event and pre-event contributions to streamflow and surface-subsurface flow interactions. The study uses flow and nitrate concentration time series data for understanding the relationship between these two variables. Unsupervised artificial neural network based learning method called self organizing map is used for the identification of clusters in the datasets. Based on the cluster results, five different pattern in the datasets are identified which correspond to (i) baseflow, (ii) subsurface flow increase, (iii) surface runoff increase, (iv) surface runoff recession, and (v) subsurface flow decrease regions. The cluster results in combination with a hydrologic model are used for discharge separation. For this purpose, a multi-objective optimization tool NSGA-II is used, where violation of cluster results is used as one of the objective functions. The results show that the use of cluster results as supplementary information for the calibration of a hydrologic model gives a plausible simulation of subsurface flow as well total runoff at the catchment outlet. The study is undertaken using data from the Weida catchment in the North-Eastern Germany, which is a sub-catchment of the Weisse Elster river in the Elbe river basin.

Shallow subsurface storm flow in a forested headwater catchment: Observations and modeling using a modified TOPMODEL

USGS Publications Warehouse

Scanlon, Todd M.; Raffensperger, Jeff P.; Hornberger, George M.; Clapp, Roger B.

2000-01-01

Transient, perched water tables in the shallow subsurface are observed at the South Fork Brokenback Run catchment in Shenandoah National Park, Virginia. Crest piezometers installed along a hillslope transect show that the development of saturated conditions in the upper 1.5 m of the subsurface is controlled by total precipitation and antecedent conditions, not precipitation intensity, although soil heterogeneities strongly influence local response. The macroporous subsurface storm flow zone provides a hydrological pathway for rapid runoff generation apart from the underlying groundwater zone, a conceptualization supported by the two‐storage system exhibited by hydrograph recession analysis. A modified version of TOPMODEL is used to simulate the observed catchment dynamics. In this model, generalized topographic index theory is applied to the subsurface storm flow zone to account for logarithmic storm flow recessions, indicative of linearly decreasing transmissivity with depth. Vertical drainage to the groundwater zone is required, and both subsurface reservoirs are considered to contribute to surface saturation.

Sequential geophysical and flow inversion to characterize fracture networks in subsurface systems

DOE PAGES

Mudunuru, Maruti Kumar; Karra, Satish; Makedonska, Nataliia; ...

2017-09-05

Subsurface applications, including geothermal, geological carbon sequestration, and oil and gas, typically involve maximizing either the extraction of energy or the storage of fluids. Fractures form the main pathways for flow in these systems, and locating these fractures is critical for predicting flow. However, fracture characterization is a highly uncertain process, and data from multiple sources, such as flow and geophysical are needed to reduce this uncertainty. We present a nonintrusive, sequential inversion framework for integrating data from geophysical and flow sources to constrain fracture networks in the subsurface. In this framework, we first estimate bounds on the statistics formore » the fracture orientations using microseismic data. These bounds are estimated through a combination of a focal mechanism (physics-based approach) and clustering analysis (statistical approach) of seismic data. Then, the fracture lengths are constrained using flow data. In conclusion, the efficacy of this inversion is demonstrated through a representative example.« less

Sequential geophysical and flow inversion to characterize fracture networks in subsurface systems

DOE Office of Scientific and Technical Information (OSTI.GOV)

Mudunuru, Maruti Kumar; Karra, Satish; Makedonska, Nataliia

Subsurface applications, including geothermal, geological carbon sequestration, and oil and gas, typically involve maximizing either the extraction of energy or the storage of fluids. Fractures form the main pathways for flow in these systems, and locating these fractures is critical for predicting flow. However, fracture characterization is a highly uncertain process, and data from multiple sources, such as flow and geophysical are needed to reduce this uncertainty. We present a nonintrusive, sequential inversion framework for integrating data from geophysical and flow sources to constrain fracture networks in the subsurface. In this framework, we first estimate bounds on the statistics formore » the fracture orientations using microseismic data. These bounds are estimated through a combination of a focal mechanism (physics-based approach) and clustering analysis (statistical approach) of seismic data. Then, the fracture lengths are constrained using flow data. In conclusion, the efficacy of this inversion is demonstrated through a representative example.« less

Subsurface conditions in hydrothermal vents inferred from diffuse flow composition, and models of reaction and transport

NASA Astrophysics Data System (ADS)

Larson, B. I.; Houghton, J. L.; Lowell, R. P.; Farough, A.; Meile, C. D.

2015-08-01

Chemical gradients in the subsurface of mid-ocean ridge hydrothermal systems create an environment where minerals precipitate and dissolve and where chemosynthetic organisms thrive. However, owing to the lack of easy access to the subsurface, robust knowledge of the nature and extent of chemical transformations remains elusive. Here, we combine measurements of vent fluid chemistry with geochemical and transport modeling to give new insights into the under-sampled subsurface. Temperature-composition relationships from a geochemical mixing model are superimposed on the subsurface temperature distribution determined using a heat flow model to estimate the spatial distribution of fluid composition. We then estimate the distribution of Gibb's free energies of reaction beneath mid oceanic ridges and by combining flow simulations with speciation calculations estimate anhydrite deposition rates. Applied to vent endmembers observed at the fast spreading ridge at the East Pacific Rise, our results suggest that sealing times due to anhydrite formation are longer than the typical time between tectonic and magmatic events. The chemical composition of the neighboring low temperature flow indicates relatively uniform energetically favorable conditions for commonly inferred microbial processes such as methanogenesis, sulfate reduction and numerous oxidation reactions, suggesting that factors other than energy availability may control subsurface microbial biomass distribution. Thus, these model simulations complement fluid-sample datasets from surface venting and help infer the chemical distribution and transformations in subsurface flow.

Groundwater Flow Through a Constructed Treatment Wetland

DTIC Science & Technology

2002-03-01

sediments or has the water found preferential flow paths? (2) Does the behavior of groundwater flow change with varying loading rates or environmental...surface of the wetland. Water flows through a subsurface flow wetland in a similar fashion as groundwater flows through an aquifer. The concept is...circuiting of the wetland media. Groundwater Flow Various physical properties influence the flow of water through soil. In wetlands, the type of soil

Application of subsurface vertical flow constructed wetlands to reject water treatment in dairy wastewater treatment plant.

PubMed

Dąbrowski, Wojciech; Karolinczak, Beata; Gajewska, Magdalena; Wojciechowska, Ewa

2017-01-01

The paper presents the effects of applying subsurface vertical flow constructed wetlands (SS VF) for the treatment of reject water generated in the process of aerobic sewage sludge stabilization in the biggest dairy wastewater treatment plant (WWTP) in Poland. Two SS VF beds were built: bed (A) with 0.65 m depth and bed (B) with 1.0 m depth, planted with reeds. Beds were fed with reject water with hydraulic load of 0.1 m d -1 in order to establish the differences in treatment efficiency. During an eight-months research period, a high removal efficiency of predominant pollutants was shown: BOD 5 88.1% (A) and 90.5% (B); COD 84.5% (A) and 87.5% (B); TSS 87.6% (A) and 91.9% (B); TKN 82.4% (A) and 76.5% (B); N-NH 4 + 89.2% (A) and 85.7% (B); TP 30.2% (A) and 40.6% (B). There were not statistically significant differences in the removal efficiencies between bed (B) with 1.0 m depth and bed (A) with 0.65 m depth. The research indicated that SS VF beds could be successfully applied to reject water treatment in dairy WWTPs. The study proved that the use of SS VF beds in full scale in dairy WWTPs would result in a significant decrease in pollutants' load in reject water. In the analyzed case, decreasing the load of ammonia nitrogen was of greatest importance, as it constituted 58% of the total load treated in dairy WWTP and posed a hazard to the stability of the treatment process.

Real rock-microfluidic flow cell: A test bed for real-time in situ analysis of flow, transport, and reaction in a subsurface reactive transport environment.

PubMed

Singh, Rajveer; Sivaguru, Mayandi; Fried, Glenn A; Fouke, Bruce W; Sanford, Robert A; Carrera, Martin; Werth, Charles J

2017-09-01

Physical, chemical, and biological interactions between groundwater and sedimentary rock directly control the fundamental subsurface properties such as porosity, permeability, and flow. This is true for a variety of subsurface scenarios, ranging from shallow groundwater aquifers to deeply buried hydrocarbon reservoirs. Microfluidic flow cells are now commonly being used to study these processes at the pore scale in simplified pore structures meant to mimic subsurface reservoirs. However, these micromodels are typically fabricated from glass, silicon, or polydimethylsiloxane (PDMS), and are therefore incapable of replicating the geochemical reactivity and complex three-dimensional pore networks present in subsurface lithologies. To address these limitations, we developed a new microfluidic experimental test bed, herein called the Real Rock-Microfluidic Flow Cell (RR-MFC). A porous 500μm-thick real rock sample of the Clair Group sandstone from a subsurface hydrocarbon reservoir of the North Sea was prepared and mounted inside a PDMS microfluidic channel, creating a dynamic flow-through experimental platform for real-time tracking of subsurface reactive transport. Transmitted and reflected microscopy, cathodoluminescence microscopy, Raman spectroscopy, and confocal laser microscopy techniques were used to (1) determine the mineralogy, geochemistry, and pore networks within the sandstone inserted in the RR-MFC, (2) analyze non-reactive tracer breakthrough in two- and (depth-limited) three-dimensions, and (3) characterize multiphase flow. The RR-MFC is the first microfluidic experimental platform that allows direct visualization of flow and transport in the pore space of a real subsurface reservoir rock sample, and holds potential to advance our understandings of reactive transport and other subsurface processes relevant to pollutant transport and cleanup in groundwater, as well as energy recovery. Copyright © 2017 Elsevier B.V. All rights reserved.

Subsurface temperatures and surface heat flow in the Michigan Basin and their relationships to regional subsurface fluid movement

USGS Publications Warehouse

Vugrinovich, R.

1989-01-01

Linear regression of 405 bottomhole temperature (BHT) measurements vs. associated depths from Michigan's Lower Peninsula results in the following equation relating BHT and depth: BHT(??C) = 14.5 + 0.0192 ?? depth(m) Temperature residuals, defined as (BHT measured)-(BHT calculated), were determined for each of the 405 BHT's. Areas of positive temperature residuals correspond to areas of regional groundwater discharge (determined from maps of equipotential surface) while areas of negative temperature residuals correspond to areas of regional groundwater recharge. These relationships are observed in the principal aquifers in rocks of Devonian and Ordovician age and in a portion of the principal aquifer in rocks of Silurian age. There is a similar correspondence between high surface heat flow (determined using the silica geothermometer) and regional groundwater discharge areas and low surface heat flow and regional groundwater recharge areas. Post-Jurassic depositional and tectonic histories suggest that the observed coupling of subsurface temperature and groundwater flow systems may have persisted since Jurassic time. Thus the higher subsurface palaeotemperatures (and palaeogeothermal gradients) indicated by recent studies most likely pre-date the Jurassic. ?? 1989.

Development of stream-subsurface flow module in sub-daily simulation of Escherichia coli using SWAT

NASA Astrophysics Data System (ADS)

Kim, Minjeong; Boithias, Laurie; Cho, Kyung Hwa; Silvera, Norbert; Thammahacksa, Chanthamousone; Latsachack, Keooudone; Rochelle-Newall, Emma; Sengtaheuanghoung, Oloth; Pierret, Alain; Pachepsky, Yakov A.; Ribolzi, Olivier

2017-04-01

Water contaminated with pathogenic bacteria poses a large threat to public health, especially in the rural areas in the tropics where sanitation and drinking water facilities are often lacking. Several studies have used the Soil and Water Assessment Tool (SWAT) to predict the export of in-stream bacteria at a watershed-scale. However, SWAT is limited to in-stream processes, such as die-off, resuspension and, deposition; and it is usually implemented on a daily time step using the SCS Curve Number method, making it difficult to explore the dynamic fate and transport of bacteria during short but intense events such as flash floods in tropical humid montane headwaters. To address these issues, this study implemented SWAT on an hourly time step using the Green-Ampt infiltration method, and tested the effects of subsurface flow (LATQ+GWQ in SWAT) on bacterial dynamics. We applied the modified SWAT model to the 60-ha Houay Pano catchment in Northern Laos, using sub-daily rainfall and discharge measurements, electric conductivity-derived fractions of overland and subsurface flows, suspended sediments concentrations, and the number of fecal indicator organism Escherichia coli monitored at the catchment outlet from 2011 to 2013. We also took into account land use change by delineating the watershed with the 3-year composite land use map. The results show that low subsurface flow of less than 1 mm recovered the underestimation of E. coli numbers during the dry season, while high subsurface flow caused an overestimation during the wet season. We also found that it is more reasonable to apply the stream-subsurface flow interaction to simulate low in-stream bacteria counts. Using fecal bacteria to identify and understand the possible interactions between overland and subsurface flows may well also provide some insight into the fate of other bacteria, such as those involved in biogeochemical fluxes both in-stream and in the adjacent soils and hyporheic zones.

Heat flow and subsurface temperature as evidence for basin-scale ground-water flow, North Slope of Alaska

USGS Publications Warehouse

Deming, D.; Sass, J.H.; Lachenbruch, A.H.; De Rito, R. F.

1992-01-01

Several high-resolution temperature logs were made in each of 21 drillholes and a total of 601 thermal conductivity measurements were made on drill cuttings and cores. Near-surface heat flow (??20%) is inversely correlated with elevation and ranges from a low of 27 mW/m2 in the foothills of the Brooks Range in the south, to a high of 90 mW/m2 near the north coast. Subsurface temperatures and thermal gradients estimated from corrected BHTs are similarly much higher on the coastal plain than in the foothills province to the south. Significant east-west variation in heat flow and subsurface temperature is also observed; higher heat flow and temperature coincide with higher basement topography. The observed thermal pattern is consistent with forced convection by a topographically driven ground-water flow system. Average ground-water (Darcy) velocity in the postulated flow system is estimated to be of the order of 0.1 m/yr; the effective basin-scale permeability is estimated to be of the order of 10-14 m2. -from Authors

Numerical simulation for horizontal subsurface flow constructed wetlands: A short review including geothermal effects and solution bounding in biodegradation procedures

NASA Astrophysics Data System (ADS)

Liolios, K.; Tsihrintzis, V.; Angelidis, P.; Georgiev, K.; Georgiev, I.

2016-10-01

Current developments on modeling of groundwater flow and contaminant transport and removal in the porous media of Horizontal Subsurface Flow Constructed Wetlands (HSF CWs) are first reviewed in a short way. The two usual environmental engineering approaches, the black-box and the process-based one, are briefly presented. Next, recent research results obtained by using these two approaches are briefly discussed as application examples, where emphasis is given to the evaluation of the optimal design and operation parameters concerning HSF CWs. For the black-box approach, the use of Artificial Neural Networks is discussed for the formulation of models, which predict the removal performance of HSF CWs. A novel mathematical prove is presented, which concerns the dependence of the first-order removal coefficient on the Temperature and the Hydraulic Residence Time. For the process-based approach, an application example is first discussed which concerns procedures to evaluate the optimal range of values for the removal coefficient, dependent on either the Temperature or the Hydraulic Residence Time. This evaluation is based on simulating available experimental results of pilot-scale units operated in Democritus University of Thrace, Xanthi, Greece. Further, in a second example, a novel enlargement of the system of Partial Differential Equations is presented, in order to include geothermal effects. Finally, in a third example, the case of parameters uncertainty concerning biodegradation procedures is considered and the use of upper and a novel approach is presented, which concerns the upper and the lower solution bound for the practical draft design of HSF CWs.

Solar-Cycle Variation of Subsurface-Flow Divergence: A Proxy of Magnetic Activity?

NASA Astrophysics Data System (ADS)

Komm, R.; Howe, R.; Hill, F.

2017-09-01

We study the solar-cycle variation of subsurface flows from the surface to a depth of 16 Mm. We have analyzed Global Oscillation Network Group (GONG) Dopplergrams with a ring-diagram analysis covering about 15 years and Helioseismic and Magnetic Imager (HMI) Dopplergrams covering more than 6 years. After subtracting the average rotation rate and meridional flow, we have calculated the divergence of the horizontal residual flows from the maximum of Solar Cycle 23 through the declining phase of Cycle 24. The subsurface flows are mainly divergent at quiet regions and convergent at locations of high magnetic activity. The relationship is essentially linear between divergence and magnetic activity at all activity levels at depths shallower than about 10 Mm. At greater depths, the relationship changes sign at locations of high activity; the flows are increasingly divergent at locations with a magnetic activity index (MAI) greater than about 24 G. The flows are more convergent by about a factor of two during the rising phase of Cycle 24 than during the declining phase of Cycle 23 at locations of medium and high activity (about 10 to 40 G MAI) from the surface to at least 10 Mm. The subsurface divergence pattern of Solar Cycle 24 first appears during the declining phase of Cycle 23 and is present during the extended minimum. It appears several years before the magnetic pattern of the new cycle is noticeable in synoptic maps. Using linear regression, we estimate the amount of magnetic activity that would be required to generate the precursor pattern and find that it should be almost twice the amount of activity that is observed.

Effect of multilayer substrate configuration in horizontal subsurface flow constructed wetlands: assessment of treatment performance, biofilm development, and solids accumulation.

PubMed

Ding, Yanli; Lyu, Tao; Bai, Shaoyuan; Li, Zhenling; Ding, Haijing; You, Shaohong; Xie, Qinglin

2018-01-01

This study investigates the influence of multilayer substrate configuration in horizontal subsurface flow constructed wetlands (HSCWs) on their treatment performance, biofilm development, and solids accumulation. Three pilot-scale HSCWs were built to treat campus sewage and have been operational for 3 years. The HSCWs included monolayer (CW1), three-layer (CW3), and six-layer (CW6) substrate configurations with hydraulic conductivity of the substrate increasing from the surface to bottom in the multilayer CWs. It was demonstrated the pollutant removal performance after a 3-year operation improved in the multilayer HSCWs (49-80%) compared to the monolayer HSCW (29-41%). Simultaneously, the multilayer HSCWs exhibited significant features that prevented clogging compared to the monolayer configuration. The amount of accumulated solids was notably higher in the monolayer CW compared to multilayer CWs. Further, multilayer HSCWs could delay clogging by providing higher biofilm development for organics removal and consequently, lesser solids accumulations. Principal component analysis strongly supported the visualization of the performance patterns in the present study and showed that multilayer substrate configuration, season, and sampling locations significantly influenced biofilm growth and solids accumulation. Finally, the present study provided important information to support the improved multilayer configured HSCW implication in the future.

Effects of step-feeding and intermittent aeration on organics and nitrogen removal in a horizontal subsurface flow constructed wetland.

PubMed

Patil, Sagar; Chakraborty, Saswati

2017-03-21

The effect of step feed strategy and intermittent aeration on removal of chemical oxygen demand (COD) and nitrogen was investigated in a laboratory scale horizontal subsurface flow constructed wetland (HSSFCW). Wetland was divided into four zones along the length (zone I to IV), and influent was introduced into first and third zones by step feeding. Continuous study was carried out in four phases. In phases I to III, 30% of influent was bypassed to zone III for denitrification along with organics removal. Intermittent aeration was provided only in zone II at 2.5 L/min for 4 h/day, during phases II, III and IV. In phase I, 87% COD and 43% NH 4 + -N (ammonia-nitrogen) removal were obtained from influents of 331 and 30 mg/L, respectively. In phase II study, external aeration resulted in 97% COD and 71% NH 4 + -N removal in the wetland. In phase IV, 40% of feed was delivered to zone III. Higher supply of organic in zone III resulted in higher denitrification, and total nitrogen removal rate increased to 70% from 56%. In the final effluent, concentration of NO 3 - -N was 9-11 mg/L in phase I to III and decreased to 4 mg/L in phase IV. Batch study showed that COD and NH 4 + -N removal followed first order kinetics in different zones of wetland.

Subsurface thermal regime to delineate the paleo-groundwater flow system in an arid area, Al Kufra, Libya

NASA Astrophysics Data System (ADS)

Salem, Zenhom El-Said

2016-12-01

The purpose of this study was to understand the groundwater flow system in Al Kufra basin, Libya, as a case study of arid areas using subsurface temperature. The temperature-depth profiles and water levels were measured in eight boreholes in the area. Well 6 is considered a recharge type profile with low geothermal gradient (0.0068 °C/m) and an estimated paleo-temperature around 19.5 °C. The other profiles are of discharge type with higher geothermal gradient (0.0133 to 0.0166 °C/m). The constructed horizontal 2D distribution maps of the hydraulic heads and the subsurface temperature measurements reveal that the main recharge area is located to the south with low temperature while the main discharge area is located to the north with higher temperature. Vertical 2D distribution maps show that location of well 4 has low hydraulic heads and higher temperature indicating that the fault defined in the area may have affected the groundwater flow system. The estimated groundwater flux ranges from 0.001 to 0.1 mm/day for the recharge area and from -0.3 to -0.7 mm/day in average in the discharge area.

The Influence of Plant Root Systems on Subsurface Flow: Implications for Slope Stability

EPA Science Inventory

Although research has explained how plant roots mechanically stabilize soils, in this article we explore how root systems create networks of preferential flow and thus influence water pressures in soils to trigger landslides. Root systems may alter subsurface flow: Hydrological m...

High performance constructed wetlands for cold climates.

PubMed

Jenssen, Petter D; Maehlum, Trend; Krogstad, Tore; Vråle, Lasse

2005-01-01

In 1991, the first subsurface flow constructed wetland for treatment of domestic wastewater was built in Norway. Today, this method is rapidly becoming a popular method for wastewater treatment in rural Norway. This is due to excellent performance even during winter and low maintenance. The systems can be constructed regardless of site conditions. The Norwegian concept for small constructed wetlands is based on the use of a septic tank followed by an aerobic vertical down-flow biofilter succeeded by a subsurface horizontal-flow constructed wetland. The aerobic biofilter, prior to the subsurface flow stage, is essential to remove BOD and achieve nitrification in a climate where the plants are dormant during the cold season. When designed according to present guidelines a consistent P-removal of > 90% can be expected for 15 years using natural iron or calcium rich sand or a new manufactured lightweight aggregate with P-sorption capacities, which exceeds most natural media. When the media is saturated with P it can be used as soil conditioner and P-fertilizer. Nitrogen removal in the range of 40-60% is achieved. Removal of indicator bacteria is high and < 1000 thermotolerant coliforms/100 ml is normally achieved.

Constructed wetland attenuation of nitrogen exported in subsurface drainage from irrigated and rain-fed dairy pastures.

PubMed

Tanner, C C; Nguyen, M L; Sukias, J P S

2005-01-01

Nitrogen removal performance is reported for constructed wetlands treating subsurface drainage from irrigated and rain-fed dairy pastures in North Island, New Zealand. Flow-proportional sampling of inflow and outflow concentrations were combined with continuous flow records to calculate mass balances for the wetlands. Drainage flows from the irrigated catchment were 2.5-4 fold higher and N exports up to 5 fold higher per unit area than for the rain-fed catchment. Hydraulic and associated N loadings to the wetlands were highly pulsed, associated with rainfall, soil water status, and irrigation events. Transient pulses of organic nitrogen were an important form of N loss from the rain-fed landscape in the first year, and were very effectively removed in the wetland (> 90%). Median nitrate concentrations of approximately 10 g m(-3) in the drainage inflows were reduced by 15-67% during passage through the wetlands and annual nitrate-N loads by 16-61% (38-31 7 g N m(-2)y(-1)). Generation in the wetlands of net ammoniacal-N and organic-N (irrigated site) partially negated reduction in nitrate-N loads. The results show that constructed wetlands comprising 1-2% of catchment area can provide moderate reductions in TN export via pastoral drainage, but performance is markedly influenced by variations in seasonal loading and establishment/maturation factors.

Nitrogen removal from wastewater by an aerated subsurface-flow constructed wetland in cold climates.

PubMed

Redmond, Eric D; Just, Craig L; Parkin, Gene F

2014-04-01

The objective of this study was to assess the role of cyclic aeration, vegetation, and temperature on nitrogen removal by subsurface-flow engineered wetlands. Aeration was shown to enhance total nitrogen and ammonia removal and to enhance removal of carbonaceous biochemical oxygen demand, chemical oxygen demand, and phosphorus. Effluent ammonia and total nitrogen concentrations were significantly lower in aerated wetland cells when compared with unaerated cells. There was no significant difference in nitrogen removal between planted and unplanted cells. Effluent total nitrogen concentrations ranged from 9 to 12 mg N/L in the aerated cells and from 23 to 24 mg N/L in unaerated cells. Effluent ammonia concentrations ranged from 3 to 7 mg N/L in aerated wetland cells and from 22 to 23 mg N/L in unaerated cells. For the conditions tested, temperature had only a minimal effect on effluent ammonia or total nitrogen concentrations. The tanks-in-series and the PkC models predicted the general trends in effluent ammonia and total nitrogen concentrations, but did not do well predicting short-term variability. Rate coefficients for aerated systems were 2 to 10 times greater than those for unaerated systems.

Understanding the rapidity of subsurface storm flow response from a fracture-oriented shallow vadose through a new perspective

NASA Astrophysics Data System (ADS)

Zhao, Peng; Zhao, Pei; Liang, Chuan; Li, Tianyang; Zhou, Baojia

2017-01-01

Velocity and celerity in hydrologic systems are controlled by different mechanisms. Efforts were made through joint sample collection and the use of hydrographs and tracers to understand the rapidity of the subsurface flow response to rainstorms on hourly time scales. Three deep subsurface flows during four natural rainstorm events were monitored. The results show that (1) deeper discharge was observed early in responding rainfall events and yielded a high hydrograph amplitude; (2) a ratio index, k, reflecting the dynamic change of the rainfall perturbation intensity in subsurface flow, might reveal inner causal relationships between the flow index and the tracer signal index. Most values of k were larger than 1 at the perturbation stage but approximated 1 at the no-perturbation stage; and (3) for statistical analysis of tracer signals in subsurface flows, the total standard deviation was 17.2, 11.9, 7.4 and 3.5 at perturbation stages and 4.4, 2.5, 1.1, and 0.95 at the non-perturbation stage for observed events. These events were 3-7 times higher in the former rather than the later, reflecting that the variation of tracer signals primarily occurred under rainfall perturbation. Thus, we affirmed that the dynamic features of rainfall have a key effect on rapid processes because, besides the gravity, mechanical waves originating from dynamic rainfall features are another driving factor for conversion between different types of rainfall mechanical energy. A conceptual model for pressure wave propagation was proposed, in which virtual subsurface flow processes in a heterogeneous vadose zone under rainfall are analogous to the water hammer phenomenon in complex conduit systems. Such an analogy can allow pressure in a shallow vadose to increase and decrease and directly influence the velocity and celerity of the flow reflecting a mechanism for rapid subsurface hydrologic response processes in the shallow vadose zone.

Enhancement of the complete autotrophic nitrogen removal over nitrite process in a modified single-stage subsurface vertical flow constructed wetland: Effect of saturated zone depth.

PubMed

Huang, Menglu; Wang, Zhen; Qi, Ran

2017-06-01

This study was conducted to explore enhancement of the complete autotrophic nitrogen removal over nitrite (CANON) process in a modified single-stage subsurface vertical flow constructed wetland (VSSF) with saturated zone, and nitrogen transformation pathways in the VSSF treating digested swine wastewater were investigated at four different saturated zone depths (SZDs). SZD significantly affected nitrogen transformation pathways in the VSSF throughout the experiment. As the SZD was 45cm, the CANON process was enhanced most effectively in the system owing to the notable enhancement of anammox. Correspondingly, the VSSF had the best TN removal performance [(76.74±7.30)%] and lower N 2 O emission flux [(3.50±0.22)mg·(m 2 ·h) - 1 ]. It could be concluded that autotrophic nitrogen removal via CANON process could become a primary route for nitrogen removal in the VSSF with optimized microenvironment that developed as a result of the appropriate SZD. Copyright © 2017 Elsevier Ltd. All rights reserved.

MODELING THREE-DIMENSIONAL SUBSURFACE FLOW, FATE AND TRANSPORT OF MICROBES AND CHEMICALS (3DFATMIC)

EPA Science Inventory

A three-dimensional model simulating the subsurface flow, microbial growth and degradation, microbial-chemical reaction, and transport of microbes and chemicals has been developed. he model is designed to solve the coupled flow and transport equations. asically, the saturated-uns...

Evaluation of SCS-CN method using a fully distributed physically based coupled surface-subsurface flow model

NASA Astrophysics Data System (ADS)

Shokri, Ali

2017-04-01

The hydrological cycle contains a wide range of linked surface and subsurface flow processes. In spite of natural connections between surface water and groundwater, historically, these processes have been studied separately. The current trend in hydrological distributed physically based model development is to combine distributed surface water models with distributed subsurface flow models. This combination results in a better estimation of the temporal and spatial variability of the interaction between surface and subsurface flow. On the other hand, simple lumped models such as the Soil Conservation Service Curve Number (SCS-CN) are still quite common because of their simplicity. In spite of the popularity of the SCS-CN method, there have always been concerns about the ambiguity of the SCS-CN method in explaining physical mechanism of rainfall-runoff processes. The aim of this study is to minimize these ambiguity by establishing a method to find an equivalence of the SCS-CN solution to the DrainFlow model, which is a fully distributed physically based coupled surface-subsurface flow model. In this paper, two hypothetical v-catchment tests are designed and the direct runoff from a storm event are calculated by both SCS-CN and DrainFlow models. To find a comparable solution to runoff prediction through the SCS-CN and DrainFlow, the variance between runoff predictions by the two models are minimized by changing Curve Number (CN) and initial abstraction (Ia) values. Results of this study have led to a set of lumped model parameters (CN and Ia) for each catchment that is comparable to a set of physically based parameters including hydraulic conductivity, Manning roughness coefficient, ground surface slope, and specific storage. Considering the lack of physical interpretation in CN and Ia is often argued as a weakness of SCS-CN method, the novel method in this paper gives a physical explanation to CN and Ia.

Storm event-scale nutrient attenuation in constructed wetlands experiencing a Mediterranean climate: A comparison of a surface flow and hybrid surface-subsurface flow system.

PubMed

Adyel, Tanveer M; Oldham, Carolyn E; Hipsey, Matthew R

2017-11-15

Among different Water Sensitive Urban Design options, constructed wetlands (CWs) are used to protect and restore downstream water quality by attenuating nutrients generated by stormwater runoff. This research compared the nutrient attenuation ability during a diverse population of storm events of two CWs: (a) a hybrid CW with multiple alternating surface flow (SF) and laterite-based subsurface flow (SSF) compartments, and (b) a single stage SF CW. Within-storm variability, nutrient concentrations were assessed at 2 to 3-h intervals at both the main inlet and outlet of each CW. Dissolved oxygen concentrations of the surface waters were also monitored at 10-min intervals using high frequency in situ sensors. Nutrient loads into the CWs were observed to be higher when a high rainfall event occurred, particularly after longer antecedent dry conditions. Longer hydraulic retention times promoted higher attenuation at both sites. However, the relative extent of nutrient attenuation differed between the CW types; the mean total nitrogen (TN) attenuation in the hybrid and SF CW was 45 and 48%, respectively. The hybrid CW attenuated 67% total phosphorus (TP) loads on average, while the SF CW acted as a net TP source. Periodic storm events transitioned the lentic CW into a lotic CW and caused riparian zone saturation; it was therefore hypothesized that such saturation of organic matter rich-riparian zones led to release of TP in the system. The hybrid CW attenuated the released TP in the downstream laterite-based SSF compartments. Diel oxygen metabolism calculated before and after the storm events was found to be strongly correlated with water temperature, solar exposure and antecedent dry condition during the pre-storm conditions. Furthermore, the SF CW showed a significant relationship between overall nutrient load attenuation and the change in oxygen metabolism during the storm perturbation, suggesting oxygen variation could be a useful proxy indicator of CW function

Lateral and subsurface flows impact arctic coastal plain lake water budgets

USGS Publications Warehouse

Koch, Joshua C.

2016-01-01

Arctic thaw lakes are an important source of water for aquatic ecosystems, wildlife, and humans. Many recent studies have observed changes in Arctic surface waters related to climate warming and permafrost thaw; however, explaining the trends and predicting future responses to warming is difficult without a stronger fundamental understanding of Arctic lake water budgets. By measuring and simulating surface and subsurface hydrologic fluxes, this work quantified the water budgets of three lakes with varying levels of seasonal drainage, and tested the hypothesis that lateral and subsurface flows are a major component of the post-snowmelt water budgets. A water budget focused only on post-snowmelt surface water fluxes (stream discharge, precipitation, and evaporation) could not close the budget for two of three lakes, even when uncertainty in input parameters was rigorously considered using a Monte Carlo approach. The water budgets indicated large, positive residuals, consistent with up to 70% of mid-summer inflows entering lakes from lateral fluxes. Lateral inflows and outflows were simulated based on three processes; supra-permafrost subsurface inflows from basin-edge polygonal ground, and exchange between seasonally drained lakes and their drained margins through runoff and evapotranspiration. Measurements and simulations indicate that rapid subsurface flow through highly conductive flowpaths in the polygonal ground can explain the majority of the inflow. Drained lakes were hydrologically connected to marshy areas on the lake margins, receiving water from runoff following precipitation and losing up to 38% of lake efflux to drained margin evapotranspiration. Lateral fluxes can be a major part of Arctic thaw lake water budgets and a major control on summertime lake water levels. Incorporating these dynamics into models will improve our ability to predict lake volume changes, solute fluxes, and habitat availability in the changing Arctic.

Advances in Multiphase Flow and Transport in the Subsurface Environment

DOE Office of Scientific and Technical Information (OSTI.GOV)

Shi, Xiaoqing; Finsterle, Stefan; Zhang, Keni

Multiphase flow and transport processes in the subsurface environment are extremely important in a number of industrial and environmental applications at various spatial and temporal scales. Thus, it is necessary to identify, understand, and predict these processes to improve the production of conventional and unconventional oil and gas, to increase the safety of geological sequestration of carbon dioxide and nuclear waste disposal, and to make remediation of contaminated aquifers more effective.

Advances in Multiphase Flow and Transport in the Subsurface Environment

DOE PAGES

Shi, Xiaoqing; Finsterle, Stefan; Zhang, Keni; ...

2018-03-04

Multiphase flow and transport processes in the subsurface environment are extremely important in a number of industrial and environmental applications at various spatial and temporal scales. Thus, it is necessary to identify, understand, and predict these processes to improve the production of conventional and unconventional oil and gas, to increase the safety of geological sequestration of carbon dioxide and nuclear waste disposal, and to make remediation of contaminated aquifers more effective.

Three-Dimensional Subsurface Flow, Fate and Transport of Microbes and Chemicals (3DFATMIC) Model

EPA Pesticide Factsheets

This model simulates subsurface flow, fate and transport of contaminants that are undergoing chemical or biological transformations. The model is applicable to transient conditions in both saturated and unsaturated zones.

Two-Dimensional Subsurface Flow, Fate and Transport of Microbes and Chemicals (2DFATMIC) Model

EPA Pesticide Factsheets

This model simulates subsurface flow, fate, and transport of contaminants that are undergoing chemical or biological transformations. This model is applicable to transient conditions in both saturated and unsaturated zones.

Removal of antibiotics and antibiotic resistance genes from domestic sewage by constructed wetlands: Effect of flow configuration and plant species.

PubMed

Chen, Jun; Ying, Guang-Guo; Wei, Xiao-Dong; Liu, You-Sheng; Liu, Shuang-Shuang; Hu, Li-Xin; He, Liang-Ying; Chen, Zhi-Feng; Chen, Fan-Rong; Yang, Yong-Qiang

2016-11-15

This study aims to investigate the removal of antibiotics and antibiotic resistance genes (ARGs) in raw domestic wastewater by various mesocosm-scale constructed wetlands (CWs) with different flow configurations or plant species including the constructed wetland with or without plant. Six mesocosm-scale CWs with three flow types (surface flow, horizontal subsurface flow and vertical subsurface flow) and two plant species (Thaliadealbata Fraser and Iris tectorum Maxim) were set up in the outdoor. 8 antibiotics including erythromycin-H2O (ETM-H2O), monensin (MON), clarithromycin (CTM), leucomycin (LCM), sulfamethoxazole (SMX), trimethoprim (TMP), sulfamethazine (SMZ) and sulfapyridine (SPD) and 12 genes including three sulfonamide resistance genes (sul1, sul2 and sul3), four tetracycline resistance genes (tetG, tetM, tetO and tetX), two macrolide resistance genes (ermB and ermC), two chloramphenicol resistance genes (cmlA and floR) and 16S rRNA (bacteria) were determined in different matrices (water, particle, substrate and plant phases) from the mesocosm-scale systems. The aqueous removal efficiencies of total antibiotics ranged from 75.8 to 98.6%, while those of total ARGs varied between 63.9 and 84.0% by the mesocosm-scale CWs. The presence of plants was beneficial to the removal of pollutants, and the subsurface flow CWs had higher pollutant removal than the surface flow CWs, especially for antibiotics. According to the mass balance analysis, the masses of all detected antibiotics during the operation period were 247,000, 4920-10,600, 0.05-0.41 and 3500-60,000μg in influent, substrate, plant and effluent of the mesocosm-scale CWs. In the CWs, biodegradation, substrate adsorption and plant uptake all played certain roles in reducing the loadings of nutrients, antibiotics and ARGs, but biodegradation was the most important process in the removal of these pollutants. Copyright © 2016 Elsevier B.V. All rights reserved.

Subsurface Flow and Moisture Dynamics in Response to Swash Motions: Effects of Beach Hydraulic Conductivity and Capillarity

NASA Astrophysics Data System (ADS)

Geng, Xiaolong; Heiss, James W.; Michael, Holly A.; Boufadel, Michel C.

2017-12-01

A combined field and numerical study was conducted to investigate dynamics of subsurface flow and moisture response to waves in the swash zone of a sandy beach located on Cape Henlopen, DE. A density-dependent variably saturated flow model MARUN was used to simulate subsurface flow beneath the swash zone. Values of hydraulic conductivity (K) and characteristic pore size (α, a capillary fringe property) were varied to evaluate their effects on subsurface flow and moisture dynamics in response to swash motions in beach aquifers. The site-specific modeling results were validated against spatiotemporal measurements of moisture and pore pressure in the beach. Sensitivity analyses indicated that the hydraulic conductivity and capillary fringe thickness of the beach greatly influenced groundwater flow pathways and associated transit times in the swash zone. A higher value of K enhanced swash-induced seawater infiltration into the beach, thereby resulting in a faster expansion of a wedge of high moisture content induced by swash cycles, and a flatter water table mound beneath the swash zone. In contrast, a thicker capillary fringe retained higher moisture content near the beach surface, and thus, significantly reduced the available pore space for infiltration of seawater. This attenuated wave effects on pore water flow in the unsaturated zone of the beach. Also, a thicker capillary fringe enhanced horizontal flow driven by the larger-scale hydraulic gradient caused by tides.

Evidence for debris flow gully formation initiated by shallow subsurface water on Mars

USGS Publications Warehouse

Lanza, N.L.; Meyer, G.A.; Okubo, C.H.; Newsom, Horton E.; Wiens, R.C.

2010-01-01

The morphologies of some martian gullies appear similar to terrestrial features associated with debris flow initiation, erosion, and deposition. On Earth, debris flows are often triggered by shallow subsurface throughflow of liquid water in slope-mantling colluvium. This flow causes increased levels of pore pressure and thus decreased shear strength, which can lead to slide failure of slope materials and subsequent debris flow. The threshold for pore pressure-induced failure creates a distinct relationship between the contributing area supplying the subsurface flow and the slope gradient. To provide initial tests of a similar debris flow initiation hypothesis for martian gullies, measurements of the contributing areas and slope gradients were made at the channel heads of martian gullies seen in three HiRISE stereo pairs. These gullies exhibit morphologies suggestive of debris flows such as leveed channels and lobate debris fans, and have well-defined channel heads and limited evidence for multiple flows. Our results show an area-slope relationship for these martian gullies that is consistent with that observed for terrestrial gullies formed by debris flow, supporting the hypothesis that these gullies formed as the result of saturation of near-surface regolith by a liquid. This model favors a source of liquid that is broadly distributed within the source area and shallow; we suggest that such liquid could be generated by melting of broadly distributed icy materials such as snow or permafrost. This interpretation is strengthened by observations of polygonal and mantled terrain in the study areas, which are both suggestive of near-surface ice. ?? 2009 Elsevier Inc.

Influence of UV radiation on chlorophyll, and antioxidant enzymes of wetland plants in different types of constructed wetland.

PubMed

Xu, Defu; Wu, Yinjuan; Li, Yingxue; Howard, Alan; Jiang, Xiaodong; Guan, Yidong; Gao, Yongxia

2014-09-01

A surface- and vertical subsurface-flow-constructed wetland were designed to study the response of chlorophyll and antioxidant enzymes to elevated UV radiation in three types of wetland plants (Canna indica, Phragmites austrail, and Typha augustifolia). Results showed that (1) chlorophyll content of C. indica, P. austrail, and T. augustifolia in the constructed wetland was significantly lower where UV radiation was increased by 10 and 20 % above ambient solar level than in treatment with ambient solar UV radiation (p < 0.05). (2) The malondialdehyde (MDA) content, guaiacol peroxidase (POD), superoxide dismutase (SOD), and catalase (CAT) activities of wetland plants increased with elevated UV radiation intensity. (3) The increased rate of MDA, SOD, POD, and CAT activities of C. indica, P. australis, and T. angustifolia by elevated UV radiation of 10 % was higher in vertical subsurface-flow-constructed wetland than in surface-flow-constructed wetland. The sensitivity of MDA, SOD, POD, and CAT activities of C. indica, P. austrail, and T. augustifolia to the elevated UV radiation was lower in surface-flow-constructed wetland than in the vertical subsurface-flow-constructed wetland, which was related to a reduction in UV radiation intensity through the dissolved organic carbon and suspended matter in the water. C. indica had the highest SOD and POD activities, which implied it is more sensitive to enhanced UV radiation. Therefore, different wetland plants had different antioxidant enzymes by elevated UV radiation, which were more sensitive in vertical subsurface-flow-constructed wetland than in surface-flow-constructed wetland.

Performance comparison and economics analysis of waste stabilization ponds and horizontal subsurface flow constructed wetlands treating domestic wastewater: a case study of the Juja sewage treatment works.

PubMed

Mburu, Njenga; Tebitendwa, Sylvie M; van Bruggen, Johan J A; Rousseau, Diederik P L; Lens, Piet N L

2013-10-15

The performance, effluent quality, land area requirement, investment and operation costs of a full-scale waste stabilization pond (WSP) and a pilot scale horizontal subsurface flow constructed wetland (HSSF-CW) at Jomo Kenyatta University of Agriculture and Technology (JKUAT) were investigated between November 2010 to January 2011. Both systems gave comparable medium to high levels of organic matter and suspended solids removal. However, the WSP showed a better removal for Total Phosphorus (TP) and Ammonium (NH4(+)-N). Based on the population equivalent calculations, the land area requirement per person equivalent of the WSP system was 3 times the area that would be required for the HSSF-CW to treat the same amount of wastewater. The total annual cost estimates consisting of capital, operation and maintenance (O&M) costs were comparable for both systems. However, the evaluation of the capital cost of either system showed that it is largely influenced by the size of the population served, local cost of land and the construction materials involved. Hence, one can select either system in terms of treatment efficiency. When land is available other factor including the volume of wastewater or the investment, and O&M costs determine the technology selection. Copyright © 2013 Elsevier Ltd. All rights reserved.

Subsurface drainage processes and management impacts

Treesearch

Elizabeth T. Keppeler; David Brown

1998-01-01

Storm-induced streamflow in forested upland watersheds is linked to rainfall by transient, variably saturated flow through several different flow paths. In the absence of exposed bedrock, shallow flow-restrictive layers, or compacted soil surfaces, virtually all of the infiltrated rainfall reaches the stream as subsurface flow. Subsurface runoff can occur within...

Subsurface flow in lowland river gravel bars

NASA Astrophysics Data System (ADS)

Bray, E. N.; Dunne, T.

2017-09-01

Geomorphic and hydraulic processes, which form gravel bars in large lowland rivers, have distinctive characteristics that control the magnitude and spatial patterns of infiltration and exfiltration between rivers and their immediate subsurface environments. We present a bedform-infiltration relation together with a set of field measurements along two reaches of the San Joaquin River, CA to illustrate the conditions required for infiltration and exfiltration of flow between a stream and its undulating bed, and a numerical model to investigate the factors that affect paths and residence times of flow through barforms at different discharges. It is shown that asymmetry of bar morphology is a first-order control on the extent and location of infiltration, which would otherwise produce equal areas of infiltration and exfiltration under the assumption of sinusoidal bedforms. Hydraulic conductivity varies by orders of magnitude due to fine sediment accumulation and downstream coarsening related to the process of bar evolution. This systematic variability not only controls the magnitude of infiltration, but also the residence time of flow through the bed. The lowest hydraulic conductivity along the reach occurred where the difference between the topographic gradient and the water-surface gradient is at a maximum and thus where infiltration would be greatest into a homogeneous bar, indicating the importance of managing sand supply to maintain the ventilation and flow through salmon spawning riffles. Numerical simulations corroborate our interpretation that infiltration patterns and rates are controlled by distinctive features of bar morphology.

A multiscale fixed stress split iterative scheme for coupled flow and poromechanics in deep subsurface reservoirs

NASA Astrophysics Data System (ADS)

Dana, Saumik; Ganis, Benjamin; Wheeler, Mary F.

2018-01-01

In coupled flow and poromechanics phenomena representing hydrocarbon production or CO2 sequestration in deep subsurface reservoirs, the spatial domain in which fluid flow occurs is usually much smaller than the spatial domain over which significant deformation occurs. The typical approach is to either impose an overburden pressure directly on the reservoir thus treating it as a coupled problem domain or to model flow on a huge domain with zero permeability cells to mimic the no flow boundary condition on the interface of the reservoir and the surrounding rock. The former approach precludes a study of land subsidence or uplift and further does not mimic the true effect of the overburden on stress sensitive reservoirs whereas the latter approach has huge computational costs. In order to address these challenges, we augment the fixed-stress split iterative scheme with upscaling and downscaling operators to enable modeling flow and mechanics on overlapping nonmatching hexahedral grids. Flow is solved on a finer mesh using a multipoint flux mixed finite element method and mechanics is solved on a coarse mesh using a conforming Galerkin method. The multiscale operators are constructed using a procedure that involves singular value decompositions, a surface intersections algorithm and Delaunay triangulations. We numerically demonstrate the convergence of the augmented scheme using the classical Mandel's problem solution.

Composting plant leachate treatment by a pilot-scale, three-stage, horizontal flow constructed wetland in central Iran.

PubMed

Bakhshoodeh, Reza; Alavi, Nadali; Paydary, Pooya

2017-10-01

Handling and treatment of composting leachate is difficult and poses major burdens on composting facilities. The main goal of this study was to evaluate usage of a three-stage, constructed wetland to treat leachate produced in Isfahan composting facility. A pilot-scale, three-stage, subsurface, horizontal flow constructed wetland, planted with vetiver with a flow rate of 24 L/day and a 15-day hydraulic retention time, was used. Removal of organic matter, ammonia, nitrate, total nitrogen, suspended solids, and several heavy metals from Isfahan composting facility leachate was monitored over a 3-month period. Constructed wetland system was capable of efficiently removing BOD 5 (87.3%), COD (74.5%), ammonia (91.5%), nitrate (87.9%), total nitrogen (87.8%), total suspended solids (85.5%), and heavy metals (ranging from 70 to 90%) from the composting leachate. High contaminant removal efficiencies were achieved, but effluent still failed to meet Iranian standards for treated wastewater. This study shows that although a three-stage horizontal flow constructed wetland planted with vetiver cannot be used alone to treat Isfahan composting facility leachate, but it has the potential to be used as a leachate pre-treatment step, along with another complementary method.

Modelling rapid subsurface flow at the hillslope scale with explicit representation of preferential flow paths

NASA Astrophysics Data System (ADS)

Wienhöfer, J.; Zehe, E.

2012-04-01

Rapid lateral flow processes via preferential flow paths are widely accepted to play a key role for rainfall-runoff response in temperate humid headwater catchments. A quantitative description of these processes, however, is still a major challenge in hydrological research, not least because detailed information about the architecture of subsurface flow paths are often impossible to obtain at a natural site without disturbing the system. Our study combines physically based modelling and field observations with the objective to better understand how flow network configurations influence the hydrological response of hillslopes. The system under investigation is a forested hillslope with a small perennial spring at the study area Heumöser, a headwater catchment of the Dornbirnerach in Vorarlberg, Austria. In-situ points measurements of field-saturated hydraulic conductivity and dye staining experiments at the plot scale revealed that shrinkage cracks and biogenic macropores function as preferential flow paths in the fine-textured soils of the study area, and these preferential flow structures were active in fast subsurface transport of artificial tracers at the hillslope scale. For modelling of water and solute transport, we followed the approach of implementing preferential flow paths as spatially explicit structures of high hydraulic conductivity and low retention within the 2D process-based model CATFLOW. Many potential configurations of the flow path network were generated as realisations of a stochastic process informed by macropore characteristics derived from the plot scale observations. Together with different realisations of soil hydraulic parameters, this approach results in a Monte Carlo study. The model setups were used for short-term simulation of a sprinkling and tracer experiment, and the results were evaluated against measured discharges and tracer breakthrough curves. Although both criteria were taken for model evaluation, still several model setups

Vulnerability assessment of groundwater-dependent ecosystems based on integrated groundwater flow modell construction

NASA Astrophysics Data System (ADS)

Tóth, Ádám; Simon, Szilvia; Galsa, Attila; Havril, Timea; Monteiro Santos, Fernando A.; Müller, Imre; Mádl-Szőnyi, Judit

2017-04-01

Groundwater-dependent ecosystems (GDEs) are highly influenced by the amount of groundwater, seasonal variation of precipitation and consequent water table fluctuation and also the anthropogenic activities. They can be regarded as natural surface manifestations of the flowing groundwater. The preservation of environment and biodiversity of these GDEs is an important issue worldwide, however, the water management policy and action plan could not be constructed in absense of proper hydrogeological knowledge. The concept of gravity-driven regional groundwater flow could aid the understanding of flow pattern and interpretation of environmental processes and conditions. Unless the required well data are available, the geological-hydrogeological numerical model of the study area cannot be constructed based only on borehole information. In this case, spatially continuous geophysical data can support groundwater flow model building: systematically combined geophysical methods can provide model input. Integration of lithostratigraphic, electrostratigraphic and hydrostratigraphic information could aid groundwater flow model construction: hydrostratigraphic units and their hydraulic behaviour, boundaries and geometry can be obtained. Groundwater-related natural manifestations, such as GDEs, can be explained with the help of the revealed flow pattern and field mapping of features. Integrated groundwater flow model construction for assessing the vulnerability of GDEs was presented via the case study of the geologically complex area of Tihany Peninsula, Hungary, with the aims of understanding the background and occurrence of groundwater-related environmental phenomena, surface water-groundwater interaction, and revealing the potential effect of anthropogenic activity and climate change. In spite of its important and protected status, fluid flow model of the area, which could support water management and natural protection policy, had not been constructed previously. The 3D

Effects of Surface and Subsurface Bed Material Composition on Gravel Transport and Flow Competence Relations—Possibilities for Prediction

NASA Astrophysics Data System (ADS)

Bunte, K.; Abt, S. R.; Swingle, K. W.; Cenderelli, D. A.; Gaeuman, D. A.

2014-12-01

Bedload transport and flow competence relations are difficult to predict in coarse-bedded steep streams where widely differing sediment supply, bed stability, and complex flow hydraulics greatly affect amounts and sizes of transported gravel particles. This study explains how properties of bed material surface and subsurface size distributions are directly related to gravel transport and may be used for prediction of gravel transport and flow competence relations. Gravel transport, flow competence, and bed material size were measured in step-pool and plane-bed streams. Power functions were fitted to gravel transport QB=aQb and flow competence Dmax=cQd relations; Q is water discharge. Frequency distributions of surface FDsurf and subsurface FDsub bed material were likewise described by power functions FDsurf=hD j and FDsub=kDm fitted over six 0.5-phi size classes within 4 to 22.4 mm. Those gravel sizes are typically mobile even in moderate floods. Study results show that steeper subsurface bed material size distributions lead to steeper gravel transport and flow competence relations, whereas larger amounts of sediment contained in those 6 size bedmaterial classes (larger h and k) flatten the relations. Similarly, steeper surface size distributions decrease the coefficients of the gravel transport and flow competence relations, whereas larger amounts of sediment within the six bed material classes increase the intercepts of gravel transport and flow competence relations. Those relations are likely causative in streams where bedload stems almost entirely from the channel bed as opposed to direct (unworked) contributions from hillslopes and tributaries. The exponent of the subsurface bed material distribution m predicted the gravel transport exponent b with r2 near 0.7 and flow competence exponent d with r2 near 0.5. The intercept of bed surface distributions h increased the intercept a of gravel transport and c of the flow competence relations with r2 near 0.6.

A locally conservative stabilized continuous Galerkin finite element method for two-phase flow in poroelastic subsurfaces

NASA Astrophysics Data System (ADS)

Deng, Q.; Ginting, V.; McCaskill, B.; Torsu, P.

2017-10-01

We study the application of a stabilized continuous Galerkin finite element method (CGFEM) in the simulation of multiphase flow in poroelastic subsurfaces. The system involves a nonlinear coupling between the fluid pressure, subsurface's deformation, and the fluid phase saturation, and as such, we represent this coupling through an iterative procedure. Spatial discretization of the poroelastic system employs the standard linear finite element in combination with a numerical diffusion term to maintain stability of the algebraic system. Furthermore, direct calculation of the normal velocities from pressure and deformation does not entail a locally conservative field. To alleviate this drawback, we propose an element based post-processing technique through which local conservation can be established. The performance of the method is validated through several examples illustrating the convergence of the method, the effectivity of the stabilization term, and the ability to achieve locally conservative normal velocities. Finally, the efficacy of the method is demonstrated through simulations of realistic multiphase flow in poroelastic subsurfaces.

Finding the best windows: An apparent environmental threshold determines which diffuse flows are dominated by subsurface microbes

NASA Astrophysics Data System (ADS)

Olins, H. C.; Rogers, D.; Scholin, C. A.; Preston, C. J.; Vidoudez, C.; Ussler, W.; Pargett, D.; Jensen, S.; Roman, B.; Birch, J. M.; Girguis, P. R.

2014-12-01

Hydrothermal vents are hotspots of microbial primary productivity often described as "windows into the subsurface biosphere." High temperature vents have received the majority of research attention, but cooler diffuse flows are as, if not more, important a source of heat and chemicals to the overlying ocean. We studied patterns of in situ gene expression and co-registered geochemistry in order to 1) describe the diversity and physiological poise of active microbial communities that span thermal and geochemical gradients from active diffuse flow to background vent field seawater, and 2) determine to what extent seawater or subsurface microbes were active throughout this environment. Analyses of multiple metatranscriptomes from 5 geochemically distinct sites (some from samples preserved in situ) show that proximate diffuse flows showed strikingly different transcription profiles. Specifically, caldera background and some diffuse flows were similar, both dominated by seawater-derived Gammaproteobacteria despite having distinct geochemistries. Intra-field community shows evidence of increased primary productivity throughout the entire vent field and not just at individual diffuse flows. In contrast, a more spatially limited, Epsilonproteobacteria-dominated transcription profile from the most hydrothermally-influenced diffuse flow appeared to be driven by the activity of vent-endemic microbes, likely reflecting subsurface microbial activity. We suggest that the microbial activity within many diffuse flow vents is primarily attributable to seawater derived Gammaproteobacterial sulfur oxidizers, while in certain other flows vent-endemic Epsilonproteobactiera are most active. These data reveal a diversity in microbial activity at diffuse flows that has not previously been recognized, and reshapes our thinking about the relative influence that different microbial communities may have on local processes (such as primary production) and potentially global biogeochemical

Vertical flow constructed wetlands: kinetics of nutrient and organic matter removal.

PubMed

Pérez, M M; Hernández, J M; Bossens, J; Jiménez, T; Rosa, E; Tack, F

2014-01-01

The kinetics of organic matter and nutrient removal in a pilot vertical subsurface wetland with red ferralitic soil as substrate were evaluated. The wetland (20 m(2)) was planted with Cyperus alternifolius. The domestic wastewater that was treated in the wetland had undergone a primary treatment consisting of a septic moat and a buffer tank. From the sixth week of operation, the performance of the wetland stabilized, and a significant reduction in pollutant concentration of the effluent wastewater was obtained. Also a significant increase of dissolved oxygen (5 mg/l) was obtained. The organic matter removal efficiency was greater than 85% and the nutrient removal efficiency was greater than 75% in the vertical subsurface wetland. Nitrogen and biochemical oxygen demand (BOD) removal could be described by a first-order model. The kinetic constants were 3.64 and 3.27 d(-1) for BOD and for total nitrogen, respectively. Data on the removal of phosphorus were adapted to a second-order model. The kinetic constant was 0.96 (mg/l)(-1) d(-1). The results demonstrated the potential of vertical flow constructed wetlands to clean treated domestic wastewater before discharge into the environment.

Emerging organic contaminant removal depending on primary treatment and operational strategy in horizontal subsurface flow constructed wetlands: influence of redox.

PubMed

Avila, Cristina; Reyes, Carolina; Bayona, Josep María; García, Joan

2013-01-01

This study aimed at assessing the influence of primary treatment (hydrolytic upflow sludge blanket (HUSB) reactor vs. conventional settling) and operational strategy (alternation of saturated/unsaturated phases vs. permanently saturated) on the removal of various emerging organic contaminants (i.e. ibuprofen, diclofenac, acetaminophen, tonalide, oxybenzone, bisphenol A) in horizontal subsurface flow constructed wetlands. For that purpose, a continuous injection experiment was carried out in an experimental treatment plant for 26 days. The plant had 3 treatment lines: a control line (settler-wetland permanently saturated), a batch line (settler-wetland operated with saturate/unsaturated phases) and an anaerobic line (HUSB reactor-wetland permanently saturated). In each line, wetlands had a surface area of 2.95 m(2), a water depth of 25 cm and a granular medium D(60) = 7.3 mm, and were planted with common reed. During the study period the wetlands were operated at a hydraulic and organic load of 25 mm/d and about 4.7 g BOD/m(2)d, respectively. The injection experiment delivered very robust results that show how the occurrence of higher redox potentials within the wetland bed promotes the elimination of conventional quality parameters as well as emerging microcontaminants. Overall, removal efficiencies were always greater for the batch line than for the control and anaerobic lines, and to this respect statistically significantly differences were found for ibuprofen, diclofenac, oxybenzone and bisphenol A. As an example, ibuprofen, whose major removal mechanism has been reported to be biodegradation under aerobic conditions, showed a higher removal in the batch line (85%) than in the control (63%) and anaerobic (52%) lines. Bisphenol A showed also a great dependence on the redox status of the wetlands, finding an 89% removal rate for the batch line, as opposed to the control and anaerobic lines (79 and 65%, respectively). Furthermore, diclofenac showed a greater

Modeling subsurface stormflow initiation in low-relief landscapes

NASA Astrophysics Data System (ADS)

Hopp, Luisa; Vaché, Kellie B.; Rhett Jackson, C.; McDonnell, Jeffrey J.

2015-04-01

Shallow lateral subsurface flow as a runoff generating mechanism at the hillslope scale has mostly been studied in steeper terrain with typical hillside angles of 10 - 45 degrees. These studies have shown that subsurface stormflow is often initiated at the interface between a permeable upper soil layer and a lower conductivity impeding layer, e.g. a B horizon or bedrock. Many studies have identified thresholds of event size and soil moisture states that need to be exceeded before subsurface stormflow is initiated. However, subsurface stormflow generation on low-relief hillslopes has been much less studied. Here we present a modeling study that investigates the initiation of subsurface stormflow on low-relief hillslopes in the Upper Coastal Plain of South Carolina, USA. Hillslopes in this region typically have slope angles of 2-5 degrees. Topsoils are sandy, underlain by a low-conductivity sandy clay loam Bt horizon. Subsurface stormflow has only been intercepted occasionally in a 120 m long trench, and often subsurface flow was not well correlated with stream signals, suggesting a disconnect between subsurface flow on the hillslopes and stream flow. We therefore used a hydrologic model to better understand which conditions promote the initiation of subsurface flow in this landscape, addressing following questions: Is there a threshold event size and soil moisture state for producing lateral subsurface flow? What role does the spatial pattern of depth to the impeding clay layer play for subsurface stormflow dynamics? We reproduced a section of a hillslope, for which high-resolution topographic data and depth to clay measurements were available, in the hydrologic model HYDRUS-3D. Soil hydraulic parameters were based on experimentally-derived data. The threshold analysis was first performed using hourly climate data records for 2009-2010 from the study site to drive the simulation. For this period also trench measurements of subsurface flow were available. In addition

Channel water balance and exchange with subsurface flow along a mountain headwater stream in Montana, United States

USGS Publications Warehouse

Payn, R.A.; Gooseff, M.N.; McGlynn, B.L.; Bencala, K.E.; Wondzell, S.M.

2009-01-01

Channel water balances of contiguous reaches along streams represent a poorly understood scale of stream-subsurface interaction. We measured reach water balances along a headwater stream in Montana, United States, during summer base flow recessions. Reach water balances were estimated from series of tracer tests in 13 consecutive reaches delineated evenly along a 2.6 km valley segment. For each reach, we estimated net change in discharge, gross hydrologic loss, and gross hydrologic gain from tracer dilution and mass recovery. Four series of tracer tests were performed during relatively high, intermediate, and low base flow conditions. The relative distribution of channel water along the stream was strongly related to a transition in valley structure, with a general increase in gross losses through the recession. During tracer tests at intermediate and low flows, there were frequent substantial losses of tracer mass (>10%) that could not be explained by net loss in flow over the reach, indicating that many of the study reaches were concurrently losing and gaining water. For example, one reach with little net change in discharge exchanged nearly 20% of upstream flow with gains and losses along the reach. These substantial bidirectional exchanges suggest that some channel interactions with subsurface flow paths were not measurable by net change in flow or transient storage of recovered tracer. Understanding bidirectional channel water balances in stream reaches along valleys is critical to an accurate assessment of stream solute fate and transport and to a full assessment of exchanges between the stream channel and surrounding subsurface.

Channel water balance and exchange with subsurface flow along a mountain headwater stream in Montana, United States

USGS Publications Warehouse

Payn, R.A.; Gooseff, M.N.; McGlynn, B.L.; Bencala, K.E.; Wondzell, S.M.

2009-01-01

Channel water balances of contiguous reaches along streams represent a poorly understood scale of stream-subsurface interaction. We measured reach water balances along a headwater stream in Montana, United States, during summer base flow recessions. Reach water balances were estimated from series of tracer tests in 13 consecutive reaches delineated evenly along a 2.6 km valley segment. For each reach, we estimated net change in discharge, gross hydrologic loss, and gross hydrologic gain from tracer dilution and mass recovery. Four series of tracer tests were performed during relatively high, intermediate, and low base flow conditions. The relative distribution of channel water along the stream was strongly related to a transition in valley structure, with a general increase in gross losses through the recession. During tracer tests at intermediate and low flows, there were frequent substantial losses of tracer mass (>10%) that could not be explained by net loss in flow over the reach, indicating that many of the study reaches were concurrently losing and gaining water. For example, one reach with little net change in discharge exchanged nearly 20% of upstream flow with gains and losses along the reach. These substantial bidirectional exchanges suggest that some channel interactions with subsurface flow paths were not measurable by net change in flow or transient storage of recovered tracer. Understanding bidirectional channel water balances in stream reaches along valleys is critical to an accurate assessment of stream solute fate and transport and to a full assessment of exchanges between the stream channel and surrounding subsurface. Copyright 2009 by the American Geophysical Union.

Comparison of simple, small, full-scale sewage treatment systems in Brazil: UASB-maturation ponds-coarse filter; UASB-horizontal subsurface-flow wetland; vertical-flow wetland (first stage of French system).

PubMed

von Sperling, M

2015-01-01

This paper presents a comparison between three simple sewage treatment lines involving natural processes: (a) upflow anaerobic sludge blanket (UASB) reactor-three maturation ponds in series-coarse rock filter; (b) UASB reactor-horizontal subsurface-flow constructed wetland; and (c) vertical-flow constructed wetlands treating raw sewage (first stage of the French system). The evaluation was based on several years of practical experience with three small full-scale plants receiving the same influent wastewater (population equivalents of 220, 60 and 100 inhabitants) in the city of Belo Horizonte, Brazil. The comparison included interpretation of concentrations and removal efficiencies based on monitoring data (organic matter, solids, nitrogen, phosphorus, coliforms and helminth eggs), together with an evaluation of practical aspects, such as land and volume requirements, sludge production and handling, plant management, clogging and others. Based on an integrated evaluation of all aspects involved, it is worth emphasizing that each system has its own specificities, and no generalization can be made on the best option. The overall conclusion is that the three lines are suitable for sewage treatment in small communities in warm-climate regions.

The thermal impact of subsurface building structures on urban groundwater resources - A paradigmatic example.

PubMed

Epting, Jannis; Scheidler, Stefan; Affolter, Annette; Borer, Paul; Mueller, Matthias H; Egli, Lukas; García-Gil, Alejandro; Huggenberger, Peter

2017-10-15

Shallow subsurface thermal regimes in urban areas are increasingly impacted by anthropogenic activities, which include infrastructure development like underground traffic lines as well as industrial and residential subsurface buildings. In combination with the progressive use of shallow geothermal energy systems, this results in the so-called subsurface urban heat island effect. This article emphasizes the importance of considering the thermal impact of subsurface structures, which commonly is underestimated due to missing information and of reliable subsurface temperature data. Based on synthetic heat-transport models different settings of the urban environment were investigated, including: (1) hydraulic gradients and conductivities, which result in different groundwater flow velocities; (2) aquifer properties like groundwater thickness to aquitard and depth to water table; and (3) constructional features, such as building depths and thermal properties of building structures. Our results demonstrate that with rising groundwater flow velocities, the heat-load from building structures increase, whereas down-gradient groundwater temperatures decrease. Thermal impacts on subsurface resources therefore have to be related to the permeability of aquifers and hydraulic boundary conditions. In regard to the urban settings of Basel, Switzerland, flow velocities of around 1 md -1 delineate a marker where either down-gradient temperature deviations or heat-loads into the subsurface are more relevant. Furthermore, no direct thermal influence on groundwater resources should be expected for aquifers with groundwater thicknesses larger 10m and when the distance of the building structure to the groundwater table is higher than around 10m. We demonstrate that measuring temperature changes down-gradient of subsurface structures is insufficient overall to assess thermal impacts, particularly in urban areas. Moreover, in areas which are densely urbanized, and where groundwater flow

Opportunities and challenges for the application of SP measurements to monitor subsurface flow (Invited)

NASA Astrophysics Data System (ADS)

Jackson, M.; Vinogradov, J.; MacAllister, D.; Butler, A. P.; Leinov, E.; Zhang, J.

2013-12-01

Measurements of self-potential (SP) have been proposed or applied to monitor flow in the shallow subsurface in numerous settings, including volcanoes, earthquake zones, geothermal fields and hydrocarbon reservoirs, to detect leaks from dams, tanks and embankments, and to characterize groundwater flow and hydraulic properties. To interpret the measurements, it is generally assumed that the SP is dominated by the streaming potential, arising from the drag of excess electrical charge in the diffuse part of the electrical double layer at the mineral-fluid interfaces. The constitutive equation relating electrical current density j to the driving forces ▽V and ▽P is then j = -σ▽V -σC▽P=-σ▽V + Qv (1) where V is the streaming potential, P is the water pressure, σ is the saturated rock conductivity, v is the Darcy velocity, C is the streaming potential coupling coefficient, and Q is the excess charge transported by the flow. Equation (1) shows that there is a close relationship between flow properties of interest, such as the pressure gradient or Darcy velocity, and the streaming potential component of the SP. Hence SP measurements are an attractive method to monitor subsurface flow. However, the problem with interpreting the measurements is that both C and Q can vary over orders of magnitude, in response to variations in pore-water salinity, temperature, rock texture, and the presence of NAPLs in the pore-space. Moreover, additional current sources may be present if there are gradients in concentration or temperature, arising from differential rates of ion migration down gradient (diffusion potentials), and because of charge exclusion from the pore-space (exclusion potentials). In general, these additional current sources are neglected. This talk suggests a potential new opportunity for the application of SP measurements to monitor subsurface flow, in which the signal of interest arises from salinity rather than pressure gradients. Saline intrusion into

Trace Element and Cu Isotopic Tracers of Subsurface Flow and Transport in Wastewater Irrigated Soils

NASA Astrophysics Data System (ADS)

Carte, J.; Fantle, M. S.

2017-12-01

An understanding of subsurface flow paths is critical for quantifying the fate of contaminants in wastewater irrigation systems. This study investigates the subsurface flow of wastewater by quantifying the distribution of trace contaminants in wastewater irrigated soils. Soil samples were collected from the upper 1m of two wetlands at Penn State University's wastewater irrigation site, at which all effluent from the University's wastewater treatment plant has been sprayed since 1983. Major and trace element and Cu isotopic composition were determined for these samples, in addition to wastewater effluent and bedrock samples. The upper 20 cm of each wetland shows an enrichment of Bi, Cd, Cr, Cu, Mo, Ni, Pb, and Zn concentrations relative to deep (>1m) soils at the site by a factor of 1.7-3.5. Each wetland also has a subsurface clay rich horizon with Bi, Cu, Li, Ni, Pb, and Zn concentrations enriched by a factor of 1.4 to 5 relative to deep soils. These subsurface horizons directly underlie intervals that could facilitate preferential effluent flow: a gravel layer in one wetland, and a silty loam with visible mottling, an indication of dynamic water saturation, in the other. Trace metal concentrations in other horizons from both wetlands fall in the range of the deep soils. Significant variability in Cu isotopic composition is present in soils from both wetlands, with δ65Cu values ranging from 0.74‰ to 5.09‰. Soil δ65Cu correlates well with Cu concentrations, with lighter δ65Cu associated with higher concentrations. The Cu isotopic composition of the zones of metal enrichment are comparable to the ostensible average wastewater effluent δ65Cu value (0.61‰), while other horizons have considerably heavier δ65Cu values. We hypothesize that wastewater is the source of the metal enrichments, as each of the enriched elements are present as contaminants in wastewater, and the enrichments are located in clay-rich horizons conducive to trace metal immobilization due

Stochastic collocation using Kronrod-Patterson-Hermite quadrature with moderate delay for subsurface flow and transport

NASA Astrophysics Data System (ADS)

Liao, Q.; Tchelepi, H.; Zhang, D.

2015-12-01

Uncertainty quantification aims at characterizing the impact of input parameters on the output responses and plays an important role in many areas including subsurface flow and transport. In this study, a sparse grid collocation approach, which uses a nested Kronrod-Patterson-Hermite quadrature rule with moderate delay for Gaussian random parameters, is proposed to quantify the uncertainty of model solutions. The conventional stochastic collocation method serves as a promising non-intrusive approach and has drawn a great deal of interests. The collocation points are usually chosen to be Gauss-Hermite quadrature nodes, which are naturally unnested. The Kronrod-Patterson-Hermite nodes are shown to be more efficient than the Gauss-Hermite nodes due to nestedness. We propose a Kronrod-Patterson-Hermite rule with moderate delay to further improve the performance. Our study demonstrates the effectiveness of the proposed method for uncertainty quantification through subsurface flow and transport examples.

On the Representation of the Porosity-Pressure Relationship in General Subsurface Flow Codes

DOE PAGES

Birdsell, Daniel Traver; Karra, Satish; Rajaram, Harihar

2018-01-11

The governing equations for subsurface flow codes in a deformable porous media are derived from the balance of fluid mass and Darcy's equation. One class of these codes, which we call general subsurface flow codes (GSFs), allow for more general constitutive relations for material properties such as porosity, permeability and density. Examples of GSFs include PFLOTRAN, FEHM, TOUGH2, STOMP, and some reservoir simulators such as BOAST. Depending on the constitutive relations used in GSFs, an inconsistency arises between the standard groundwater flow equation and the governing equation of GSFs, and we clarify that the reason for this inconsistency is becausemore » the Darcy's equation used in the GSFs should account for the velocity of fluid with respect to solid. Due to lack of awareness of this inconsistency, users of the GSFs tend to use a porosity-pressure relationship that comes from the standard groundwater flow equation and assumes that the relative velocity is already accounted for. For the Theis problem, we show that using this traditional relationship in the GSFs leads to significantly large errors. We propose an alternate porosity-pressure relationship that is consistent with the derivation of the governing equations in the GSFs where the solid velocity is not tracked, and show that, with this relationship, the results are more accurate for the Theis problem. In conclusion, the purpose of this note is to make the users and developers of these GSFs aware of this inconsistency and to advocate that the alternate porosity model derived here should be incorporated in GSFs.« less

On the Representation of the Porosity-Pressure Relationship in General Subsurface Flow Codes

DOE Office of Scientific and Technical Information (OSTI.GOV)

Birdsell, Daniel Traver; Karra, Satish; Rajaram, Harihar

The governing equations for subsurface flow codes in a deformable porous media are derived from the balance of fluid mass and Darcy's equation. One class of these codes, which we call general subsurface flow codes (GSFs), allow for more general constitutive relations for material properties such as porosity, permeability and density. Examples of GSFs include PFLOTRAN, FEHM, TOUGH2, STOMP, and some reservoir simulators such as BOAST. Depending on the constitutive relations used in GSFs, an inconsistency arises between the standard groundwater flow equation and the governing equation of GSFs, and we clarify that the reason for this inconsistency is becausemore » the Darcy's equation used in the GSFs should account for the velocity of fluid with respect to solid. Due to lack of awareness of this inconsistency, users of the GSFs tend to use a porosity-pressure relationship that comes from the standard groundwater flow equation and assumes that the relative velocity is already accounted for. For the Theis problem, we show that using this traditional relationship in the GSFs leads to significantly large errors. We propose an alternate porosity-pressure relationship that is consistent with the derivation of the governing equations in the GSFs where the solid velocity is not tracked, and show that, with this relationship, the results are more accurate for the Theis problem. In conclusion, the purpose of this note is to make the users and developers of these GSFs aware of this inconsistency and to advocate that the alternate porosity model derived here should be incorporated in GSFs.« less

On the Representation of the Porosity-Pressure Relationship in General Subsurface Flow Codes

NASA Astrophysics Data System (ADS)

Birdsell, Daniel T.; Karra, Satish; Rajaram, Harihar

2018-02-01

The governing equations for subsurface flow codes in a deformable porous media are derived from the balance of fluid mass and Darcy's equation. One class of these codes, which we call general subsurface flow codes (GSFs), allow for more general constitutive relations for material properties such as porosity, permeability and density. Examples of GSFs include PFLOTRAN, FEHM, TOUGH2, STOMP, and some reservoir simulators such as BOAST. Depending on the constitutive relations used in GSFs, an inconsistency arises between the standard groundwater flow equation and the governing equation of GSFs, and we clarify that the reason for this inconsistency is because the Darcy's equation used in the GSFs should account for the velocity of fluid with respect to solid. Due to lack of awareness of this inconsistency, users of the GSFs tend to use a porosity-pressure relationship that comes from the standard groundwater flow equation and assumes that the relative velocity is already accounted for. For the Theis problem, we show that using this traditional relationship in the GSFs leads to significantly large errors. We propose an alternate porosity-pressure relationship that is consistent with the derivation of the governing equations in the GSFs where the solid velocity is not tracked, and show that, with this relationship, the results are more accurate for the Theis problem. The purpose of this note is to make the users and developers of these GSFs aware of this inconsistency and to advocate that the alternate porosity model derived here should be incorporated in GSFs.

On the Representation of Aquifer Compressibility in General Subsurface Flow Codes: How an Alternate Definition of Aquifer Compressibility Matches Results from the Groundwater Flow Equation

NASA Astrophysics Data System (ADS)

Birdsell, D.; Karra, S.; Rajaram, H.

2017-12-01

The governing equations for subsurface flow codes in deformable porous media are derived from the fluid mass balance equation. One class of these codes, which we call general subsurface flow (GSF) codes, does not explicitly track the motion of the solid porous media but does accept general constitutive relations for porosity, density, and fluid flux. Examples of GSF codes include PFLOTRAN, FEHM, STOMP, and TOUGH2. Meanwhile, analytical and numerical solutions based on the groundwater flow equation have assumed forms for porosity, density, and fluid flux. We review the derivation of the groundwater flow equation, which uses the form of Darcy's equation that accounts for the velocity of fluids with respect to solids and defines the soil matrix compressibility accordingly. We then show how GSF codes have a different governing equation if they use the form of Darcy's equation that is written only in terms of fluid velocity. The difference is seen in the porosity change, which is part of the specific storage term in the groundwater flow equation. We propose an alternative definition of soil matrix compressibility to correct for the untracked solid velocity. Simulation results show significantly less error for our new compressibility definition than the traditional compressibility when compared to analytical solutions from the groundwater literature. For example, the error in one calculation for a pumped sandstone aquifer goes from 940 to <70 Pa when the new compressibility is used. Code users and developers need to be aware of assumptions in the governing equations and constitutive relations in subsurface flow codes, and our newly-proposed compressibility function should be incorporated into GSF codes.

On the Representation of Aquifer Compressibility in General Subsurface Flow Codes: How an Alternate Definition of Aquifer Compressibility Matches Results from the Groundwater Flow Equation

NASA Astrophysics Data System (ADS)

Birdsell, D.; Karra, S.; Rajaram, H.

2016-12-01

The governing equations for subsurface flow codes in deformable porous media are derived from the fluid mass balance equation. One class of these codes, which we call general subsurface flow (GSF) codes, does not explicitly track the motion of the solid porous media but does accept general constitutive relations for porosity, density, and fluid flux. Examples of GSF codes include PFLOTRAN, FEHM, STOMP, and TOUGH2. Meanwhile, analytical and numerical solutions based on the groundwater flow equation have assumed forms for porosity, density, and fluid flux. We review the derivation of the groundwater flow equation, which uses the form of Darcy's equation that accounts for the velocity of fluids with respect to solids and defines the soil matrix compressibility accordingly. We then show how GSF codes have a different governing equation if they use the form of Darcy's equation that is written only in terms of fluid velocity. The difference is seen in the porosity change, which is part of the specific storage term in the groundwater flow equation. We propose an alternative definition of soil matrix compressibility to correct for the untracked solid velocity. Simulation results show significantly less error for our new compressibility definition than the traditional compressibility when compared to analytical solutions from the groundwater literature. For example, the error in one calculation for a pumped sandstone aquifer goes from 940 to <70 Pa when the new compressibility is used. Code users and developers need to be aware of assumptions in the governing equations and constitutive relations in subsurface flow codes, and our newly-proposed compressibility function should be incorporated into GSF codes.

Higher and lowest order mixed finite element approximation of subsurface flow problems with solutions of low regularity

NASA Astrophysics Data System (ADS)

Bause, Markus

2008-02-01

In this work we study mixed finite element approximations of Richards' equation for simulating variably saturated subsurface flow and simultaneous reactive solute transport. Whereas higher order schemes have proved their ability to approximate reliably reactive solute transport (cf., e.g. [Bause M, Knabner P. Numerical simulation of contaminant biodegradation by higher order methods and adaptive time stepping. Comput Visual Sci 7;2004:61-78]), the Raviart- Thomas mixed finite element method ( RT0) with a first order accurate flux approximation is popular for computing the underlying water flow field (cf. [Bause M, Knabner P. Computation of variably saturated subsurface flow by adaptive mixed hybrid finite element methods. Adv Water Resour 27;2004:565-581, Farthing MW, Kees CE, Miller CT. Mixed finite element methods and higher order temporal approximations for variably saturated groundwater flow. Adv Water Resour 26;2003:373-394, Starke G. Least-squares mixed finite element solution of variably saturated subsurface flow problems. SIAM J Sci Comput 21;2000:1869-1885, Younes A, Mosé R, Ackerer P, Chavent G. A new formulation of the mixed finite element method for solving elliptic and parabolic PDE with triangular elements. J Comp Phys 149;1999:148-167, Woodward CS, Dawson CN. Analysis of expanded mixed finite element methods for a nonlinear parabolic equation modeling flow into variably saturated porous media. SIAM J Numer Anal 37;2000:701-724]). This combination might be non-optimal. Higher order techniques could increase the accuracy of the flow field calculation and thereby improve the prediction of the solute transport. Here, we analyse the application of the Brezzi- Douglas- Marini element ( BDM1) with a second order accurate flux approximation to elliptic, parabolic and degenerate problems whose solutions lack the regularity that is assumed in optimal order error analyses. For the flow field calculation a superiority of the BDM1 approach to the RT0 one is

Vertical Subsurface Flow Mixing and Horizontal Anisotropy in Coarse Fluvial Aquifers: Structural Aspects

NASA Astrophysics Data System (ADS)

Huggenberger, P.; Huber, E.

2014-12-01

Detailed descriptions of the subsurface heterogeneities in coarse fluvial aquifer gravel often lack in concepts to distinguish between the essence and the noise of a permeability structure and the ability to extrapolate site specific hydraulic information at the tens to several hundred meters scale. At this scale the heterogeneity strongly influences the anisotropies of the flow field and the mixing processes in groundwater. However, in many hydrogeological models the complexity of natural systems is oversimplified. Understanding the link between the dynamics of the surface processes of braided-river systems and the resulting subsurface sedimentary structures is the key to characterizing the complexity of horizontal and vertical mixing processes in groundwater. From the different depositional elements of coarse braided-river systems, the largest permeability contrasts can be observed in the scour-fills. Other elements (e.g. different types of gravel sheets) show much smaller variabilities and could be considered as a kind of matrix. Field experiments on the river Tagliamento (Northeast Italy) based on morphological observation and ground-penetrating radar (GPR) surveys, as well as outcrop analyses of gravel pit exposures (Switzerland) allowed us to define the shape, sizes, spatial distribution and preservation potential of scour-fills. In vertical sections (e.g. 2D GPR data, vertical outcrop), the spatial density of remnant erosional bounding surfaces of scours is an indicator for the dynamics of the braided-river system (lateral mobility of the active floodplain, rate of sediment net deposition and spatial distribution of the confluence scours). In case of combined low aggradation rate and low lateral mobility the deposits may be dominated by a complex overprinting of scour-fills. The delineation of the erosional bounding surfaces, that are coherent over the survey area, is based on the identification of angular discontinuities of the reflectors. Fence diagrams

DIVERGENT HORIZONTAL SUB-SURFACE FLOWS WITHIN ACTIVE REGION 11158

DOE Office of Scientific and Technical Information (OSTI.GOV)

Jain, Kiran; Tripathy, S. C.; Hill, F., E-mail: kjain@nso.edu, E-mail: stripathy@nso.edu, E-mail: fhill@nso.edu

We measure the horizontal subsurface flow in a fast emerging active region (AR; NOAA 11158) using the ring-diagram technique and the Helioseismic and Magnetic Imager high spatial resolution Dopplergrams. This AR had a complex magnetic structure and displayed significant changes in morphology during its disk passage. Over a period of six days from 2011 February 11 to 16, the temporal variation in the magnitude of the total velocity is found to follow the trend of magnetic field strength. We further analyze regions of individual magnetic polarity within AR 11158 and find that the horizontal velocity components in these sub-regions havemore » significant variation with time and depth. The leading and trailing polarity regions move faster than the mixed-polarity region. Furthermore, both zonal and meridional components have opposite signs for trailing and leading polarity regions at all depths showing divergent flows within the AR. We also find a sharp decrease in the magnitude of total horizontal velocity in deeper layers around major flares. It is suggested that the re-organization of magnetic fields during flares, combined with the sunspot rotation, decreases the magnitude of horizontal flows or that the flow kinetic energy has been converted into the energy released by flares. After the decline in flare activity and sunspot rotation, the flows tend to follow the pattern of magnetic activity. We also observe less variation in the velocity components near the surface but these tend to increase with depth, further demonstrating that the deeper layers are more affected by the topology of ARs.« less

System analysis to estimate subsurface flow: from global level to the State of Minnesota

NASA Astrophysics Data System (ADS)

Shmagin, Boris A.; Kanivetsky, Roman

2002-06-01

Stream runoff data globally and in the state of Minnesota were used to estimate subsurface water flow. This system approach is based, in principal, on unity of groundwater and surface water systems, and it is in stark contrast to the traditional deterministic approach based on modeling. In coordination with methodology of system analysis, two levels of study were used to estimate subsurface flow. First, the global stream runoff data were assessed to estimate the temporal-spatial variability of surface water runoff. Factor analysis was used to study the temporal-spatial variability of global runoff for the period from 1918 to 1967. Results of these analysis demonstrate that the variability of global runoff could be represented by seven major components (factor scores) that could be grouped into seven distinct independent grouping from the total of 18 continental slopes on the Earth. Computed variance value in this analysis is 76% and supports such analysis. The global stream runoff for this period is stationary, and is more closely connected with the stream flow of Asia to the Pacific Ocean as well as with the stream runoff of North America towards the Arctic and Pacific Oceans. The second level examines the distribution of river runoff (annual and for February) for various landscapes and the hydrogeological conditions in the State of Minnesota (218,000 km2). The annual and minimal monthly rate of stream runoff for 115 gauging stations with a period of observation of 47 years (1935-1981) were used to characterize the spatio-temporal distribution of stream runoff in Minnesota. Results of this analysis demonstrate that the annual stream runoff rate changes from 6.3, towards 3.95, and then to 2.09 l s-1 km-2 (the difference is significant based on Student's criteria). These values in Minnesota correspond to ecological provinces from a mixed forest province towards the broadleaf forest and to prairie province, respectively. The distribution of minimal monthly stream

System analysis to estimate subsurface flow: From global level to the State of Minnesota

USGS Publications Warehouse

Shmagin, B.A.; Kanivetsky, R.

2002-01-01

Stream runoff data globally and in the state of Minnesota were used to estimate subsurface water flow. This system approach is based, in principal, on unity of groundwater and surface water systems, and it is in stark contrast to the traditional deterministic approach based on modeling. In coordination with methodology of system analysis, two levels of study were used to estimate subsurface flow. First, the global stream runoff data were assessed to estimate the temporal-spatial variability of surface water runoff. Factor analysis was used to study the temporal-spatial variability of global runoff for the period from 1918 to 1967. Results of these analysis demonstrate that the variability of global runoff could be represented by seven major components (factor scores) that could be grouped into seven distinct independent grouping from the total of 18 continental slopes on the Earth. Computed variance value in this analysis is 76% and supports such analysis. The global stream runoff for this period is stationary, and is more closely connected with the stream flow of Asia to the Pacific Ocean as well as with the stream runoff of North America towards the Arctic and Pacific Oceans. The second level examines the distribution of river runoff (annual and for February) for various landscapes and the hydrogeological conditions in the State of Minnesota (218,000 km2). The annual and minimal monthly rate of stream runoff for 115 gauging stations with a period of observation of 47 years (1935-1981) were used to characterize the spatio-temporal distribution of stream runoff in Minnesota. Results of this analysis demonstrate that the annual stream runoff rate changes from 6.3, towards 3.95, and then to 2.09 1 s-1 km-2 (the difference is significant based on Student's criteria). These values in Minnesota correspond to ecological provinces from a mixed forest province towards the broadleaf forest and to prairie province, respectively. The distribution of minimal monthly stream

Effects of basin size on low-flow stream chemistry and subsurface contact time in the neversink river watershed, New York

USGS Publications Warehouse

Wolock, D.M.; Fan, J.; Lawrence, G.B.

1997-01-01

The effects of basin size on low-flow stream chemistry and subsurface contact time were examined for a part of the Neversink River watershed in southern New York State. Acid neutralizing capacity (ANC), the sum of base cation concentrations (SBC), pH and concentrations of total aluminum (Al), dissolved organic carbon (DOC) and silicon (Si) were measured during low stream flow at the outlets of nested basins ranging in size from 0.2 to 166.3 km2. ANC, SBC, pH, Al and DOC showed pronounced changes as basin size increased from 0.2 to 3 km2, but relatively small variations were observed as basin size increased beyond 3 km2. An index of subsurface contact time computed from basin topography and soil hydraulic conductivity also showed pronounced changes as basin size increased from 0.2 to 3 km2 and smaller changes as basin size increased beyond 3 km2. These results suggest that basin size affects low-flow stream chemistry because of the effects of basin size on subsurface contact time. ?? 1997 by John Wiley & Sons, Ltd.

Performance of hybrid subsurface constructed wetland system for piggery wastewater treatment.

PubMed

Zhang, X; Inoue, T; Kato, K; Harada, J; Izumoto, H; Wu, D; Sakuragi, H; Ietsugu, H; Sugawara, Y

2016-01-01

The objective of this study was to evaluate performance of a hybrid constructed wetland (CW) built for high organic content piggery wastewater treatment in a cold region. The system consists of four vertical and one horizontal flow subsurface CWs. The wetland was built in 2009 and water quality was monitored from the outset. Average purification efficiency of this system was 95±5, 91±7, 89±8, 70±10, 84±15, 90±6, 99±2, and 93±16% for biochemical oxygen demand (BOD5), chemical oxygen demand (COD), total carbon (TC), total nitrogen (TN), ammonium-N (NH4-N), total phosphorus (TP), total coliform (T. Coliform), and suspended solids (SS), respectively during August 2010-December 2013. Pollutant removal rate was 15±18 g m(-2) d(-1), 49±52 g m(-2) d(-1), 6±4 g m(-2) d(-1), 7±5 g m(-2) d(-1), and 1±1 g m(-2) d(-1) for BOD5, COD, TN, NH4-N, and TP, respectively. The removal efficiency of BOD5, COD, NH4-N, and SS improved yearly since the start of operation. With respect to removal of TN and TP, efficiency improved in the first three years but slightly declined in the fourth year. The system performed well during both warm and cold periods, but was more efficient in the warm period. The nitrate increase may be attributed to a low C/N ratio, due to limited availability of carbon required for denitrification.

Subsurface flow pathway dynamics in the active layer of coupled permafrost-hydrogeological systems under seasonal and annual temperature variability.

NASA Astrophysics Data System (ADS)

Frampton, Andrew

2017-04-01

There is a need for improved understanding of the mechanisms controlling subsurface solute transport in the active layer in order to better understand permafrost-hydrological-carbon feedbacks, in particular with regards to how dissolved carbon is transported in coupled surface and subsurface terrestrial arctic water systems under climate change. Studying solute transport in arctic systems is also relevant in the context of anthropogenic pollution which may increase due to increased activity in cold region environments. In this contribution subsurface solute transport subject to ground surface warming causing permafrost thaw and active layer change is studied using a physically based model of coupled cryotic and hydrogeological flow processes combined with a particle tracking method. Changes in subsurface water flows and solute transport travel times are analysed for different modelled geological configurations during a 100-year warming period. Results show that for all simulated cases, the minimum and mean travel times increase non-linearly with warming irrespective of geological configuration and heterogeneity structure. The timing of the start of increase in travel time depends on heterogeneity structure, combined with the rate of permafrost degradation that also depends on material thermal and hydrogeological properties. These travel time changes are shown to depend on combined warming effects of increase in pathway length due to deepening of the active layer, reduced transport velocities due to a shift from horizontal saturated groundwater flow near the surface to vertical water percolation deeper into the subsurface, and pathway length increase and temporary immobilization caused by cryosuction-induced seasonal freeze cycles. The impact these change mechanisms have on solute and dissolved substance transport is further analysed by integrating pathway analysis with a Lagrangian approach, incorporating considerations for both dissolved organic and inorganic

Multiyear nutrient removal performance of three constructed wetlands intercepting tile drain flows from grazed pastures.

PubMed

Tanner, Chris C; Sukias, James P S

2011-01-01

Subsurface tile drain flows can be a major s ource of nurient loss from agricultural landscapes. This study quantifies flows and nitrogen and phosphorus yields from tile drains at three intensively grazed dairy pasture sites over 3- to 5-yr periods and evaluates the capacity of constructed wetlands occupying 0.66 to 1.6% of the drained catchments too reduce nutrient loads. Continuous flow records are combined with automated flow-proportional sampling of nutrient concentrations to calculate tile drain nutrient yields and wetland mass removal rates. Annual drainage water yields rangedfrom 193 to 564 mm (16-51% of rainfall) at two rain-fed sites and from 827 to 853 mm (43-51% of rainfall + irrigation) at an irrigated site. Annually, the tile drains exported 14 to 109 kg ha(-1) of total N (TN), of which 58 to 90% was nitrate-N. Constructed wetlands intercepting these flows removed 30 to 369 gTN m(-2) (7-63%) of influent loadings annually. Seasonal percentage nitrate-N and TN removal were negatively associated with wetland N mass loadings. Wetland P removal was poor in all wetlands, with 12 to 115% more total P exported annually overall than received. Annually, the tile drains exported 0.12 to 1.38 kg ha of total P, of which 15 to 93% was dissolved reactive P. Additional measures are required to reduce these losses or provide supplementary P removal. Wetland N removal performance could be improved by modifying drainage systems to release flows more gradually and improving irrigation practices to reduce drainage losses.

Paleomagnetic correlation of surface and subsurface basaltic lava flows and flow groups in the southern part of the Idaho National Laboratory, Idaho, with paleomagnetic data tables for drill cores

USGS Publications Warehouse

Champion, Duane E.; Hodges, Mary K.V.; Davis, Linda C.; Lanphere, Marvin A.

2011-01-01

Paleomagnetic inclination and polarity studies have been conducted on thousands of subcore samples from 51 coreholes located at and near the Idaho National Laboratory. These studies are used to paleomagnetically characterize and correlate successive stratigraphic intervals in each corehole to similar depth intervals in adjacent coreholes. Paleomagnetic results from 83 surface paleomagnetic sites, within and near the INL, are used to correlate these buried lava flow groups to basaltic shield volcanoes still exposed on the surface of the eastern Snake River Plain. Sample handling and demagnetization protocols are described as well as the paleomagnetic data averaging process. Paleomagnetic inclination comparisons between coreholes located only kilometers apart show comparable stratigraphic successions of mean inclination values over tens of meters of depth. At greater distance between coreholes, comparable correlation of mean inclination values is less consistent because flow groups may be missing or additional flow groups may be present and found at different depth intervals. Two shallow intersecting cross-sections, A-A- and B-B- (oriented southwest-northeast and northwest-southeast, respectively), drawn through southwest Idaho National Laboratory coreholes show the corehole to corehole or surface to corehole correlations derived from the paleomagnetic inclination data. From stratigraphic top to bottom, key results included the (1) Quaking Aspen Butte flow group, which erupted from Quaking Aspen Butte southwest of the Idaho National Laboratory, flowed northeast, and has been found in the subsurface in corehole USGS 132; (2) Vent 5206 flow group, which erupted near the southwestern border of the Idaho National Laboratory, flowed north and east, and has been found in the subsurface in coreholes USGS 132, USGS 129, USGS 131, USGS 127, USGS 130, USGS 128, and STF-AQ-01; and (3) Mid Butte flow group, which erupted north of U.S. Highway 20, flowed northwest, and has been

The Subsurface Flow and Transport Laboratory: A New Department of Energy User's Facility for Intermediate-Scale Experimentation

NASA Astrophysics Data System (ADS)

Wietsma, T. W.; Oostrom, M.; Foster, N. S.

2003-12-01

Intermediate-scale experiments (ISEs) for flow and transport are a valuable tool for simulating subsurface features and conditions encountered in the field at government and private sites. ISEs offer the ability to study, under controlled laboratory conditions, complicated processes characteristic of mixed wastes and heterogeneous subsurface environments, in multiple dimensions and at different scales. ISEs may, therefore, result in major cost savings if employed prior to field studies. A distinct advantage of ISEs is that researchers can design physical and/or chemical heterogeneities in the porous media matrix that better approximate natural field conditions and therefore address research questions that contain the additional complexity of processes often encountered in the natural environment. A new Subsurface Flow and Transport Laboratory (SFTL) has been developed for ISE users in the Environmental Spectroscopy & Biogeochemistry Facility in the Environmental Molecular Sciences Laboratory (EMSL) at Pacific Northwest National Laboratory (PNNL). The SFTL offers a variety of columns and flow cells, a new state-of-the-art dual-energy gamma system, a fully automated saturation-pressure apparatus, and analytical equipment for sample processing. The new facility, including qualified staff, is available for scientists interested in collaboration on conducting high-quality flow and transport experiments, including contaminant remediation. Close linkages exist between the SFTL and numerical modelers to aid in experimental design and interpretation. This presentation will discuss the facility and outline the procedures required to submit a proposal to use this unique facility for research purposes. The W. R. Wiley Environmental Molecular Sciences Laboratory, a national scientific user facility, is sponsored by the U.S. Department of Energy's Office of Biological and Environmental Research and located at Pacific Northwest National Laboratory.

Deformation band clusters on Mars and implications for subsurface fluid flow

USGS Publications Warehouse

Okubo, C.H.; Schultz, R.A.; Chan, M.A.; Komatsu, G.

2009-01-01

High-resolution imagery reveals unprecedented lines of evidence for the presence of deformation band clusters in layered sedimentary deposits in the equatorial region of Mars. Deformation bands are a class of geologic structural discontinuity that is a precursor to faults in clastic rocks and soils. Clusters of deformation bands, consisting of many hundreds of individual subparallel bands, can act as important structural controls on subsurface fluid flow in terrestrial reservoirs, and evidence of diagenetic processes is often preserved along them. Deformation band clusters are identified on Mars based on characteristic meter-scale architectures and geologic context as observed in data from the High-Resolution Imaging Science Experiment (HiRISE) camera. The identification of deformation band clusters on Mars is a key to investigating the migration of fluids between surface and subsurface reservoirs in the planet's vast sedimentary deposits. Similar to terrestrial examples, evidence of diagenesis in the form of light- and dark-toned discoloration and wall-rock induration is recorded along many of the deformation band clusters on Mars. Therefore, these structures are important sites for future exploration and investigations into the geologic history of water and water-related processes on Mars. ?? 2008 Geological Society of America.

A trench study to assess transfer of pesticides in subsurface lateral flow for a soil with contrasting texture on a sloping vineyard in Beaujolais.

PubMed

Peyrard, X; Liger, L; Guillemain, C; Gouy, V

2016-01-01

Subsurface lateral flow in both texture-contrast soils and catchments with shallow bedrock is suspected to be a non-point source of contamination of watercourses by pesticides used in agriculture. As a case study, the north of the Beaujolais region (eastern France) provides a favorable environment for such contamination due to its agro-pedo-climatic conditions. Environments seen in the Beaujolais region include intense viticulture, permeable and shallow soils, steep hillslopes, and storms that occur during the periods of pesticide application. Watercourse contamination by pesticides has been widely observed in this region, and offsite pesticide transport by subsurface lateral flow is suspected to be involved in diffuse and chronic presence of pesticides in surface water. In order to confirm and quantify the potential role of such processes in pesticide transfer, an automated trench system has been designed. The trench was set up on a steep farmed hillslope in a texture-contrast soil. It was equipped with a tipping bucket flow meter and an automatic sampler to monitor pesticide concentrations in lateral flow at fine resolution, by means of a flow-dependent sampling strategy. Four pesticides currently used in vine growing were studied to provide a range of mobility properties: one insecticide (chlorpyrifos-methyl) and three fungicides (spiroxamine, tebuconazole, and dimethomorph). With this system, it was possible to study pesticide concentration dynamics in the subsurface lateral flow, generated by substantial rainfall events following pesticide applications. The experimental design ascertained to be a suitable method in which to monitor subsurface lateral flow and related transfer of pesticides.

Analysis of chemical reaction kinetics of depredating organic pollutants from secondary effluent of wastewater treatment plant in constructed wetlands.

PubMed

Wang, Hao; Jiang, Dengling; Yang, Yong; Cao, Guoping

2013-01-01

Four subsurface constructed wetlands were built to treat the secondary effluent of a wastewater treatment plant in Tangshan, China. The chemical pollutant indexes of chemical oxygen demand (COD) were analyzed to evaluate the removal efficiency of organic pollutants from the secondary effluent of the wastewater treatment plant. In all cases, the subsurface constructed wetlands were efficient in treating organic pollutants. Under the same hydraulic loading condition, the horizontal flow wetlands exhibited better efficiency of COD removal than vertical flow wetlands: the removal rates in horizontal flow wetlands could be maintained at 68.4 ± 2.42% to 92.2 ± 1.61%, compared with 63.8 ± 1.19% to 85.0 ± 1.25% in the vertical flow wetlands. Meanwhile, the chemical reaction kinetics of organic pollutants was analyzed, and the results showed that the degradation courses of the four subsurface wetlands all corresponded with the first order reaction kinetics to a large extent.

Paleomagnetic correlation of basalt flows in selected coreholes near the Advanced Test Reactor Complex, the Idaho Nuclear Technology and Engineering Center, and along the southern boundary, Idaho National Laboratory, Idaho

USGS Publications Warehouse

Hodges, Mary K.V.; Champion, Duane E.

2016-10-03

The U.S. Geological Survey, in cooperation with the U.S. Department of Energy, used paleomagnetic data from 18 coreholes to construct three cross sections of subsurface basalt flows in the southern part of the Idaho National Laboratory (INL). These cross sections, containing descriptions of the subsurface horizontal and vertical distribution of basalt flows and sediment layers, will be used in geological studies, and to construct numerical models of groundwater flow and contaminant transport.Subsurface cross sections were used to correlate surface vents to their subsurface flows intersected by coreholes, to correlate subsurface flows between coreholes, and to identify possible subsurface vent locations of subsurface flows. Correlations were identified by average paleomagnetic inclinations of flows, and depth from land surface in coreholes, normalized to the North American Datum of 1927. Paleomagnetic data were combined, in some cases, with other data, such as radiometric ages of flows. Possible vent locations of buried basalt flows were identified by determining the location of the maximum thickness of flows penetrated by more than one corehole.Flows from the surface volcanic vents Quaking Aspen Butte, Vent 5206, Mid Butte, Lavatoo Butte, Crater Butte, Pond Butte, Vent 5350, Vent 5252, Tin Cup Butte, Vent 4959, Vent 5119, and AEC Butte are found in coreholes, and were correlated to the surface vents by matching their paleomagnetic inclinations, and in some cases, their stratigraphic positions.Some subsurface basalt flows that do not correlate to surface vents, do correlate over several coreholes, and may correlate to buried vents. Subsurface flows which correlate across several coreholes, but not to a surface vent include the D3 flow, the Big Lost flow, the CFA buried vent flow, the Early, Middle, and Late Basal Brunhes flows, the South Late Matuyama flow, the Matuyama flow, and the Jaramillo flow. The location of vents buried in the subsurface by younger basalt

Sequential Nitrification/Denitrification in Subsurface Flow Constructed Wetlands. A Literature Review

DTIC Science & Technology

1992-12-01

I I I I I iii I I I 3 TABLE OF CONTENTS SIG NATURE PAG E...3 SECTION I1. CONSTRUCTED WETLANDS: AN OVERVIEW ......................................... 4 A. Natural W...pollutants from secondary effluent in natural w etlands

dfnWorks: A discrete fracture network framework for modeling subsurface flow and transport

DOE Office of Scientific and Technical Information (OSTI.GOV)

Hyman, Jeffrey D.; Karra, Satish; Makedonska, Nataliia

DFNWORKS is a parallelized computational suite to generate three-dimensional discrete fracture networks (DFN) and simulate flow and transport. Developed at Los Alamos National Laboratory over the past five years, it has been used to study flow and transport in fractured media at scales ranging from millimeters to kilometers. The networks are created and meshed using DFNGEN, which combines FRAM (the feature rejection algorithm for meshing) methodology to stochastically generate three-dimensional DFNs with the LaGriT meshing toolbox to create a high-quality computational mesh representation. The representation produces a conforming Delaunay triangulation suitable for high performance computing finite volume solvers in anmore » intrinsically parallel fashion. Flow through the network is simulated in dfnFlow, which utilizes the massively parallel subsurface flow and reactive transport finite volume code PFLOTRAN. A Lagrangian approach to simulating transport through the DFN is adopted within DFNTRANS to determine pathlines and solute transport through the DFN. Example applications of this suite in the areas of nuclear waste repository science, hydraulic fracturing and CO 2 sequestration are also included.« less

dfnWorks: A discrete fracture network framework for modeling subsurface flow and transport

DOE PAGES

Hyman, Jeffrey D.; Karra, Satish; Makedonska, Nataliia; ...

2015-11-01

DFNWORKS is a parallelized computational suite to generate three-dimensional discrete fracture networks (DFN) and simulate flow and transport. Developed at Los Alamos National Laboratory over the past five years, it has been used to study flow and transport in fractured media at scales ranging from millimeters to kilometers. The networks are created and meshed using DFNGEN, which combines FRAM (the feature rejection algorithm for meshing) methodology to stochastically generate three-dimensional DFNs with the LaGriT meshing toolbox to create a high-quality computational mesh representation. The representation produces a conforming Delaunay triangulation suitable for high performance computing finite volume solvers in anmore » intrinsically parallel fashion. Flow through the network is simulated in dfnFlow, which utilizes the massively parallel subsurface flow and reactive transport finite volume code PFLOTRAN. A Lagrangian approach to simulating transport through the DFN is adopted within DFNTRANS to determine pathlines and solute transport through the DFN. Example applications of this suite in the areas of nuclear waste repository science, hydraulic fracturing and CO 2 sequestration are also included.« less

Effect of recirculation on organic matter removal in a hybrid constructed wetland system.

PubMed

Ayaz, S C; Findik, N; Akça, L; Erdoğan, N; Kinaci, C

2011-01-01

This research project aimed to determine the technologically feasible and applicable wastewater treatment systems which will be constructed to solve environmental problems caused by small communities in Turkey. Pilot-scale treatment of a small community's wastewater was performed over a period of more than 2 years in order to show applicability of these systems. The present study involves removal of organic matter and suspended solids in serially operated horizontal (HFCW) and vertical (VFCW) sub-surface flow constructed wetlands. The pilot-scale wetland was constructed downstream of anaerobic reactors at the campus of TUBITAK-MRC. Anaerobically pretreated wastewater was introduced into this hybrid two-stage sub-surface flow wetland system (TSCW). Wastewater was first introduced into the horizontal sub-surface flow system and then the vertical flow system before being discharged. Recirculation of the effluent was tested in the system. When the recirculation ratio was 100%, average removal efficiencies for TSCW were 91 +/- 4% for COD, 83 +/- 10% for BOD and 96 +/- 3% for suspended solids with average effluent concentrations of 9 +/- 5 mg/L COD, 6 +/- 3 mg/L BOD and 1 mg/L for suspended solids. Comparing non-recirculation and recirculation periods, the lowest effluent concentrations were obtained with a 100% recirculation ratio. The effluent concentrations met the Turkish regulations for discharge limits of COD, BOD and TSS in each case. The study showed that a hybrid constructed wetland system with recirculation is a very effective method of obtaining very low effluent organic matter and suspended solids concentrations downstream of anaerobic pretreatment of domestic wastewaters in small communities.

Moment Analysis Characterizing Water Flow in Repellent Soils from On- and Sub-Surface Point Sources

NASA Astrophysics Data System (ADS)

Xiong, Yunwu; Furman, Alex; Wallach, Rony

2010-05-01

Water repellency has a significant impact on water flow patterns in the soil profile. Flow tends to become unstable in such soils, which affects the water availability to plants and subsurface hydrology. In this paper, water flow in repellent soils was experimentally studied using the light reflection method. The transient 2D moisture profiles were monitored by CCD camera for tested soils packed in a transparent flow chamber. Water infiltration experiments and subsequent redistribution from on-surface and subsurface point sources with different flow rates were conducted for two soils of different repellency degrees as well as for wettable soil. We used spatio-statistical analysis (moments) to characterize the flow patterns. The zeroth moment is related to the total volume of water inside the moisture plume, and the first and second moments are affinitive to the center of mass and spatial variances of the moisture plume, respectively. The experimental results demonstrate that both the general shape and size of the wetting plume and the moisture distribution within the plume for the repellent soils are significantly different from that for the wettable soil. The wetting plume of the repellent soils is smaller, narrower, and longer (finger-like) than that of the wettable soil compared with that for the wettable soil that tended to roundness. Compared to the wettable soil, where the soil water content decreases radially from the source, moisture content for the water-repellent soils is higher, relatively uniform horizontally and gradually increases with depth (saturation overshoot), indicating that flow tends to become unstable. Ellipses, defined around the mass center and whose semi-axes represented a particular number of spatial variances, were successfully used to simulate the spatial and temporal variation of the moisture distribution in the soil profiles. Cumulative probability functions were defined for the water enclosed in these ellipses. Practically identical

Coarsening of three-dimensional structured and unstructured grids for subsurface flow

NASA Astrophysics Data System (ADS)

Aarnes, Jørg Espen; Hauge, Vera Louise; Efendiev, Yalchin

2007-11-01

We present a generic, semi-automated algorithm for generating non-uniform coarse grids for modeling subsurface flow. The method is applicable to arbitrary grids and does not impose smoothness constraints on the coarse grid. One therefore avoids conventional smoothing procedures that are commonly used to ensure that the grids obtained with standard coarsening procedures are not too rough. The coarsening algorithm is very simple and essentially involves only two parameters that specify the level of coarsening. Consequently the algorithm allows the user to specify the simulation grid dynamically to fit available computer resources, and, e.g., use the original geomodel as input for flow simulations. This is of great importance since coarse grid-generation is normally the most time-consuming part of an upscaling phase, and therefore the main obstacle that has prevented simulation workflows with user-defined resolution. We apply the coarsening algorithm to a series of two-phase flow problems on both structured (Cartesian) and unstructured grids. The numerical results demonstrate that one consistently obtains significantly more accurate results using the proposed non-uniform coarsening strategy than with corresponding uniform coarse grids with roughly the same number of cells.

3DFATMIC: THREE DIMENSIONAL SUBSURFACE FLOW, FATE AND TRANSPORT OF MICROBES AND CHEMICALS MODEL - USER'S MANUAL VERSION 1.0

EPA Science Inventory

This document is the user's manual of 3DFATMIC, a 3-Dimensional Subsurface Flow, Fate and Transport of Microbes and Chemicals Model using a Lagrangian-Eulerian adapted zooming and peak capturing (LEZOOMPC) algorithm.

Hydraulic characterization and optimization of total nitrogen removal in an aerated vertical subsurface flow treatment wetland.

PubMed

Boog, Johannes; Nivala, Jaime; Aubron, Thomas; Wallace, Scott; van Afferden, Manfred; Müller, Roland Arno

2014-06-01

In this study, a side-by-side comparison of two pilot-scale vertical subsurface flow constructed wetlands (6.2 m(2)×0.85 m, q(i)=95 L/m(2) d, τ(n)=3.5 d) handling primary treated domestic sewage was conducted. One system (VA-i) was set to intermittent aeration while the other was aerated continuously (VAp-c). Intermittent aeration was provided to VA-i in an 8 h on/4 h off pattern. The intermittently aerated wetland, VA-i, was observed to have 70% less nitrate nitrogen mass outflow than the continuously aerated wetland, VAp-c. Intermittent aeration was shown to increase treatment performance for TN while saving 33% of running energy cost for aeration. Parallel tracer experiments in the two wetlands showed hydraulic characteristics similar to one Continuously Stirred Tank Reactor (CSTR). Intermittent aeration did not significantly affect the hydraulic functioning of the system. Hydraulic efficiencies were 78% for VAp-c and 76% for VA-i. Copyright © 2014 Elsevier Ltd. All rights reserved.

Subsurface Controls on Stream Intermittency in a Semi-Arid Landscape

NASA Astrophysics Data System (ADS)

Dohman, J.; Godsey, S.; Thackray, G. D.; Hale, R. L.; Wright, K.; Martinez, D.

2017-12-01

Intermittent streams currently constitute 30% to greater than 50% of the global river network. In addition, the number of intermittent streams is expected to increase due to changes in land use and climate. These streams provide important ecosystem services, such as water for irrigation, increased biodiversity, and high rates of nutrient cycling. Many hydrological studies have focused on mapping current intermittent flow regimes or evaluating long-term flow records, but very few have investigated the underlying causes of stream intermittency. The disconnection and reconnection of surface flow reflects the capacity of the subsurface to accommodate flow, so characterizing subsurface flow is key to understanding stream drying. We assess how subsurface flow paths control local surface flows during low-flow periods, including intermittency. Water table dynamics were monitored in an intermittent reach of Gibson Jack Creek in southeastern Idaho. Four transects were delineated with a groundwater well located in the hillslope, riparian zone, and in the stream, for a total of 12 groundwater wells. The presence or absence of surface flow was determined by frequent visual observations as well as in situ loggers every 30m along the 200m study reach. The rate of surface water drying was measured in conjunction with temperature, precipitation, subsurface hydraulic conductivity, hillslope-riparian-stream connectivity and subsurface travel time. Initial results during an unusually wet year suggest different responses in reaches that were previously observed to occasionally cease flowing. Flows in the intermittent reaches had less coherent and lower amplitude diel variations during base flow periods than reaches that had never been observed to dry out. Our findings will help contribute to our understanding of mechanisms driving expansion and contraction cycles in intermittent streams, increase our ability to predict how land use and climate change will affect flow regimes, and

Groundwater flow pattern and related environmental phenomena in complex geologic setting based on integrated model construction

NASA Astrophysics Data System (ADS)

Tóth, Ádám; Havril, Tímea; Simon, Szilvia; Galsa, Attila; Monteiro Santos, Fernando A.; Müller, Imre; Mádl-Szőnyi, Judit

2016-08-01

Groundwater flow, driven, controlled and determined by topography, geology and climate, is responsible for several natural surface manifestations and affected by anthropogenic processes. Therefore, flowing groundwater can be regarded as an environmental agent. Numerical simulation of groundwater flow could reveal the flow pattern and explain the observed features. In complex geologic framework, where the geologic-hydrogeologic knowledge is limited, the groundwater flow model could not be constructed based solely on borehole data, but geophysical information could aid the model building. The integrated model construction was presented via the case study of the Tihany Peninsula, Hungary, with the aims of understanding the background and occurrence of groundwater-related environmental phenomena, such as wetlands, surface water-groundwater interaction, slope instability, and revealing the potential effect of anthropogenic activity and climate change. The hydrogeologic model was prepared on the basis of the compiled archive geophysical database and the results of recently performed geophysical measurements complemented with geologic-hydrogeologic data. Derivation of different electrostratigraphic units, revealing fracturing and detecting tectonic elements was achieved by systematically combined electromagnetic geophysical methods. The deduced information can be used as model input for groundwater flow simulation concerning hydrostratigraphy, geometry and boundary conditions. The results of numerical modelling were interpreted on the basis of gravity-driven regional groundwater flow concept and validated by field mapping of groundwater-related phenomena. The 3D model clarified the hydraulic behaviour of the formations, revealed the subsurface hydraulic connection between groundwater and wetlands and displayed the groundwater discharge pattern, as well. The position of wetlands, their vegetation type, discharge features and induced landslides were explained as

Reduced-Order Models Based on POD-Tpwl for Compositional Subsurface Flow Simulation

NASA Astrophysics Data System (ADS)

Durlofsky, L. J.; He, J.; Jin, L. Z.

2014-12-01

A reduced-order modeling procedure applicable for compositional subsurface flow simulation will be described and applied. The technique combines trajectory piecewise linearization (TPWL) and proper orthogonal decomposition (POD) to provide highly efficient surrogate models. The method is based on a molar formulation (which uses pressure and overall component mole fractions as the primary variables) and is applicable for two-phase, multicomponent systems. The POD-TPWL procedure expresses new solutions in terms of linearizations around solution states generated and saved during previously simulated 'training' runs. High-dimensional states are projected into a low-dimensional subspace using POD. Thus, at each time step, only a low-dimensional linear system needs to be solved. Results will be presented for heterogeneous three-dimensional simulation models involving CO2 injection. Both enhanced oil recovery and carbon storage applications (with horizontal CO2 injectors) will be considered. Reasonably close agreement between full-order reference solutions and compositional POD-TPWL simulations will be demonstrated for 'test' runs in which the well controls differ from those used for training. Construction of the POD-TPWL model requires preprocessing overhead computations equivalent to about 3-4 full-order runs. Runtime speedups using POD-TPWL are, however, very significant - typically O(100-1000). The use of POD-TPWL for well control optimization will also be illustrated. For this application, some amount of retraining during the course of the optimization is required, which leads to smaller, but still significant, speedup factors.

Fate of estrone in laboratory-scale constructed wetlands

USDA-ARS?s Scientific Manuscript database

A horizontal, subsurface, laboratory-scale constructed wetland (CW) consisting of four cells in series was used to determine the attenuation of the steroid hormone estrone (E1) present in animal wastewater. Liquid swine manure diluted 1:80 with farm pond water and dosed with [14C]E1 flowed through ...

Wastewater treatment in tsunami affected areas of Thailand by constructed wetlands.

PubMed

Brix, H; Koottatep, T; Laugesen, C H

2007-01-01

The tsunami of December 2004 destroyed infrastructure in many coastal areas in South-East Asia. In January 2005, the Danish Government gave a tsunami relief grant to Thailand to re-establish the wastewater management services in some of the areas affected by the tsunami. This paper describes the systems which have been built at three locations: (a) Baan Pru Teau: A newly-built township for tsunami victims which was constructed with the contribution of the Thai Red Cross. Conventional septic tanks were installed for the treatment of blackwater from each household and its effluent and grey water (40 m3/day) are collected and treated at a 220 m2 subsurface flow constructed wetland. (b) Koh Phi Phi Don island: A wastewater collection system for the main business and hotel area of the island, a pumping station and a pressure pipe to the treatment facility, a multi-stage constructed wetland system and a system for reuse of treated wastewater. The constructed wetland system (capacity 400 m3/day) consists of vertical flow, horizontal subsurface flow, free water surface flow and pond units. Because the treatment plant is surrounded by resorts, restaurants and shops, the constructed wetland systems are designed with terrains as scenic landscaping. (c) Patong: A 5,000 m2 constructed wetland system has been established to treat polluted water from drainage canals which collect overflow from septic tanks and grey water from residential areas. It is envisaged that these three systems will serve as prototype demonstration systems for appropriate wastewater management in Thailand and other tropical countries.

The impact of urbanization on subsurface flow paths - A paired-catchment isotopic study

NASA Astrophysics Data System (ADS)

Bonneau, Jeremie; Burns, Matthew J.; Fletcher, Tim D.; Witt, Roman; Drysdale, Russell N.; Costelloe, Justin F.

2018-06-01

Urbanization disturbs groundwater flow through the sealing of native soils with impervious surfaces and through modifications to the subsoil by constructed drainage and other infrastructure (trenches and excavations, e.g. water supply). The impact of these disturbances on groundwater contributions to urban streams (i.e. baseflow) is poorly understood. While high flows caused by impervious runoff to streams are a common focus of urban studies, the route taken by groundwater to become streamflow in urban landscapes is not generally considered. To assess the impact of urbanization on groundwater sources to streams, both rainfall and baseflow were sampled weekly for stable isotopes of water in two nearby streams-one draining a peri-urban catchment and the other draining a forested, natural catchment. In addition, to study the rate of groundwater discharge to the stream, monthly baseflow recession behavior was investigated. We found that baseflow in the forested catchment was constant in stable isotope values (δ18O = -5.73‰ ± 0.14‰) throughout the year. Monthly baseflow recession constants were close to 1 and had little variation (ranging 0.951-0.992), indicating a well-mixed groundwater store and long residence times. In contrast, the urban baseflow isotopic composition featured distinct seasonal variations (δ18O = -3.35‰ ± 1.20‰ from October to March and δ18O = -4.54‰ ± 0.43‰ from April to September) and high week-to-week variability in summer, reflecting a contribution of recent rainfall to baseflow. Recession constants were lower (ranging 0.727-0.955) with pronounced seasonal variations, suggesting shorter residence times and the likely presence of a variety of stores and pathways. These results provide evidence that the urban catchment has diversified groundwater pathways, and its groundwater storage is drained faster than that of the forested catchment. It highlights some of the subsurface hydrological consequences of urbanization. Restoring low-flow

Distribution and mass balance of hexavalent and trivalent chromium in a subsurface, horizontal flow (SF-h) constructed wetland operating as post-treatment of textile wastewater for water reuse.

PubMed

Fibbi, Donatella; Doumett, Saer; Lepri, Luciano; Checchini, Leonardo; Gonnelli, Cristina; Coppini, Ester; Del Bubba, Massimo

2012-01-15

In this study, during a two-year period, we investigated the fate of hexavalent and trivalent chromium in a full-scale subsurface horizontal flow constructed wetland planted with Phragmites australis. The reed bed operated as post-treatment of the effluent wastewater from an activated sludge plant serving the textile industrial district and the city of Prato (Italy). Chromium speciation was performed in influent and effluent wastewater and in water-suspended solids, at different depths and distances from the inlet; plants were also analyzed for total chromium along the same longitudinal profile. Removals of hexavalent and trivalent chromium equal to 72% and 26%, respectively were achieved. The mean hexavalent chromium outlet concentration was 1.6 ± 0.9 μg l(-1) and complied with the Italian legal limits for water reuse. Chromium in water-suspended solids was in the trivalent form, thus indicating that its removal from wastewater was obtained by the reduction of hexavalent chromium to the trivalent form, followed by accumulation of the latter inside the reed bed. Chromium in water-suspended solids was significantly affected by the distance from the inlet. Chromium concentrations in the different plant organs followed the same trend of suspended solids along the longitudinal profile and were much lower than those found in the solid material, evidencing a low metal accumulation in P. australis. Copyright © 2011 Elsevier B.V. All rights reserved.

Application of a fully integrated surface-subsurface physically based flow model for evaluating groundwater recharge from a flash flood event

NASA Astrophysics Data System (ADS)

Pino, Cristian; Herrera, Paulo; Therrien, René

2017-04-01

In many arid regions around the world groundwater recharge occurs during flash floods. This transient spatially and temporally concentrated flood-recharge process takes place through the variably saturated zone between surface and usually the deep groundwater table. These flood events are characterized by rapid and extreme changes in surface flow depth and velocity and soil moisture conditions. Infiltration rates change over time controlled by the hydraulic gradients and the unsaturated hydraulic conductivity at the surface-subsurface interface. Today is a challenge to assess the spatial and temporal distribution of groundwater recharge from flash flood events under real field conditions at different scales in arid areas. We apply an integrated surface-subsurface variably saturated physically-based flow model at the watershed scale to assess the recharge process during and after a flash flood event registered in an arid fluvial valley in Northern Chile. We are able to reproduce reasonably well observed groundwater levels and surface flow discharges during and after the flood with a calibrated model. We also investigate the magnitude and spatio-temporal distribution of recharge and the response of the system to variations of different surface and subsurface parameters, initial soil moisture content and groundwater table depths and surface flow conditions. We demonstrate how an integrated physically based model allows the exploration of different spatial and temporal system states, and that the analysis of the results of the simulations help us to improve our understanding of the recharge processes in similar type of systems that are common to many arid areas around the world.

BPA and NP removal from municipal wastewater by tropical horizontal subsurface constructed wetlands.

PubMed

Toro-Vélez, A F; Madera-Parra, C A; Peña-Varón, M R; Lee, W Y; Bezares-Cruz, J C; Walker, W S; Cárdenas-Henao, H; Quesada-Calderón, S; García-Hernández, H; Lens, P N L

2016-01-15

It has been recognized that numerous synthetic compounds like Bisphenol A (BPA) and nonylphenols (NP) are present in effluents from wastewater treatment plants (WWTP) at levels of parts per billion (μg L(-1)) or even parts per trillion (ng L(-1)) with a high potential to cause endocrine disruption in the aquatic environment. Constructed wetlands (CW) are a cost-effective wastewater treatment alternative with promising performance to treat these afore mentioned compounds. This research was aimed to evaluate the efficacy of CW treatment of WWTP effluent for mitigating the effects endocrine disrupting compounds (EDCs). This research goal was accomplished by (1) quantifying the removal of BPA and NP in CWs; (2) isolating CW fungal strains and testing for laccase production; and (3) performing endocrine disruption (reproduction) bioassays using the fruit fly Drosophila melanogaster. Three pilot scale horizontal subsurface flow constructed wetlands (HSSF-CW) were operated for eight weeks: one planted with Phragmites australis; one planted with Heliconia psitacorum; and one unplanted. The Heliconia CW showed a removal efficiency of 73.3(± 19%) and 62.8(± 20.1%) for BPA and NP, respectively; while the Phragmites CW demonstrated a similar removal for BPA (70.2 ± 27%) and lower removal efficiency for NP 52.1(± 37.1%).The unplanted CW achieved 62.2 (± 33%) removal for BPA and 25.3(± 37%) removal for NP. Four of the eleven fungal strains isolated from the Heliconia-CW showed the capacity to produce laccase. Even though complete removal of EDCs was not achieved by the CWs, the bioassay confirmed a significant improvement (p < 0.05) in fly viability for all CWs, with Heliconia sp. being the most effective at mitigating adverse effects on first and second generational reproduction. This study showed that a CW planted with a native Heliconia sp. CW demonstrated a higher removal of endocrine disrupting compounds and better mitigation of reproductive disruption in the

Channel water balance and exchange with subsurface flow along a mountain headwater stream in Montana, United States

Treesearch

R.A. Payn; M.N. Gooseff; B.L. McGlynn; K.E. Bencala; S.M. Wondzell

2009-01-01

Channel water balances of contiguous reaches along streams represent a poorly understood scale of stream-subsurface interaction. We measured reach water balances along a headwater stream in Montana, United States, during summer base flow recessions. Reach water balances were estimated from series of tracer tests in 13 consecutive reaches delineated evenly along a 2.6-...

Subsurface Water Flow and its Subsequent Impact on Chemical Behavior

USDA-ARS?s Scientific Manuscript database

The impact of the subsurface stratigraphy on crop growth and agrichemical behavior has been studied for several years at the OPE3 research site located at the USDA-ARS Beltsville Agricultural Research Center, in Beltsville Maryland. This site contains subsurface restricting layers that have been id...

Robust Representation of Integrated Surface-subsurface Hydrology at Watershed Scales

NASA Astrophysics Data System (ADS)

Painter, S. L.; Tang, G.; Collier, N.; Jan, A.; Karra, S.

2015-12-01

A representation of integrated surface-subsurface hydrology is the central component to process-rich watershed models that are emerging as alternatives to traditional reduced complexity models. These physically based systems are important for assessing potential impacts of climate change and human activities on groundwater-dependent ecosystems and water supply and quality. Integrated surface-subsurface models typically couple three-dimensional solutions for variably saturated flow in the subsurface with the kinematic- or diffusion-wave equation for surface flows. The computational scheme for coupling the surface and subsurface systems is key to the robustness, computational performance, and ease-of-implementation of the integrated system. A new, robust approach for coupling the subsurface and surface systems is developed from the assumption that the vertical gradient in head is negligible at the surface. This tight-coupling assumption allows the surface flow system to be incorporated directly into the subsurface system; effects of surface flow and surface water accumulation are represented as modifications to the subsurface flow and accumulation terms but are not triggered until the subsurface pressure reaches a threshold value corresponding to the appearance of water on the surface. The new approach has been implemented in the highly parallel PFLOTRAN (www.pflotran.org) code. Several synthetic examples and three-dimensional examples from the Walker Branch Watershed in Oak Ridge TN demonstrate the utility and robustness of the new approach using unstructured computational meshes. Representation of solute transport in the new approach is also discussed. Notice: This manuscript has been authored by UT-Battelle, LLC, under Contract No. DE-AC0500OR22725 with the U.S. Department of Energy. The United States Government retains and the publisher, by accepting the article for publication, acknowledges that the United States Government retains a non-exclusive, paid

The integration of constructed wetlands into a treatment system for airport runoff.

PubMed

Revitt, D M; Worral, P; Brewer, D

2001-01-01

A new surface runoff treatment system has been designed for London Heathrow Airport, which incorporates separate floating constructed wetlands or reedbeds and sub-surface flow constructed wetlands as major pollutant removal systems. The primary requirement of the newly developed treatment system is to control the concentrations of glycols following their use as de-icers and anti-icers within the airport. The ability of reedbeds to contribute to this treatment role was fully tested through pilot scale, on-site experiments over a 2 year period. The average reductions in runoff BOD concentrations achieved by pilot scale surface flow and sub-surface flow reedbeds were 30.9% and 32.9%, respectively. The corresponding average glycol removal efficiencies were 54.2% and 78.3%, following shock dosing inputs. These treatment performances are used to predict the required full scale constructed wetland surface areas needed to attain the desired effluent water quality. The treatment system also incorporates aeration, storage and, combined with reedbed technology, has been designed to reduce a mixed inlet BOD concentration of 240 mg/l to less than 40 mg/l for water temperatures varying between 6 degrees C and 20 degrees C.

Reactive transport codes for subsurface environmental simulation

DOE PAGES

Steefel, C. I.; Appelo, C. A. J.; Arora, B.; ...

2014-09-26

A general description of the mathematical and numerical formulations used in modern numerical reactive transport codes relevant for subsurface environmental simulations is presented. The formulations are followed by short descriptions of commonly used and available subsurface simulators that consider continuum representations of flow, transport, and reactions in porous media. These formulations are applicable to most of the subsurface environmental benchmark problems included in this special issue. The list of codes described briefly here includes PHREEQC, HPx, PHT3D, OpenGeoSys (OGS), HYTEC, ORCHESTRA, TOUGHREACT, eSTOMP, HYDROGEOCHEM, CrunchFlow, MIN3P, and PFLOTRAN. The descriptions include a high-level list of capabilities for each of themore » codes, along with a selective list of applications that highlight their capabilities and historical development.« less

Enhanced P, N and C removal from domestic wastewater using constructed wetland employing construction solid waste (CSW) as main substrate.

PubMed

Yang, Y; Wang, Z M; Liu, C; Guo, X C

2012-01-01

Construction solid waste (CSW), an inescapable by-product of the construction and demolition process, was used as main substrate in a four-stage vertical subsurface flow constructed wetland system to improve phosphorus P removal from domestic wastewater. A 'tidal flow' operation was also employed in the treatment system. Under a hydraulic loading rate (HLR) of 0.76 m3/m2 d for 1st and 3rd stage and HLR of 0.04 m3/m2 d for 2nd and 4th stage of the constructed wetland system respectively and tidal flow operation strategy, average removal efficiencies of 99.4% for P, 95.4% for ammoniacal-nitrogen, 56.5% for total nitrogen and 84.5% for total chemical oxygen demand were achieved during the operation period. The CSW-based constructed wetland system presents excellent P removal performance. The adoption of tidal flow strategy creates the aerobic/anoxic condition intermittently in the treatment system. This can achieve better oxygen transfer and hence lead to more complete nitrification and organic matter removal and enhanced denitrification. Overall, the CSW-based tidal flow constructed wetland system holds great promise for enabling high rate removal of P, ammoniacal-nitrogen and organic matter from domestic wastewater, and transforms CSW from a waste into a useful material.

The Revival of a Failed Constructed Wetland Treating of a High Fe Load AMD

Treesearch

A.D. Karathanasis; C.D. Barton

1999-01-01

Acid mine drainage (AMD) from abandoned mines has significantly impaired water quality in eastern Kentucky. A small surface flow wetland constructed in 1989 to reduce AMD effects and subsequently failed after six months of operation was renovated by incorporating anoxic limestone drains (ALDs) and anaerobic subsurface drains promoting vertical flow through successive...

Evaluation of constructed wetland treatment performance for winery wastewater.

PubMed

Grismer, Mark E; Carr, Melanie A; Shepherd, Heather L

2003-01-01

Rapid expansion of wineries in rural California during the past three decades has created contamination problems related to winery wastewater treatment and disposal; however, little information is available about performance of on-site treatment systems. Here, the project objective was to determine full-scale, subsurface-flow constructed wetland retention times and treatment performance through assessment of water quality by daily sampling of total dissolved solids, pH, total suspended solids, chemical oxygen demand (COD), tannins, nitrate, ammonium, total Kjeldahl nitrogen, phosphate, sulfate, and sulfide across operating systems for winery wastewater treatment. Measurements were conducted during both the fall crush season of heavy loading and the spring following bottling and racking operations at the winery. Simple decay model coefficients for these constituents as well as COD and tannin removal efficiencies from winery wastewater in bench-scale reactors are also determined. The bench-scale study used upward-flow, inoculated attached-growth (pea-gravel substrate) reactors fed synthetic winery wastewater. Inlet and outlet tracer studies for determination of actual retention times were essential to analyses of treatment performance from an operational subsurface-flow constructed wetland that had been overloaded due to failure to install a pretreatment system for suspended solids removal. Less intensive sampling conducted at a smaller operational winery wastewater constructed wetland that had used pretreatment suspended solids removal and aeration indicated that the constructed wetlands were capable of complete organic load removal from the winery wastewater.

Study on nitrogen removal enhanced by shunt distributing wastewater in a constructed subsurface infiltration system under intermittent operation mode.

PubMed

Li, Yinghua; Li, Haibo; Sun, Tieheng; Wang, Xin

2011-05-15

Subsurface wastewater infiltration system is an efficient and economic technology in treating small scattered sewage. The removal rates are generally satisfactory in terms of COD, BOD(5), TP and SS removal; while nitrogen removal is deficient in most of the present operating SWIS due to the different requirements for the presence of oxygen for nitrification and denitrification processes. To study the enhanced nitrogen removal technologies, two pilot subsurface wastewater infiltration systems were constructed in a village in Shenyang, China. The filled matrix was a mixture of 5% activated sludge, 65% brown soil and 30% coal slag in volume ratio for both systems. Intermittent operation mode was applied in to supply sufficient oxygen to accomplish the nitrification; meanwhile sewage was supplemented as the carbon source to the lower part in to denitrify. The constructed subsurface wastewater infiltration systems worked successfully under wetting-drying ratio of 1:1 with hydraulic loading of 0.081 m(3)/(m(2)d) for over 4 months. Carbon source was supplemented with shunt ratio of 1:1 and shunt position at the depth of 0.5m. The experimental results showed that intermittent operation mode and carbon source supplementation could significantly enhance the nitrogen removal efficiency with little influence on COD and TP removal. The average removal efficiencies for NH(3)-N and TN were 87.7 ± 1.4 and 70.1 ± 1.0%, increased by 12.5 ± 1.0 and 8.6 ± 0.7%, respectively. Copyright © 2011 Elsevier B.V. All rights reserved.

Probing Subsurface Flows in NOAA Active Region 12192: Comparison with NOAA 10486

NASA Astrophysics Data System (ADS)

Jain, Kiran; Tripathy, S. C.; Hill, F.

2017-11-01

NOAA Active Region (AR) 12192 is the biggest AR observed in solar cycle 24 so far. This was a long-lived AR that survived for four Carrington rotations (CRs) and exhibited several unusual phenomena. We measure the horizontal subsurface flows in this AR in multiple rotations using the ring-diagram technique of local helioseismology and the Global Oscillation Network Group (GONG+) Dopplergrams, and we investigate how different was the plasma flow in AR 12192 from that in AR 10486. Both regions produced several high M- and X-class flares, but they had different coronal mass ejection (CME) productivity. Our analysis suggests that these ARs had unusually large horizontal flow amplitude with distinctly different directions. While meridional flow in AR 12192 was poleward that supports the flux transport to poles, it was equatorward in AR 10486. Furthermore, there was a sudden increase in the magnitude of estimated zonal flow in shallow layers in AR 12192 during the X3.1 flare; however, it reversed direction in AR 10486 with the X17.2 flare. These flow patterns produced strong twists in horizontal velocity with depth in AR 10486 that persisted throughout the disk passage, as opposed to AR 12192, which produced a twist only after the eruption of the X3.1 flare that disappeared soon after. Our study indicates that the sunspot rotation combined with the reorganization of magnetic field in AR 10486 was not sufficient to decrease the flow energy even after several large flares that might have triggered CMEs. Furthermore, in the absence of sunspot rotation in AR 12192, this reorganization of magnetic field contributed significantly to the substantial release of flow energy after the X3.1 flare.

Stability of Gas Hydrates on Continental Margins: Implications of Subsurface Fluid Flow

NASA Astrophysics Data System (ADS)

Nunn, J. A.

2008-12-01

Gas hydrates are found at or just below the sediment-ocean interface in continental margins settings throughout the world. They are also found on land in high latitude regions such as the north slope of Alaska. While gas hydrate occurrence is common, gas hydrates are stable under a fairly restricted range of temperatures and pressures. In a purely conductive thermal regime, near surface temperatures depend on basal heat flow, thermal conductivity of sediments, and temperature at the sediment-water or sediment-air interface. Thermal conductivity depends on porosity and sediment composition. Gas hydrates are most stable in areas of low heat flow and high thermal conductivity which produce low temperature gradients. Older margins with thin continental crust and coarse grained sediments would tend to be colder. Another potentially important control on subsurface temperatures is advective heat transport by recharge/discharge of groundwater. Upward fluid flow depresses temperature gradients over a purely conductive regime with the same heat flow which would make gas hydrates more stable. Downward fluid flow would have the opposite effect. However, regional scale fluid flow may substantially increase heat flow in discharge areas which would destabilize gas hydrates. For example, discharge of topographically driven groundwater along the coast in the Central North Slope of Alaska has increased surface heat flow in some areas by more than 50% over a purely conductive thermal regime. Fluid flow also alters the pressure regime which can affect gas hydrate stability. Modeling results suggest a positive feedback between gas hydrate formation/disassociation and fluid flow. Disassociation of gas hydrates or permafrost due to global warming could increase permeability. This could enhance fluid flow and associated heat transport causing a more rapid and/or more spatially extensive gas hydrate disassociation than predicted solely from conductive propagation of temporal changes in

Enhanced organics and nitrogen removal in batch-operated vertical flow constructed wetlands by combination of intermittent aeration and step feeding strategy.

PubMed

Fan, Jinlin; Liang, Shuang; Zhang, Bo; Zhang, Jian

2013-04-01

Oxygen and carbon source supply are usually insufficient in subsurface flow constructed wetlands. Simultaneous removal of organic pollutants and nitrogen in five batch-operated vertical flow constructed wetlands under different operating conditions was investigated. Alternate aerobic and anaerobic regions were created well with intermittent aeration. Four-month experiments showed that the wetland-applied intermittent aeration combined with step feeding strategy (reactor E) greatly improved the removal of organics, ammonium nitrogen (NH4-N), and total nitrogen (TN) simultaneously, which were 97, 96, and 82%, respectively. It was much better than non-aerated reactors A and B and outperformed intermittently aerated reactor D without step feeding. Continuous aeration (reactor C) significantly enhanced the organics removal and nitrification, but it limited the TN removal (29%) seriously as a result of low denitrification level, and the high operation cost remained a question. The effect of plants was confirmed in this study, and the monitoring data showed that the plants could grow normally. Intermittent aeration as well as step feeding had no obvious influence on the growth of wetland plants in this study.

Exploring the role of mixing between subsurface flow paths on transit time distributions using a Lagrangian model

NASA Astrophysics Data System (ADS)

Zehe, Erwin; Jackisch, Conrad; Rodriguez, Nicolas; Klaus, Julian

2017-04-01

Only a minute amount of global fresh water is stored in the unsaturated zone. Yet this tiny compartment controls soil microbial activity and associated trace gas emissions, transport and transformations of contaminants, plant productivity, runoff generation and groundwater recharge. To date, the processes controlling renewal and age of different fractions of the soil water stock are far from being understood. Current theories and process concepts were largely inferred either from over-simplified laboratory experiments, or non-exhaustive point observations and tracer data in the field. Tracer data provide key but yet integrated information about the distribution of travel times of the tracer molecules to a certain depth or on their travel depth distribution within a given time. We hence are able to observe the "effect" of soil structure i.e. partitioning of infiltrating water between fast preferential and slow flow paths and imperfect subsequent mixing between these flow paths in the subsurface and the related plant water uptake. However, we are not able to study the "cause" - because technologies for in-situ observations of flow, flow path topology and exchange processes at relevant interfaces have up to now not been at hand. In the present study we will make use of a Lagrangian model for subsurface water dynamics to explore how subsurface heterogeneity and mixing among different storage fractions affects residence time distribution in the unsaturated zone in a forward approach. Soil water is represented by particles of constant mass, which travel according to the Itô form of the Fokker Planck equation. The model concept builds on established soil physics by estimating the drift velocity and the diffusion term based on the soil water characteristics. The model has been shown to simulate capillary driven soil moisture dynamics in good accordance with a) the Richards equation and b) observed soil moisture data in different soil. The particle model may furthermore

Recovery and fractionation of phosphorus retained by lightweight expanded shale and masonry sand used as media in subsurface flow treatment wetlands.

PubMed

Forbes, Margaret G; Dickson, Kenneth L; Saleh, Farida; Waller, William T; Doyle, Robert D; Hudak, Paul

2005-06-15

Most subsurface flow treatment wetlands, also known as reed bed or root zone systems, use sand or gravel substrates to reduce organics, solids, and nutrients in septic tank effluents. Phosphorus (P) retention in these systems is highly variable and few studies have identified the fate of retained P. In this study, two substrates, expanded shale and masonry sand, were used as filter media in five subsurface flow pilot-scale wetlands (2.7 m3). After 1 year of operation, we estimated the annual rate of P sorption by taking the difference between total P (TP) of substrate in the pilot cells and TP of substrate not exposed to wastewater (control). Means and standard deviations of TP retained by expanded shale were 349 +/- 171 mg kg(-1), respectively. For a substrate depth of 0.9 m, aerial P retention by shale was 201 +/- 98.6 g of P m(-2) year(-1), respectively. Masonry sand retained an insignificant quantity of wastewater P (11.9 +/- 21.8 mg kg(-1)) and on occasion exported P. Substrate samples were also sequentially fractionated into labile P, microbial P, (Fe + Al) P, humic P, (Ca + Mg) P, and residual P. In expanded shale samples, the greatest increase in P was in the relatively permanent form of (Fe + Al) P (108 mg kg(-1)), followed by labile P (46.7 mg kg(-1)) and humic P (39.8 mg kg(-1)). In masonry sand, there was an increase in labile P (9.71 mg kg(-1)). Results suggest that sand is a poor candidate for long-term P storage, but its efficiency is similar to that reported for many sand, gravel, and rock systems. By contrast, expanded shale and similar products with high hydraulic conductivity and P sorption capacity could greatly improve performance of P retention in constructed wetlands.

Subsurface lateral preferential flow network revealed by time-lapse ground-penetrating radar in a hillslope

NASA Astrophysics Data System (ADS)

Guo, Li; Chen, Jin; Lin, Henry

2014-12-01

Subsurface lateral preferential flow (LPF) has been observed to contribute substantially to hillslope and catchment runoff. However, the complex nature of LPF and the lack of an appropriate investigation method have hindered direct LPF observation in the field. Thus, the initiation, persistence, and dynamics of LPF networks remain poorly understood. This study explored the application of time-lapse ground-penetrating radar (GPR) together with an artificial infiltration to shed light on the nature of LPF and its dynamics in a hillslope. Based on our enhanced field experimental setup and carefully refined GPR data postprocessing algorithms, we developed a new protocol to reconstruct LPF networks with centimeter resolution. This is the first time that a detailed LPF network and its dynamics have been revealed noninvasively along a hillslope. Real-time soil water monitoring and field soil investigation confirmed the locations of LPF mapped by time-lapse GPR surveys. Our results indicated the following: (1) Increased spatial variations of radar signals after infiltration suggested heterogeneous soil water changes within the studied soil, which reflected the generation and dynamics of LPF; (2) Two types of LPF networks were identified, the network at the location of soil permeability contrasts and that formed via a series of connected preferential flow paths; and (3) The formation and distribution of LPF networks were influenced by antecedent soil water condition. Overall, this study demonstrates clearly that carefully designed time-lapse GPR surveys with enhanced data postprocessing offer a practical and nondestructive way of mapping LPF networks in the field, thereby providing a potentially significant enhancement in our ability to study complex subsurface flow processes across the landscape.

Information systems for material flow management in construction processes

NASA Astrophysics Data System (ADS)

Mesároš, P.; Mandičák, T.

2015-01-01

The article describes the options for the management of material flows in the construction process. Management and resource planning is one of the key factors influencing the effectiveness of construction project. It is very difficult to set these flows correctly. The current period offers several options and tools to do this. Information systems and their modules can be used just for the management of materials in the construction process.

Deep subsurface microbial processes

USGS Publications Warehouse

Lovley, D.R.; Chapelle, F.H.

1995-01-01

Information on the microbiology of the deep subsurface is necessary in order to understand the factors controlling the rate and extent of the microbially catalyzed redox reactions that influence the geophysical properties of these environments. Furthermore, there is an increasing threat that deep aquifers, an important drinking water resource, may be contaminated by man's activities, and there is a need to predict the extent to which microbial activity may remediate such contamination. Metabolically active microorganisms can be recovered from a diversity of deep subsurface environments. The available evidence suggests that these microorganisms are responsible for catalyzing the oxidation of organic matter coupled to a variety of electron acceptors just as microorganisms do in surface sediments, but at much slower rates. The technical difficulties in aseptically sampling deep subsurface sediments and the fact that microbial processes in laboratory incubations of deep subsurface material often do not mimic in situ processes frequently necessitate that microbial activity in the deep subsurface be inferred through nonmicrobiological analyses of ground water. These approaches include measurements of dissolved H2, which can predict the predominant microbially catalyzed redox reactions in aquifers, as well as geochemical and groundwater flow modeling, which can be used to estimate the rates of microbial processes. Microorganisms recovered from the deep subsurface have the potential to affect the fate of toxic organics and inorganic contaminants in groundwater. Microbial activity also greatly influences 1 the chemistry of many pristine groundwaters and contributes to such phenomena as porosity development in carbonate aquifers, accumulation of undesirably high concentrations of dissolved iron, and production of methane and hydrogen sulfide. Although the last decade has seen a dramatic increase in interest in deep subsurface microbiology, in comparison with the study of

Groundwater Flow Through a Constructed Treatment Wetland

DTIC Science & Technology

2003-03-01

the treatment wetland is to biodegrade perchloroethylene, which is present in the groundwater as a contaminant. Contaminated water enters the...characterizing groundwater flow through a constructed treatment wetland, one can visualize the flow paths of water through various types of soil. With...flowing groundwater and are now appearing in drinking water wells. Since contamination originated from government practices at many of these sites

Stochastic porous media modeling and high-resolution schemes for numerical simulation of subsurface immiscible fluid flow transport

NASA Astrophysics Data System (ADS)

Brantson, Eric Thompson; Ju, Binshan; Wu, Dan; Gyan, Patricia Semwaah

2018-04-01

This paper proposes stochastic petroleum porous media modeling for immiscible fluid flow simulation using Dykstra-Parson coefficient (V DP) and autocorrelation lengths to generate 2D stochastic permeability values which were also used to generate porosity fields through a linear interpolation technique based on Carman-Kozeny equation. The proposed method of permeability field generation in this study was compared to turning bands method (TBM) and uniform sampling randomization method (USRM). On the other hand, many studies have also reported that, upstream mobility weighting schemes, commonly used in conventional numerical reservoir simulators do not accurately capture immiscible displacement shocks and discontinuities through stochastically generated porous media. This can be attributed to high level of numerical smearing in first-order schemes, oftentimes misinterpreted as subsurface geological features. Therefore, this work employs high-resolution schemes of SUPERBEE flux limiter, weighted essentially non-oscillatory scheme (WENO), and monotone upstream-centered schemes for conservation laws (MUSCL) to accurately capture immiscible fluid flow transport in stochastic porous media. The high-order schemes results match well with Buckley Leverett (BL) analytical solution without any non-oscillatory solutions. The governing fluid flow equations were solved numerically using simultaneous solution (SS) technique, sequential solution (SEQ) technique and iterative implicit pressure and explicit saturation (IMPES) technique which produce acceptable numerical stability and convergence rate. A comparative and numerical examples study of flow transport through the proposed method, TBM and USRM permeability fields revealed detailed subsurface instabilities with their corresponding ultimate recovery factors. Also, the impact of autocorrelation lengths on immiscible fluid flow transport were analyzed and quantified. A finite number of lines used in the TBM resulted into visual

Linking Chaotic Advection with Subsurface Biogeochemical Processes

NASA Astrophysics Data System (ADS)

Mays, D. C.; Freedman, V. L.; White, S. K.; Fang, Y.; Neupauer, R.

2017-12-01

This work investigates the extent to which groundwater flow kinematics drive subsurface biogeochemical processes. In terms of groundwater flow kinematics, we consider chaotic advection, whose essential ingredient is stretching and folding of plumes. Chaotic advection is appealing within the context of groundwater remediation because it has been shown to optimize plume spreading in the laminar flows characteristic of aquifers. In terms of subsurface biogeochemical processes, we consider an existing model for microbially-mediated reduction of relatively mobile uranium(VI) to relatively immobile uranium(IV) following injection of acetate into a floodplain aquifer beneath a former uranium mill in Rifle, Colorado. This model has been implemented in the reactive transport code eSTOMP, the massively parallel version of STOMP (Subsurface Transport Over Multiple Phases). This presentation will report preliminary numerical simulations in which the hydraulic boundary conditions in the eSTOMP model are manipulated to simulate chaotic advection resulting from engineered injection and extraction of water through a manifold of wells surrounding the plume of injected acetate. This approach provides an avenue to simulate the impact of chaotic advection within the existing framework of the eSTOMP code.

Subsurface clade of Geobacteraceae that predominates in a diversity of Fe(III)-reducing subsurface environments

USGS Publications Warehouse

Holmes, Dawn E.; O'Neil, Regina A.; Vrionis, Helen A.; N'Guessan, Lucie A.; Ortiz-Bernad, Irene; Larrahondo, Maria J.; Adams, Lorrie A.; Ward, Joy A.; Nicoll , Julie S.; Nevin, Kelly P.; Chavan, Milind A.; Johnson, Jessica P.; Long, Philip E.; Lovely, Derek R.

2007-01-01

There are distinct differences in the physiology of Geobacter species available in pure culture. Therefore, to understand the ecology of Geobacter species in subsurface environments, it is important to know which species predominate. Clone libraries were assembled with 16S rRNA genes and transcripts amplified from three subsurface environments in which Geobacter species are known to be important members of the microbial community: (1) a uranium-contaminated aquifer located in Rifle, CO, USA undergoing in situ bioremediation; (2) an acetate-impacted aquifer that serves as an analog for the long-term acetate amendments proposed for in situ uranium bioremediation and (3) a petroleum-contaminated aquifer in which Geobacter species play a role in the oxidation of aromatic hydrocarbons coupled with the reduction of Fe(III). The majority of Geobacteraceae 16S rRNA sequences found in these environments clustered in a phylogenetically coherent subsurface clade, which also contains a number of Geobacter species isolated from subsurface environments. Concatamers constructed with 43 Geobacter genes amplified from these sites also clustered within this subsurface clade. 16S rRNA transcript and gene sequences in the sediments and groundwater at the Rifle site were highly similar, suggesting that sampling groundwater via monitoring wells can recover the most active Geobacter species. These results suggest that further study of Geobacter species in the subsurface clade is necessary to accurately model the behavior of Geobacter species during subsurface bioremediation of metal and organic contaminants.

Ebullition, Plant-Mediated Transport, and Subsurface Horizontal Water Flow Dominate Methane Transport in an Arctic Sphagnum Bog

NASA Astrophysics Data System (ADS)

Wehr, R. A.; McCalley, C. K.; Logan, T. A.; Chanton, J.; Crill, P. M.; Rich, V. I.; Saleska, S. R.

2017-12-01

Emission of the greenhouse gas methane from wetlands is of prime concern in the prediction of climate change - especially emission associated with thawing permafrost, which may drive a positive feedback loop of emission and warming. In addition to the biochemistry of methane production and consumption, wetland methane emission depends critically on the transport mechanisms by which methane moves through and out of the ecosystem. We therefore developed a model of methane biochemistry and transport for a sphagnum bog representing an intermediate permafrost thaw stage in Stordalen Mire, Sweden. In order to simultaneously reproduce measured profiles of both the concentrations and isotopic compositions of both methane and carbon dioxide in the peat pore water (Fig. 1) - as well as the surface methane emission - it was necessary for the model to include ebullition, plant-mediated transport via aerenchyma, and subsurface horizontal water flow. Diffusion of gas through the pore water was relatively unimportant. As a result, 90% of the produced methane escaped the wetland rather than being consumed by methanotrophic organisms in the near-surface pore water. Our model provides a comprehensive picture of methane emission from this bog site by quantifying the vertical profiles of: acetoclastic methanogenesis, hydrogenotrophic methanogenesis, methane oxidation, aerobic respiration, ebullition, plant-mediated transport, subsurface horizontal water flow, and diffusion.

Surface-subsurface turbulent interaction at the interface of a permeable bed: influence of the wall permeability

NASA Astrophysics Data System (ADS)

Kim, T.; Blois, G.; Best, J.; Christensen, K. T.

2017-12-01

Coarse-gravel river beds possess a high degree of permeability. Flow interactions between surface and subsurface flow across the bed interface is key to a number of natural processes occurring in the hyporheic zone. In fact, it is increasingly recognized that these interactions drive mass, momentum and energy transport across the interface, and consequently control biochemical processes as well as stability of sediments. The current study explores the role of the wall permeability in surface and subsurface flow interaction under controlled experimental conditions on a physical model of a gravel bed. The present wall model was constructed by five layers of cubically arranged spheres (d=25.4mm, where d is a diameter) providing 48% of porosity. Surface topography was removed by cutting half of a diameter on the top layer of spheres to render the flow surface smooth and highlight the impact of the permeability on the overlying flow. An impermeable smooth wall was also considered as a baseline of comparison for the permeable wall flow. To obtain basic flow statistics, low-frame-rate high-resolution PIV measurements were performed first in the streamwise-wall-normal (x-y) plane and refractive-index matching was employed to optically access the flow within the permeable wall. Time-resolved PIV experiments in the same facility were followed to investigate the flow interaction across the wall interface in sptaio-temporal domain. In this paper, a detailed analysis of the first and second order velocity statistics as well as the amplitude modulation for the flow overlying the permeable smooth wall will be presented.

Modeling of ground-water flow in subsurface Austin Chalk and Taylor marl in Ellis County, Texas, near the superconducting super collider site

DOE Office of Scientific and Technical Information (OSTI.GOV)

Mace, R.E.

1993-02-01

Numerical models are useful tools for developing an understanding of ground-water flow in sparsely characterized low-permeability aquifers. Finite-difference, cross-sectional models of Cretaceous chalk and marl formations near the Superconducting Super Collider (SSC) were constructed using MODFLOW to evaluate ground-water circulation paths and travel times. Weathered and fractured zones with enhanced permeability were included to assess the effect these features had on flow paths and times. Pump tests, slug tests, packer tests, core tests, and estimates were used to define hydraulic properties for model input. The model was calibrated with water-level data from monitor wells and from wire-line piezometers near amore » test shaft excavated by the SSC project. A ratio of vertical-to-horizontal permeability of 0.0085 was estimated through model calibration. A chalk-to-marl permeability ratio of 18 was needed to reproduce artesian head in a well completed in chalk beneath marl. Hydraulic head distributions and ground-water flow paths reflected local, intermediate, and regional flow systems with recharge beneath upland surface-water divides and discharge in valleys. Most of the flow (99%) occurred in the weathered zone, with average residence times of 5 to 10 years. Residence time in unweathered chalk bedrock was substantially longer, at an average of 1.7 Ma. As expected, the model demonstrated that deep and rapid ground-water circulation might occur in fracture zones. Particle paths calculated using MODPATH showed that ground-water travel times from recharge areas to the SSC subsurface facilities might be 20 to 60 years where flow is through fracture zones.« less

SUB-SURFACE MERIDIONAL FLOW, VORTICITY, AND THE LIFETIME OF SOLAR ACTIVE REGIONS

DOE Office of Scientific and Technical Information (OSTI.GOV)

Maurya, R. A.; Ambastha, A., E-mail: ramajor@prl.res.i, E-mail: ambastha@prl.res.i

Solar sub-surface fluid topology provides an indirect approach to examine the internal characteristics of active regions (ARs). Earlier studies have revealed the prevalence of strong flows in the interior of ARs having complex magnetic fields. Using the Doppler data obtained by the Global Oscillation Network Group project for a sample of 74 ARs, we have discovered the presence of steep gradients in meridional velocity at depths ranging from 1.5 to 5 Mm in flare productive ARs. The sample of these ARs is taken from the Carrington rotations 1980-2052 covering the period 2001 August-2007 January. The gradients showed an interesting hemisphericmore » trend of negative (positive) signs in the northern (southern) hemisphere, i.e., directed toward the equator. We have discovered three sheared layers in the depth range of 0-10 Mm, providing evidence of complex flow structures in several ARs. An important inference derived from our analysis is that the location of the deepest zero vertical vorticity is correlated with the remaining lifetime of ARs. This new finding may be employed as a tool for predicting the life expectancy of an AR.« less

Time-lapse 3-D seismic imaging of shallow subsurface contaminant flow.

PubMed

McKenna, J; Sherlock, D; Evans, B

2001-12-01

This paper presents a physical modelling study outlining a technique whereby buoyant contaminant flow within water-saturated unconsolidated sand was remotely monitored utilizing the time-lapse 3-D (TL3-D) seismic response. The controlled temperature and pressure conditions, along with the high level of acquisition repeatability attainable using sandbox physical models, allow the TL3-D seismic response to pore fluid movement to be distinguished from all other effects. TL3-D seismic techniques are currently being developed to monitor hydrocarbon reserves within producing reservoirs in an endeavour to improve overall recovery. However, in many ways, sandbox models under atmospheric conditions more accurately simulate the shallow subsurface than petroleum reservoirs. For this reason, perhaps the greatest application for analogue sandbox modelling is to improve our understanding of shallow groundwater and environmental flow mechanisms. Two fluid flow simulations were conducted whereby air and kerosene were injected into separate water-saturated unconsolidated sand models. In both experiments, a base 3-D seismic volume was recorded and compared with six later monitor surveys recorded while the injection program was conducted. Normal incidence amplitude and P-wave velocity information were extracted from the TL3-D seismic data to provide visualization of contaminant migration. Reflection amplitudes displayed qualitative areal distribution of fluids when a suitable impedance contrast existed between pore fluids. TL3-D seismic reflection tomography can potentially monitor the change in areal distribution of fluid contaminants over time, indicating flow patterns. However, other research and this current work have not established a quantifiable relationship between either normal reflection amplitudes and attenuation and fluid saturation. Generally, different pore fluids will have unique seismic velocities due to differences in compressibility and density. The predictable

An assessment of the performance of municipal constructed wetlands in Ireland.

PubMed

Hickey, Anthony; Arnscheidt, Joerg; Joyce, Eadaoin; O'Toole, James; Galvin, Gerry; O' Callaghan, Mark; Conroy, Ken; Killian, Darran; Shryane, Tommy; Hughes, Francis; Walsh, Katherine; Kavanagh, Emily

2018-03-15

While performance assessments of constructed wetlands sites around the world have appraised their capacity for effective removal of organics, a large variance remains in these sites' reported ability to retain nutrients, which appears to depend on differences in design, operation and climate factors. Nutrient retention is a very important objective for constructed wetlands, to avoid eutrophication of aquatic environments receiving their effluents. This study assessed the performance of constructed wetlands in terms of nutrient retention and associated parameters under the humid conditions of Ireland's temperate maritime climate. A review of the performance of 52 constructed wetland sites from 17 local authorities aimed to identify the best performing types of constructed wetlands and the treatment factors determining successful compliance with environmental standards. Data analysis compared effluent results from constructed wetlands with secondary free surface flow or tertiary horizontal subsurface flow, hybrid systems and integrated constructed wetlands with those from small-scale mechanical wastewater treatment plants of the same size class. Nutrient concentrations in effluents of constructed wetlands were negatively correlated (p < .01) with specific area, i.e. the ratio of surface area and population equivalents. The latest generation of integrated constructed wetlands, which had applied design guidelines issued by the Department of the Environment, performed best. Storm management design features improved treatment performance of constructed wetlands significantly (p < .05) for total suspended solids concentrations and exceedance frequency of limit values for total nitrogen. Mechanical wastewater treatment plants, secondary free surface water and tertiary horizontal subsurface flow wetlands showed a very large variance in effluent concentrations for organic and nutrient parameters. E. coli numbers in effluents were lowest for integrated constructed

THE HYDROLOGIC SYSTEM: GEOMORPHIC AND HYDROGEOLOGIC CONTROLS ON SURFACE AND SUBSURFACE FLOW REGIMES IN RIPARIAN MEADOW ECOSYSTEMS IN THE CENTRAL GREAT BASIN

EPA Science Inventory

Riparian corridors in upland watersheds in the Great Basin of central Nevada contain the majority of the region's biodiversity. Water, in both surface and subsurface flow regimes, is an important resource sustaining these sensitive ecosystems and other similar riparian ecosystem...

Deep thermal disturbances related to the sub-surface groundwater flow (Western Alps, France)

NASA Astrophysics Data System (ADS)

Mommessin, Grégoire; Dzikowski, Marc; Menard, Gilles; Monin, Nathalie

2013-04-01

In mountain area, the bedrock of the valley side is affected by a thickness of decompressed rock in subsurface (decompressed zone). Groundwater flowing in this zone disrupts the depth geothermal gradients. The evolution of thermal gradients under the decompressed zone depends of groundwater temperature changes into the decompressed zone. In this study, the phenomenon is studied from data acquired in exploration drilling prior to the construction of the France - Italy transalpine tunnel (High Speed Line project between Lyon and Turin). The study area is located in the Vanoise siliceous series between Modane and Avrieux (Western Alps, France). Of 31 boreholes, we selected 14 wells showing a natural thermal disturbance (not due to the drilling) linked to the groundwater flow in decompressed zone. The drill holes have a length between 200 and 1380m and well logs were carried out (gamma log, acoustic log, temperature log, flowmeter log). The rocks are constituted mainly by quartzite with high thermal conductivity or by schist and gneiss with low thermal conductivity. The decompressed zone concerns the quartzite with thicknesses ranging from 50m to 750m where groundwater flow imposes a constant temperature throughout the rock thickness. In the very low permeability rocks under the decompressed zone, the thermal gradient shows variations with depth. These variations suggest a water temperature change in the decompressed zone probably due to a paleoclimate event. We used the derived of the equation describing the propagation of a temperature in a 1D semi-infinite, in response to a sudden temperature disturbance at the boundary of the medium, to estimate the age and the amplitude of temperature change in the decompressed zone. The medium under the decompressed zone is supposed to be initially in a steady state and only conductive. Numerical tests assess that the 1D model is applicable in the slope context. The results obtained from 13 wells data show a few warming degrees

Tidal Response of Europa's Subsurface Ocean

NASA Astrophysics Data System (ADS)

Karatekin, O.; Comblen, R.; Deleersnijder, E.; Dehant, V. M.

2010-12-01

Time-variable tides in the subsurface oceans of icy satellites cause large periodic surface displacements and tidal dissipation can become a major energy source that can affect long-term orbital and internal evolution. In the present study, we investigate the response of the subsurface ocean of Europa to a time-varibale tidal potential. Two-dimensional nonlinear shallow water equations are solved on a sphere by means of a finite element code. The resulting ocean tidal flow velocities,dissipation and surface displacements will be presented.

GEOCHEMISTRY OF SUBSURFACE REACTIVE BARRIERS FOR REMEDIATION OF CONTAMINATED GROUND WATER

EPA Science Inventory

Reactive barriers that couple subsurface fluid flow with a passive chemical treatment zone are emerging, cost effective approaches for in-situ remediation of contaminated groundwater. Factors such as the build-up of surface precipitates, bio-fouling, and changes in subsurface tr...

Surface-subsurface flow modeling: an example of large-scale research at the new NEON user facility

NASA Astrophysics Data System (ADS)

Powell, H.; McKnight, D. M.

2009-12-01

Climate change is predicted to alter surface-subsurface interactions in freshwater ecosystems. These interactions are hypothesized to control nutrient release at diel and seasonal time scales, which may then exert control over epilithic algal growth rates. The mechanisms underlying shifts in complex physical-chemical-biological patterns can be elucidated by long-term observations at sites that span hydrologic and climate gradients across the continent. Development of the National Ecological Observatory Network (NEON) will provide researchers the opportunity to investigate continental-scale patterns by combining investigator-driven measurements with Observatory data. NEON is a national-scale research platform for analyzing and understanding the impacts of climate change, land-use change, and invasive species on ecology. NEON features sensor networks and experiments, linked by advanced cyberinfrastructure to record and archive ecological data for at least 30 years. NEON partitions the United States into 20 ecoclimatic domains. Each domain hosts one fully instrumented Core Aquatic site in a wildland area and one Relocatable site, which aims to capture ecologically significant gradients (e.g. landuse, nitrogen deposition, urbanization). In the current definition of NEON there are 36 Aquatic sites: 30 streams/rivers and 6 ponds/lakes. Each site includes automated, in-situ sensors for groundwater elevation and temperature; stream flow (discharge and stage); pond water elevation; atmospheric chemistry (Tair, barometric pressure, PAR, radiation); and surface water chemistry (DO, Twater, conductivity, pH, turbidity, cDOM, nutrients). Groundwater and surface water sites shall be regularly sampled for selected chemical and isotopic parameters. The hydrologic and geochemical monitoring design provides basic information on water and chemical fluxes in streams and ponds and between groundwater and surface water, which is intended to support investigator-driven modeling studies

3-D Reconstructions of Subsurface Pleistocene Basalt Flows from Paleomagnetic Inclination Data and 40Ar/39Ar Ages in the Southern Part of the Idaho National Laboratory (INL), Idaho (USA)

NASA Astrophysics Data System (ADS)

Hodges, M. K.; Champion, D. E.; Turrin, B. D.; Swisher, C. C.

2012-12-01

The U. S. Geological Survey, in cooperation with the U.S. Department of Energy, is mapping the distribution of basalt flows and sedimentary interbeds at the Idaho National Laboratory in three dimensions to provide data for refining numerical models of groundwater flow and contaminant transport in the eastern Snake River Plain aquifer. Paleomagnetic inclination and polarity data from basalt samples from 47 coreholes are being used to create a three-dimensional (3-D) model of the subsurface of the southern part of the INL. Surface and sub-surface basalt flows can be identified in individual cores and traced in three dimensions on the surface and in the subsurface for distances of more than 20 km using a combination of paleomagnetic, stratigraphic, and 40Ar/39Ar data. Eastern Snake River Plain olivine tholeiite basalts have K2O contents of 0.2 to 1.0 weight per cent. In spite of the low-K content, high-precision 40Ar/39Ar ages were obtained by applying a protocol that employs short irradiation times (minimizing interferences from Ca derived 36Ar), frequent measurement of various size atmospheric Ar pipettes to monitor and correct for temporal variation, and signal size dependent nonlinearity in spectrometer mass bias, resulting in age dates with resolution generally between 2 to 10% of the age. 3-D models of subsurface basalt flows are being used to: (1) Estimate eruption volumes; (2) locate the approximate vent areas and extent of sub-surface flows; and (3) Help locate high and low transmissivity zones. Results indicate that large basalt eruptions (>3 km3) occurred at and near the Central Facilities Area between 637 ka and 360 ka; at and near the Radioactive Waste Management Complex before 540 ka; and north of the Naval Reactors Facility at about 580 ka. Since about 360 ka, large basalt flows have erupted along the Arco-Big Southern Butte Volcanic Rift Zone and the Axial Volcanic Zone, and flowed northerly towards the Central Facilities Area. Basalt eruptions shifted

3-D reconstructions of subsurface Pleistocene basalt flows from paleomagnetic inclination data and 40Ar/39Ar ages in the southern part of the Idaho National Laboratory (INL), Idaho (USA)

USGS Publications Warehouse

Hodges, Mary K. V.; Champion, Duane E.; Turrin, B.D.; Swisher, C. C.

2012-01-01

The U. S. Geological Survey, in cooperation with the U.S. Department of Energy, is mapping the distribution of basalt flows and sedimentary interbeds at the Idaho National Laboratory in three dimensions to provide data for refining numerical models of groundwater flow and contaminant transport in the eastern Snake River Plain aquifer. Paleomagnetic inclination and polarity data from basalt samples from 47 coreholes are being used to create a three-dimensional (3-D) model of the subsurface of the southern part of the INL. Surface and sub-surface basalt flows can be identified in individual cores and traced in three dimensions on the surface and in the subsurface for distances of more than 20 km using a combination of paleomagnetic, stratigraphic, and 40Ar/39Ar data. Eastern Snake River Plain olivine tholeiite basalts have K2O contents of 0.2 to 1.0 weight per cent. In spite of the low-K content, high-precision 40Ar/39Ar ages were obtained by applying a protocol that employs short irradiation times (minimizing interferences from Ca derived 36Ar), frequent measurement of various size atmospheric Ar pipettes to monitor and correct for temporal variation, and signal size dependent nonlinearity in spectrometer mass bias, resulting in age dates with resolution generally between 2 to 10% of the age. 3-D models of subsurface basalt flows are being used to: (1) Estimate eruption volumes; (2) locate the approximate vent areas and extent of sub-surface flows; and (3) Help locate high and low transmissivity zones. Results indicate that large basalt eruptions (>3 km3) occurred at and near the Central Facilities Area between 637 ka and 360 ka; at and near the Radioactive Waste Management Complex before 540 ka; and north of the Naval Reactors Facility at about 580 ka. Since about 360 ka, large basalt flows have erupted along the Arco-Big Southern Butte Volcanic Rift Zone and the Axial Volcanic Zone, and flowed northerly towards the Central Facilities Area. Basalt eruptions shifted

Untangling the effects of urban development on subsurface storage in Baltimore

NASA Astrophysics Data System (ADS)

Bhaskar, Aditi S.; Welty, Claire; Maxwell, Reed M.; Miller, Andrew J.

2015-02-01

The impact of urban development on surface flow has been studied extensively over the last half century, but effects on groundwater systems are still poorly understood. Previous studies of the influence of urban development on subsurface storage have not revealed any consistent pattern, with results showing increases, decreases, and negligible change in groundwater levels. In this paper, we investigated the effects of four key features that impact subsurface storage in urban landscapes. These include reduced vegetative cover, impervious surface cover, infiltration and inflow (I&I) of groundwater and storm water into wastewater pipes, and other anthropogenic recharge and discharge fluxes including water supply pipe leakage and well and reservoir withdrawals. We applied the integrated groundwater-surface water-land surface model ParFlow.CLM to the Baltimore metropolitan area. We compared the base case (all four features) to simulations in which an individual urban feature was removed. For the Baltimore region, the effect of infiltration of groundwater into wastewater pipes had the greatest effect on subsurface storage (I&I decreased subsurface storage 11.1% relative to precipitation minus evapotranspiration after 1 year), followed by the impact of water supply pipe leakage and lawn irrigation (combined anthropogenic discharges and recharges led to a 7.4% decrease) and reduced vegetation (1.9% increase). Impervious surface cover led to a small increase in subsurface storage (0.56% increase) associated with decreased groundwater discharge as base flow. The change in subsurface storage due to infiltration of groundwater into wastewater pipes was largest despite the smaller spatial extent of surface flux modifications, compared to other features.

A multiple-point geostatistical method for characterizing uncertainty of subsurface alluvial units and its effects on flow and transport

USGS Publications Warehouse

Cronkite-Ratcliff, C.; Phelps, G.A.; Boucher, A.

2012-01-01

This report provides a proof-of-concept to demonstrate the potential application of multiple-point geostatistics for characterizing geologic heterogeneity and its effect on flow and transport simulation. The study presented in this report is the result of collaboration between the U.S. Geological Survey (USGS) and Stanford University. This collaboration focused on improving the characterization of alluvial deposits by incorporating prior knowledge of geologic structure and estimating the uncertainty of the modeled geologic units. In this study, geologic heterogeneity of alluvial units is characterized as a set of stochastic realizations, and uncertainty is indicated by variability in the results of flow and transport simulations for this set of realizations. This approach is tested on a hypothetical geologic scenario developed using data from the alluvial deposits in Yucca Flat, Nevada. Yucca Flat was chosen as a data source for this test case because it includes both complex geologic and hydrologic characteristics and also contains a substantial amount of both surface and subsurface geologic data. Multiple-point geostatistics is used to model geologic heterogeneity in the subsurface. A three-dimensional (3D) model of spatial variability is developed by integrating alluvial units mapped at the surface with vertical drill-hole data. The SNESIM (Single Normal Equation Simulation) algorithm is used to represent geologic heterogeneity stochastically by generating 20 realizations, each of which represents an equally probable geologic scenario. A 3D numerical model is used to simulate groundwater flow and contaminant transport for each realization, producing a distribution of flow and transport responses to the geologic heterogeneity. From this distribution of flow and transport responses, the frequency of exceeding a given contaminant concentration threshold can be used as an indicator of uncertainty about the location of the contaminant plume boundary.

Removal of nutrients from combined sewer overflows and lake water in a vertical-flow constructed wetland system.

PubMed

Gervin, L; Brix, H

2001-01-01

Lake Utterslev is situated in a densely built-up area of Copenhagen, and is heavily eutrophicated from combined sewer overflows. At the same time the lake suffers from lack of water. Therefore, a 5,000 m2 vertical flow wetland system was constructed in 1998 to reduce the phosphorus discharge from combined sewer overflows without reducing the water supply to the lake. During dry periods the constructed wetland is used to remove phosphorus from the lake water. The system is designed as a 90 m diameter circular bed with a bed depth of c. 2 m. The system is isolated from the surroundings by a polyethylene membrane. The bed medium consists of a mixture of gravel and crushed marble, which has a high binding capacity for phosphorus. The bed is located within the natural littoral zone of the lake and is planted with common reed (Phragmites australis). The constructed wetland is intermittently loaded with combined sewer overflow water or lake water and, after percolation through the bed medium, the water is collected in a network of drainage pipes at the bottom of the bed and pumped to the lake. The fully automated loading cycle results in alternating wet and dry periods. During the initial two years of operation, the phosphorus removal for combined sewer overflows has been consistently high (94-99% of inflow concentrations). When loaded with lake water, the phosphorus removal has been high during summer (71-97%) and lower during winter (53-75%) partly because of lower inlet concentrations. Effluent phosphorus concentrations are consistently low (0.03-0.04 mg/L). Ammonium nitrogen is nitrified in the constructed wetland, and total suspended solids and COD are generally reduced to concentrations below 5 mg/L and 25 mg/L, respectively. The study documents that a subsurface flow constructed wetland system can be designed and operated to effectively remove phosphorus and other pollutants from combined sewer overflows and eutrophicated lake water.

Whole-stream metabolism of a perennial spring-fed aufeis field in Alaska, with coincident surface and subsurface flow

NASA Astrophysics Data System (ADS)

Hendrickson, P. J.; Gooseff, M. N.; Huryn, A. D.

2017-12-01

Aufeis (icings or naleds) are seasonal arctic and sub-arctic features that accumulate through repeated overflow and freeze events of river or spring discharge. Aufeis fields, defined as the substrate on which aufeis form and the overlaying ice, have been studied to mitigate impacts on engineering structures; however, ecological characteristics and functions of aufeis fields are poorly understood. The perennial springs that supply warm water to aufeis fields create unique fluvial habitats, and are thought to act as winter and summer oases for biota. To investigate ecosystem function, we measured whole-stream metabolism at the Kuparuk River Aufeis (North Slope, AK), a large ( 5 km2) field composed of cobble substrate and predominately subsurface flow dynamics. The single-station open channel diel oxygen method was utilized at several dissolved oxygen (DO) stations located within and downstream of the aufeis field. DO loggers were installed in August 2016, and data downloaded summer 2017. Daily ecosystem respiration (ER), gross primary production (GPP) and reaeration rates were modeled using BASE, a package freely available in the open-source software R. Preliminary results support net heterotrophy during a two-week period of DO measurements in the fall season when minimum ice extent is observed. GPP, ER, and net metabolism are greater at the upstream reach near the spring source (P/R = 0.53), and decrease as flow moves downstream. As flow exits the aufeis field, surface and subsurface flow are incorporated into the metabolism model, and indicate the stream system becomes dependent on autochthonous production (P/R = 0.91). Current work is directed towards spring and summer discharge and metabolic parameter estimation, which is associated with maximum ice extent and rapid melting of the aufeis feature.

Evaluation of Subsurface Flow and Free-water Surface Wetlands Treating NPR-3 Produced Water - Year No. 1

DOE Office of Scientific and Technical Information (OSTI.GOV)

Myers, J. E.; Jackson, L. M.

2001-10-13

This paper is a summary of some of the activities conducted during the first year of a three-year cooperative research and development agreement (CRADA) between the Department of Energy (DOE) Rocky Mountain Oilfield Testing Center (RMOTC) and Texaco relating to the treatment of produced water by constructed wetlands. The first year of the CRADA is for design, construction and acclimation of the wetland pilot units. The second and third years of the CRADA are for tracking performance of pilot wetlands as the plant and microbial communities mature. A treatment wetland is a proven technology for the secondary and tertiary treatmentmore » of produced water, storm water and other wastewaters. Treatment wetlands are typically classified as either free-water surface (FWS) or subsurface flow (SSF). Both FWS and SSF wetlands work well when properly designed and operated. This paper presents a collection of kinetic data gathered from pilot units fed a slipstream of Wyoming (NPR-3) produced water. The pilot units are set up outdoors to test climatic influences on treatment. Monitoring parameters include evapotranspiration, plant growth, temperature, and NPDES discharge limits. The pilot wetlands (FWS and SSF) consist of a series of 100-gal plastic tubs filled with local soils, gravel, sharp sand and native wetland plants (cattail (Typha spp.), bulrush (Scirpus spp.), dwarf spikerush (Eleocharis)). Feed pumps control hydraulic retention time (HRT) and simple water control structures control the depth of water. The treated water is returned to the existing produced water treatment system. All NPDES discharge limits are met. Observations are included on training RMOTC summer students to do environmental work.« less

Form and function in hillslope hydrology: characterization of subsurface flow based on response observations

NASA Astrophysics Data System (ADS)

Angermann, Lisa; Jackisch, Conrad; Allroggen, Niklas; Sprenger, Matthias; Zehe, Erwin; Tronicke, Jens; Weiler, Markus; Blume, Theresa

2017-07-01

The phrase form and function was established in architecture and biology and refers to the idea that form and functionality are closely correlated, influence each other, and co-evolve. We suggest transferring this idea to hydrological systems to separate and analyze their two main characteristics: their form, which is equivalent to the spatial structure and static properties, and their function, equivalent to internal responses and hydrological behavior. While this approach is not particularly new to hydrological field research, we want to employ this concept to explicitly pursue the question of what information is most advantageous to understand a hydrological system. We applied this concept to subsurface flow within a hillslope, with a methodological focus on function: we conducted observations during a natural storm event and followed this with a hillslope-scale irrigation experiment. The results are used to infer hydrological processes of the monitored system. Based on these findings, the explanatory power and conclusiveness of the data are discussed. The measurements included basic hydrological monitoring methods, like piezometers, soil moisture, and discharge measurements. These were accompanied by isotope sampling and a novel application of 2-D time-lapse GPR (ground-penetrating radar). The main finding regarding the processes in the hillslope was that preferential flow paths were established quickly, despite unsaturated conditions. These flow paths also caused a detectable signal in the catchment response following a natural rainfall event, showing that these processes are relevant also at the catchment scale. Thus, we conclude that response observations (dynamics and patterns, i.e., indicators of function) were well suited to describing processes at the observational scale. Especially the use of 2-D time-lapse GPR measurements, providing detailed subsurface response patterns, as well as the combination of stream-centered and hillslope-centered approaches

Urban heat fluxes in the subsurface of Cologne, Germany

NASA Astrophysics Data System (ADS)

Zhu, K.; Bayer, P.; Blum, P.

2012-04-01

Urbanization during the last hundred years has led to both environmental and thermal impacts on the subsurface. The urban heat island (UHI) effect is mostly described as an atmospheric phenomenon, where the measured aboveground temperatures in cities are elevated in comparison to undisturbed rural regions. However, UHIs can be found below, as well as above ground. A large amount of anthropogenic heat migrates into the urban subsurface, which also raises the ground temperature and permanently changes the thermal conditions in shallow aquifers. The main objective of our work is to study and determine the urban heat fluxes in Cologne, Germany, and to improve our understanding of the dynamics of subsurface energy fluxes in UHIs. Ideally, our findings will contribute to strategic and more sustainable geothermal use in cities. For a quantitative analysis of the energy fluxes within the subsurface and across the atmospheric boundary, two and three-dimensional coupled numerical flow and heat transport models were developed. The simulation results indicate that during the past hundred years, an average vertical urban heat flux that ranges between 80 and 375 mW m-2 can be deduced. Thermal anomalies have migrated into the local urban aquifer system and they reach a depth of about 150 m. In this context, the influence of the regional groundwater flow on the subsurface heat transport and temperature development is comprehensively discussed.

Subsurface flow recently triggered the development of taliks under a sub-Arctic road embankment: a prelude to the collapse of the Alaska Highway?

NASA Astrophysics Data System (ADS)

Chen, L.; Fortier, D.; Sliger, M.; McKenzie, J. M.; Murchison, P.

2017-12-01

The Alaska Highway extends over 2200 km between central Alaska, U.S.A. and northern British-Columbia, Canada. This transportation corridor is crucial for the economy of Alaska as it is the only terrestrial link between mainland Alaska and the contiguous United States. Northern British Columbia and southwestern Yukon also greatly benefit from this highway for the transportation of goods and people across this remote corner of Canada. About a quarter of the Alaska Highway is built on permafrost, which is typically ice-rich and at a temperature near the point of thawing. Degradation of the permafrost under the embankment has led to severe structural damages to the highway such as deep longitudinal cracks, extended depressions, potholes and sinkholes. Here we present thermal data from the Beaver Creek experimental road test section in southwestern Yukon. Our study investigates convective heat transfers linked to subsurface water flow under the road embankment based on seven years (2009 to 2016) of thermal monitoring. Observation results demonstrate that snowmelt water infiltration in the spring causes rapid temperature increase of the upper portion of the embankment. Later in the summer, subsurface flow under the highway embankment can lead to step temperature-increase rates, which can be 200 times larger than those via conductive heat transfers. In the fall water trapped under the road significantly delays freeze back of the active layer and contributes to higher permafrost temperature. During the monitoring period, we observed the initiation and growth of taliks along sub-surface flow paths. Positive feedback mechanisms related to water flow through the taliks significantly increased permafrost degradation. Such taliks represent an un-precedent and presumably irreversible thermal state of the highway. Similar terrain conditions which severely threaten the structural integrity of the infrastructure on the short term are numerous along the Alaska Highway corridor.

Coupling a three-dimensional subsurface flow and transport model with a land surface model to simulate stream-aquifer-land interactions (CP v1.0)

NASA Astrophysics Data System (ADS)

Bisht, Gautam; Huang, Maoyi; Zhou, Tian; Chen, Xingyuan; Dai, Heng; Hammond, Glenn E.; Riley, William J.; Downs, Janelle L.; Liu, Ying; Zachara, John M.

2017-12-01

A fully coupled three-dimensional surface and subsurface land model is developed and applied to a site along the Columbia River to simulate three-way interactions among river water, groundwater, and land surface processes. The model features the coupling of the Community Land Model version 4.5 (CLM4.5) and a massively parallel multiphysics reactive transport model (PFLOTRAN). The coupled model, named CP v1.0, is applied to a 400 m × 400 m study domain instrumented with groundwater monitoring wells along the Columbia River shoreline. CP v1.0 simulations are performed at three spatial resolutions (i.e., 2, 10, and 20 m) over a 5-year period to evaluate the impact of hydroclimatic conditions and spatial resolution on simulated variables. Results show that the coupled model is capable of simulating groundwater-river-water interactions driven by river stage variability along managed river reaches, which are of global significance as a result of over 30 000 dams constructed worldwide during the past half-century. Our numerical experiments suggest that the land-surface energy partitioning is strongly modulated by groundwater-river-water interactions through expanding the periodically inundated fraction of the riparian zone, and enhancing moisture availability in the vadose zone via capillary rise in response to the river stage change. Meanwhile, CLM4.5 fails to capture the key hydrologic process (i.e., groundwater-river-water exchange) at the site, and consequently simulates drastically different water and energy budgets. Furthermore, spatial resolution is found to significantly impact the accuracy of estimated the mass exchange rates at the boundaries of the aquifer, and it becomes critical when surface and subsurface become more tightly coupled with groundwater table within 6 to 7 meters below the surface. Inclusion of lateral subsurface flow influenced both the surface energy budget and subsurface transport processes as a result of river-water intrusion into the

Coupling a three-dimensional subsurface flow and transport model with a land surface model to simulate stream–aquifer–land interactions (CP v1.0)

DOE Office of Scientific and Technical Information (OSTI.GOV)

Bisht, Gautam; Huang, Maoyi; Zhou, Tian

A fully coupled three-dimensional surface and subsurface land model is developed and applied to a site along the Columbia River to simulate three-way interactions among river water, groundwater, and land surface processes. The model features the coupling of the Community Land Model version 4.5 (CLM4.5) and a massively parallel multiphysics reactive transport model (PFLOTRAN). The coupled model, named CP v1.0, is applied to a 400 m × 400 m study domain instrumented with groundwater monitoring wells along the Columbia River shoreline. CP v1.0 simulations are performed at three spatial resolutions (i.e., 2, 10, and 20 m) over a 5-year period to evaluate themore » impact of hydroclimatic conditions and spatial resolution on simulated variables. Results show that the coupled model is capable of simulating groundwater–river-water interactions driven by river stage variability along managed river reaches, which are of global significance as a result of over 30 000 dams constructed worldwide during the past half-century. Our numerical experiments suggest that the land-surface energy partitioning is strongly modulated by groundwater–river-water interactions through expanding the periodically inundated fraction of the riparian zone, and enhancing moisture availability in the vadose zone via capillary rise in response to the river stage change. Meanwhile, CLM4.5 fails to capture the key hydrologic process (i.e., groundwater–river-water exchange) at the site, and consequently simulates drastically different water and energy budgets. Furthermore, spatial resolution is found to significantly impact the accuracy of estimated the mass exchange rates at the boundaries of the aquifer, and it becomes critical when surface and subsurface become more tightly coupled with groundwater table within 6 to 7 meters below the surface. Inclusion of lateral subsurface flow influenced both the surface energy budget and subsurface transport processes as a result of river

Coupling a three-dimensional subsurface flow and transport model with a land surface model to simulate stream–aquifer–land interactions (CP v1.0)

DOE PAGES

Bisht, Gautam; Huang, Maoyi; Zhou, Tian; ...

2017-01-01

A fully coupled three-dimensional surface and subsurface land model is developed and applied to a site along the Columbia River to simulate three-way interactions among river water, groundwater, and land surface processes. The model features the coupling of the Community Land Model version 4.5 (CLM4.5) and a massively parallel multiphysics reactive transport model (PFLOTRAN). The coupled model, named CP v1.0, is applied to a 400 m × 400 m study domain instrumented with groundwater monitoring wells along the Columbia River shoreline. CP v1.0 simulations are performed at three spatial resolutions (i.e., 2, 10, and 20 m) over a 5-year period to evaluate themore » impact of hydroclimatic conditions and spatial resolution on simulated variables. Results show that the coupled model is capable of simulating groundwater–river-water interactions driven by river stage variability along managed river reaches, which are of global significance as a result of over 30 000 dams constructed worldwide during the past half-century. Our numerical experiments suggest that the land-surface energy partitioning is strongly modulated by groundwater–river-water interactions through expanding the periodically inundated fraction of the riparian zone, and enhancing moisture availability in the vadose zone via capillary rise in response to the river stage change. Meanwhile, CLM4.5 fails to capture the key hydrologic process (i.e., groundwater–river-water exchange) at the site, and consequently simulates drastically different water and energy budgets. Furthermore, spatial resolution is found to significantly impact the accuracy of estimated the mass exchange rates at the boundaries of the aquifer, and it becomes critical when surface and subsurface become more tightly coupled with groundwater table within 6 to 7 meters below the surface. Inclusion of lateral subsurface flow influenced both the surface energy budget and subsurface transport processes as a result of river

Coupling a three-dimensional subsurface flow and transport model with a land surface model to simulate stream–aquifer–land interactions (CP v1.0)

DOE Office of Scientific and Technical Information (OSTI.GOV)

Bisht, Gautam; Huang, Maoyi; Zhou, Tian

A fully coupled three-dimensional surface and subsurface land model is developed and applied to a site along the Columbia River to simulate three-way interactions among river water, groundwater, and land surface processes. The model features the coupling of the Community Land Model version 4.5 (CLM4.5) and a massively parallel multiphysics reactive transport model (PFLOTRAN). The coupled model, named CP v1.0, is applied to a 400 m × 400 m study domain instrumented with groundwater monitoring wells along the Columbia River shoreline. CP v1.0 simulations are performed at three spatial resolutions (i.e., 2, 10, and 20 m) over a 5-year periodmore » to evaluate the impact of hydroclimatic conditions and spatial resolution on simulated variables. Results show that the coupled model is capable of simulating groundwater–river-water interactions driven by river stage variability along managed river reaches, which are of global significance as a result of over 30 000 dams constructed worldwide during the past half-century. Our numerical experiments suggest that the land-surface energy partitioning is strongly modulated by groundwater–river-water interactions through expanding the periodically inundated fraction of the riparian zone, and enhancing moisture availability in the vadose zone via capillary rise in response to the river stage change. Meanwhile, CLM4.5 fails to capture the key hydrologic process (i.e., groundwater–river-water exchange) at the site, and consequently simulates drastically different water and energy budgets. Furthermore, spatial resolution is found to significantly impact the accuracy of estimated the mass exchange rates at the boundaries of the aquifer, and it becomes critical when surface and subsurface become more tightly coupled with groundwater table within 6 to 7 meters below the surface. Inclusion of lateral subsurface flow influenced both the surface energy budget and subsurface transport processes as a result of river-water intrusion

Coupling a three-dimensional subsurface flow and transport model with a land surface model to simulate stream–aquifer–land interactions (CP v1.0)

DOE PAGES

Bisht, Gautam; Huang, Maoyi; Zhou, Tian; ...

2017-12-12

A fully coupled three-dimensional surface and subsurface land model is developed and applied to a site along the Columbia River to simulate three-way interactions among river water, groundwater, and land surface processes. The model features the coupling of the Community Land Model version 4.5 (CLM4.5) and a massively parallel multiphysics reactive transport model (PFLOTRAN). The coupled model, named CP v1.0, is applied to a 400 m × 400 m study domain instrumented with groundwater monitoring wells along the Columbia River shoreline. CP v1.0 simulations are performed at three spatial resolutions (i.e., 2, 10, and 20 m) over a 5-year periodmore » to evaluate the impact of hydroclimatic conditions and spatial resolution on simulated variables. Results show that the coupled model is capable of simulating groundwater–river-water interactions driven by river stage variability along managed river reaches, which are of global significance as a result of over 30 000 dams constructed worldwide during the past half-century. Our numerical experiments suggest that the land-surface energy partitioning is strongly modulated by groundwater–river-water interactions through expanding the periodically inundated fraction of the riparian zone, and enhancing moisture availability in the vadose zone via capillary rise in response to the river stage change. Meanwhile, CLM4.5 fails to capture the key hydrologic process (i.e., groundwater–river-water exchange) at the site, and consequently simulates drastically different water and energy budgets. Furthermore, spatial resolution is found to significantly impact the accuracy of estimated the mass exchange rates at the boundaries of the aquifer, and it becomes critical when surface and subsurface become more tightly coupled with groundwater table within 6 to 7 meters below the surface. Inclusion of lateral subsurface flow influenced both the surface energy budget and subsurface transport processes as a result of river-water intrusion

Induction heaters used to heat subsurface formations

DOEpatents

Nguyen, Scott Vinh [Houston, TX; Bass, Ronald M [Houston, TX

2012-04-24

A heating system for a subsurface formation includes an elongated electrical conductor located in the subsurface formation. The electrical conductor extends between at least a first electrical contact and a second electrical contact. A ferromagnetic conductor at least partially surrounds and at least partially extends lengthwise around the electrical conductor. The electrical conductor, when energized with time-varying electrical current, induces sufficient electrical current flow in the ferromagnetic conductor such that the ferromagnetic conductor resistively heats to a temperature of at least about 300.degree. C.

On the characterization of subsurface flow and hydraulic conductivity from surface SP measurements: correcting for electrical heterogeneities.

NASA Astrophysics Data System (ADS)

Sailhac, P.; Marquis, G.; Darnet, M.; Szalai, S.

2003-04-01

Surface self potential measurements (SP) are useful to characterize underground fluid flow or chemical reactions (as redox) and can be used in addition to NMR and electrical prospecting in hydrological investigations. Assuming that the SP anomalies have an electrokinetic origin, the source of SP data is the divergence of underground fluid flow; one important problem with surface SP data is then its interpretation in terms of fluid flow geometry. Some integral transform techniques have been shown to be powerful for SP interpretation (e.g. Fournier 1989, Patella, 1997; Sailhac &Marquis 2001). All these techniques are based upon Green’{ }s functions to characterize underground water flow, but they assume a constant electrical conductivity in the subsurface. This unrealistic approximation results in the appearance of non-electrokinetic sources at strong lateral electrical conductivity contrasts. We present here new Green’{ }s functions suitable for media of heterogeneous electrical conductivity. This new approach allows the joint interpretation of electrical resistivity tomography and SP measurements to detect electrokinetic sources caused by fluid flow. Tests on synthetic examples show that it gives more realistic results that when a constant electrical conductivity is assumed.

A field evaluation of subsurface and surface runoff. II. Runoff processes

USGS Publications Warehouse

Pilgrim, D.H.; Huff, D.D.; Steele, T.D.

1978-01-01

Combined use of radioisotope tracer, flow rate, specific conductance and suspended-sediment measurements on a large field plot near Stanford, California, has provided more detailed information on surface and subsurface storm runoff processes than would be possible from any single approach used in isolation. Although the plot was surficially uniform, the runoff processes were shown to be grossly nonuniform, both spatially over the plot, and laterally and vertically within the soil. The three types of processes that have been suggested as sources of storm runoff (Horton-type surface runoff, saturated overland flow, and rapid subsurface throughflow) all occurred on the plot. The nonuniformity of the processes supports the partial- and variable-source area concepts. Subsurface storm runoff occurred in a saturated layer above the subsoil horizon, and short travel times resulted from flow through macropores rather than the soil matrix. Consideration of these observations would be necessary for physically realistic modeling of the storm runoff process. ?? 1978.

Integrating non-colocated well and geophysical data to capture subsurface heterogeneity at an aquifer recharge and recovery site

NASA Astrophysics Data System (ADS)

Gottschalk, Ian P.; Hermans, Thomas; Knight, Rosemary; Caers, Jef; Cameron, David A.; Regnery, Julia; McCray, John E.

2017-12-01

Geophysical data have proven to be very useful for lithological characterization. However, quantitatively integrating the information gained from acquiring geophysical data generally requires colocated lithological and geophysical data for constructing a rock-physics relationship. In this contribution, the issue of integrating noncolocated geophysical and lithological data is addressed, and the results are applied to simulate groundwater flow in a heterogeneous aquifer in the Prairie Waters Project North Campus aquifer recharge site, Colorado. Two methods of constructing a rock-physics transform between electrical resistivity tomography (ERT) data and lithology measurements are assessed. In the first approach, a maximum likelihood estimation (MLE) is used to fit a bimodal lognormal distribution to horizontal crosssections of the ERT resistivity histogram. In the second approach, a spatial bootstrap is applied to approximate the rock-physics relationship. The rock-physics transforms provide soft data for multiple point statistics (MPS) simulations. Subsurface models are used to run groundwater flow and tracer test simulations. Each model's uncalibrated, predicted breakthrough time is evaluated based on its agreement with measured subsurface travel time values from infiltration basins to selected groundwater recovery wells. We find that incorporating geophysical information into uncalibrated flow models reduces the difference with observed values, as compared to flow models without geophysical information incorporated. The integration of geophysical data also narrows the variance of predicted tracer breakthrough times substantially. Accuracy is highest and variance is lowest in breakthrough predictions generated by the MLE-based rock-physics transform. Calibrating the ensemble of geophysically constrained models would help produce a suite of realistic flow models for predictive purposes at the site. We find that the success of breakthrough predictions is highly

Effect of subsurface heterogeneity on free-product recovery from unconfined aquifers

NASA Astrophysics Data System (ADS)

Kaluarachchi, Jagath J.

1996-03-01

Free-product record system designs for light-hydrocarbon-contaminated sites were investigated to evaluate the effects of subsurface heterogeneity using a vertically integrated three-phase flow model. The input stochastic variable of the areal flow analysis was the log-intrinsic permeability and it was generated using the Turning Band method. The results of a series of hypothetical field-scale simulations showed that subsurface heterogeneity has a substantial effect on free-product recovery predictions. As the heterogeneity increased, the recoverable oil volume decreased and the residual trapped oil volume increased. As the subsurface anisotropy increased, these effects together with free- and total-oil contaminated areas were further enhanced. The use of multiple-stage water pumping was found to be insignificant compared to steady uniform pumping due to reduced recovery efficiency and increased residual oil volume. This observation was opposite to that produced under homogeneous scenarios. The effect of subsurface heterogeneity was enhanced at relatively low water pumping rates. The difference in results produced by homogeneous and heterogeneous simulations was substantial, indicating greater attention should be paid in modeling free-product recovery systems with appropriate subsurface heterogeneity.

The Effectiveness of Organic Pollutants Removal in Constructed Wetland with Horizontal Sub-Surface Flow / Efektywność Usuwania Zanieczyszczeń Organicznych W Oczyszczalni Hydrofitowej

NASA Astrophysics Data System (ADS)

Jakubaszek, Anita; Sadecka, Zofia

2015-03-01

This paper presents the results of the research work related to the removal efficiency from wastewater organic pollutants and suspended solids at HSSF (horizontal subsurface flow) constructed wetland. The average effectiveness defined as loss of value COD in wastewater has reached 77%, for BOD5 - 80% and TOC - 82%. The effect of seasonal temperature changes and the period of plant vegetation and rest on the effectiveness of wastewater treatment were also analyzed. The results of the presented research showed a decrease in the efficiency of removing organic pollutants from wastewater and suspended solids in the autumn and winter. During the vegetation the object in Małyszyn has been characterized by the effectiveness of wastewater treatment at the level of 78% for COD, 82% for BOD5, and in the non-vegetation period the effectiveness has decreased up to 75% for COD and 74% for BOD5. During the plants growth the total suspension was removed in 88%, whereas during the plants rest efficiency of removing lowered to 69%. W pracy przedstawiono wyniki badań dotyczące efektywności usuwania ze ścieków zanieczyszczeń organicznych w oczyszczalni hydrofitowej. Średnia skuteczność oczyszczania wyrażona jako obniżenie wartości ChZT w ściekach była na poziomie 77%, dla BZT5 80%, a dla OWO 82%. Analizowano również wpływ sezonowych zmian temperatury oraz okresu wegetacji i spoczynku roślin na skuteczność oczyszczania ścieków. Wyniki badań wykazały obniżenie efektywności usuwania zanieczyszczeń organicznych ze ścieków wyrażonych przez ChZT i BZT5 oraz zawiesiny ogólnej w okresie jesienno-zimowym. W okresie wegetacyjnym obiekt w Małyszynie charakteryzował się efektywnością oczyszczania ścieków na poziomie: 78% dla ChZT, 82% dla BZT5, a w sezonie pozawegetacyjnym skuteczność uległa obniżeniu do 75% w przypadku ChZT oraz 74% dla BZT5. Zawiesina ogólna w okresie wegetacji trzciny usuwana była w 88%, a w okresie powegetacyjnym w 69%.

Coupled surface and subsurface flow modeling of natural hillslopes in the Aburrá Valley (Medellín, Colombia)

NASA Astrophysics Data System (ADS)

Blessent, Daniela; Barco, Janet; Temgoua, André Guy Tranquille; Echeverrri-Ramirez, Oscar

2017-03-01

Numerical results are presented of surface-subsurface water modeling of a natural hillslope located in the Aburrá Valley, in the city of Medellín (Antioquia, Colombia). The integrated finite-element hydrogeological simulator HydroGeoSphere is used to conduct transient variably saturated simulations. The objective is to analyze pore-water pressure and saturation variation at shallow depths, as well as volumes of water infiltrated in the porous medium. These aspects are important in the region of study, which is highly affected by soil movements, especially during the high-rain seasons that occur twice a year. The modeling exercise considers rainfall events that occurred between October and December 2014 and a hillslope that is currently monitored because of soil instability problems. Simulation results show that rainfall temporal variability, mesh resolution, coupling length, and the conceptual model chosen to represent the heterogeneous soil, have a noticeable influence on results, particularly for high rainfall intensities. Results also indicate that surface-subsurface coupled modeling is required to avoid unrealistic increase in hydraulic heads when high rainfall intensities cause top-down saturation of soil. This work is a first effort towards fostering hydrogeological modeling expertise that may support the development of monitoring systems and early landslide warning in a country where the rainy season is often the cause of hydrogeological tragedies associated with landslides, mud flow or debris flow.

Removal and factors influencing removal of sulfonamides and trimethoprim from domestic sewage in constructed wetlands.

PubMed

Dan A; Yang, Yang; Dai, Yu-Nv; Chen, Chun-Xing; Wang, Su-Yu; Tao, Ran

2013-10-01

Twelve pilot-scale constructed wetlands with different configurations were set up in the field to evaluate the removal and factors that influence removal of sulfonamides (sulfadiazine, sulfapyridine, sulfacetamide, sulfamethazine and sulfamethoxazole) and trimethoprim from domestic sewage. The treatments included four flow types, three substrates, two plants and three hydraulic loading rates across two seasons (summer and winter). Most target antibiotics were efficiently removed by specific constructed wetlands; in particular, all types of constructed wetlands performed well for the degradation of sulfapyridine. Flow types were the most important influencing factor in this study, and the best removal of sulfonamides was achieved in vertical subsurface-flow constructed wetlands; however, the opposite phenomenon was found with trimethoprim. Significant relationships were observed between antibiotic degradation and higher temperature and redox potential, which indicated that microbiological pathways were the most probable degradation route for sulfonamides and trimethoprim in constructed wetlands. Copyright © 2013 Elsevier Ltd. All rights reserved.

The Lakshmi Plateau structure as an indicator of asthenosphere horizontal flows on Venus

NASA Technical Reports Server (NTRS)

Pronin, A. A.

1986-01-01

The structure of Lakshmi Planum in the western part of Ishtar Terra in a fold-fault setting which conforms to the basic massif of the plateau with eruptive centers is constructed concentrically and is interpreted from the point of view of the subsurface flow of materials in the form of horizontally diverging asthenospheric flows and gravitational creep. The surrounding structures are formed by the deformation of the more rigid lithosphere as it breaks away from the asthenospheric flow.

Flow diagram analysis of electrical fatalities in construction industry.

PubMed

Chi, Chia-Fen; Lin, Yuan-Yuan; Ikhwan, Mohamad

2012-01-01

The current study reanalyzed 250 electrical fatalities in the construction industry from 1996 to 2002 into seven patterns based on source of electricity (power line, energized equipment, improperly installed or damaged equipment), direct contact or indirect contact through some source of injury (boom vehicle, metal bar or pipe, and other conductive material). Each fatality was coded in terms of age, company size, experience, performing tasks, source of injury, accident cause and hazard pattern. The Chi-square Automatic Interaction Detector (CHAID) was applied to the coded data of the fatal electrocution to find a subset of predictors that might derive meaningful classifications or accidents scenarios. A series of Flow Diagrams was constructed based on CHAID result to illustrate the flow of electricity travelling from electrical source to human body. Each of the flow diagrams can be directly linked with feasible prevention strategies by cutting the flow of electricity.

Subsurface architecture of Las Bombas volcano circular structure (Southern Mendoza, Argentina) from geophysical studies

NASA Astrophysics Data System (ADS)

Prezzi, Claudia; Risso, Corina; Orgeira, María Julia; Nullo, Francisco; Sigismondi, Mario E.; Margonari, Liliana

2017-08-01

The Plio-Pleistocene Llancanelo volcanic field is located in the south-eastern region of the province of Mendoza, Argentina. This wide back-arc lava plateau, with hundreds of monogenetic pyroclastic cones, covers a large area behind the active Andean volcanic arc. Here we focus on the northern Llancanelo volcanic field, particularly in Las Bombas volcano. Las Bombas volcano is an eroded, but still recognizable, scoria cone located in a circular depression surrounded by a basaltic lava flow, suggesting that Las Bombas volcano was there when the lava flow field formed and, therefore, the lava flow engulfed it completely. While this explanation seems reasonable, the common presence of similar landforms in this part of the field justifies the need to establish correctly the stratigraphic relationship between lava flow fields and these circular depressions. The main purpose of this research is to investigate Las Bombas volcano 3D subsurface architecture by means of geophysical methods. We carried out a paleomagnetic study and detailed topographic, magnetic and gravimetric land surveys. Magnetic anomalies of normal and reverse polarity and paleomagnetic results point to the occurrence of two different volcanic episodes. A circular low Bouguer anomaly was detected beneath Las Bombas scoria cone indicating the existence of a mass deficit. A 3D forward gravity model was constructed, which suggests that the mass deficit would be related to the presence of fracture zones below Las Bombas volcano cone, due to sudden degassing of younger magma beneath it, or to a single phreatomagmatic explosion. Our results provide new and detailed information about Las Bombas volcano subsurface architecture.

The removal efficiency of constructed wetlands filled with the zeolite-slag hybrid substrate for the rural landfill leachate treatment.

PubMed

He, Hailing; Duan, Zhiwei; Wang, Zhenqing; Yue, Bo

2017-07-01

The removal efficiencies of two horizontal subsurface flow constructed wetlands (HSSF CWs, down-flow (F1) and up-flow (F2)) filled with the zeolite-slag hybrid substrate for the rural landfill leachate treatment were investigated. The adsorption experiment was conducted to evaluate the potential of zeolite and slag as the wetland substrate. The effects of distance variations along the longitudinal profile of wetland bed on pollutant removal were assessed by sampling at four locations (inlet, outlet, 0.55 m, and 1.10 m from the inlet). During the operation time, the influent and effluent concentrations of chemical oxygen demand (COD), ammonia nitrogen (NH 3 -N), total nitrogen (TN), heavy metals, and polycyclic aromatic hydrocarbon (PAH) were measured. The results showed that the constructed wetlands were capable of removing COD, 20.5-48.2% (F1) and 18.6-61.2% (F2); NH 3 -N, 84.0-99.9% (F1) and 93.5-99.2% (F2); TN, 80.3-92.1% (F1) and 80.3-91.2% (F2); and heavy metals, about 90% (F1 and F2). The zeolite-slag hybrid substrate performed excellent removal efficiency for the nitrogen and heavy metals. The inlet area was the most active region of leachate removal. The up-flow constructed wetland (F2) has a higher removal efficiency for the PAH compounds. The significant removal efficiency illustrated that the rural landfill leachate can be treated using the horizontal subsurface flow constructed wetland filled with the zeolite-slag hybrid substrate.

Effects of interspecific competition on the growth of macrophytes and nutrient removal in constructed wetlands: A comparative assessment of free water surface and horizontal subsurface flow systems.

PubMed

Zheng, Yucong; Wang, Xiaochang; Dzakpasu, Mawuli; Zhao, Yaqian; Ngo, Huu Hao; Guo, Wenshan; Ge, Yuan; Xiong, Jiaqing

2016-05-01

The outcome of competition between adjoining interspecific colonies of Phragmites and Typha in two large field pilot-scale free water surface (FWS) and subsurface flow (SSF) CWs is evaluated. According to findings, the effect of interspecific competition was notable for Phragmites australis, whereby it showed the highest growth performance in both FWS and SSF wetland. In a mixed-culture, P. australis demonstrates superiority in terms of competitive interactions for space between plants. Furthermore, the interspecific competition among planted species seemed to cause different ecological responses of plant species in the two CWs. For example, while relatively high density and shoot height determined the high aboveground dry weight of P. australis in the FWS wetland, this association was not evident in the SSF. Additionally, while plants nutrients uptake accounts for a higher proportion of the nitrogen removal in FWS, that in the SSF accounts for a higher proportion of the phosphorous removal. Copyright © 2016 Elsevier Ltd. All rights reserved.

Construction and geometric stability of physiological flow rate wall-less stenosis phantoms.

PubMed

Ramnarine, K V; Anderson, T; Hoskins, P R

2001-02-01

Wall-less flow phantoms are preferred for ultrasound (US) because tissue-mimicking material (TMM) with good acoustical properties can be made and cast to form anatomical models. The construction and geometrical stability of wall-less TMM flow phantoms is described using a novel method of sealing to prevent leakage of the blood-mimicking fluid (BMF). Wall-less stenosis flow models were constructed using a robust agar-based TMM and sealed using reticulated foam at the inlet and outlet tubes. There was no BMF leakage at the highest flow rate of 2.8 L/min in 0%, 35% and 57% diameter reduction stenoses models. Failure of the 75% stenosis model, due to TMM fracture, occurred at maximum flow rate of 2 L/min (mean velocity 10 m/s within the stenosis). No change of stenosis geometry was measured over 4 days. The construction is simple and effective and extends the possibility for high flow rate studies using robust TMM wall-less phantoms.

Crushed cement concrete substitution for construction aggregates; a materials flow analysis

USGS Publications Warehouse

Kelly, Thomas

1998-01-01

An analysis of the substitution of crushed cement concrete for natural construction aggregates is performed by using a materials flow diagram that tracks all material flows into and out of the cement concrete portion of the products made with cement concrete: highways, roads, and buildings. Crushed cement concrete is only one of the materials flowing into these products, and the amount of crushed cement concrete substituted influences the amount of other materials in the flow. Factors such as availability and transportation costs, as well as physical properties, that can affect stability and finishability, influence whether crushed cement concrete or construction aggregates should be used or predominate for a particular end use.

Transport of Chemical Vapors from Subsurface Sources to Atmosphere as Affected by Shallow Subsurface and Atmospheric Conditions

NASA Astrophysics Data System (ADS)

Rice, A. K.; Smits, K. M.; Hosken, K.; Schulte, P.; Illangasekare, T. H.

2012-12-01

Understanding the movement and modeling of chemical vapor through unsaturated soil in the shallow subsurface when subjected to natural atmospheric thermal and mass flux boundary conditions at the land surface is of importance to applications such as landmine detection and vapor intrusion into subsurface structures. New, advanced technologies exist to sense chemical signatures at the land/atmosphere interface, but interpretation of these sensor signals to make assessment of source conditions remains a challenge. Chemical signatures are subject to numerous interactions while migrating through the unsaturated soil environment, attenuating signal strength and masking contaminant source conditions. The dominant process governing movement of gases through porous media is often assumed to be Fickian diffusion through the air phase with minimal or no quantification of other processes contributing to vapor migration, such as thermal diffusion, convective gas flow due to the displacement of air, expansion/contraction of air due to temperature changes, temporal and spatial variations of soil moisture and fluctuations in atmospheric pressure. Soil water evaporation and interfacial mass transfer add to the complexity of the system. The goal of this work is to perform controlled experiments under transient conditions of soil moisture, temperature and wind at the land/atmosphere interface and use the resulting dataset to test existing theories on subsurface gas flow and iterate between numerical modeling efforts and experimental data. Ultimately, we aim to update conceptual models of shallow subsurface vapor transport to include conditionally significant transport processes and inform placement of mobile sensors and/or networks. We have developed a two-dimensional tank apparatus equipped with a network of sensors and a flow-through head space for simulation of the atmospheric interface. A detailed matrix of realistic atmospheric boundary conditions was applied in a series of

Apparatus for passive removal of subsurface contaminants

DOEpatents

Pemberton, Bradley E.; May, Christopher P.; Rossabi, Joseph

1997-01-01

An apparatus is provided which passively removes contaminated gases from a subsurface. The apparatus includes a riser pipe extending into a subsurface which has an exterior end in fluid communication with a valve. When well pressure is greater than atmospheric pressure, the valve opens to release contaminants into the atmosphere, and when well pressure is less than atmospheric pressure, the valve closes to prevent flow of air into the well. The valve assembly of the invention comprises a lightweight ball which is lifted from its valve seat with a slight pressure drop between the well and the atmosphere.

Apparatus for passive removal of subsurface contaminants

DOEpatents

Pemberton, B.E.; May, C.P.; Rossabi, J.

1997-06-24

An apparatus is provided which passively removes contaminated gases from a subsurface. The apparatus includes a riser pipe extending into a subsurface which has an exterior end in fluid communication with a valve. When well pressure is greater than atmospheric pressure, the valve opens to release contaminants into the atmosphere, and when well pressure is less than atmospheric pressure, the valve closes to prevent flow of air into the well. The valve assembly of the invention comprises a lightweight ball which is lifted from its valve seat with a slight pressure drop between the well and the atmosphere. 7 figs.

On the effects of subsurface parameters on evaporite dissolution (Switzerland)

NASA Astrophysics Data System (ADS)

Zidane, Ali; Zechner, Eric; Huggenberger, Peter; Younes, Anis

2014-05-01

Uncontrolled subsurface evaporite dissolution could lead to hazards such as land subsidence. Observed subsidences in a study area of Northwestern Switzerland were mainly due to subsurface dissolution (subrosion) of evaporites such as halite and gypsum. A set of 2D density driven flow simulations were evaluated along 1000 m long and 150 m deep 2D cross sections within the study area that is characterized by tectonic horst and graben structures. The simulations were conducted to study the effect of the different subsurface parameters that could affect the dissolution process. The heterogeneity of normal faults and its impact on the dissolution of evaporites is studied by considering several permeable faults that include non-permeable areas. The mixed finite element method (MFE) is used to solve the flow equation, coupled with the multipoint flux approximation (MPFA) and the discontinuous Galerkin method (DG) to solve the diffusion and the advection parts of the transport equation.

Nested sparse grid collocation method with delay and transformation for subsurface flow and transport problems

NASA Astrophysics Data System (ADS)

Liao, Qinzhuo; Zhang, Dongxiao; Tchelepi, Hamdi

2017-06-01

In numerical modeling of subsurface flow and transport problems, formation properties may not be deterministically characterized, which leads to uncertainty in simulation results. In this study, we propose a sparse grid collocation method, which adopts nested quadrature rules with delay and transformation to quantify the uncertainty of model solutions. We show that the nested Kronrod-Patterson-Hermite quadrature is more efficient than the unnested Gauss-Hermite quadrature. We compare the convergence rates of various quadrature rules including the domain truncation and domain mapping approaches. To further improve accuracy and efficiency, we present a delayed process in selecting quadrature nodes and a transformed process for approximating unsmooth or discontinuous solutions. The proposed method is tested by an analytical function and in one-dimensional single-phase and two-phase flow problems with different spatial variances and correlation lengths. An additional example is given to demonstrate its applicability to three-dimensional black-oil models. It is found from these examples that the proposed method provides a promising approach for obtaining satisfactory estimation of the solution statistics and is much more efficient than the Monte-Carlo simulations.

Debates - Stochastic subsurface hydrology from theory to practice: Introduction

NASA Astrophysics Data System (ADS)

Rajaram, Harihar

2016-12-01

This paper introduces the papers in the "Debates - Stochastic Subsurface Hydrology from Theory to Practice" series. Beginning in the 1970s, the field of stochastic subsurface hydrology has been an active field of research, with over 3500 journal publications, of which over 850 have appeared in Water Resources Research. We are fortunate to have insightful contributions from four groups of distinguished authors who discuss the reasons why the advanced research framework established in stochastic subsurface hydrology has not impacted the practice of groundwater flow and transport modeling and design significantly. There is reasonable consensus that a community effort aimed at developing "toolboxes" for applications of stochastic methods will make them more accessible and encourage practical applications.

Data on subsurface storage of liquid waste near Pensacola, Florida, 1963-1980

USGS Publications Warehouse

Hull, R.W.; Martin, J.B.

1982-01-01

Since 1963, when industrial waste was first injected into the subsurface in northwest Florida, considerable data have been collected relating to the geochemistry of subsurface waste storage. This report presents hydrogeologic data on two subsurface waste storage. This report presents hydrogeologic data on two subsurface storage systems near Pensacola, Fla., which inject liquid industrial waste through deep wells into a saline aquifer. Injection sites are described giving a history of well construction, injection, and testing; geologic data from cores and grab samples; hydrographs of injection rates, volume, pressure, and water levels; and chemical and physical data from water-quality samples collected from injection and monitor wells. (USGS)

Translucent Radiosity: Efficiently Combining Diffuse Inter-Reflection and Subsurface Scattering.

PubMed

Sheng, Yu; Shi, Yulong; Wang, Lili; Narasimhan, Srinivasa G

2014-07-01

It is hard to efficiently model the light transport in scenes with translucent objects for interactive applications. The inter-reflection between objects and their environments and the subsurface scattering through the materials intertwine to produce visual effects like color bleeding, light glows, and soft shading. Monte-Carlo based approaches have demonstrated impressive results but are computationally expensive, and faster approaches model either only inter-reflection or only subsurface scattering. In this paper, we present a simple analytic model that combines diffuse inter-reflection and isotropic subsurface scattering. Our approach extends the classical work in radiosity by including a subsurface scattering matrix that operates in conjunction with the traditional form factor matrix. This subsurface scattering matrix can be constructed using analytic, measurement-based or simulation-based models and can capture both homogeneous and heterogeneous translucencies. Using a fast iterative solution to radiosity, we demonstrate scene relighting and dynamically varying object translucencies at near interactive rates.

Subsurface Biodegradation in a Fractured Basement Reservoir, Shropshire, UK

NASA Astrophysics Data System (ADS)

Parnell, John; Baba, Mas'ud; Bowden, Stephen; Muirhead, David

2017-04-01

Subsurface Biodegradation in a Fractured Basement Reservoir, Shropshire, UK. John Parnell, Mas'ud Baba, Stephen Bowden, David Muirhead Subsurface biodegradation in current oil reservoirs is well established, but there are few examples of fossil subsurface degradation. Biomarker compositions of viscous and solid oil residues ('bitumen') in fractured Precambrian and other basement rocks below the Carboniferous cover in Shropshire, UK, show that they are variably biodegraded. High levels of 25-norhopanes imply that degradation occurred in the subsurface. Lower levels of 25-norhopanes occur in active seepages. Liquid oil trapped in fluid inclusions in mineral veins in the fractured basement confirm that the oil was emplaced fresh before subsurface degradation. A Triassic age for the veins implies a 200 million year history of hydrocarbon migration in the basement rocks. The data record microbial colonization of a fractured basement reservoir, and add to evidence in modern basement aquifers for microbial activity in deep fracture systems. Buried basement highs may be especially favourable to colonization, through channelling fluid flow to shallow depths and relatively low temperatures

Stochastic constructions of flows of rank 1

NASA Astrophysics Data System (ADS)

Prikhod'ko, A. A.

2001-12-01

Automorphisms of rank 1 appeared in the well-known papers of Chacon (1965), who constructed an example of a weakly mixing automorphism not having the strong mixing property, and Ornstein (1970), who proved the existence of mixing automorphisms without a square root. Ornstein's construction is essentially stochastic, since its parameters are chosen in a "sufficiently random manner" according to a certain random law.In the present article it is shown that mixing flows of rank 1 exist. The construction given is also stochastic and is based to a large extent on ideas in Ornstein's paper. At the same time it complements Ornstein's paper and makes it more transparent. The construction can be used also to obtain automorphisms with various approximation and statistical properties. It is established that the new examples of dynamical systems are not isomorphic to Ornstein automorphisms, that is, they are qualitatively new.

Stochastic constructions of flows of rank 1

DOE Office of Scientific and Technical Information (OSTI.GOV)

Prikhod'ko, A A

2001-12-31

Automorphisms of rank 1 appeared in the well-known papers of Chacon (1965), who constructed an example of a weakly mixing automorphism not having the strong mixing property, and Ornstein (1970), who proved the existence of mixing automorphisms without a square root. Ornstein's construction is essentially stochastic, since its parameters are chosen in a 'sufficiently random manner' according to a certain random law. In the present article it is shown that mixing flows of rank 1 exist. The construction given is also stochastic and is based to a large extent on ideas in Ornstein's paper. At the same time it complementsmore » Ornstein's paper and makes it more transparent. The construction can be used also to obtain automorphisms with various approximation and statistical properties. It is established that the new examples of dynamical systems are not isomorphic to Ornstein automorphisms, that is, they are qualitatively new.« less

A Method for Partitioning Surface and Subsurface Flow Using Rainfall Simulaton and Two-Dimensional Surface Electrical Resistivity Imaging

NASA Astrophysics Data System (ADS)

Carey, A. M.; Paige, G. B.; Miller, S. N.; Carr, B. J.; Holbrook, W. S.

2014-12-01

In semi-arid rangeland environments understanding how surface and subsurface flow processes and their interactions are influenced by watershed and rainfall characteristics is critical. However, it is difficult to resolve the temporal variations between mechanisms controlling these processes and challenging to obtain field measurements that document their interactions. Better insight into how these complex systems respond hydrologically is necessary in order to refine hydrologic models and decision support tools. We are conducting field studies integrating high resolution, two-dimensional surface electrical resistivity imaging (ERI) with variable intensity rainfall simulation, to quantify real-time partitioning of rainfall into surface and subsurface response. These studies are being conducted at the hillslope scale on long-term runoff plots on four different ecological sites in the Upper Crow Creek Watershed in southeastern Wyoming. Variable intensity rainfall rates were applied using the Walnut Gulch Rainfall Simulator in which intensities were increased incrementally from 49 to 180 mm hr-1 and steady-state runoff rates for each intensity were measured. Two 13.5 m electrode arrays at 0.5 m spacing were positioned on the surface perpendicular to each plot and potentials were measured at given time intervals prior to, during and following simulations using a dipole-dipole array configuration. The configuration allows for a 2.47 m depth of investigation in which magnitude and direction of subsurface flux can be determined. We used the calculated steady state infiltration rates to quantify the variability in the partial area runoff response on the ecological sites. Coupling this information with time-lapse difference inversions of ERI data, we are able to track areas of increasing and decreasing resistivity in the subsurface related to localized areas of infiltration during and following rainfall events. We anticipate implementing this method across a variety of

Nutrient transport through a Vegetative Filter Strip with subsurface drainage.

PubMed

Bhattarai, Rabin; Kalita, Prasanta Kumar; Patel, Mita Kanu

2009-04-01

The transport of nutrients and soil sediments in runoff has been recognized as a noteworthy environmental issue. Vegetative Filter Strips (VFS) have been used as one of the best management practices (BMPs) for retaining nutrients and sediments from surface runoff, thus preventing the pollutants from reaching receiving waters. However, the effectiveness of a VFS when combined with a subsurface drainage system has not been investigated previously. This study was undertaken to monitor the retention and transport of nutrients within a VFS that had a subsurface drainage system installed at a depth of 1.2 m below the soil surface. Nutrient concentrations of NO(3)-N (Nitrate Nitrogen), PO(-)(4) (Orthophosphorus), and TP (Total Phosphorus) were measured in surface water samples (entering and leaving the VFS), and subsurface outflow. Soil samples were collected and analyzed for plant available Phosphorus (Bray P1) and NO(3)-N concentrations. Results showed that PO(-)(4), NO(3)-N, and TP concentrations decreased in surface flow through the VFS. Many surface outflow water samples from the VFS showed concentration reductions of as much as 75% for PO(-)(4) and 70% for TP. For subsurface outflow water samples through the drainage system, concentrations of PO(-)(4) and TP decreased but NO(3)-N concentrations increased in comparison to concentrations in surface inflow samples. Soil samples that were collected from various depths in the VFS showed a minimal buildup of nutrients in the top soil profile but indicated a gradual buildup of nutrients at the depth of the subsurface drain. Results demonstrate that although a VFS can be very effective in reducing runoff and nutrients from surface flow, the presence of a subsurface drain underneath the VFS may not be environmentally beneficial. Such a combination may increase NO(3)-N transport from the VFS, thus invalidating the purpose of the BMP.

Recirculation or artificial aeration in vertical flow constructed wetlands: a comparative study for treating high load wastewater.

PubMed

Foladori, Paola; Ruaben, Jenny; Ortigara, Angela R C

2013-12-01

Vertical subsurface-flow constructed wetlands at pilot-scale have been applied to treat high hydraulic and organic loads by implementing the following configurations: (1) intermittent recirculation of the treated wastewater from the bottom to the top of the bed, (2) intermittent artificial aeration supplied at the bottom of the bed and (3) the combination of both. These configurations were operated with a saturated bottom layer for a 6h-treatment phase, followed by a free drainage phase prior to a new feeding. COD removal efficiency was 85-90% in all the configurations and removed loads were 54-70 gCOD m(-2)d(-1). The aerated and recirculated wetland resulted in a higher total nitrogen removal (8.6 gN m(-2)d(-1)) due to simultaneous nitrification/denitrification, even in the presence of intermittent aeration (6.8 Nm(3)m(-2)d(-1)). The extra investment needed for implementing aeration/recirculation would be compensated for by a reduction of the surface area per population equivalent, which decreased to 1.5m(2)/PE. Copyright © 2013 Elsevier Ltd. All rights reserved.

Characterizing subsurface water flow to artificial drain lines using fiber-optic distributed temperature sensing

NASA Astrophysics Data System (ADS)

Shults, D.; Brooks, E. S.; Heinse, R.; Keller, C. K.

2017-12-01

Over the last several years growers have experienced increasingly wet spring conditions in the Palouse Region located in North Idaho, Eastern Washington and Eastern Oregon. As a result more artificial drain lines are being installed so growers can access their fields earlier in the growing season. Additionally there has been increasing adoption of no-tillage practices among growers in order minimize erosion and runoff in the region. There is a growing body of evidence that suggests long-term no-tillage may lead to the establishment of large macropore networks through increased earthworm activity and the preservation of root channels. These macropore networks, in conjunctions with the presence of artificial drains lines, may create connected preferential flow paths from agricultural fields to receiving streams. This connectivity of flow paths from agricultural fields to receiving water bodies may increase the loading of nutrients and agricultural chemicals as some flow paths may largely bypass soil matrix interaction where materials can be sequestered. Our primary objective for this study was to characterize subsurface flow to two artificial drain lines, one under conventional tillage and the other under no-tillage, using distributed temperature sensing (DTS) technology. During the study (November 2016-April 2017) the near surface soil-water temperature was consistently colder than that of deeper depths. Temperature was thus used as a tracer as snow melt and soil-water moved from the near surface to the drain lines during snowmelt and precipitation events. The spatial and temporal variability of the temperature along the artificial drain line under no-tillage practices was found to be greater than that of the conventional tilled field. It is hypothesized that preferential flow paths are responsible for the increased variability of temperature seen in the drain line under long term no-till management. The temperature along the conventional till drain line showed a

Modeling Phosphorus Losses through Surface Runoff and Subsurface Drainage Using ICECREAM.

PubMed

Qi, Hongkai; Qi, Zhiming; Zhang, T Q; Tan, C S; Sadhukhan, Debasis

2018-03-01

Modeling soil phosphorus (P) losses by surface and subsurface flow pathways is essential in developing successful strategies for P pollution control. We used the ICECREAM model to simultaneously simulate P losses in surface and subsurface flow, as well as to assess effectiveness of field practices in reducing P losses. Monitoring data from a mineral-P-fertilized clay loam field in southwestern Ontario, Canada, were used for calibration and validation. After careful adjustment of model parameters, ICECREAM was shown to satisfactorily simulate all major processes of surface and subsurface P losses. When the calibrated model was used to assess tillage and fertilizer management scenarios, results point to a 10% reduction in total P losses by shifting autumn tillage to spring, and a 25.4% reduction in total P losses by injecting fertilizer rather than broadcasting. Although the ICECREAM model was effective in simulating surface and subsurface P losses when thoroughly calibrated, further testing is needed to confirm these results with manure P application. As illustrated here, successful use of simulation models requires careful verification of model routines and comprehensive calibration to ensure that site-specific processes are accurately represented. Copyright © by the American Society of Agronomy, Crop Science Society of America, and Soil Science Society of America, Inc.

EXPOSURE ASSESSMENT MODELING FOR HYDROCARBON SPILLS INTO THE SUBSURFACE

EPA Science Inventory

Hydrocarbons which enter the subsurface through spills or leaks may create serious, long-lived ground-water contamination problems. onventional finite difference and finite element models of multiphase, multicomponent flow often have extreme requirements for both computer time an...

Integrated surface/subsurface permafrost thermal hydrology: Model formulation and proof-of-concept simulations

DOE PAGES

Painter, Scott L.; Coon, Ethan T.; Atchley, Adam L.; ...

2016-08-11

The need to understand potential climate impacts and feedbacks in Arctic regions has prompted recent interest in modeling of permafrost dynamics in a warming climate. A new fine-scale integrated surface/subsurface thermal hydrology modeling capability is described and demonstrated in proof-of-concept simulations. The new modeling capability combines a surface energy balance model with recently developed three-dimensional subsurface thermal hydrology models and new models for nonisothermal surface water flows and snow distribution in the microtopography. Surface water flows are modeled using the diffusion wave equation extended to include energy transport and phase change of ponded water. Variation of snow depth in themore » microtopography, physically the result of wind scour, is also modeled heuristically with a diffusion wave equation. The multiple surface and subsurface processes are implemented by leveraging highly parallel community software. Fully integrated thermal hydrology simulations on the tilted open book catchment, an important test case for integrated surface/subsurface flow modeling, are presented. Fine-scale 100-year projections of the integrated permafrost thermal hydrological system on an ice wedge polygon at Barrow Alaska in a warming climate are also presented. Finally, these simulations demonstrate the feasibility of microtopography-resolving, process-rich simulations as a tool to help understand possible future evolution of the carbon-rich Arctic tundra in a warming climate.« less

Triaxial thermopile array geo-heat-flow sensor

DOEpatents

Carrigan, C.R.; Hardee, H.C.; Reynolds, G.D.; Steinfort, T.D.

1990-01-01

A triaxial thermopile array geothermal heat flow sensor is designed to measure heat flow in three dimensions in a reconstituted or unperturbed subsurface regime. Heat flow can be measured in conductive or permeable convective media. The sensor may be encased in protective pvc tubing and includes a plurality of thermistors and an array of heat flow transducers produce voltage proportional to heat flux along the subsurface regime and permit direct measurement of heat flow in the subsurface regime. The presence of the thermistor array permits a comparison to be made between the heat flow estimates obtained from the transducers and heat flow calculated using temperature differences and Fourier's Law. The device is extremely sensitive with an accuracy of less than 0.1 Heat Flow Units (HFU) and may be used for long term readings. 6 figs.

Triaxial thermopile array geo-heat-flow sensor

DOEpatents

Carrigan, Charles R.; Hardee, Harry C.; Reynolds, Gerald D.; Steinfort, Terry D.

1992-01-01

A triaxial thermopile array geothermal heat flow sensor is designed to measure heat flow in three dimensions in a reconstituted or unperturbed subsurface regime. Heat flow can be measured in conductive or permeable convective media. The sensor may be encased in protective pvc tubing and includes a plurality of thermistors and an array of heat flow transducers arranged in a vertical string. The transducers produce voltage proportional to heat flux along the subsurface regime and permit direct measurement of heat flow in the subsurface regime. The presence of the thermistor array permits a comparison to be made between the heat flow estimates obtained from the transducers and heat flow calculated using temperature differences and Fourier's Law. The device is extremely sensitive with an accuracy of less than 0.1 Heat Flow Units (HFU) and may be used for long term readings.

The influence of subsurface hydrodynamics on convective precipitation

NASA Astrophysics Data System (ADS)

Rahman, A. S. M. M.; Sulis, M.; Kollet, S. J.

2014-12-01

The terrestrial hydrological cycle comprises complex processes in the subsurface, land surface, and atmosphere, which are connected via complex non-linear feedback mechanisms. The influence of subsurface hydrodynamics on land surface mass and energy fluxes has been the subject of previous studies. Several studies have also investigated the soil moisture-precipitation feedback, neglecting however the connection with groundwater dynamics. The objective of this study is to examine the impact of subsurface hydrodynamics on convective precipitation events via shallow soil moisture and land surface processes. A scale-consistent Terrestrial System Modeling Platform (TerrSysMP) that consists of an atmospheric model (COSMO), a land surface model (CLM), and a three-dimensional variably saturated groundwater-surface water flow model (ParFlow), is used to simulate hourly mass and energy fluxes over days with convective rainfall events over the Rur catchment, Germany. In order to isolate the effect of groundwater dynamics on convective precipitation, two different model configurations with identical initial conditions are considered. The first configuration allows the groundwater table to evolve through time, while a spatially distributed, temporally constant groundwater table is prescribed as a lower boundary condition in the second configuration. The simulation results suggest that groundwater dynamics influence land surface soil moisture, which in turn affects the atmospheric boundary layer (ABL) height by modifying atmospheric thermals. It is demonstrated that because of this sensitivity of ABL height to soil moisture-temperature feedback, the onset and magnitude of convective precipitation is influenced by subsurface hydrodynamics. Thus, the results provide insight into the soil moisture-precipitation feedback including groundwater dynamics in a physically consistent manner by closing the water cycle from aquifers to the atmosphere.

Quantification of the effect of temperature gradients in soils on subsurface radon signal

NASA Astrophysics Data System (ADS)

Haquin, Gustavo; Ilzycer, Danielle; Kamai, Tamir; Zafrir, Hovav; Weisbrod, Noam

2017-04-01

Temperature gradients that develop in soils due to atmospheric temperature cycles are factors of primary importance in determining the rates and directions of subsurface gas flow. Models including mechanisms of thermal convection and thermal diffusion partially explain the impact of temperature gradients on subsurface radon transport. However, the overall impact of temperature gradients on subsurface radon transport is still not well understood. A laboratory setup was designed and built to experimentally investigate the influence of temperature gradients on radon transport under well controlled conditions. A 60 cm diameter and 120 cm tall column was thermally insulated except from the atmosphere-soil interface, such that it was constructed to simulate field conditions where temperature gradients in soils are developed following atmospheric temperature cycles. The column was filled with fine grinded phosphate rock which provided the porous media with radon source. Radon in soil-air was continuously monitored using NaI gamma detectors positioned at different heights along the column. Soil temperature, differential pressure, and relative humidity were monitored along the column. Experiments based on steep and gradual stepwise changes in ambient temperature were conducted. Absolute changes on radon levels in the order of 10-30% were measured at temperature gradients of up to ±20oC/m. Results showed a non-linear correlation between the temperature gradient and the subsurface radon concentration. An asymmetric relationship between the radon concentration and the temperature gradients for ΔT>0 and ΔT<0 was also observed. Laboratory simulations of the time- and depth-dependent temperature wave functions with frequencies ranged from a daily cycle to few days were performed. In response to the harmonic temperature behaviour radon oscillations at similar frequencies were detected correspondingly. In this work a quantitative relationship between radon and temperature

Understanding leachate flow in municipal solid waste landfills by combining time-lapse ERT and subsurface flow modelling - Part II: Constraint methodology of hydrodynamic models.

PubMed

Audebert, M; Oxarango, L; Duquennoi, C; Touze-Foltz, N; Forquet, N; Clément, R

2016-09-01

Leachate recirculation is a key process in the operation of municipal solid waste landfills as bioreactors. To ensure optimal water content distribution, bioreactor operators need tools to design leachate injection systems. Prediction of leachate flow by subsurface flow modelling could provide useful information for the design of such systems. However, hydrodynamic models require additional data to constrain them and to assess hydrodynamic parameters. Electrical resistivity tomography (ERT) is a suitable method to study leachate infiltration at the landfill scale. It can provide spatially distributed information which is useful for constraining hydrodynamic models. However, this geophysical method does not allow ERT users to directly measure water content in waste. The MICS (multiple inversions and clustering strategy) methodology was proposed to delineate the infiltration area precisely during time-lapse ERT survey in order to avoid the use of empirical petrophysical relationships, which are not adapted to a heterogeneous medium such as waste. The infiltration shapes and hydrodynamic information extracted with MICS were used to constrain hydrodynamic models in assessing parameters. The constraint methodology developed in this paper was tested on two hydrodynamic models: an equilibrium model where, flow within the waste medium is estimated using a single continuum approach and a non-equilibrium model where flow is estimated using a dual continuum approach. The latter represents leachate flows into fractures. Finally, this methodology provides insight to identify the advantages and limitations of hydrodynamic models. Furthermore, we suggest an explanation for the large volume detected by MICS when a small volume of leachate is injected. Copyright © 2016 Elsevier Ltd. All rights reserved.

Textural evidence for jamming and dewatering of a sub-surface, fluid-saturated granular flow

NASA Astrophysics Data System (ADS)

Sherry, T. J.; Rowe, C. D.; Kirkpatrick, J. D.; Brodsky, E. E.

2011-12-01

Sand injectites are spectacular examples of large-scale granular flows involving migration of hundreds of cubic meters of sand slurry over hundreds of meters to kilometers in the sub-surface. By studying the macro- and microstructural textures of a kilometer-scale sand injectite, we interpret the fluid flow regimes during emplacement and define the timing of formation of specific textures in the injected material. Fluidized sand sourced from the Santa Margarita Fm., was injected upward into the Santa Cruz Mudstone, Santa Cruz County, California. The sand injectite exposed at Yellow Bank Beach records emplacement of both hydrocarbon and aqueous sand slurries. Elongate, angular mudstone clasts were ripped from the wall rock during sand migration, providing evidence for high velocity, turbid flow. However, clast long axis orientations are consistently sub-horizontal suggesting the slurry transitioned to a laminar flow as the flow velocity decreased in the sill-like intrusion. Millimeter to centimeter scale laminations are ubiquitous throughout the sand body and are locally parallel to the mudstone clast long axes. The laminations are distinct in exposure because alternating layers are preferentially cemented with limonite sourced from later groundwater infiltration. Quantitative microstructural analyses show that the laminations are defined by subtle oscillations in grain alignment between limonite and non-limonite stained layers. Grain packing, size and shape distributions do not vary. The presence of limonite in alternating layers results from differential infiltration of groundwater, indicating permeability changes between the layers despite minimal grain scale differences. Convolute dewatering structures deform the laminations. Dolomite-cemented sand, a signature of hydrocarbon saturation, forms irregular bodies that cross-cut the laminations and dewatering structures. Laminations are not formed in the dolomite-cemented sand. The relative viscosity difference

Assessing coastal plain risk indices for subsurface phosphorus loss

USDA-ARS?s Scientific Manuscript database

Phosphorus (P) Indices are important tools for nutrient management planning in the U.S. whose evaluation often has been deemphasized in favor of research and development. Assessing P Indices in artificially drained agroecosystems is especially important, as subsurface flow is the predominant mode of...

Occurrence and removal of estrogens, progesterone and testosterone in three constructed wetlands treating municipal sewage in the Czech Republic.

PubMed

Vymazal, Jan; Březinová, Tereza; Koželuh, Milan

2015-12-01

Estrogenic hormones, progesterone and testosterone are endocrine-disrupting chemicals and their presence in aquatic environments represents a potentially adverse environmental and public health impact. There is a considerable amount of information about removal of estrogens, progesterone and testosterone in conventional wastewater treatment plants, namely activated sludge systems. However, the information about removal of these compounds in constructed wetlands is very limited. Three constructed wetlands with horizontal subsurface flow in the Czech Republic have been selected to evaluate removal of estrogens (estrone, estriol, 17β-estradiol, 17α-ethinylestradiol), testosterone and progesterone. Monitored constructed wetlands for 100, 150 and 200 PE have been in operation for more than 10 years and all systems exhibit very high treatment efficiency for organics and suspended solids. The results indicate that removal of all estrogens, progesterone and testosterone was high and only estrone was found in the outflow from one constructed wetland in concentrations above the limit of quantification 1 ng l(-1). The limits of quantification for other estrogens, i.e., 10 ng l(-1) for estriol, 1 ng l(-1) for 17β-estradiol and 2 ng l(-1) for 17α-ethinylestradiol were not exceeded in the outflow of all monitored constructed wetlands. Also, for progesterone and testosterone, all outflow concentrations were below the LOQ of 0.5 ng l(-1). The results indicated that constructed wetlands with horizontal subsurface flow are a promising technology for elimination of estrogens, progesterone and testosterone from municipal sewage but more information is needed to confirm this finding. Copyright © 2015 Elsevier B.V. All rights reserved.

Machine learning to construct reduced-order models and scaling laws for reactive-transport applications

NASA Astrophysics Data System (ADS)

Mudunuru, M. K.; Karra, S.; Vesselinov, V. V.

2017-12-01

The efficiency of many hydrogeological applications such as reactive-transport and contaminant remediation vastly depends on the macroscopic mixing occurring in the aquifer. In the case of remediation activities, it is fundamental to enhancement and control of the mixing through impact of the structure of flow field which is impacted by groundwater pumping/extraction, heterogeneity, and anisotropy of the flow medium. However, the relative importance of these hydrogeological parameters to understand mixing process is not well studied. This is partially because to understand and quantify mixing, one needs to perform multiple runs of high-fidelity numerical simulations for various subsurface model inputs. Typically, high-fidelity simulations of existing subsurface models take hours to complete on several thousands of processors. As a result, they may not be feasible to study the importance and impact of model inputs on mixing. Hence, there is a pressing need to develop computationally efficient models to accurately predict the desired QoIs for remediation and reactive-transport applications. An attractive way to construct computationally efficient models is through reduced-order modeling using machine learning. These approaches can substantially improve our capabilities to model and predict remediation process. Reduced-Order Models (ROMs) are similar to analytical solutions or lookup tables. However, the method in which ROMs are constructed is different. Here, we present a physics-informed ML framework to construct ROMs based on high-fidelity numerical simulations. First, random forests, F-test, and mutual information are used to evaluate the importance of model inputs. Second, SVMs are used to construct ROMs based on these inputs. These ROMs are then used to understand mixing under perturbed vortex flows. Finally, we construct scaling laws for certain important QoIs such as degree of mixing and product yield. Scaling law parameters dependence on model inputs are

Facilitated strontium transport by remobilization of strontium-containing secondary precipitates in Hanford Site subsurface.

PubMed

Wang, Guohui; Um, Wooyong

2013-03-15

Significantly enhanced immobilization of radionuclides (such as (90)Sr and (137)Cs) due to adsorption and coprecipitation with neo-formed colloid-sized secondary precipitates has been reported at the U.S. Department of Energy's Hanford Site. However, the stability of these secondary precipitates containing radionuclides in the subsurface under changeable field conditions is not clear. Here, the authors tested the remobilization possibility of Sr-containing secondary precipitates (nitrate-cancrinite) in the subsurface using saturated column experiments under different geochemical and flow conditions. The columns were packed with quartz sand that contained secondary precipitates (nitrate-cancrinite containing Sr), and leached using colloid-free solutions under different flow rates, varying pH, and ionic strength conditions. The results indicate remobilization of the neo-formed secondary precipitates could be possible given a change of pH of ionic strength and flow rate conditions. The remobility of the neo-formed precipitates increased with the rise in the leaching solution flow rate and pH (in a range of pH 4-11), as well as with decreasing solution ionic strength. The increased mobility of Sr-containing secondary precipitates with changing background conditions can be a potential source for additional radionuclide transport in Hanford Site subsurface environments. Published by Elsevier B.V.

TOUGH2Biot - A simulator for coupled thermal-hydrodynamic-mechanical processes in subsurface flow systems: Application to CO2 geological storage and geothermal development

NASA Astrophysics Data System (ADS)

Lei, Hongwu; Xu, Tianfu; Jin, Guangrong

2015-04-01

Coupled thermal-hydrodynamic-mechanical processes have become increasingly important in studying the issues affecting subsurface flow systems, such as CO2 sequestration in deep saline aquifers and geothermal development. In this study, a mechanical module based on the extended Biot consolidation model was developed and incorporated into the well-established thermal-hydrodynamic simulator TOUGH2, resulting in an integrated numerical THM simulation program TOUGH2Biot. A finite element method was employed to discretize space for rock mechanical calculation and the Mohr-Coulomb failure criterion was used to determine if the rock undergoes shear-slip failure. Mechanics is partly coupled with the thermal-hydrodynamic processes and gives feedback to flow through stress-dependent porosity and permeability. TOUGH2Biot was verified against analytical solutions for the 1D Terzaghi consolidation and cooling-induced subsidence. TOUGH2Biot was applied to evaluate the thermal, hydrodynamic, and mechanical responses of CO2 geological sequestration at the Ordos CCS Demonstration Project, China and geothermal exploitation at the Geysers geothermal field, California. The results demonstrate that TOUGH2Biot is capable of analyzing change in pressure and temperature, displacement, stress, and potential shear-slip failure caused by large scale underground man-made activity in subsurface flow systems. TOUGH2Biot can also be easily extended for complex coupled process problems in fractured media and be conveniently updated to parallel versions on different platforms to take advantage of high-performance computing.

Subsurface phosphorus transport through a no-till field in the semi arid Palouse region

NASA Astrophysics Data System (ADS)

Norby, J. C.; Brooks, E. S.; Strawn, D. G.

2017-12-01

Excess application of fertilizers containing nitrogen and phosphorus for farming use has led to ongoing water quality issues in the United States. When these nutrients leave agronomic systems, and enter water bodies in large quantities, algal bloom and eutrophication can occur. Extensive studies focusing on phosphorus as a pollutant from agronomic systems have been conducted in the many regions of the United States; however, there has been a lack of studies completed in the semiarid Palouse region of eastern Washington and western Idaho. The goal of this research study was to better understand how no-till farm management has altered soil P temporally and the current availability for off-site transport of P throughout an artificially drained catchment at the Cook Agronomy Farm in Pullman, WA. We also attempted to determine the processes responsible for subsurface flow of phosphorus, specifically through preferential flow pathways. Dissolved reactive P (DRP)concentrations of subsurface drainage from a artificial drain exceeded TMDL threshold concentrations during numerous seasonal high flow events over the two-year study time frame. Soil analyses show a highly variable distribution of water-extractable P across the sub-catchment area and initial results suggest a translocation of P species deeper into the soil profile after implementing no-till practices in 1998. We hypothesized that a greater network of macropores from lack of soil disturbance allow for preferential flow of nutrient-laden water deeper into the subsurface and to the artificial drain system. Simulated flow experiments on soil cores from the study site showed large-scale macropore development, extreme variability in soil conductivity, and high P adsorption potential for the soils, suggesting a disconnect between P movement through macropore soil and subsurface drainage water rich in DRP at the artificial drain line outlet.

Identifying future directions for subsurface hydrocarbon migration research

NASA Astrophysics Data System (ADS)

Leifer, I.; Clark, J. F.; Luyendyk, B.; Valentine, D.

Subsurface hydrocarbon migration is important for understanding the input and impacts of natural hydrocarbon seepage on the environment. Great uncertainties remain in most aspects of hydrocarbon migration, including some basic mechanisms of this four-phase flow of tar, oil, water, and gas through the complex fracture-network geometry particularly since the phases span a wide range of properties. Academic, government, and industry representatives recently attended a workshop to identify the areas of greatest need for future research in shallow hydrocarbon migration.Novel approaches such as studying temporal and spatial seepage variations and analogous geofluid systems (e.g., geysers and trickle beds) allow deductions of subsurface processes and structures that remain largely unclear. Unique complexities exist in hydrocarbon migration due to its multiphase flow and complex geometry, including in-situ biological weathering. Furthermore, many aspects of the role of hydrocarbons (positive and negative) in the environment are poorly understood, including how they enter the food chain (respiration, consumption, etc.) and “percolate” to higher trophic levels. But understanding these ecological impacts requires knowledge of the emissions' temporal and spatial variability and trajectories.

Intercomparison of Multiscale Modeling Approaches in Simulating Subsurface Flow and Transport

NASA Astrophysics Data System (ADS)

Yang, X.; Mehmani, Y.; Barajas-Solano, D. A.; Song, H. S.; Balhoff, M.; Tartakovsky, A. M.; Scheibe, T. D.

2016-12-01

Hybrid multiscale simulations that couple models across scales are critical to advance predictions of the larger system behavior using understanding of fundamental processes. In the current study, three hybrid multiscale methods are intercompared: multiscale loose-coupling method, multiscale finite volume (MsFV) method and multiscale mortar method. The loose-coupling method enables a parallel workflow structure based on the Swift scripting environment that manages the complex process of executing coupled micro- and macro-scale models without being intrusive to the at-scale simulators. The MsFV method applies microscale and macroscale models over overlapping subdomains of the modeling domain and enforces continuity of concentration and transport fluxes between models via restriction and prolongation operators. The mortar method is a non-overlapping domain decomposition approach capable of coupling all permutations of pore- and continuum-scale models with each other. In doing so, Lagrange multipliers are used at interfaces shared between the subdomains so as to establish continuity of species/fluid mass flux. Subdomain computations can be performed either concurrently or non-concurrently depending on the algorithm used. All the above methods have been proven to be accurate and efficient in studying flow and transport in porous media. However, there has not been any field-scale applications and benchmarking among various hybrid multiscale approaches. To address this challenge, we apply all three hybrid multiscale methods to simulate water flow and transport in a conceptualized 2D modeling domain of the hyporheic zone, where strong interactions between groundwater and surface water exist across multiple scales. In all three multiscale methods, fine-scale simulations are applied to a thin layer of riverbed alluvial sediments while the macroscopic simulations are used for the larger subsurface aquifer domain. Different numerical coupling methods are then applied between

Efficient uncertainty quantification in fully-integrated surface and subsurface hydrologic simulations

NASA Astrophysics Data System (ADS)

Miller, K. L.; Berg, S. J.; Davison, J. H.; Sudicky, E. A.; Forsyth, P. A.

2018-01-01

Although high performance computers and advanced numerical methods have made the application of fully-integrated surface and subsurface flow and transport models such as HydroGeoSphere common place, run times for large complex basin models can still be on the order of days to weeks, thus, limiting the usefulness of traditional workhorse algorithms for uncertainty quantification (UQ) such as Latin Hypercube simulation (LHS) or Monte Carlo simulation (MCS), which generally require thousands of simulations to achieve an acceptable level of accuracy. In this paper we investigate non-intrusive polynomial chaos for uncertainty quantification, which in contrast to random sampling methods (e.g., LHS and MCS), represents a model response of interest as a weighted sum of polynomials over the random inputs. Once a chaos expansion has been constructed, approximating the mean, covariance, probability density function, cumulative distribution function, and other common statistics as well as local and global sensitivity measures is straightforward and computationally inexpensive, thus making PCE an attractive UQ method for hydrologic models with long run times. Our polynomial chaos implementation was validated through comparison with analytical solutions as well as solutions obtained via LHS for simple numerical problems. It was then used to quantify parametric uncertainty in a series of numerical problems with increasing complexity, including a two-dimensional fully-saturated, steady flow and transient transport problem with six uncertain parameters and one quantity of interest; a one-dimensional variably-saturated column test involving transient flow and transport, four uncertain parameters, and two quantities of interest at 101 spatial locations and five different times each (1010 total); and a three-dimensional fully-integrated surface and subsurface flow and transport problem for a small test catchment involving seven uncertain parameters and three quantities of interest at

Pattern-Based Inverse Modeling for Characterization of Subsurface Flow Models with Complex Geologic Heterogeneity

NASA Astrophysics Data System (ADS)

Golmohammadi, A.; Jafarpour, B.; M Khaninezhad, M. R.

2017-12-01

Calibration of heterogeneous subsurface flow models leads to ill-posed nonlinear inverse problems, where too many unknown parameters are estimated from limited response measurements. When the underlying parameters form complex (non-Gaussian) structured spatial connectivity patterns, classical variogram-based geostatistical techniques cannot describe the underlying connectivity patterns. Modern pattern-based geostatistical methods that incorporate higher-order spatial statistics are more suitable for describing such complex spatial patterns. Moreover, when the underlying unknown parameters are discrete (geologic facies distribution), conventional model calibration techniques that are designed for continuous parameters cannot be applied directly. In this paper, we introduce a novel pattern-based model calibration method to reconstruct discrete and spatially complex facies distributions from dynamic flow response data. To reproduce complex connectivity patterns during model calibration, we impose a feasibility constraint to ensure that the solution follows the expected higher-order spatial statistics. For model calibration, we adopt a regularized least-squares formulation, involving data mismatch, pattern connectivity, and feasibility constraint terms. Using an alternating directions optimization algorithm, the regularized objective function is divided into a continuous model calibration problem, followed by mapping the solution onto the feasible set. The feasibility constraint to honor the expected spatial statistics is implemented using a supervised machine learning algorithm. The two steps of the model calibration formulation are repeated until the convergence criterion is met. Several numerical examples are used to evaluate the performance of the developed method.

Subsurface Hydrologic Processes Revealed by Time-lapse GPR in Two Contrasting Soils in the Shale Hills CZO

NASA Astrophysics Data System (ADS)

Guo, L.; Lin, H.; Nyquist, J.; Toran, L.; Mount, G.

2017-12-01

Linking subsurface structures to their functions in determining hydrologic processes, such as soil moisture dynamics, subsurface flow patterns, and discharge behaviours, is a key to understanding and modelling hydrological systems. Geophysical techniques provide a non-invasive approach to investigate this form-function dualism of subsurface hydrology at the field scale, because they are effective in visualizing subsurface structure and monitoring the distribution of water. In this study, we used time-lapse ground-penetrating radar (GPR) to compare the hydrologic responses of two contrasting soils in the Shale Hills Critical Zone Observatory. By integrating time-lapse GPR with artificial water injection, we observed distinct flow patterns in the two soils: 1) in the deep Rushtown soil (over 1.5 m depth to bedrock) located in a concave hillslope, a lateral preferential flow network extending as far as 2 m downslope was identified above a less permeable layer and via a series of connected macropores; whereas 2) in the shallow Weikert soil ( 0.3 m depth to saprock) located in a planar hillslope, vertical infiltration into the permeable fractured shale dominated the flow field, while the development of lateral preferential flow along the hillslope was restrained. At the Weikert soil site, the addition of brilliant blue dye to the water injection followed by in situ excavation supported GPR interpretation that only limited lateral preferential flow formed along the soil-saprock interface. Moreover, seasonally repeated GPR surveys indicated different patterns of profile moisture distribution in the two soils that in comparison with the dry season, a dense layer within the BC horizon in the deep Rushtown soil prevented vertical infiltration in the wet season, leading to the accumulation of soil moisture above this layer; whereas, in the shallow Weikert soil, water infiltrated into saprock in wet seasons, building up water storage within the fractured bedrock (i.e., the

Inputs and losses by surface runoff and subsurface leaching for pastures managed by continuous or rotational stocking.

PubMed

Owens, L B; Barker, D J; Loerch, S C; Shipitalo, M J; Bonta, J V; Sulc, R M

2012-01-01

Pasture management practices can affect forage quality and production, animal health and production, and surface and groundwater quality. In a 5-yr study conducted at the North Appalachian Experimental Watershed near Coshocton, Ohio, we compared the effects of two contrasting grazing methods on surface and subsurface water quantity and quality. Four pastures, each including a small, instrumented watershed (0.51-1.09 ha) for surface runoff measurements and a developed spring for subsurface flow collection, received 112 kg N ha(-1) yr(-1) and were grazed at similar stocking rates (1.8-1.9 cows ha(-1)). Two pastures were continuously stocked; two were subdivided so that they were grazed with frequent rotational stocking (5-6 times weekly). In the preceding 5 yr, these pastures received 112 kg N ha(-1) yr(-1) after several years of 0 N fertilizer and were grazed with weekly rotational stocking. Surface runoff losses of N were minimal. During these two periods, some years had precipitation up to 50% greater than the long-term average, which increased subsurface flow and NO(3)-N transport. Average annual NO(3)-N transported in subsurface flow from the four watersheds during the two 5-yr periods ranged from 11.3 to 22.7 kg N ha(-1), which was similar to or less than the mineral-N received in precipitation. Flow and transport variations were greater among seasons than among watersheds. Flow-weighted seasonal NO(3)-N concentrations in subsurface flow did not exceed 7 mg L(-1). Variations in NO(3)-N leached from pastures were primarily due to variable precipitation rather than the effects of continuous, weekly rotational, or frequent rotational stocking practices. This suggests that there was no difference among these grazing practices in terms of NO(3)-N leaching. Copyright © by the American Society of Agronomy, Crop Science Society of America, and Soil Science Society of America, Inc.

Analytical solution for tension-saturated and unsaturated flow from wicking porous pipes in subsurface irrigation: The Kornev-Philip legacies revisited

NASA Astrophysics Data System (ADS)

Kacimov, A. R.; Obnosov, Yu. V.

2017-03-01

The Russian engineer Kornev in his 1935 book raised perspectives of subsurface "negative pressure" irrigation, which have been overlooked in modern soil science. Kornev's autoirrigation utilizes wicking of a vacuumed water from a porous pipe into a dry adjacent soil. We link Kornev's technology with a slightly modified Philip (1984)'s analytical solutions for unsaturated flow from a 2-D cylindrical pipe in an infinite domain. Two Darcian flows are considered and connected through continuity of pressure along the pipe-soil contact. The first fragment is a thin porous pipe wall in which water seeps at tension saturation; the hydraulic head is a harmonic function varying purely radially across the wall. The Thiem solution in this fragment gives the boundary condition for azimuthally varying suction pressure in the second fragment, ambient soil, making the exterior of the pipe. The constant head, rather than Philip's isobaricity boundary condition, along the external wall slightly modifies Philip's formulae for the Kirchhoff potential and pressure head in the soil fragment. Flow characteristics (magnitudes of the Darcian velocity, total flow rate, and flow net) are explicitly expressed through series of Macdonald's functions. For a given pipe's external diameter, wall thickness, position of the pipe above a free water datum in the supply tank, saturated conductivities of the wall and soil, and soil's sorptive number, a nonlinear equation with respect to the total discharge from the pipe is obtained and solved by a computer algebra routine. Efficiency of irrigation is evaluated by computation of the moisture content within selected zones surrounding the porous pipe.