Search USGS⌕ Search

SEARCH · Search USGS

Results for “Journal of Maps”

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 793 records · Page 44Linked to original sources
USGS2009

Growth trishear model and its application to the Gilbertown graben system, southwest Alabama

Fault-propagation folding associated with an upward propagating fault in the Gilbertown graben system is revealed by well-based 3-D subsurface mapping and dipmeter analysis. The fold is developed in the Selma chalk, which is an oil reservoir along the southern margin of the graben. Area-depth-strain analysis suggests that the Cretaceous strata were growth units, the Jurassic strata were pregrowth units, and the graben system is detached in the Louann Salt. The growth trishear model has been applied in this paper to study the evolution and kinematics of extensional fault-propagation folding. Models indicate that the propagation to slip (p/s) ratio of the underlying fault plays an important role in governing the geometry of the resulting extensional fault-propagation fold. With a greater p/s ratio, the fold is more localized in the vicinity of the propagating fault. The extensional fault-propagation fold in the Gilbertown graben is modeled by both a compactional and a non-compactional growth trishear model. Both models predict a similar geometry of the extensional fault-propagation fold. The trishear model with compaction best predicts the fold geometry. ?? 2008 Elsevier Ltd. All rights reserved.

Journal of Structural Geology↗
USGS2018

Rapid, remote assessment of Hurricane Matthew impacts using four-dimensional structure-from-motion photogrammetry

Timely assessment of coastal landforms and structures after storms is important for evaluating storm impacts, aiding emergency response and restoration, and initializing and assessing morphological models. Four-dimensional multiview photogrammetry, also known as structure from motion (4D SfM), provides a method for generating three-dimensional reconstructions of landscapes at two times (before and after events) using only photos and existing information for ground control points. Here, these techniques were applied using National Oceanic and Atmospheric Administration (NOAA)-obtained oblique aerial photos taken before (2015) and immediately after Hurricane Matthew (2016) to assess coastal changes near Matanzas, Florida. This work demonstrated that 3D digital elevation models can be constructed within 48 hours of postevent photo collection without on-site ground control measurements. One advantage of timely SfM elevation-change assessments is that they avoid confusion of storm impacts with changes that occur after the event but before LIDAR surveys can be performed. The accuracy and precision of the 4D SfM maps were assessed a posteriori using the first-available LIDAR data, which were collected more than a month after the hurricane, and 11 independent ground-truth survey points measured a week after the hurricane. Horizontal coordinates of the 4D SfM reconstruction were biased by an average of 0.79 m (0.83 m root-mean-square difference; RMSD) compared with the ground-truth points, but vertical elevations were more accurate. They were biased from the LIDAR by −0.09 to −0.25 m, with ∼0.20 m RMSD from both the LIDAR data and five ground-truth points with good vertical positioning and 0.25 m RMSD from LIDAR data along a 60-m stretch of pavement. This level of precision was sufficient to quantify geomorphological change that was often in excess of 1 m. The methodology is conducive for rapid assessment of changes along short stretches (tens of kilometers) of coast with modest resources and could be scaled up for larger regions.

Florida↗
USGS2017

Spatial heterogeneity of within-stream methane concentrations

Streams, rivers, and other freshwater features may be significant sources of CH 4 to the atmosphere. However, high spatial and temporal variabilities hinder our ability to understand the underlying processes of CH 4 production and delivery to streams and also challenge the use of scaling approaches across large areas. We studied a stream having high geomorphic variability to assess the underlying scale of CH 4 spatial variability and to examine whether the physical structure of a stream can explain the variation in surface CH 4 . A combination of high-resolution CH 4 mapping, a survey of groundwater CH 4 concentrations, quantitative analysis of methanogen DNA, and sediment CH 4 production potentials illustrates the spatial and geomorphic controls on CH 4 emissions to the atmosphere. We observed significant spatial clustering with high CH 4 concentrations in organic-rich stream reaches and lake transitions. These sites were also enriched in the methane-producing mcrA gene and had highest CH 4 production rates in the laboratory. In contrast, mineral-rich reaches had significantly lower concentrations and had lesser abundances of mcrA . Strong relationships between CH 4 and the physical structure of this aquatic system, along with high spatial variability, suggest that future investigations will benefit from viewing streams as landscapes, as opposed to ecosystems simply embedded in larger terrestrial mosaics. In light of such high spatial variability, we recommend that future workers evaluate stream networks first by using similar spatial tools in order to build effective sampling programs.

