SEARCH · Geology Reports

Results for “Cells”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 487 records · Page 27Linked to original sources

USGS2011

The high life: Transport of microbes in the atmosphere

Microbes (bacteria, fungi, algae, and viruses) are the most successful types of life on Earth because of their ability to adapt to new environments, reproduce quickly, and disperse globally. Dispersal occurs through a number of vectors, such as migrating animals or the hydrological cycle, but transport by wind may be the most common way microbes spread. General awareness of airborne microbes predates the science of microbiology. People took advantage of wild airborne yeasts to cultivate lighter, more desirable bread as far back as ancient Egypt by simply leaving a mixture of grain and liquids near an open window. In 1862, Louis Pasteur's quest to disprove spontaneous generation resulted in the discovery that microbes were actually single-celled, living creatures, prevalent in the environment and easily killed with heat (pasteurization). His rudimentary experiments determined that any nutrient medium left open to the air would eventually teem with microbial life because of free-floating, colonizing cells. The same can happen in a kitchen: Opportunistic fungal and bacterial cells cause food items exposed to the air to eventually spoil.

Eos, Transactions, American Geophysical Union↗

USGS2020

Sediment connectivity: A framework for analyzing coastal sediment transport pathways

Connectivity provides a framework for analyzing coastal sediment transport pathways, building on conceptual advances in graph theory from other scientific disciplines. Connectivity schematizes sediment pathways as a directed graph (i.e., a set of nodes and links). This study presents a novel application of graph theory and connectivity metrics like modularity and centrality to coastal sediment dynamics, exemplified here using Ameland Inlet in the Netherlands. We divide the study site into geomorphic cells (i.e., nodes) and then quantify sediment transport between these cells (i.e., links) using a numerical model. The system of cells and fluxes between them is then schematized in a network described by an adjacency matrix. Network metrics like link density, asymmetry, and modularity quantify system‐wide connectivity. The degree, strength, and centrality of individual nodes identify key locations and pathways throughout the system. For instance, these metrics indicate that under strictly tidal forcing, sand originating near shore predominantly bypasses Ameland Inlet via the inlet channels, whereas sand on the deeper foreshore mainly bypasses the inlet via the outer delta shoals. Connectivity analysis can also inform practical management decisions about where to place sand nourishments, the fate of nourishment sand, or how to monitor locations vulnerable to perturbations. There are still open challenges associated with quantifying connectivity at varying space and time scales and the development of connectivity metrics specific to coastal systems. Nonetheless, connectivity provides a promising technique for predicting the response of our coasts to climate change and the human adaptations it provokes.

Ameland Inlet↗

USGS2025

Physicochemical properties and bioreactivity of sub-10 µm geogenic particles: Comparison of volcanic ash and desert dust

Exposure to ambient particulate matter (PM) with an aerodynamic diameter of <10 μm (PM 10 ) is a well-established health hazard. There is increasing evidence that geogenic (Earth-derived) particles can induce adverse biological effects upon inhalation, though there is high variability in particle bioreactivity that is associated with particle source and physicochemical properties. In this study, we investigated physicochemical properties and biological reactivity of volcanic ash from the April 2021 eruption of La Soufrière volcano, St. Vincent, and two desert dust samples: a standardized test dust from Arizona and an aeolian Gobi Desert dust sampled in China. We determined particle size, morphology, mineralogy, surface texture and chemistry in sub-10 μm material to investigate associations between particle physicochemical properties and observed bioreactivity. We assessed cellular responses (cytotoxic and pro-inflammatory effects) to acute particle exposures (24 hr) in monocultures at the air-liquid interface using two types of cells of the human airways: BEAS-2B bronchial epithelial cells and A549 alveolar type II epithelial cells. In acellular assays, we also assessed particle oxidative potential and the presence of microorganisms. The results showed that volcanic ash and desert dust exhibit intrinsically different particle morphology, surface textures and chemistry, and variable mineralogical content. We found that Gobi Desert dust is more bioreactive than freshly erupted volcanic ash and Arizona test dust, which is possibly linked to the presence of microorganisms (bacteria) and/or nanoscale elongated silicate minerals (potentially clay such as illite or vermiculite) on particle surfaces.

