Search USGSSearch

SEARCH · Search USGS

Results for “Structure”

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

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

At least 523 records · Page 29Linked to original sources

IASPEI workshop: Seismic modeling of laterally varying structures

During the past 10 years, significant progress has been made in the methods of collection and analysis of seismic reflection and refraction data. This progress has led to the development of new models for the structure and composition of the earth's crust, based on sophisticated analysis of numerous profiles in many areas of geologic importance. The third triannual meeting of the IASPEI (International Association of Seismology and Physics of the Earth's Interior) Commission of Controlled Source Seismology was convened in Park City, Utah, on August 11–17, 1980, to bring together seismologists and geologists to explore and assess the progress of controlled source techniques (controlled sources include explosions, air guns, and Vibroseis‐type sources), and to evaluate its significance in terms of current models of the seismic velocity structure and composition of the crust and upper mantle. Particular attention was paid to the progress and problems in the modeling of two‐ and three‐dimensional structures.

Eos, Transactions, American Geophysical Union

Electrical structure of Newberry Volcano, Oregon

From the interpretation of magnetotelluric, transient electromagnetic, and Schlumberger resistivity soundings, the electrical structure of Newberry Volcano in central Oregon is found to consist of four units. From the surface downward, the geoelectrical units are (1) very resistive, young, unaltered volcanic rock, (2) a conductive layer of older volcanic material composed of altered tuffs, (3) a thick resistive layer thought to be in part intrusive rocks, and (4) a lower-crustal conductor. This model is similar to the regional geoelectrical structure found throughout the Cascade Range. Inside the caldera, the conductive second layer corresponds to the steep temperature gradient and alteration minerals observed in the USGS Newberry 2 test hole. Drill hole information on the south and north flanks of the volcano (test holes GEO N-1 and GEO N-3, respectively) indicates that outside the caldera the conductor is due to alteration minerals (primarily smectite) and not high-temperature pore fluids. On the flanks of Newberry the conductor is generally deeper than inside the caldera, and it deepens with distance from the summit. A notable exception to this pattern is seen just west of the caldera rim, where the conductive zone is shallower than at other flank locations. The volcano sits atop a rise in the resistive layer, interpreted to be due to intrusive rocks. The intrusive material has served as a heat source to produce enhanced hydrothermal alteration and, perhaps in the case of the west-flank anomaly, elevated fluid temperatures. While no public drill hole information is available to confirm this hypothesis, the west-flank anomaly appears to be a good geothermal target. In addition to the possibility that a region on the west side of the volcano could be favorable for prospecting, part of the resistive structure under the center of the volcano could be due to a vapor-dominated environment with temperatures above 300°C. In other parts of the Cascades, pervasive alteration has produced mixed layer clays and zeolites, resulting in low-resistivity anomalies. Low resistivities cannot be assumed to indicate high-temperature pore fluids. The use of electrical methods that measure resistivity as a function of excitation frequency, such as spectral induced polarization, may provide a way of obtaining information about the type and extent of alteration.

Journal of Geophysical Research Solid Earth

Does bird community structure vary with landscape patchiness? A Chihuahuan Desert perspective

During the springs of 1995-1997, we studied birds and landscapes at 70 sites in the Chihuahuan Desert to assess relations between bird community structure and landscape patchiness. Within each of two spatial extents (1-kin and 2-km-radius areas centered on each site), we measured the number of patches of individual land-cover types and the total number of patches of all land-cover types. Mean bird richness, and the mean abundance and probability of occurrence of most bird species were significantly correlated with one or more of these variables. Contrary to evidence from other systems, positive association with landscape patchiness did not increase with the degree to which species were habitat generalists, was not negatively related to body size, and did not differ between neotropical migrants and nonmigrants. For the communities' primary constituent species as a group, the strength of positive and negative associations with patchiness did not differ between landscape extents. Within the 1-km but not the 2-km extent, habitat specialists were more positively and negatively associated with patchiness than were habitat generalists. In general, however, neither habitat breadth, body size, nor migratory status seemed to be responsible for associations with landscape patchiness. Mean richness, and the mean abundance and probability of occurrence of most species were significantly correlated with patchiness within one or both extents, and patchiness of all of the most extensive land-cover types was influential. The simplest explanation for most of the bird-patchiness relations we found is that the associations reflected species-specific habitat needs. Through effects on avian richness, abundance, and occurrence, landscape patchiness affected bird community structure. A more complete understanding of the effects of landscape patchiness on bird community structure is likely to emerge when ecologists study the patchiness of major land-cover types at various spatial extents.