Journal of Geophysical Research G: Biogeosciences↗
USGS2012

Gravity survey of the Mount Toondina impact structure, South Australia

Gravity and seismic reflection data, together with geologic mapping, indicate that the Mount Toondina feature in South Australia is best interpreted as an eroded 4-km-diameter impact structure consisting of a ring structural depression surrounding a pronounced central uplift. Beds at the center of the structure within the central uplift have been raised as much as 200 m from depth and deformed by convergent flow. Seismic reflection data indicate that deformation extends to depths of only ≈800 m; at greater depths the reflectors are nearly flat lying, indicating little or no deformation. Gravity data show residual anomalies of +1.0 mGal coincident with the central uplift and a −0.5 Mgal low associated with the ring structural depression. Modeling of the gravity data indicates that relatively high-density material occurs within the central uplift, whereas the ring depression is filled with low-density material. The deformation at Mount Toondina is typical of a complex impact crater; the 4-km diameter is consistent with the expected threshold size for complex craters formed in weak to moderate strength sedimentary rocks.

South Australia↗
USGS1973

Critically refracted waves in a spherically symmetric radially heterogeneous Earth model

A theoretical analysis of acoustic waves refracted by a spherical boundary across which velocity and density increase abruptly and below which velocity and density may either increase or decrease continuously with depth is formulated in terms of waves generated at a harmonic point source and scattered by a radially heterogeneous spherical body. Through the application of an Earth-flattening transformation on the radial solution and the Watson transform on the sum over eigenfunctions, the solution to the spherical problem for high frequencies is expressed as an integral for the corresponding half-space problem in which the effect of boundary curvature maps into an effective positive velocity gradient with depth. The results of both analytical and numerical evaluation of this integral can be summarized as follows for body waves in the crust and upper mantle: (1) In the special case of a critical velocity gradient (a gradient equal and opposite to the effective curvature gradient), waves interacting with the boundary at the critical angle of incidence have the same form as the classical head wave for flat, homogeneous layers. (2) For gradients more negative than critical, the amplitude of waves incident at the critical angle decay more rapidly with distance than the classical head wave. (3) For gradients that are positive, null, and less negative than critical, the amplitude of waves near the critical angle decays less rapidly with distance than the classical head wave, and at sufficiently large distances, the refracted wave field can be adequately described in terms of ray-theoretical diving waves. At intermediate distances from the critical point, the spectral amplitude of the refracted wave is scalloped due to multiple diving wave interference.

Geophysical Journal of the Royal Astronomical Soci↗
USGS2013

Survival and home-range size of Northern Spotted Owls in southwestern Oregon

In the Klamath province of southwestern Oregon, Northern Spotted Owls (Strix occidentalis caurina) occur in complex, productive forests that historically supported frequent fires of variable severity. However, little is known about the relationships between Spotted Owl survival and home-range size and the characteristics of fire-prone, mixed-conifer forests of the Klamath province. Thus, the objectives of this study were to estimate monthly survival rates and home-range size in relation to habitat characteristics for Northern Spotted Owls in southwestern Oregon. Home-range size and survival of 15 Northern Spotted Owls was monitored using radiotelemetry in the Ashland Ranger District of the Rogue River–Siskiyou National Forest from September 2006 to October 2008. Habitat classes within Spotted Owl home ranges were characterized using a remote-sensed vegetation map of the study area. Estimates of monthly survival ranged from 0.89 to 1.0 and were positively correlated with the number of late-seral habitat patches and the amount of edge, and negatively correlated with the mean nearest neighbor distance between late-seral habitats. Annual home-range size varied from to 189 to 894 ha ( x = 576; SE = 75), with little difference between breeding and nonbreeding home ranges. Breeding-season home-range size increased with the amount of hard edge, and the amount of old and mature forest combined. Core area, annual and nonbreeding season home-range sizes all increased with increased amounts of hard edge, suggesting that increased fragmentation is associated with larger core and home-range sizes. Although no effect of the amount of late-seral stage forest on either survival or home-range size was detected, these results are the first to concurrently demonstrate increased forest fragmentation with decreased survival and increased home-range size of Northern Spotted Owls.

Oregon↗
USGS2019

Insights into pāhoehoe lava emplacement using visible and thermal structure-from-motion photogrammetry