GeoHealth↗

USGS1998

Nonlinear refraction and reflection travel time tomography

We develop a rapid nonlinear travel time tomography method that simultaneously inverts refraction and reflection travel times on a regular velocity grid. For travel time and ray path calculations, we apply a wave front method employing graph theory. The first-arrival refraction travel times are calculated on the basis of cell velocities, and the later refraction and reflection travel times are computed using both cell velocities and given interfaces. We solve a regularized nonlinear inverse problem. A Laplacian operator is applied to regularize the model parameters (cell slownesses and reflector geometry) so that the inverse problem is valid for a continuum. The travel times are also regularized such that we invert travel time curves rather than travel time points. A conjugate gradient method is applied to minimize the nonlinear objective function. After obtaining a solution, we perform nonlinear Monte Carlo inversions for uncertainty analysis and compute the posterior model covariance. In numerical experiments, we demonstrate that combining the first arrival refraction travel times with later reflection travel times can better reconstruct the velocity field as well as the reflector geometry. This combination is particularly important for modeling crustal structures where large velocity variations occur in the upper crust. We apply this approach to model the crustal structure of the California Borderland using ocean bottom seismometer and land data collected during the Los Angeles Region Seismic Experiment along two marine survey lines. Details of our image include a high-velocity zone under the Catalina Ridge, but a smooth gradient zone between. Catalina Ridge and San Clemente Ridge. The Moho depth is about 22 km with lateral variations. Copyright 1998 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

USGS2001

Large carbon isotope fractionation associated with oxidation of methyl halides by methylotrophic bacteria