Oikos

New links between the Chicxulub impact structure and the Cretaceous/Tertiary boundary

The 200-km-diameter Chicxulub structure 1–3 in northern Yucatan, Mexico has emerged as the prime candidate for the Cretaceous/Tertiary (K/T) boundary impact crater 3–6 . Concentric geophysical anomalies associated with enigmatic occurrences of Upper Cretaceous breccias and andesitic rocks led Penfield and Camargo 1 to suspect that this structure was a buried impact basin. More recently, the discovery of shocked quartz grains in a Chicxulub breccia 3 , and chemical similarities between Chicxulub rocks and K/T tektite-like glasses 3–6 have been advanced as evidence that the Chicxulub structure is a K/T impact site. Here we present evidence from core samples that Chicxulub is indeed a K/T source crater, and can apparently account for all the evidence of impact distributed globally at the K/T boundary without the need for simultaneous multiple impacts or comet showers. Shocked breccia clasts found in the cores are similar to shocked lithic fragments found worldwide in the K/T boundary ejecta layer 7,8 . The Chicxulub melt rocks that we studied contain anomalously high levels of iridium (up to 13.5 parts per 10 9 ), also consistent with the iridium-enriched K/T boundary layer 9 . Our best estimate of the crystallization age of these melt rocks, as determined by 40 Ar/< 39 Ar analyses, is 65.2 ±0.4 (1σ) Myr, in good agreement with the mean plateau age of 64.98 ± 0.05 Myr recently reported 10 . Furthermore, these melt rocks acquired a remanent magnetization indicating that they cooled during an episode of reversed geomagnetic polarity. The only such episode consistent with 40 Ar/< 39 Ar constraints is chron 29R, which includes the K/T boundary.

northern Yucatan

Trans-crustal structural control of CO2-rich extensional magmatic systems revealed at Mount Erebus Antarctica

Erebus volcano, Antarctica, with its persistent phonolite lava lake, is a classic example of an evolved, CO 2 -rich rift volcano. Seismic studies provide limited images of the magmatic system. Here we show using magnetotelluric data that a steep, melt-related conduit of low electrical resistivity originating in the upper mantle undergoes pronounced lateral re-orientation in the deep crust before reaching shallower magmatic storage and the summit lava lake. The lateral turn represents a structural fault-valve controlling episodic flow of magma and CO 2 vapour, which replenish and heat the high level phonolite differentiation zone. This magmatic valve lies within an inferred, east-west structural trend forming part of an accommodation zone across the southern termination of the Terror Rift, providing a dilatant magma pathway. Unlike H 2 O-rich subduction arc volcanoes, CO 2 -dominated Erebus geophysically shows continuous magmatic structure to shallow crustal depths of < 1 km, as the melt does not experience decompression-related volatile supersaturation and viscous stalling.

Nature Communications

Crystal structure of the heptamolybdate(VI)(paramolybdate) ion, [Mo7O24]6–, in the ammonium and potassium tetrahydrate salts

