Low‐Resolution Modeling of Dense Drainage Networks in Confining Layers

Groundwater‐surface water (GW‐SW) interaction in numerical groundwater flow models is generally simulated using a Cauchy boundary condition, which relates the flow between the surface water and the groundwater to the product of the head difference between the node and the surface water level, and a coefficient, often referred to as the “conductance.” Previous studies have shown that in models with a low grid resolution, the resistance to GW‐SW interaction below the surface water bed should often be accounted for in the parameterization of the conductance, in addition to the resistance across the surface water bed. Three conductance expressions that take this resistance into account were investigated: two that were presented by Mehl and Hill (2010) and the one that was presented by De Lange (1999). Their accuracy in low‐resolution models regarding salt and water fluxes to a dense drainage network in a confined aquifer system was determined. For a wide range of hydrogeological conditions, the influence of (1) variable groundwater density; (2) vertical grid discretization; and (3) simulation of both ditches and tile drains in a single model cell was investigated. The results indicate that the conductance expression of De Lange (1999) should be used in similar hydrogeological conditions as considered in this paper, as it is better taking into account the resistance to flow below the surface water bed. For the cases that were considered, the influence of variable groundwater density and vertical grid discretization on the accuracy of the conductance expression of De Lange (1999) is small.     

Pesticides Removal Using Conventional and Low‐Cost Adsorbents: A Review

Worldwide pesticide usage has increased dramatically during the last three decades, coinciding with changing practices and increasing by intensive agriculture. This widespread use of pesticides for agricultural and non‐agricultural purposes has resulted in the presence of their residues in various environmental matrices. The occurrence of pesticides and their metabolite transported in rivers, channels, lakes, sea, air, soils, groundwater, and even drinking water, proves the high risk of these chemicals to human health and the environment. Therefore, pesticide removal is of an increasing concern. In this study, a review of the published literature dealing with pesticides removal process is presented. Firstly, pesticide removal by conventional means is briefly considered. Secondly, the use of the low‐cost sorbent through biosorption process is discussed comprehensively. The effect of factors such as pH, contact time, sorbent dosage, initial pesticide concentration, and optimization of biosorption conditions is also discussed. Kinetic, thermodynamic, and mechanism studies are also given. This study shows that both microorganisms and other materials with biological origin like agricultural by‐products may be used to this end. There is a significant potential for pesticide uptake by the use of various pristine and especially modified biosorbents. In the case of living organisms used as removal agents, degradation may also play a role in the total removal observed.     

Biofilm Detection in a Model Well‐Bore Environment Using Low‐Field NMR

This research addresses the challenges of the lack of non‐invasive methods and poor spatiotemporal resolution associated with monitoring biogeochemical activity central to bioremediation of subsurface contaminants. Remediation efforts often include growth of biofilm to contain or degrade chemical contaminants, such as nitrates, hydrocarbons, heavy metals, and some chlorinated solvents. Previous research indicates that nuclear magnetic resonance (NMR) is sensitive to the biogeochemical processes of biofilm accumulation. The current research focuses on developing methods to use low‐cost NMR technology to support in situ monitoring of biofilm growth and geochemical remediation processes in the subsurface. Biofilm was grown in a lab‐scale radial flow bioreactor designed to model the near wellbore subsurface environment. The Vista Clara Javelin NMR logging device, a slim down‐the‐borehole probe, collected NMR measurements over the course of eight days while biofilm was cultivated in the sand‐packed reactor. Measured NMR mean log T₂ relaxation times decreased from approximately 710 to 389 ms, indicating that the pore environment and bulk fluid properties were changing due to biofilm growth. Destructive sampling employing drop plate microbial population analysis and scanning electron and stereoscopic microscopy confirmed biofilm formation. Our findings demonstrate that the NMR logging tool can detect small to moderate changes in T₂ distribution associated with environmentally relevant quantities of biofilm in quartz sand.     

Tide‐Induced Air Pressure Fluctuations in a Coastal Unsaturated Zone: Effects of Thin Low‐Permeability Pavements