We present the evolution over 3 months of a 2016–2017 pāhoehoe flow at Kīlauea as it changed from a narrow sheet flow into a compound lava field fed by a stable system of tubes. The portion of the flow located on Kīlauea's coastal plain was characterized using helicopter‐based visible and thermal structure‐from‐motion photogrammetry to construct a series of georeferenced digital surface models and thermal maps on eight different days. Results reveal key influences on the emplacement and evolution of such long‐lived pāhoehoe flows. This region of the flow grew by ~12 × 10 6 m 3 with a near‐constant time‐average discharge rate of 1.2–2.7 m 3 /s. The development of two tube systems is captured and shows an initial nascent tube enhanced by a narrow topographic confinement, which later inflated and created a topographic inversion that modulated the emplacement of a second flow lobe with its own tube system. The analysis of breakouts at various stages of the field's life suggests that the evolution of the thermal and morphological properties of the flow surface reflect its maturity. Thermal properties of breakouts were used to expand the empirical relationship of breakout cooling to longer timescales. This study contributes to the long‐term development and validation of more accurate predictive models for pāhoehoe, required during the management of long‐lasting lava flow crises in Hawai'i and elsewhere.

Hawaii↗
USGS2014

Precise determination of δ 88 Sr in rocks, minerals, and waters by double-spike TIMS: A powerful tool in the study of chemical, geologic, hydrologic and biologic processes

We present strontium isotopic ( 88 Sr/ 86 Sr and 87 Sr/ 86 Sr) results obtained by 87 Sr– 84 Sr double spike thermal ionization mass-spectrometry (DS-TIMS) for several standards as well as natural water samples and mineral samples of abiogenic and biogenic origin. The detailed data reduction algorithm and a user-friendly Sr-specific stand-alone computer program used for the spike calibration and the data reduction are also presented. Accuracy and precision of our δ 88 Sr measurements, calculated as permil (‰) deviations from the NIST SRM-987 standard, were evaluated by analyzing the NASS-6 seawater standard, which yielded δ 88 Sr = 0.378 ± 0.009‰. The first DS-TIMS data for the NIST SRM-607 potassium feldspar standard and for several US Geological Survey carbonate, phosphate, and silicate standards (EN-1, MAPS-4, MAPS-5, G-3, BCR-2, and BHVO-2) are also reported. Data obtained during this work for Sr-bearing solids and natural waters show a range of δ 88 Sr values of about 2.4‰, the widest observed so far in terrestrial materials. This range is easily resolvable analytically because the demonstrated external error (±SD, standard deviation) for measured δ 88 Sr values is typically ≤0.02‰. It is shown that the “true” 87 Sr/ 86 Sr value obtained by the DS-TIMS or any other external normalization method combines radiogenic and mass-dependent mass-fractionation effects, which cannot be separated. Therefore, the “true” 87 Sr/ 86 Sr and the δ 87 Sr parameter derived from it are not useful isotope tracers. Data presented in this paper for a wide range of naturally occurring sample types demonstrate the potential of the δ 88 Sr isotope tracer in combination with the traditional radiogenic 87 Sr/ 86 Sr tracer for studying a variety of biological, hydrological, and geological processes.

Journal of Analytical Atomic Spectrometry↗
USGS2007

Evaluation of an experimental LiDAR for surveying a shallow, braided, sand-bedded river

Reaches of a shallow (<1.0m), braided, sand-bedded river were surveyed in 2002 and 2005 with the National Aeronautics and Space Administration's Experimental Advanced Airborne Research LiDAR (EAARL) and concurrently with conventional survey-grade, real-time kinematic, global positioning system technology. The laser pulses transmitted by the EAARL instrument and the return backscatter waveforms from exposed sand and submerged sand targets in the river were completely digitized and stored for postflight processing. The vertical mapping accuracy of the EAARL was evaluated by comparing the ellipsoidal heights computed from ranging measurements made using an EAARL terrestrial algorithm to nearby (<0.5m apart) ground-truth ellipsoidal heights. After correcting for apparent systematic bias in the surveys, the root mean square error of these heights with the terrestrial algorithm in the 2002 survey was 0.11m for the 26 measurements taken on exposed sand and 0.18m for the 59 measurements taken on submerged sand. In the 2005 survey, the root mean square error was 0.18m for 92 measurements taken on exposed sand and 0.24m for 434 measurements on submerged sand. In submerged areas the waveforms were complicated by reflections from the surface, water column entrained turbidity, and potentially the riverbed. When applied to these waveforms, especially in depths greater than 0.4m, the terrestrial algorithm calculated the range above the riverbed. A bathymetric algorithm has been developed to approximate the position of the riverbed in these convolved waveforms and preliminary results are encouraging. ?? 2007 ASCE.