The largest biological fractionations of stable carbon isotopes observed in nature occur during production of methane by methanogenic archaea. These fractionations result in substantial (as much as ≈70‰) shifts in δ 13 C relative to the initial substrate. We now report that a stable carbon isotopic fractionation of comparable magnitude (up to 70‰) occurs during oxidation of methyl halides by methylotrophic bacteria. We have demonstrated biological fractionation with whole cells of three methylotrophs (strain IMB-1, strain CC495, and strain MB2) and, to a lesser extent, with the purified cobalamin-dependent methyltransferase enzyme obtained from strain CC495. Thus, the genetic similarities recently reported between methylotrophs, and methanogens with respect to their pathways for C 1 -unit metabolism are also reflected in the carbon isotopic fractionations achieved by these organisms. We found that only part of the observed fractionation of carbon isotopes could be accounted for by the activity of the corrinoid methyltransferase enzyme, suggesting fractionation by enzymes further along the degradation pathway. These observations are of potential biogeochemical significance in the application of stable carbon isotope ratios to constrain the tropospheric budgets for the ozone-depleting halocarbons, methyl bromide and methyl chloride. Methyl bromide (MeBr) and methyl chloride (MeCl) are, respectively, the most abundant volatile bromo- and chlorocarbons in the troposphere and are major contributors to stratospheric ozone destruction (1). Both compounds have natural and human-influenced sources and a predominant sink by reaction with OH in the troposphere (2–4). MeBr also has a bacterial soil sink (5) that represents about 20% of the estimated total removal from the troposphere, and it is likely that a soil sink of similar magnitude exists for MeCl (6). Hence, if an isotopic fractionation is associated with the soil sink, it will influence the isotopic compositions of MeBr and MeCl in the lower atmosphere (7). The δ 13 C value of industrially produced MeBr ranges between −43.5‰ and −66.4‰ (7), but δ 13 C values of tropospheric MeBr and natural sources are not yet known. The δ 13 C of atmospheric MeCl has been measured from −22‰ to −45‰ (8, 9). If carbon isotope ratios are to be used to constrain the budgets of these methyl halides, it is essential to determine the extent of carbon isotope fractionation that occurs during biological degradation of these compounds. Methylotrophic bacteria use C 1 compounds, which are simple organic molecules that contain no carbon–carbon bonds. Strains IMB-1, CC495, and MB2 are as-yet-unnamed facultative methylotrophs isolated from agricultural soil, woodland leaf litter, and coastal seawater, respectively (10–13), environments where methyl halides are produced. They are members of the α subgroup of the Proteobacteria . On the basis of 16S rRNA gene sequences, strains IMB-1 and CC495 show some phylogenetic alignment with the genus Rhizobium (10, 11) and are very closely related to the new genus Pseudoaminobacte r (I. McDonald, personal communication). Strain MB2 aligns within the Ruegeria clade [J. K. Schaefer, K. D. Goodwin, I. R. McDonald, J. C. Murrell and R.S.O., unpublished work]. All of these aerobic bacteria are methylotrophs in that they can grow by using MeBr or MeCl as their sole carbon source, but they do not metabolize methane. They oxidize MeBr, MeCl, and methyl iodide (MeI) to CO 2 . Soil bacteria are known to consume MeBr at the ambient tropospheric mixing ratio of around 10 parts per trillion by volume (5). Preliminary experiments with strain IMB-1 indicate that it can oxidize MeBr at these mixing ratios** and is therefore likely to be characteristic of bacteria associated with MeBr uptake by soils. We examined δ 13 C of MeCl, MeBr, and MeI during oxidation by whole-cell suspensions of IMB-1 and CC495 and also the change in δ 13 C values of the three methyl halides during oxidation by the marine strain MB2. In addition, we measured the fractionation of carbon isotopes during formation of methane thiol (MeSH) from MeCl by the purified cobalamin-dependent enzyme, halomethane:bisulphide/halide ion methyltransferase (11) from CC495, to determine whether this initial step in MeCl degradation could account for the observed fractionation by whole cells. Finally, we determined the fractionation associated with the degradation of MeBr during field studies with agricultural soil by monitoring MeBr concentration and δ 13 C of MeBr in the headspace of flux chambers under fumigation conditions.

Proceedings of the National Academy of Sciences of↗

USGS2000

Transformations in methane hydrates

Detailed study of pure methane hydrate in a diamond cell with in situ optical, Raman, and x-ray microprobe techniques reveals two previously unknown structures, structure II and structure H, at high pressures. The structure II methane hydrate at 250 MPa has a cubic unit cell of a = 17.158(2) A?? and volume V = 5051.3(13) A??3; structure H at 600 MPa has a hexagonal unit cell of a = 11.980(2) A??, c = 9.992(3) A??, and V = 1241.9(5) A??3. The compositions of these two investigated phases are still not known. With the effects of pressure and the presence of other gases in the structure, the structure II phase is likely to dominate over the known structure I methane hydrate within deep hydrate-bearing sediments underlying continental margins.

Proceedings of the National Academy of Sciences of↗

USGS2017

Using interviews and biological sign surveys to infer seasonal use of forested and agricultural portions of a human-dominated landscape by Asian elephants in Nepal

Understanding how wide-ranging animals use landscapes in which human use is highly heterogeneous is important for determining patterns of human–wildlife conflict and designing mitigation strategies. Here, we show how biological sign surveys in forested components of a human-dominated landscape can be combined with human interviews in agricultural portions of a landscape to provide a full picture of seasonal use of different landscape components by wide-ranging animals and resulting human–wildlife conflict. We selected Asian elephants ( Elephas maximus ) in Nepal to illustrate this approach. Asian elephants are threatened throughout their geographic range, and there are large gaps in our understanding of their landscape-scale habitat use. We identified all potential elephant habitat in Nepal and divided the potential habitat into sampling units based on a 10 km by 10 km grid. Forested areas within grids were surveyed for signs of elephant use, and local villagers were interviewed regarding elephant use of agricultural areas and instances of conflict. Data were analyzed using single-season and multi-season (dynamic) occupancy models. A single-season occupancy model applied to data from 139 partially or wholly forested grid cells estimated that 0.57 of grid cells were used by elephants. Dynamic occupancy models fit to data from interviews across 158 grid cells estimated that monthly use of non-forested, human-dominated areas over the preceding year varied between 0.43 and 0.82 with a minimum in February and maximum in October. Seasonal patterns of crop raiding by elephants coincided with monthly elephant use of human-dominated areas, and serious instances of human–wildlife conflict were common. Efforts to mitigate human–elephant conflict in Nepal are likely to be most effective if they are concentrated during August through December when elephant use of human-dominated landscapes and human–elephant conflict are most common.