The crystal structures of the isomorphous salts M I 6 [Mo 7 O 24 ],4H 2 O (M = NH 4 or K) have been refined by three-dimensional X -ray diffraction methods. Unit cell dimensions of these monoclinic compounds, space group P 2 1 / C with Z = 4, are, ammonium salt: a = 8·3934 ± 0·0008, b = 36·1703 ± 0·0045, c = 10·4715 ± 0·0011 Å, β= 115·958°± 0·008°; and potassium salt: a = 8·15 ± 0·02, b = 35·68 ± 0·1, c = 10·30 ± 0·02 Å, β= 115·2°± 02°. By use of multiple Weissenberg patterns, 8197 intensity data (Mo- K α radiation) for the ammonium compound and 2178 (Cu- K α radiation) for the potassium compound were estimated visually and used to test and refine Lindqvist's proposed structure in the space group P 2 1 / c . Lindqvist's structure was confirmed and the full matrix least-squares isotropic refinement led to R 0·076 (ammonium) 0·120 (potassium), with direct unambiguous location of the cations and water molecules in the potassium compound.

Journal of the Chemical Society, Dalton Transactio

Founding events influence genetic population structure of sockeye salmon (Oncorhynchus nerka) in Lake Clark, Alaska

Bottlenecks can have lasting effects on genetic population structure that obscure patterns of contemporary gene flow and drift. Sockeye salmon are vulnerable to bottleneck effects because they are a highly structured species with excellent colonizing abilities and often occupy geologically young habitats. We describe genetic divergence among and genetic variation within spawning populations of sockeye salmon throughout the Lake Clark area of Alaska. Fin tissue was collected from sockeye salmon representing 15 spawning populations of Lake Clark, Six-mile Lake, and Lake Iliamna. Allele frequencies differed significantly at 11 microsatellite loci in 96 of 105 pairwise population comparisons. Pairwise estimates of FST ranged from zero to 0.089. Six-mile Lake and Lake Clark populations have historically been grouped together for management purposes and are geographically proximate. However, Six-mile Lake populations are genetically similar to Lake Iliamna populations and are divergent from Lake Clark populations. The reduced allelic diversity and strong divergence of Lake Clark populations relative to Six-mile Lake and Lake Iliamna populations suggest a bottleneck associated with the colonization of Lake Clark by sockeye salmon. Geographic distance and spawning habitat differences apparently do not contribute to isolation and divergence among populations. However, temporal isolation based on spawning time and founder effects associated with ongoing glacial retreat and colonization of new spawning habitats contribute to the genetic population structure of Lake Clark sock-eye salmon. Nonequilibrium conditions and the strong influence of genetic drift caution against using estimates of divergence to estimate gene flow among populations of Lake Clark sockeye salmon.

Molecular Ecology

Spatial structuring within a reservoir fish population: implications for management

Spatial structuring in reservoir fish populations can exist because of environmental gradients, species-specific behaviour, or even localised fishing effort. The present study investigated whether white crappie exhibited evidence of improved population structure where the northern more productive half of a lake is closed to fishing to provide waterfowl hunting opportunities. Population response to angling was modelled for each substock of white crappie (north (protected) and south (unprotected) areas), the entire lake (single-stock model) and by combining simulations of the two independent substock models (additive model). White crappie in the protected area were more abundant, consisting of larger, older individuals, and exhibited a lower total annual mortality rate than in the unprotected area. Population modelling found that fishing mortality rates between 0.1 and 0.3 resulted in sustainable populations (spawning potential ratios (SPR) >0.30). The population in the unprotected area appeared to be more resilient (SPR > 0.30) at the higher fishing intensities (0.35&ndash;0.55). Considered additively, the whole-lake fishery appeared more resilient than when modelled as a single-panmictic stock. These results provided evidence of spatial structuring in reservoir fish populations, and we recommend model assessments used to guide management decisions should consider those spatial differences in other populations where they exist.

Marine and Freshwater Research

A multi-scaled approach to evaluating the fish assemblage structure within southern Appalachian streams USA.