Journal of Hydraulic Engineering↗
USGS2013

Geologic controls on regional and local erosion rates of three northern Gulf of Mexico barrier-island systems

The stratigraphy of sections of three barrier island systems in the northeastern Gulf of Mexico (Apalachicola, Mississippi, and Chandeleur) have been mapped using geophysical and coring techniques to assess the influence of geologic variations in barrier lithosomes and adjoining inner shelf deposits on long-term rates of shoreline change at regional and local scales. Regional scale was addressed by comparing average geologic characteristics of the three areas with mean shoreline-change rates for each area. Regionally, differences in sand volume contained within the part of the barrier lithosome above sea level, sand volume on the inner shelf, and to a lesser extent, sediment grain size correlate with shoreline change rates. Larger sand volumes and coarser grain sizes are found where erosion rates are lower. Local scale was addressed by comparing alongshore variations in barrier island and inner shelf geology with alongshore variations in shoreline change. Locally, long-term shoreline change rates are highest directly shoreward of paleovalleys exposed on the inner shelf. While geology is not the sole explanation for observed differences in shoreline change along these three coastal regions, it is a significant contributor to change variability.

Journal of Coastal Research↗
USGS2016

Regional patterns of total nitrogen concentrations in the National Rivers and Streams Assessment

Patterns of nitrogen (N) concentrations in streams sampled by the National Rivers and Streams Assessment (NRSA) were examined semiquantitatively to identify regional differences in stream N levels. The data were categorized and analyzed by watershed size classes to reveal patterns of the concentrations that are consistent with the spatial homogeneity in natural and anthropogenic characteristics associated with regional differences in N levels. Ecoregions and mapped information on human activities including agricultural practices were used to determine the resultant regions. Marked differences in N levels were found among the nine aggregations of ecoregions used to report the results of the NRSA. We identified distinct regional patterns of stream N concentrations within the reporting regions that are associated with the characteristics of specific Level III ecoregions, groups of Level III ecoregions, groups of Level IV ecoregions, certain geographic characteristics within ecoregions, and/or particular watershed size classes. We described each of these regions and illustrated their areal extent and median and range in N concentrations. Understanding the spatial variability of nutrient concentrations in flowing waters and the apparent contributions that human and nonhuman factors have on different sizes of streams and rivers is critical to the development of effective water quality assessment and management plans. This semi-quantitative analysis is also intended to identify areas within which more detailed quantitative work can be conducted to determine specific regional factors associated with variations in stream N concentrations.

Journal of Soil and Water Conservation↗
USGS1977

Land use and land cover changes: A framework for monitoring

Information on the rate and kind of change in the use of land resources is essential to-the proper planning, management, and regulation of the use of such resources. Starting in 1975, the U.S. Geological Survey has been engaged in mapping and inventory of land use and land cover at a scale of 1:250000 with plans to change to new 1:100000 base maps as they become available. Such scales are-appropriate for acquiring benchmark or baseline data on land use and land cover on a nationwide basis with a reasonable degree of standardization and accuracy, within a practical time frame of 6 to 7 years, and at a reasonable level of detail in categorization. Changes in land use and land cover occur at highly variable rates from place to place. Thus variations, in the intervals between updates of land use and land cover data should be carefully evaluated in order to provide an appropriate frame of reference for the systematic updating of such data. A methodological framework is needed for monitoring changes in land use and land cover that will be timely, relatively inexpensive, and appropriate for widely varying needs at national, interstate-regional, state, multicounty, county, and city levels. In this article a partial framework for monitoring land use and land cover changes is provided with opportunity to extend its use to greater levels of detail in categorization and shorter intervals of frequency when needed for selected areas.

Journal of Research of the U.S. Geological Survey↗
USGS2015

Glaciological and marine geological controls on terminus dynamics of Hubbard Glacier, southeast Alaska