This paper presents an analytical case study to explore one‐dimensional subsurface air pressure variation in a coastal three‐layered unsaturated zone. The upper layer is thin and much less permeable than the middle layer, and water table is located in the very permeable lower layer. An analytical solution was derived to describe the air pressure variation caused by tide‐induced water table fluctuations. We revisited the case study at Hong Kong International Airport conducted by Jiao and Li (2004) who used a two‐dimensional numerical model. The analytical prediction using the parameter values equivalent to the two‐dimensional numerical model agreed very well with the observed air pressure, indicating the validity and applicability of our one‐dimensional model in approximating the actual situation in this coastal zone with adequate accuracy. The analysis revealed that the asphalt pavement played an important role in causing air pressure fluctuations below it. Abnormally high air pressure can be caused beneath the surface pavement when the air permeability decreases due to rainfall infiltration, which may lead to heaving problems during rising tides.     

Experimental and Economic Assessment of Two Surfactant Formulations for Source Zone Remediation at a Former Dry Cleaning Facility

Treatability tests and cost analyses were conducted to provide objective criteria for selection of a surfactant formulation to be used for surfactant enhanced aquifer remediation (SEAR) of a tetrachloroethene (PCE)-contaminated site in Oscoda, Michigan. Two surfactant formulations, 4% Tween 80 + 500 mg/L CaCl₂ and 8% Aerosol MA/IPA +15,000 mg/L NaCl + 1000 mg/L CaCl₂, were considered based on their capacity to solubilize PCE and prior use in SEAR applications. Results of a two-dimensional aquifer cell experiment indicated that 53% of the released PCE was recovered after flushing with approximately 8 pore volumes of 4% Tween 80. In contrast, only 3 pore volumes of 8% Aerosol MA/IPA solution were required to recover 78% of the PCE from the two-dimensional aquifer cell, although the greater recovery of PCE was attributed, in large part, to the higher concentration of Aerosol MA. However, mobilization of PCE as free product was observed during the 8% Aerosol MA/IPA flood, which was consistent with total trapping number (N_T) calculations. At the pilot-scale, SEAR treatment costs were estimated to be $222,000 and $244,000 for 4% Tween 80 and 8% Aerosol MA/IPA, respectively, which compared favorably to the estimated pump-and-treat cost of $316,000. Projected full-scale costs, based on a line-drive flushing system, were $382,000 for 4% Tween 80 and $443,000 for 8% Aerosol MA/IPA. In contrast, full-scale pump-and-treat costs were estimated to be $1,167,000. Surfactant recycling was shown to be logistically and economically infeasible at the pilot scale, and provided only a minimal cost benefit for 4% Tween 80 at the full scale. Based on the similarities in solubilization capacity and treatment cost, but substantially lower risk of PCE displacement, Tween 80 was recommended over Aerosol MA/IPA for pilot-scale testing of SEAR.     

Low‐Flow Hydraulic Conductivity Tests at Wells that Cross the Water Table

Wells with screens and sand packs that cross the water table represent a challenging problem for determining hydraulic conductivity by slug testing due to sand pack drainage and resaturation. Sand pack drainage results in a multisegmented recovery curve. One must then subjectively pick a portion of the curve to analyze. Sand pack drainage also results in a change in the effective radius of the well which requires a guess at the porosity or specific yield in analyzing the test. In the study of Robbins et al. (2009) , a method was introduced to obtain hydraulic conductivity in monitoring wells using the steady‐state drawdown and flow rate obtained during low‐flow sampling. The method was tested in this study in wells whose screens cross the water table and shown to avoid sand pack drainage problems that complicate analyzing slug tests. In applying the method to low‐flow sampling, only a single pair of steady‐state flow rate and drawdown are needed; hence, to derive meaningful results, an accurate determination of these parameters is required.     

Using Soil and Meteorologic Data Bases in Unsaturated Zone Modeling of Pesticides

Various types of models are being used to evaluate pesticide transport and transformation in the unsaturated zone. Model predictions can be used, for example, to develop alternative agricultural management strategies for pesticide use. However, intensive data requirements for transient models sometimes deter their use. Site-specific measurements are preferred, but existing data bases can be used as a source of required model parameters, especially weather and soil characteristics. These existing data bases make possible the use of models to predict leaching potential in a wide variety of environments.     