There is considerable uncertainty about the relative roles of stream habitat and landscape characteristics in structuring stream-fish assemblages. We evaluated the relative importance of environmental characteristics on fish occupancy at the local and landscape scales within the upper Little Tennessee River basin of Georgia and North Carolina. Fishes were sampled using a quadrat sample design at 525 channel units within 48 study reaches during two consecutive years. We evaluated species&ndash;habitat relationships (local and landscape factors) by developing hierarchical, multispecies occupancy models. Modeling results suggested that fish occupancy within the Little Tennessee River basin was primarily influenced by stream topology and topography, urban land coverage, and channel unit types. Landscape scale factors (e.g., urban land coverage and elevation) largely controlled the fish assemblage structure at a stream-reach level, and local-scale factors (i.e., channel unit types) influenced fish distribution within stream reaches. Our study demonstrates the utility of a multi-scaled approach and the need to account for hierarchy and the interscale interactions of factors influencing assemblage structure prior to monitoring fish assemblages, developing biological management plans, or allocating management resources throughout a stream system.

Georgia, North Carolina, Tennessee

Modeling White Sucker (Catostomus commersonii) populations to assess commercial harvest influence on age structure

Commercial harvest of White Suckers Catostomus commersonii for bait in the American Lobster Homarus americanus industry is minimally regulated in Maine and there are concerns as to the influence of increasing harvest. We built a population model using parameters (i.e. age-specific mortality, age at maturity, and size-specific fecundity) from literature and field studies to investigate the theoretical effects of harvest mortality on age structure. Because stock-recruitment relations are poorly characterized for this species, we explored the influence of both Beverton-Holt and Ricker recruitment processes. Our base model closely resembled the empirical age structures reported from three unfished lakes in Maine, with four percent of fish in the modeled spawning run being age-10 or older. We assessed the additive effects of harvest mortality on age structure using the full range of possible mortalities. As expected, increased harvest mortality in the model resulted in a decline and disappearance of older age-classes such that few fish greater than age-10 remained in the population under a realistic harvest mortality scenario. This age-truncation was qualitatively comparable to data from aggregate age distributions reported from three commercially harvested lakes in Maine. Because the loss of older fish may compromise population viability, this model is a valuable guidance tool for managers to craft regulation of this growing fishery.

Maine

Quantifying the contribution of habitats and pathways to a spatially structured population facing environmental change

The consequences of environmental disturbance and management are difficult to quantify for spatially structured populations because changes in one location carry through to other areas as a result of species movement. We develop a metric, G , for measuring the contribution of a habitat or pathway to network-wide population growth rate in the face of environmental change. This metric is different from other contribution metrics, as it quantifies effects of modifying vital rates for habitats and pathways in perturbation experiments. Perturbation treatments may range from small degradation or enhancement to complete habitat or pathway removal. We demonstrate the metric using a simple metapopulation example and a case study of eastern monarch butterflies. For the monarch case study, the magnitude of environmental change influences the ordering of node contribution. We find that habitats within which all individuals reside during one season are the most important to short-term network growth under complete removal scenarios, whereas the central breeding region contributes most to population growth over all but the strongest disturbances. The metric G provides for more efficient management interventions that proactively mitigate impacts of expected disturbances to spatially structured populations. use changes in one location carry through to other areas due to species movement. We develop a metric, G, for measuring the contribution of a habitat or pathway to network-wide population growth rate in the face of environmental change. This metric is different than other contribution metrics as it quantifies effects of modifying vital rates for habitats and pathways in perturbation experiments. Perturbation treatments may range from small degradation or enhancement to complete habitat or pathway removal. We demonstrate the metric using a simple metapopulation example and a case study of eastern monarch butterflies. For the monarch case study, the magnitude of environmental change influences ordering of node contribution. We find that habitats through which all migrants flow are the most important to short-term network growth under complete-removal scenarios. The metric G provides for more efficient management interventions that proactively mitigate impacts of expected disturbances to spatially structured populations.