Ethology Ecology and Evolution↗

USGS2012

Kinetics of viral load and erythrocytic inclusion body formation in pacific herring artificially infected with erythrocytic necrosis virus

Viral erythrocytic necrosis (VEN) is a condition that affects marine and anadromous fish species, including herrings and salmonids, in the Atlantic and Pacific oceans. Infection is frequently associated with severe anemia and causes episodic mortality among wild and hatchery fish when accompanied by additional stressors; VEN can be presumptively diagnosed by (1) light microscopic identification of a single characteristic—a round, magenta-colored, 0.8-μm-diameter inclusion body (IB) within the cytoplasm of erythrocytes and their precursors on Giemsa-stained blood films; or (2) observation (via transmission electron microscopy [TEM]) of the causative iridovirus, erythrocytic necrosis virus (ENV), within erythrocytes or their precursors. To better understand the kinetics of VEN, specific-pathogen-free Pacific herring Clupea pallasii were infected with ENV by intraperitoneal injection. At 1, 4, 7, 10, 14, 21, and 28 d postexposure, samples of blood, spleen, and kidney were collected and assessed (1) via light microscopy for the number of intracytoplasmic IBs in blood smears and (2) via TEM for the number of virions within erythrocytes. The mean prevalence of intracytoplasmic IBs in the blood cells increased from 0% at 0–4 d postexposure to 94% at 28 d postexposure. Viral load within circulating red blood cells peaked at 7 d postexposure, fell slightly, and then reached a plateau. However, blood cells observed within the kidney and spleen tissues demonstrated high levels of ENV between 14 and 28 d postexposure. The results indicate that the viral load within erythrocytes does not correlate well with IB prevalence and that the virus can persist in infected fish for more than 28 d.

Journal of Aquatic Animal Health↗

USGS2017

An accurate method for measuring triploidy of larval fish spawns

A standard flow cytometric protocol was developed for estimating triploid induction in batches of larval fish. Polyploid induction treatments are not guaranteed to be 100% efficient, thus the ability to quantify the proportion of triploid larvae generated by a particular treatment helps managers to stock high-percentage spawns and researchers to select treatments for efficient triploid induction. At 3 d posthatch, individual Grass Carp Ctenopharyngodon idella were mechanically dissociated into single-cell suspensions; nuclear DNA was stained with propidium iodide then analyzed by flow cytometry. Following ploidy identification of individuals, aliquots of diploid and triploid cell suspensions were mixed to generate 15 levels (0–100%) of known triploidy ( n = 10). Using either 20 or 50 larvae per level, the observed triploid percentages were lower than the known, actual values. Using nonlinear regression analyses, quadratic equations solved for triploid proportions in mixed samples and corresponding estimation reference plots allowed for predicting triploidy. Thus, an accurate prediction of the proportion of triploids in a spawn can be made by following a standard larval processing and analysis protocol with either 20 or 50 larvae from a single spawn, coupled with applying the quadratic equations or reference plots to observed flow cytometry results. Due to the universality of triploid DNA content being 1.5 times the diploid level and because triploid fish consist of fewer cells than diploids, this method should be applicable to other produced triploid fish species, and it may be adapted for use with bivalves or other species where batch analysis is appropriate.