Evaluation of Two Low‐Cost–High‐Performance Adsorbent Materials in the Waste‐to‐Product Approach for the Removal of Pesticides from Drinking Water

This study evaluates the performance of two low cost and high performance adsorption materials, i.e., activated carbon produced from two natural waste products: Bamboo and coconut shell, in the removal of three pesticides from drinking water sources. Due to the fact that bamboo and coconut shell are abundant and inexpensive materials in many parts of the world, they respond to the “low‐cost” aspect. The adsorption capacities of two local adsorbents have been compared with commercial activated carbon to explore their potential to respond to the “high quality” aspect. Two pesticides were selected, namely dieldrin and chlorpyrifos, because they are commonly used in agriculture activities, and may remain in high concentrations in surface water used as drinking water sources. The results indicate that the adsorption of pesticides on activated carbons is influenced by physico‐chemical properties of the activated carbon and the pesticides such as the presence of an aromatic ring, and their molar mass. The activated carbon produced from bamboo can be employed as low‐cost and high performance adsorbent, alternative to commercial activated carbon for the removal of pesticides during drinking water production. The performance of activated carbon from bamboo was better due to its relatively large macroporosity and planar surface. The effect of adsorbent and pesticide characteristics on the performance was derived from batch experiments in which the adsorption behavior was studied on the basis of Freundlich isotherms.     

Removal of Organic Pollutants from Wastewater Using Wood Fly Ash as a Low‐Cost Sorbent

In this study, untreated and treated wood fly ash (WA) was used as a low‐cost sorbent in batch sorption tests to investigate the removal of organic pollutants from a real wastewater generated by cleaning/washing of machinery in a wood‐laminate floor industry in Sweden. The experiments focused on the effect of the WA dosage and particle size on the removal efficiency for organic compounds. With a WA dosage of 160 g L^?1 and a particle size less than 1 mm, the reductions of chemical oxygen demand (COD), biologic oxygen demand, and total organic carbon were 37 ± 0.4, 24 ± 0.4, and 30 ± 0.3%, respectively. Pre‐treatment of WA with hot water improved the COD removal efficiency by absorption from 37 ± 0.4 to 42 ± 1.6% when the same dosage (160 g L^?1) was applied. Sorption isotherm and sorption kinetics for COD using untreated WA can be explained by Freundlich isotherm and pseudo‐second‐order kinetic models. Intra‐particle diffusion model indicates that pore diffusion is not the rate‐limiting step for COD removal. Based on the experimental data, WA could be used as an alternative low‐cost sorption media/filter for removal of organic compounds from real industrial wastewater.     

Hydraulic Tomography: Continuity and Discontinuity of High‐K and Low‐K Zones

Hydraulic tomography is an emerging field and modeling method that provides a continuous hydraulic conductivity (K) distribution for an investigated region. Characterization approaches that rely on interpolation between one‐dimensional (1D) profiles have limited ability to accurately identify high‐K channels, juxtapositions of lenses with high K contrast, and breaches in layers or channels between such profiles. However, locating these features is especially important for groundwater flow and transport modeling, and for design and operation of in situ remediation in complex hydrogeologic environments. We use transient hydraulic tomography to estimate 3D K in a volume of 15‐m diameter by 20‐m saturated thickness in a highly heterogeneous unconfined alluvial (clay to sand‐and‐gravel) aquifer with a K range of approximately seven orders of magnitude at an active industrial site in Assemini, Sardinia, Italy. A modified Levenberg‐Marquardt algorithm was used for geostatistical inversion to deal with the nonlinear nature of the highly heterogeneous system. The imaging results are validated with pumping tests not used in the tomographic inversion. These tests were conducted from three of five clusters of continuous multichannel tubing (CMTs) installed for observation in the tomographic testing. Locations of high‐K continuity and discontinuity, juxtaposition of very high‐K and very low‐K lenses, and low‐K “plugs” are evident in regions of the investigated volume where they likely would not have been identified with interpolation from 1D profiles at the positions of the pumping well and five CMT clusters. Quality assessment methods identified a suspect high‐K feature between the tested volume and a lateral boundary of the model.     