Hubbard Glacier, located in southeast Alaska, is the world's largest non-polar tidewater glacier. It has been steadily advancing since it was first mapped in 1895; occasionally, the advance creates an ice or sediment dam that blocks a tributary fjord (Russell Fiord). The sustained advance raises the probability of long-term closure in the near-future, which will strongly impact the ecosystem of Russell Fiord and the nearby community of Yakutat. Here, we examine a 43-year record of flow speeds and terminus position to understand the large-scale dynamics of Hubbard Glacier. Our long-term record shows that the rate of terminus advance has increased slightly since 1895, with the exception of a slowed advance between approximately 1972 and 1984. The short-lived closure events in 1986 and 2002 were not initiated by perturbations in ice velocity or environmental forcings, but were likely due to fluctuations in sedimentation patterns at the terminus. This study points to the significance of a coupled system where short-term velocity fluctuations and morainal shoal development control tidewater glacier terminus position.

Alaska↗
USGS1992

Precipitation estimation in mountainous terrain using multivariate geostatistics. Part I: structural analysis

Values of average annual precipitation (AAP) are desired for hydrologic studies within a watershed containing Yucca Mountain, Nevada, a potential site for a high-level nuclear-waste repository. Reliable values of AAP are not yet available for most areas within this watershed because of a sparsity of precipitation measurements and the need to obtain measurements over a sufficient length of time. To estimate AAP over the entire watershed, historical precipitation data and station elevations were obtained from a network of 62 stations in southern Nevada and southeastern California. Multivariate geostatistics (cokriging) was selected as an estimation method because of a significant (p = 0.05) correlation of r = .75 between the natural log of AAP and station elevation. A sample direct variogram for the transformed variable, TAAP = ln [(AAP) 1000], was fitted with an isotropic, spherical model defined by a small nugget value of 5000, a range of 190 000 ft, and a sill value equal to the sample variance of 163 151. Elevations for 1531 additional locations were obtained from topographic maps to improve the accuracy of cokriged estimates. A sample direct variogram for elevation was fitted with an isotropic model consisting of a nugget value of 5500 and three nested transition structures: a Gaussian structure with a range of 61 000 ft, a spherical structure with a range of 70 000 ft, and a quasi-stationary, linear structure. The use of an isotropic, stationary model for elevation was considered valid within a sliding-neighborhood radius of 120 000 ft. The problem of fitting a positive-definite, nonlinear model of coregionalization to an inconsistent sample cross variogram for TAAP and elevation was solved by a modified use of the Cauchy-Schwarz inequality. A selected cross-variogram model consisted of two nested structures: a Gaussian structure with a range of 61 000 ft and a spherical structure with a range of 190 000 ft. Cross validation was used for model selection and for comparing the geostatistical model with six alternate estimation methods. Multivariate geostatistics provided the best cross-validation results.

Nevada↗
USGS2025

Multi-lithofacies alluvial characterization via airborne electromagnetic-borehole fusion using ordinary interval kriging and geologic constraints

Airborne electromagnetic (AEM) survey provides extensive spatial coverage and detailed resolution at the near surface and can be used to develop hydrogeological models. However, utilization of AEM data is not straightforward because AEM resistivity is an indirect measurement for inferring sediment types. This study develops an ordinary interval kriging (OIK) algorithm and a resistivity-to-multi-lithofacies (R2ML) data fusion workflow for multi-lithofacies alluvial characterization. OIK utilizes irregular interval data to construct three-dimensional (3D) resistivity fields from one-dimensional inverted AEM resistivity models. The R2ML workflow maps the resistivity field generated from OIK into a multi-facies lithological model, incorporating geologic constraints derived from well logs and geological observations. The numerical and real-world cases demonstrate that OIK is computationally efficient, accounts for 3D anisotropy, and minimizes the smoothing effect, thereby preserving resistivity contrasts and reducing interpolation uncertainty. The methodology is applied to lithologic characterization of the Mississippi River Valley alluvial aquifer (MRVA) in the Shellmound area, Mississippi, U.S. A frequency-domain AEM survey was conducted to support groundwater studies for the managed aquifer recharge (MAR) to the MRVA. The resulting lithological model, including four types of lithofacies—clay, very fine sand, fine-medium sands, and graveliferous sands, illustrates the geomorphological processes of the MRVA and implies potential MAR. The alignment between the lithological model and existing geological and hydrogeological investigations demonstrates that OIK and R2ML workflow effectively capture the subsurface architecture of the MRVA. The methods have broad applicability for characterizing alluvial aquifers through AEM-borehole data fusion, supporting sustainable groundwater management.