North American Journal of Aquaculture↗

USGS2007

Evaluating nephrotoxicity of high-molecular-weight organic compounds in drinking water from lignite aquifers

High-molecular-weight organic compounds such as humic acids and/or fulvic acids that are naturally mobilized from lignite beds into untreated drinking-water supplies were suggested as one possible cause of Balkan endemic nephropathy (BEN) and cancer of the renal pelvis. A lab investigation was undertaken in order to assess the nephrotoxic potential of such organic compounds using an in vitro tissue culture model. Because of the infeasibility of exposing kidney tissue to low concentrations of organics for years in the lab, tangential flow ultrafiltration was employed to hyperconcentrate samples suitable for discerning effects in the short time frames necessitated by tissue culture systems. Effects on HK-2 kidney cells were measured using two different cell proliferation assays (MTT and alamarBlue). Results demonstrated that exposure of kidney tissue to high-molecular-weight organics produced excess cell death or proliferation depending on concentration and duration of exposure. Copyright ?? Taylor & Francis Group, LLC.

Journal of Toxicology and Environmental Health, Pa↗

USGS2008

Exposure of xenopus laevis tadpoles to cadmium reveals concentration-dependent bimodal effects on growth and monotonic effects on development and thyroid gland activity

Xenopus laevis were exposed to 0-855 ??g cadmium (Cd)/l (measured concentrations) in FETAX medium from fertilization to 47 days postfertilization. Measurements included embryonic survival and, at 47 days, tadpole survival, snout-vent length, tail length, total length, hindlimb length, weight, Nieuwkoop-Faber (NF) stage of development, initiation of metamorphic climax (??? NF 58), and thyroid follicle cell height. Embryonic and larval survival were unaffected by Cd. Relative to control tadpoles, reduced tail and total length were observed at 0.1- 8 and at 855 ??g Cd/l; and reduced snout-vent length, hindlimb length, and weight were observed at 0.1-1 and at 855 ??g Cd/l. Mean stage of development and rate of initiation of climax were unaffected by Cd at 0-84 ??g/l; however, none of the tadpoles exposed to 855 ??g Cd/l progressed beyond mid-premetamorphosis (NF 51). Thyroid glands with fully formed follicles were observed in all tadpoles ??? NF 49 examined. Follicle cell height was unaffected by Cd at 0-84 ??g/l but it was reduced at 855 ??g/l; in the latter, cell height was reduced even when compared with NF 49-51 tadpoles pooled from the 0 to 84 ??g Cd/l groups. In conclusion, (1) Cd affected tadpole growth in a bimodal pattern with the first and second inhibitory modes at concentrations below and above 84 ??g Cd/l, respectively; (2) exposure to high Cd concentrations (855 ??g/l) reduced thyroid activity and arrested tadpole development at mid-premetamorphosis; and (3) unlike its effect on growth, Cd inhibited tadpole development and thyroid function in a seemingly monotonic pattern.

Toxicological Sciences↗

USGS2004

Isolation and characterization of a rhabdovirus from starry flounder (Platichthys stellatus) collected from the northern portion of Puget Sound, Washington, USA

The initial characterization of a rhabdovirus isolated from a single, asymptomatic starry flounder (Platichthys stellatus) collected during a viral survey of marine fishes from the northern portion of Puget Sound, Washington, USA, is reported. Virions were bullet-shaped and approximately 100 nm long and 50 nm wide, contained a lipid envelope, remained stable for at least 14 days at temperatures ranging from -80 to 5 degrees C and grew optimally at 15 degrees C in cultures of epithelioma papulosum cyprini (EPC) cells. The cytopathic effect on EPC cell monolayers was characterized by raised foci containing rounded masses of cells. Pyknotic and dark-staining nuclei that also showed signs of karyorrhexis were observed following haematoxylin and eosin, May-Grunwald Giemsa and acridine orange staining. PAGE of the structural proteins and PCR assays using primers specific for other known fish rhabdoviruses, including Infectious hematopoietic necrosis virus, Viral hemorrhagic septicemia virus, Spring viremia of carp virus, and Hirame rhabdovirus, indicated that the new virus, tentatively termed starry flounder rhabdovirus (SFRV), was previously undescribed in marine fishes from this region. In addition, sequence analysis of 2678 nt of the amino portion of the viral polymerase gene indicated that SFRV was genetically distinct from other members of the family Rhabdoviridae for which sequence data are available. Detection of this virus during a limited viral survey of wild fishes emphasizes the void of knowledge regarding the diversity of viruses that naturally infect marine fish species in the North Pacific Ocean.