American Naturalist

Stress- and structure-controlled anisotropy in a region of complex faulting—Yuha Desert, California

We examine shear velocity anisotropy in the Yuha Desert, California using aftershocks of the 2010 M7.2 El Mayor-Cucapah earthquake. The Yuha Desert is underlain by a complex network of right- and left-lateral conjugate faults, some of which experienced triggered slip during the El Mayor-Cucapah earthquake. An automated method that implements multiple measurement windows and a range of bandpass filters is used to estimate the fast direction ( ϕ ) and delay time ( &delta;t ) of the split shear waves. We find an average ϕ oriented approximately north&ndash;south suggesting it is primarily controlled by the regional maximum compressive stress direction. However, the spatial variability in ϕ reveals that the fault structures that underlie the Yuha Desert also influence the measured splitting parameters. We infer that the northeast- and northwest-oriented ϕ reflect shear fabric subparallel to the conjugate fault structures. We do not observe a simple correlation between &delta;t and hypocentral distance. Instead, the observed spatial variation in &delta;t suggests that near-source variation in anisotropic strength may be equal to or more important than effects local to the station. No temporal variation in splitting parameters is observed during the 70-day period following the main shock. In this region of complex faulting, we observe a spatially variable pattern of anisotropy that is both stress- and structure-controlled. This study suggests that shear fabric can form even along short, discontinuous fault strands with minimal offset.

California

Range-wide genetic analysis of an endangered bumble bee (Bombus affinis, Hymenoptera: Apidae) reveals population structure, isolation by distance, and low colony abundance

Declines in bumble bee species range and abundances are documented across multiple continents and have prompted the need for research to aid species recovery and conservation. The rusty patched bumble bee ( Bombus affinis ) is the first federally listed bumble bee species in North America. We conducted a range-wide population genetics study of B. affinis from across all extant conservation units to inform conservation efforts. To understand the species’ vulnerability and help establish recovery targets, we examined population structure, patterns of genetic diversity, and population differentiation. Additionally, we conducted a site-level analysis of colony abundance to inform prioritizing areas for conservation, translocation, and other recovery actions. We find substantial evidence of population structuring along an east-to-west gradient. Putative populations show evidence of isolation by distance, high inbreeding coefficients, and a range-wide male diploidy rate of ~15%. Our results suggest the Appalachians represent a genetically distinct cluster with high levels of private alleles and substantial differentiation from the rest of the extant range. Site-level analyses suggest low colony abundance estimates for B. affinis compared to similar datasets of stable, co-occurring species. These results lend genetic support to trends from observational studies, suggesting that B. affinis has undergone a recent decline and exhibit substantial spatial structure. The low colony abundances observed here suggest caution in overinterpreting the stability of populations even where B. affinis is reliably detected interannually. These results help delineate informed management units, provide context for the potential risks of translocation programs, and help set clear recovery targets for this and other threatened bumble bee species.

Iowa, Illinois, Indiana, Kentucky, Michigan, Minne

Article navigation zonal distribution of variations in structural state of alkali feldspar within the Rader Creek pluton, Boulder Batholith, Montana