Mississippi↗
USGS2004

Granular avalanches across irregular three-dimensional terrain: 2. Experimental tests

Scaling considerations indicate that miniature experiments can be used to test models of granular avalanches in which the effects of intergranular fluid and cohesion are negligible. To test predictions of a granular avalanche model described in a companion paper, we performed bench top experiments involving avalanches of dry sand across irregular basal topography that mimicked the complexity of natural terrain. The experiments employed a novel method of laser-assisted cartography to map the three-dimensional morphology of rapidly moving avalanches, thereby providing high-resolution data for comparison with model output. Model input consisted of two material properties (angles of internal and basal Coulomb friction of the sand), which were measured in independent tests, and of initial and boundary conditions that characterized the geometry of the experimental apparatus. Experimental results demonstrate that the model accurately predicts not only the gross behavior but also many details of avalanche motion, from initiation to deposition. We attribute this accuracy to a mathematical and computational formulation that conserves mass and momentum in three-dimensional physical space and satisfies the Coulomb equation in three-dimensional stress space. Our results support the hypothesis that a Coulomb proportionality between shear and normal stresses applies in moderately rapid granular flows and that complicated constitutive postulates are unnecessary if momentum conservation is strictly enforced in continuum avalanche models. Furthermore, predictions of our Coulomb continuum model contrast with those of a Coulomb point mass model, illustrating the importance of multidimensional modeling and model testing.

Journal of Geophysical Research Earth Surface↗
USGS2014

The role of defensible space for residential structure protection during wildfires

With the potential for worsening fire conditions, discussion is escalating over how to best reduce effects on urban communities. A widely supported strategy is the creation of defensible space immediately surrounding homes and other structures. Although state and local governments publish specific guidelines and requirements, there is little empirical evidence to suggest how much vegetation modification is needed to provide significant benefits. We analysed the role of defensible space by mapping and measuring a suite of variables on modern pre-fire aerial photography for 1000 destroyed and 1000 surviving structures for all fires where homes burned from 2001 to 2010 in San Diego County, CA, USA. Structures were more likely to survive a fire with defensible space immediately adjacent to them. The most effective treatment distance varied between 5 and 20 m (16&ndash;58 ft) from the structure, but distances larger than 30 m (100 ft) did not provide additional protection, even for structures located on steep slopes. The most effective actions were reducing woody cover up to 40% immediately adjacent to structures and ensuring that vegetation does not overhang or touch the structure. Multiple-regression models showed landscape-scale factors, including low housing density and distances to major roads, were more important in explaining structure destruction. The best long-term solution will involve a suite of prevention measures that include defensible space as well as building design approach, community education and proactive land use planning that limits exposure to fire.

International Journal of Wildland Fire↗
USGS2013

Sequential Gaussian co-simulation of rate decline parameters of longwall gob gas ventholes

Gob gas ventholes (GGVs) are used to control methane inflows into a longwall mining operation by capturing the gas within the overlying fractured strata before it enters the work environment. Using geostatistical co-simulation techniques, this paper maps the parameters of their rate decline behaviors across the study area, a longwall mine in the Northern Appalachian basin. Geostatistical gas-in-place (GIP) simulations were performed, using data from 64 exploration boreholes, and GIP data were mapped within the fractured zone of the study area. In addition, methane flowrates monitored from 10 GGVs were analyzed using decline curve analyses (DCA) techniques to determine parameters of decline rates. Surface elevation showed the most influence on methane production from GGVs and thus was used to investigate its relation with DCA parameters using correlation techniques on normal-scored data. Geostatistical analysis was pursued using sequential Gaussian co-simulation with surface elevation as the secondary variable and with DCA parameters as the primary variables. The primary DCA variables were effective percentage decline rate, rate at production start, rate at the beginning of forecast period, and production end duration. Co-simulation results were presented to visualize decline parameters at an area-wide scale. Wells located at lower elevations, i.e., at the bottom of valleys, tend to perform better in terms of their rate declines compared to those at higher elevations. These results were used to calculate drainage radii of GGVs using GIP realizations. The calculated drainage radii are close to ones predicted by pressure transient tests.

Pennsylvania↗