Modeling River Stage‐Discharge Relationships Using Different Neural Network Computing Techniques

One of the most important problems in hydrology is the establishment of rating curves. The statistical tools that are commonly used for river stage‐discharge relationships are regression and curve fitting. However, these techniques are not adequate in view of the complexity of the problems involved. Three different neural network techniques, i. e., multi‐layer perceptron neural network with Levenberg‐Marquardt and quasi‐Newton algorithms and radial basis neural networks, are used for the development of river stage‐discharge relationships by constructing nonlinear relationships between stage and discharge. Daily stage and flow data from three stations, Yamula, Tuzkoy and Sogutluhan, on the Kizilirmak River in Turkey were used. Regression techniques are also applied to the same data. Different input combinations including the previous stages and discharges are used. The models' results are compared using three criteria, i. e., root mean square errors, mean absolute error and the determination coefficient. The results of the comparison reveal that the neural network techniques are much more suitable for setting up stage‐discharge relationships than the regression techniques. Among the neural network methods, the radial basis neural network is found to be slightly better than the others.     

An In‐Depth Characterization of Urban Aerosols Using Electron Microscopy and Energy Dispersive X‐Ray Analysis

The physical and chemical characterization of aerosols from three large cities, Karachi and Islamabad, Pakistan, and New York City (NYC), USA, was investigated. A scanning electron microscope equipped with an energy dispersive spectrometer (EDS) was used to determine particle morphology and elemental composition of the samples. A Bruker Spirit system in combination with a Sahara detector provided both computer controlled Automated Chemical Classification (ACC) and digital mapping features for analysis purposes. The use of these two features to characterize the elemental composition, particle size, and to determine specific classes or source types is described in this paper. Filters were analyzed for the following elements; Na, Mg, Al, Si, P, S, Cl, K, Ca, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, Zn, Ga, As, Se, Br, Sr, Y, Zr, Mo, Pd, Ag, Cd, In, Sn, Sb, Ba, La, Au, Hg, Tl, Pb, and U. The EDS work was qualitative not quantitative. Seven source types (mobile, cement, soil, steel mill, fossil fuel, sea salt and biological) contributed to particulate matter in the ambient air of both cities in Pakistan, whereas there were eight source types (diesel, road dust, automobile, iron, residual fuel oil, sea salt, ammonium sulfate and ammonium nitrate) for NYC. For all three urban centers, vehicular traffic emissions were the main contributor to particle number. Diesel emissions from trucks and buses were relatively more prominent in NYC aerosols, while gasoline emission from automotive exhaust was dominant in the two cities of Pakistan. The cement/limestone component from local cement industries was very evident in both particle surface characteristics and number for both Karachi and Islamabad, but not in NYC air. Sea salt aerosols were seen in the two coastal cities, Karachi and NYC. It was also witnessed in Islamabad aerosols and was attributable to the mining of rock salt at the world's richest salt mine, Khewra, situated upwind from the city.     

Modeling Water and Heat Regimes of Arid Areas at Varying Atmospheric Precipitation

Water and heat regimes of an arid area are simulated with the use of differential equations of the water and heat balances. The equations are developed, analyzed, and used to calculate the dynamics of the main hydrological and climatic characteristics, given the linear increase or variations in the amount of precipitation. Calculations are made for the area of the Caspian Depression and West Khazakhstan.     

Adsorption Kinetic Modeling of Safranin onto Rice Husk Biomatrix Using Pseudo‐first‐ and Pseudo‐second‐order Kinetic Models: Comparison of Linear and Non‐linear Methods