The granodioritic Rader Creek pluton of the composite Boulder batholith contains microperthitic alkali feldspar of bulk composition Or 65 to Or 86 with a structurally variable potassic phase. Complete cell parameters, 2V measurements, and bulk composition are given for 11 feldspar samples. The 131 and 131 reflections for these and 58 additional samples show the following structural types in the potassic phase: orthoclase only; orthoclase with subordinate maximum or near-maximum microcline (obliquity = 0.75–1.00); orthoclase with subordinate intermediate microcline (obliquity = 0.64–0.71); and intermediate microcline (obliquity = 0.56–0.77) with subordinate orthoclase. Within the pluton different feldspar structural types occur in zones whose boundaries are approximately parallel to contacts with younger intrusive rocks cutting the Rader Creek pluton but are, in places, nearly perpendicular to zonation within the pluton defined by rock composition. In general, the orthoclase zone is closest to the contact with younger intrusives; the intermediate microcline zone is the most distant. Bulk compositions of alkali feldspar are more potassic in the orthoclase zone than elsewhere. The data suggest a complex history for the alkali feldspar, involving at least two stages: 1. Exsolution and partial inversion of orthoclase to intermediate microcline during cooling of the Rader Creek pluton; 2. Transformation of the intermediate-microcline assemblage to orthoclase during reheating of the pluton at the time of intrusion of younger plutons of the batholith. The transitional stage in this transformation is characterized by orthoclase co-existing with subordinate microcline, whose obliquity usually approaches that of maximum microcline.

Montana

Know what you don't know: Embracing state uncertainty in disease-structured multistate models