Journal of General Virology↗

USGS2015

Extending the MODPATH algorithm to rectangular unstructured grids

The recent release of MODFLOW-USG, which allows model grids to have irregular, unstructured connections, requires a modification of the particle-tracking algorithm used by MODPATH. This paper describes a modification of the semi-analytical particle-tracking algorithm used by MODPATH that allows it to be extended to rectangular-based unstructured grids by dividing grid cells with multi-cell face connections into sub-cells. The new method will be incorporated in the next version of MODPATH which is currently under development.

Ground Water↗

USGS2011

Time-lapse three-dimensional inversion of complex conductivity data using an active time constrained (ATC) approach

Induced polarization (more precisely the magnitude and phase of impedance of the subsurface) is measured using a network of electrodes located at the ground surface or in boreholes. This method yields important information related to the distribution of permeability and contaminants in the shallow subsurface. We propose a new time-lapse 3-D modelling and inversion algorithm to image the evolution of complex conductivity over time. We discretize the subsurface using hexahedron cells. Each cell is assigned a complex resistivity or conductivity value. Using the finite-element approach, we model the in-phase and out-of-phase (quadrature) electrical potentials on the 3-D grid, which are then transformed into apparent complex resistivity. Inhomogeneous Dirichlet boundary conditions are used at the boundary of the domain. The calculation of the Jacobian matrix is based on the principles of reciprocity. The goal of time-lapse inversion is to determine the change in the complex resistivity of each cell of the spatial grid as a function of time. Each model along the time axis is called a ‘reference space model’. This approach can be simplified into an inverse problem looking for the optimum of several reference space models using the approximation that the material properties vary linearly in time between two subsequent reference models. Regularizations in both space domain and time domain reduce inversion artefacts and improve the stability of the inversion problem. In addition, the use of the time-lapse equations allows the simultaneous inversion of data obtained at different times in just one inversion step (4-D inversion). The advantages of this new inversion algorithm are demonstrated on synthetic time-lapse data resulting from the simulation of a salt tracer test in a heterogeneous random material described by an anisotropic semi-variogram.

Geophysical Journal International↗

USGS2012

Bacterial dominance in subseafloor sediments characterized by methane hydrates

The degradation of organic carbon in subseafloor sediments on continental margins contributes to the largest reservoir of methane on Earth. Sediments in the Andaman Sea are composed of ~ 1% marine-derived organic carbon and biogenic methane is present. Our objective was to determine microbial abundance and diversity in sediments that transition the gas hydrate occurrence zone (GHOZ) in the Andaman Sea. Microscopic cell enumeration revealed that most sediment layers harbored relatively low microbial abundance (10 3 –10 5 cells cm −3 ). Archaea were never detected despite the use of both DNA- and lipid-based methods. Statistical analysis of terminal restriction fragment length polymorphisms revealed distinct microbial communities from above, within, and below the GHOZ, and GHOZ samples were correlated with a decrease in organic carbon. Primer-tagged pyrosequences of bacterial 16S rRNA genes showed that members of the phylum Firmicutes are predominant in all zones. Compared with other seafloor settings that contain biogenic methane, this deep subseafloor habitat has a unique microbial community and the low cell abundance detected can help to refine global subseafloor microbial abundance.