Batch kinetic studies were carried out for the removal of safranin from aqueous solution using a biomatrix prepared from rice husk. The adsorption kinetic data were modeled using the pseudo‐first‐order and pseudo‐second‐order kinetic equations. The linear and non‐linear forms of these two widely used kinetic models were compared in this study. In order to determine the best‐fitting equation, the coefficient of determination (r²), the sum of the squares of the errors (SSE), sum of the absolute errors (SAE), average relative error (ARE), hybrid fractional error function (HYBRID), Marquardt's percent standard deviation (MPSD), and the Chi‐squared test (χ²) were used as error analysis methods. Results showed that the non‐linear forms of pseudo‐first‐order and pseudo‐second‐order models were more suitable than the linear forms for fitting the experimental data. Non‐linear method is thus more appropriate for estimating the kinetic parameters and should primarily be used to describe adsorption kinetics.     

A Low‐Cost,In Situ Resistivity and Temperature Monitoring System

We present a low‐cost, reliable method for long‐term in situ autonomous monitoring of subsurface resistivity and temperature in a shallow, moderately heterogeneous subsurface. Probes, to be left in situ, were constructed at relatively low cost with an electrode spacing of 5 cm. Once installed, these were wired to the CR‐1000 Campbell Scientific Inc. datalogger at the surface to electrically image infiltration fronts in the shallow subsurface. This system was constructed and installed in June 2005 to collect apparent resistivity and temperature data from 96 subsurface electrodes set to a pole‐pole resistivity array pattern and 14 thermistors at regular intervals of 30 cm through May of 2008. From these data, a temperature and resistivity relationship was determined within the vadose zone (to a depth of ~1 m) and within the saturated zone (at depths between 1 and 2 m). The high vertical resolution of the data with resistivity measurements on a scale of 5‐cm spacing coupled with surface precipitation measurements taken at 3‐min intervals for a period of roughly 3 years allowed unique observations of infiltration related to seasonal changes. Both the vertical resistivity instrument probes and the data logger system functioned well for the duration of the test period and demonstrated the capability of this low‐cost monitoring system.     

Calibration of the Regional Crustal Waveguide and the Retrieval of Source Parameters Using Waveform Modeling

v--vRegional crustal waveguide calibration is essential to the retrieval of source parameters and the location of smaller (M < 4.8) seismic events. This path calibration of regional seismic phases is strongly dependent on the accuracy of hypocentral locations of calibration (or master) events. This information can be difficult to obtain, especially for smaller events. Generally, explosion or quarry blast generated travel-time data with known locations and origin times are useful for developing the path calibration parameters, but in many regions such data sets are scanty or do not exist. We present a method which is useful for regional path calibration independent of such data, i.e. with earthquakes, which is applicable for events down to M_w = 4 and which has successfully been applied in India, central Asia, western Mediterranean, North Africa, Tibet and the former Soviet Union. These studies suggest that reliably determining depth is essential to establishing accurate epicentral location and origin time for events. We find that the error in source depth does not necessarily trade-off only with the origin time for events with poor azimuthal coverage, but with the horizontal location as well, thus resulting in poor epicentral locations. For example, hypocenters for some events in central Asia were found to move from their fixed-depth locations by about 20 km. Such errors in location and depth will propagate into path calibration parameters, particularly with respect to travel times. The modeling of teleseismic depth phases (pP, sP) yields accurate depths for earthquakes down to magnitude M_w = 4.7. This M_w threshold can be lowered to four if regional seismograms are used in conjunction with a calibrated velocity structure model to determine depth, with the relative amplitude of the P_nl waves to the surface waves and the interaction of regional sPmP and pPmP phases being good indicators of event depths. We also found that for deep events a seismic phase which follows an S-wave path to the surface and becomes critical, developing a head wave by S to P conversion is also indicative of depth. The detailed characteristic of this phase is controlled by the crustal waveguide. The key to calibrating regionalized crustal velocity structure is to determine depths for a set of master events by applying the above methods and then by modeling characteristic features that are recorded on the regional waveforms. The regionalization scheme can also incorporate mixed-path crustal waveguide models for cases in which seismic waves traverse two or more distinctly different crustal structures. We also demonstrate that once depths are established, we need only two-stations travel-time data to obtain reliable epicentral locations using a new adaptive grid-search technique which yields locations similar to those determined using travel-time data from local seismic networks with better azimuthal coverage.