Hidden Markov models (HMMs) are broadly applicable hierarchical models that derive their utility from separating state processes from observation processes yielding the data. Multistate models such as mark–recapture and dynamic multistate occupancy models are HMMs frequently used in ecology. In their early formulations, states, such as pathogen infection status, were assumed to be perfectly observed without ambiguity. However, state uncertainty is a pervasive feature of many ecological studies, and multievent models were developed to explicitly account for it. We developed a novel extended multievent mark–recapture model that incorporates state uncertainty at multiple levels of detection. Using a disease-structured example, both false negative and false positive state assignment errors are modelled at two levels of state assignment—the pathogen sampling process and the diagnostic process that samples are subjected to. We additionally describe methods to jointly model infection intensity to integrate heterogeneity in ecological parameters, such as mortality and infection dynamics, and the pathogen detection processes. We provide code to simulate and analyse datasets with various underlying ecological processes and fit our model to a mark–recapture dataset of Mixophyes fleayi (Fleay's barred frog) infected with the amphibian chytrid fungus ( Batrachochytrium dendrobatidis , Bd ). In our case study, we found evidence for various state assignment errors: the sampling protocol performed poorly in detecting Bd , pathogen detection was highly dependent on infection intensity and false positives were non-negligible. Incorporating state uncertainty yielded significantly higher estimates of infection prevalence and 4–5 times lower rates of infection state transitions compared to those obtained from a traditional multistate model. Our results highlight that incorporating state assignment errors improves inference on the ecological process, especially when sensitivity and specificity of the state assignment processes are low. The general model structure can be applied to other HMMs, providing a foundation for modelling state uncertainty in related models. For disease-structured multistate models, we recommend conducting robust design surveys and collecting samples during each capture event to facilitate incorporating pathogen detection errors.

Methods in Ecology and Evolution

Disease and community structure: white-nose syndrome alters spatial and temporal niche partitioning in sympatric bat species

Aim Emerging infectious diseases present a major perturbation with apparent direct effects such as reduced population density, extirpation and/or extinction. Comparatively less is known about the potential indirect effects of disease that likely alter community structure and larger ecosystem function. Since 2006, white-nose syndrome (WNS) has resulted in the loss of over 6 million hibernating bats in eastern North America. Considerable evidence exists concerning niche partitioning in sympatric bat species in this region, and the unprecedented, rapid decline in multiple species following WNS may provide an opportunity to observe a dramatic restructuring of the bat community. Location We conducted our study at Fort Drum Army Installation in Jefferson and Lewis counties, New York, USA, where WNS first impacted extant bat species in winter 2007–2008. Methods Acoustical monitoring during 2003–2011 allowed us to test the hypothesis that spatial and temporal niche partitioning by bats was relaxed post-WNS. Results We detected nine bat species pre- and post-WNS. Activity for most bat species declined post-WNS. Dramatic post-WNS declines in activity of little brown bat ( Myotis lucifugus , MYLU), formerly the most abundant bat species in the region, were associated with complex, often species-specific responses by other species that generally favoured increased spatial and temporal overlap with MYLU. Main conclusions In addition to the obvious direct effects of disease on bat populations and activity levels, our results provide evidence that disease can have cascading indirect effects on community structure. Recent occurrence of WNS in North America, combined with multiple existing stressors, is resulting in dramatic shifts in temporal and spatial niche partitioning within bat communities. These changes might influence long-term population viability of some bat species as well as broader scale ecosystem structure and function.

New York

Mixed-stock analysis in the age of genomics: Rapture genotyping enables evaluation of stock-specific exploitation in a freshwater fish population with weak genetic structure

Mixed-stock analyses using genetic markers have informed fisheries management in cases where strong genetic differentiation occurs among local spawning populations, yet many fisheries are supported by multiple spawning stocks that are weakly differentiated. Freshwater fisheries exemplify this problem, with many harvested populations supported by multiple stocks of young evolutionary age and that are isolated across small spatial scales. As a result, attempts to conduct genetic mixed-stock analyses of inland fisheries have often been unsuccessful. Advances in genomic sequencing now offer the ability to discriminate among populations with weak population structure, by providing the necessary resolution to conduct mixed-stock assignment among previously indistinguishable stocks. We demonstrate the use of genomic data to conduct a mixed-stock analysis of Lake Erie's commercial and recreational walleye (Sander vitreus) fisheries and estimate the relative harvest of weakly differentiated stocks (pairwise FST < 0.01). We used RAD-capture (Rapture) to sequence and genotype individuals at 12,081 loci that had been previously determined to be capable of discriminating between western and eastern basin stocks with 95% reassignment accuracy. An outcome not possible in the past with microsatellite markers. Genetic assignment of 1,075 fish harvested from recreational and commercial fisheries in the eastern basin indicated that western basin stocks constituted the majority of individuals harvested during peak walleye fishing season (July – September). Composition of harvest changed seasonally, with eastern basin fish comprising much of the early season harvest (May – June). Clear spatial structure in stock-specific harvest existed; more easterly sites contained more individuals of east basin origin than did westerly sites. Our study provides important stock contribution estimates for Lake Erie fishery management and demonstrates the power of genomic data to facilitate mixed-stock analysis in exploited fish populations with weak population structure or limited existing genetic resources.

Michigan, New York, Ohio, Pennsylvania

A semi-structured MODFLOW-USG model to evaluate local water sources to wells for decision support

In order to better represent the configuration of the stream network and simulate local groundwater-surface water interactions, a version of MODFLOW with refined spacing in the topmost layer was applied to a Lake Michigan Basin (LMB) regional groundwater-flow model developed by the U.S. Geological. Regional MODFLOW models commonly use coarse grids over large areas; this coarse spacing precludes model application to local management issues (e.g., surface-water depletion by wells) without recourse to labor-intensive inset models. Implementation of an unstructured formulation within the MODFLOW framework (MODFLOW-USG) allows application of regional models to address local problems. A &ldquo;semi-structured&rdquo; approach (uniform lateral spacing within layers, different lateral spacing among layers) was tested using the LMB regional model. The parent 20-layer model with uniform 5000-foot (1524-m) lateral spacing was converted to 4 layers with 500-foot (152-m) spacing in the top glacial (Quaternary) layer, where surface water features are located, overlying coarser resolution layers representing deeper deposits. This semi-structured version of the LMB model reproduces regional flow conditions, whereas the finer resolution in the top layer improves the accuracy of the simulated response of surface water to shallow wells. One application of the semi-structured LMB model is to provide statistical measures of the correlation between modeled inputs and the simulated amount of water that wells derive from local surface water. The relations identified in this paper serve as the basis for metamodels to predict (with uncertainty) surface-water depletion in response to shallow pumping within and potentially beyond the modeled area, see Fienen et al. (2015a).

Illinois, Indiana, Michigan, Wisconsin