FEMS Microbiology Ecology↗

USGS2005

Stochastic Ground Water Flow Simulation with a Fracture Zone Continuum Model

A method is presented for incorporating the hydraulic effects of vertical fracture zones into two-dimensional cell-based continuum models of ground water flow and particle tracking. High hydraulic conductivity features are used in the model to represent fracture zones. For fracture zones that are not coincident with model rows or columns, an adjustment is required for the hydraulic conductivity value entered into the model cells to compensate for the longer flowpath through the model grid. A similar adjustment is also required for simulated travel times through model cells. A travel time error of less than 8% can occur for particles moving through fractures with certain orientations. The fracture zone continuum model uses stochastically generated fracture zone networks and Monte Carlo analysis to quantify uncertainties with simulated advective travel times. An approach is also presented for converting an equivalent continuum model into a fracture zone continuum model by establishing the contribution of matrix block transmissivity to the bulk transmissivity of the aquifer. The methods are used for a case study in west-central Florida to quantify advective travel times from a potential wetland rehydration site to a municipal supply wellfield. Uncertainties in advective travel times are assumed to result from the presence of vertical fracture zones, commonly observed on aerial photographs as photolineaments.

Ground Water↗

USGS2010

Simulation and analysis of conjunctive use with MODFLOW's farm process

The extension of MODFLOW onto the landscape with the Farm Process (MF-FMP) facilitates fully coupled simulation of the use and movement of water from precipitation, streamflow and runoff, groundwater flow, and consumption by natural and agricultural vegetation throughout the hydrologic system at all times. This allows for more complete analysis of conjunctive use water-resource systems than previously possible with MODFLOW by combining relevant aspects of the landscape with the groundwater and surface water components. This analysis is accomplished using distributed cell-by-cell supply-constrained and demand-driven components across the landscape within " water-balance subregions" comprised of one or more model cells that can represent a single farm, a group of farms, or other hydrologic or geopolitical entities. Simulation of micro-agriculture in the Pajaro Valley and macro-agriculture in the Central Valley are used to demonstrate the utility of MF-FMP. For Pajaro Valley, the simulation of an aquifer storage and recovery system and related coastal water distribution system to supplant coastal pumpage was analyzed subject to climate variations and additional supplemental sources such as local runoff. For the Central Valley, analysis of conjunctive use from different hydrologic settings of northern and southern subregions shows how and when precipitation, surface water, and groundwater are important to conjunctive use. The examples show that through MF-FMP's ability to simulate natural and anthropogenic components of the hydrologic cycle, the distribution and dynamics of supply and demand can be analyzed, understood, and managed. This analysis of conjunctive use would be difficult without embedding them in the simulation and are difficult to estimate a priori. Journal compilation ?? 2010 National Ground Water Association. No claim to original US government works.

Ground Water↗

USGS2011

Approaches to the simulation of unconfined flow and perched groundwater flow in MODFLOW

Various approaches have been proposed to manage the nonlinearities associated with the unconfined flow equation and to simulate perched groundwater conditions using the MODFLOW family of codes. The approaches comprise a variety of numerical techniques to prevent dry cells from becoming inactive and to achieve a stable solution focused on formulations of the unconfined, partially-saturated, groundwater flow equation. Keeping dry cells active avoids a discontinuous head solution which in turn improves the effectiveness of parameter estimation software that relies on continuous derivatives. Most approaches implement an upstream weighting of intercell conductance and Newton-Raphson linearization to obtain robust convergence. In this study, several published approaches were implemented in a stepwise manner into MODFLOW for comparative analysis. First, a comparative analysis of the methods is presented using synthetic examples that create convergence issues or difficulty in handling perched conditions with the more common dry-cell simulation capabilities of MODFLOW. Next, a field-scale three-dimensional simulation is presented to examine the stability and performance of the discussed approaches in larger, practical, simulation settings.

Ground Water↗