USGS Science⌕ Search

SEARCH · USGS Science

Results for “Near Surface Geophysics”

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 739 records · Page 41Linked to original sources

TES premapping data: Slab ice and snow flurries in the Martian north polar night

In the 1970s, Mariner and Viking spacecraft observations of the north polar region of Mars revealed polar brightness temperatures that were significantly below the expected kinetic temperatures for CO 2 sublimation. For the past few decades, the scientific community has speculated as to the nature of these Martian polar cold spots. Thermal Emission Spectrometer (TES) thermal spectral data have shown these cold spots to result largely from fine‐grained CO 2 and have constrained most of these cold spots to the surface (or near‐surface). Cold spot formation is strongly dependent on topography, forming preferentially near craters and on polar slopes. TES data, combined with Mars Orbiter Laser Altimeter (MOLA) cloud data, suggest atmospheric condensates form a small fraction of the observed cold spots. TES observations of spectra close to a blackbody indicate that another major component of the polar cap is slab CO 2 ice; these spectrally bland regions commonly have a low albedo. The cause is uncertain but may result from most of the light being reflected toward the specular direction, from the slab ice being intrinsically dark, or from it being transparent. Regions of the cap where the difference between the brightness temperatures at 18 μm ( T 18 ) and 25 μm ( T 25 ) is less than 5° are taken to indicate deposits of slab ice. Slab ice is the dominant component of the polar cap at latitudes outside of the polar night.

Journal of Geophysical Research E: Planets↗

Regional quaternary submarine geomorphology in the Florida Keys

High-quality seismic reflection profiles fill a major gap in geophysical data along the south Florida shelf, allowing updated interpretations of the history of the Quaternary coral reef system. Incorporation of the new and existing data sets provides the basis for detailed color maps of the Pleistocene surface and thickness of overlying Holocene accretions. The maps cover the Florida Keys to a margin-wide upper-slope terrace (30 to 40 m deep) and extend from The Elbow Reef (north Key Largo) to Rebecca Shoal (Gulf of Mexico). The data indicate that Pleistocene bedrock is several meters deeper to the southwest than to the northeast, yet in general, Holocene sediments are ~3 to 4 m thick shelf-wide. The Pleistocene map demonstrates the significance of a westward-dipping bedrock surface to Holocene flooding history and coral reef evolution. Seismic facies show evidence for two possible Holocene stillstands. Aerial photographs provide information on the seabed surface, much of which is below seismic resolution. The photographs define a prominent, regional nearshore rock ledge that extends ~2.5 km seaward from the keys' shoreline. They show that bands of rock ridges exist along the outer shelf and on the upper-slope terrace. The photographs also reveal four tracts of outlier reefs on the terrace, one more than had been documented seismically. Seismic and photographic data indicate the tracts are >200 km long, nearly four times longer than previously thought. New interpretations provide insights into a youngest possible terrace age (ca. 175 ka?) and the likelihood that precise ages of oxygen isotope substage 5e ooid tidal-bar and coral reef components may differ. The tidal-bar/reef complex forms the Florida Keys.

Geological Society of America Bulletin↗

Fluctuations in groundwater levels related to regional and local withdrawals in the fractured-bedrock groundwater system in northern Wake County, North Carolina, March 2008-February 2009

A study of dewatering of the fractured-bedrock aquifer in a localized area of east-central North Carolina was conducted from March 2008 through February 2009 to gain an understanding of why some privately owned wells and monitoring wells were intermittently dry. Although the study itself was localized in nature, the resulting water-resources data and information produced from the study will help enable resource managers to make sound water-supply and water-use decisions in similar crystalline-rock aquifer setting in parts of the Piedmont and Blue Ridge Physiographic Provinces. In June 2005, homeowners in a subdivision of approximately 11 homes on lots approximately 1 to 2 acres in size in an unincorporated area of Wake County, North Carolina, reported extremely low water pressure and temporarily dry wells during a brief period. This area of the State, which is in the Piedmont Physiographic Province, is undergoing rapid growth and development. Similar well conditions were reported again in July 2007. In an effort to evaluate aquifer conditions in the area of intermittent water loss, a study was begun in March 2008 to measure and monitor water levels and groundwater use. During the study period from March 2008 through February 2009, regular dewatering of the fractured-bedrock aquifer was documented with water levels in many wells ranging between 100 and 200 feet below land surface. Prior to this period, water levels from the 1980s through the late 1990s were reported to range from 15 to 50 feet below land surface. The study area includes three community wells and more than 30 private wells within a 2,000-foot radius of the dewatered private wells. Although groundwater levels were low, recovery was observed during periods of heavy rainfall, most likely a result of decreased withdrawals owing to less demand for irrigation purposes. Similar areal patterns of low groundwater levels were delineated during nine water-level measurement periods from March 2008 through February 2009. Correlation of groundwater-level distribution patterns with orientations of geologic structures obtained from surficial mapping, borehole geophysical measurements, and interpretation of fracture traces suggests two dominant trends striking north-south and N. 65 degrees W. A variation in overall response to groundwater withdrawals was noted in the continuous groundwater-level records for the monitored observation wells and dewatered private wells. The largest overall declines during the study period were observed in an observation well in which the water-level declined as much as 247 feet from mid-July through early August 2008, during a period of heavy usage. A private well had a water-level decline of about 94 feet during the same monitoring period. The large declines recorded in the observation well and the private well indicated a substantial temporary loss of storage in the fractured-bedrock aquifer near the wells, thus reducing the amount of water available to shallow wells in the area (those wells with total depths of about 300 feet), and resulting in temporary well failures until such time as the aquifer recovered.

North Carolina↗

The Alabama, U.S.A., seismic event and strata collapse of May 7, 1986

On May 7, 1986, the residents of Tuscaloosa, Alabama, felt a seismic event of local magnitude 3.6 that occurred at the same time as a rock burst and roof collapse in an active longwall coal mine. Visual inspection of the seismograms reveals a deficiency in energy at frequencies above 20 Hz compared to tectonic earthquakes or surface blasts. The predominance of energy below 5 Hz may explain reports of body wave magnitudes (mb) greater than 4.2. Also, 1.0 Hz surface waves were more strongly excited than body waves and may explain local felt effects more typically associated with greater epicentral distances. All recorded first motions were dilatational. The concentration of stations in the northern hemisphere allows reverse motion on an east-trending near-vertical plane or strike-slip motion on northwest or southeast trending planes. The reverse focal mechanism is preferred, because the area of roof collapse and the area of active longwall mining are located between two east-striking loose vertical fracture zones. The characteristics of the seismic event suggest that it might have been sudden shear failure resulting from accumulated strain energy in overlying strata behind an active longwall. Although an alternate interpretation of the focal mechanism as an implosion or shear failure in the strata above previously mined out areas is also allowed by the first motion data, this alternate intepretation is not supported by geological data. ?? 1989 Birkha??user Verlag.

Pure and Applied Geophysics PAGEOPH↗

Resurfacing history of the northern plains of Mars based on geologic mapping of Mars Global Surveyor data

Geologic mapping of the northern plains of Mars, based on Mars Orbiter Laser Altimeter topography and Viking and Mars Orbiter Camera images, reveals new insights into geologic processes and events in this region during the Hesperian and Amazonian Periods. We propose four successive stages of lowland resurfacing likely related to the activity of near-surface volatiles commencing at the highland-lowland boundary (HLB) and progressing to lower topographic levels as follows (highest elevations indicated): Stage 1, upper boundary plains, Early Hesperian, <-2.0 to -2.9 km; Stage 2, lower boundary plains and outflow channel dissection, Late Hesperian, <-2.7 to -4.0 km; Stage 3, Vastitas Borealis Formation (VBF) surface, Late Hesperian to Early Amazonian, <-3.1 to -4.1 km; and Stage 4, local chaos zones, Early Amazonian, <-3.8 to -5.0 km. At Acidalia Mensa, Stage 2 and 3 levels may be lower (<-4.4 and -4.8 km, respectively). Contractional ridges form the dominant structure in the plains and developed from near the end of the Early Hesperian to the Early Amazonian. Geomorphic evidence for a northern-plains-filling ocean during Stage 2 is absent because one did not form or its evidence was destroyed by Stage 3 resurfacing. Remnants of possible Amazonian dust mantles occur on top of the VBF. The north polar layered deposits appear to be made up of an up to kilometer-thick lower sequence of sandy layers Early to Middle Amazonian in age overlain by Late Amazonian ice-rich dust layers; both units appear to have outliers, suggesting that they once were more extensive.

Journal of Geophysical Research E: Planets↗

Magma intrusion and volatile ascent beneath Norris Geyser Basin, Yellowstone National Park

Recent activity has provided new insights into the causes of surface deformation in and around the Yellowstone Caldera, a topic that has been debated since the discovery of caldera floor uplift more than four decades ago. An episode of unusually rapid uplift (>15 cm/yr) centered near Norris Geyser Basin along the north caldera rim began in late 2013 and continued until a M w 4.9 earthquake on 30 March 2014; thereafter, uplift abruptly switched to subsidence. Uplift at rates of several centimeters per year resumed in 2016 and continued at least through the end of 2018. Modeling of Global Positioning System and interferometric synthetic aperture radar data suggests an evolving process of deep magma intrusion during 1996–2001 followed by volatile ascent and accumulation at shallow levels, perhaps as shallow as a few hundred meters depth. The depth of shallow volatile accumulation appears to have shallowed from the 2014 to the 2016 deformation episode, and frequent eruptions of Steamboat Geyser since March 2018 are likely a surface manifestation of this ongoing process. Hydrothermal explosion features are prominent in the Norris Geyser Basin area, and the apparent shallow nature of the volatile accumulation implies an increased risk of hydrothermal explosions.

Wyoming↗

Landfill mapping using multi-disciplinary geophysical techniques at the U.S. Air Force Academy, Colorado Springs, CO

This paper describes a multi-disciplinary geophysical survey conducted over a landfill on the U.S. Air Force Academy grounds near Colorado Springs, Colorado. The landfill is known to contain waste generated during the construction of the Academy and reportedly contains buried steel drums. The purpose of the geophysical surveys was to determine the subsurface distribution of buried metallic objects within the landfill. Different geophysical techniques were evaluated along a test line to determine their relative effectiveness at this site. The geophysical methods included total magnetic field, vertical magnetic gradient, VLF, horizontal and vertical coplanar electromagnetic, GPR and seismic refraction. Magnetic and coplanar electromagnetic (EM) methods were chosen to survey the entire landfill because they easily detected magnetic and conductive sources and have better anomaly resolution than other methods evaluated, as demonstrated by the test line results. In addition, these methods are rapid and cost effective for surveys involving a large number of measurements. Surveys of the landfill identified numerous magnetic and conductive anomalies indicating the presence of buried metallic objects. The vertical gradient and EM measurements indicate that several of the large total field anomalies are produced by groups of smaller objects rather than by single, large buried sources. Many of the smaller anomalies are associated with the position of a recently dismantled railroad track and result from iron and steel parts buried along the abandoned grade. Two long, narrow conductive anomalies were identified by the electromagnetic surveys. These conductive features have no surface expression and apparently run the length of the landfill. The EM data indicates these conductors are narrow and relatively shallow. One conductor is relatively magnetic, the other conductor has no magnetic signature suggesting a different composition. The geophysical surveys determined that large areas of the landfill are relatively free of buried metal due to the lack of observed magnetic or conductive anomalies. The geophysical data also suggests the landfill may be larger than originally thought. Numerous magnetic and conductive responses were observed beyond the eastern edge of the present landfill in an area thought to be natural terrain.

Colorado↗

Evaluation of bridge-scour data at selected sites in Ohio

Scour data collected during 1989-94 were evaluated to determine whether pier scour and contraction scour occurred at 22 bridge sites in Ohio. Pier-scour depths computed from selected pier-scour prediction equations were compared with measured pier-scour depths, and the accuracy of the prediction equations were evaluated. Observed pier-scour relations were compared to relations developed through laboratory research. Mean streambed elevations were evaluated to determine the depth of contraction scour. Channel stability was assessed by use of mean streambed elevations at the approach section. Ground-penetrating radar was used at all sites to investigate the presence of historical scour. Pier scour was observed in 45 of 47 scour measurements made during floods; 84 cases of pier scour were documented, 83 at solid-wall piers and 1 at a capped-pile type pier. Estimated recurrence intervals for 27 of the 35 measured streamflows, all on unregulated streams, were less than 2 years. Seventeen pier-scour prediction equations were evaluated. The Froehlich Design equation was found to most closely meet the 'best design equation' criteria for all 84 cases of the observed data. The Larras equation was found to be the best design equation for the observed data where approach-flow attack angles were 10 degrees or less. Observed pier-scour depths and flow depths ranged from 0.5 to 6.1 feet and 3.0 to 19.8 feet, respectively. All pier-scour depths were less than 2.4 times the corresponding pier width. Selected factors were normalized by dividing by effective pier width. LOWESS curves were developed using the 84 cases of observed pier scour. Normalized scour depth increased with normalized flow depth; however, the rate of increase appeared to lessen as normalized flow depth exceeded 2.5. Normalized scour depths increased rapidly as flow intensity approached the threshold value of 1 and then decreased as flow intensities exceeded this threshold. Normalized scour depth was found to increase with Froude number, and a steeper slope was evident for Froude numbers exceeding 0.2. Normalized scour depth was found to increase with median grain size up to about 10 millimeters for bed material near the pier, then decrease for median grain sizes greater than 10 millimeters. Normalized scour depth was also found to decrease as sediment gradation of bed material near the pier increased. The observed pier-scour relations determined from the field measurements tend to support conclusions by previous researchers of streambed scour, except for the previous finding that normalized scour depth decreases consistently with increasing median grain size. Possible factors that may have influenced the observed trends in the relation between normalized scour depth and median grain size in this study are cohesion and scour measurements made at nonequilibrium conditions. LOWESS curves were developed for 45 of 84 cases of observed pier scour where approach-flow attack angles were less than or equal to 10 degrees. These curves were visually compared to LOWESS curves developed from all observations of pier scour. For three relations, differences in the trends of the LOWESS curves were of sufficient magnitude to warrant discussion. Contraction scour was observed in 4 of the 47 scour measurements and ranged from 0.8 to2.3 feet in depth. Analysis of annual mean streambed approach-section elevations indicated that approach sections were generally stable at 18 of the 22 sites. Ground-penetrating radar, a geophysical method that enables subsurface exploration of the streambed when conditions are favorable, was used at all sites to determine whether historical scour had occurred. Results of the ground-penetrating radar surveys at 20 sites in 1990 indicated the presence of historical scour surfaces at 5 sites. At four of the five sites showing evidence of possible historical scour, differences between the estimated depth of historical scour and the maximum observed scour were w

Water-Resources Investigations Report↗

HiRISE observations of new impact craters exposing Martian ground ice

Twenty small new impact craters or clusters have been observed to excavate bright material inferred to be ice at mid and high latitudes on Mars. In the northern hemisphere, the craters are widely distributed geographically and occur at latitudes as low as 39°N. Stability modeling suggests that this ice distribution requires a long-term average atmospheric water vapor content around 25 precipitable microns, more than double the present value, which is consistent with the expected effect of recent orbital variations. Alternatively, near-surface humidity could be higher than expected for current column abundances if water vapor is not well-mixed with atmospheric CO 2 , or the vapor pressure at the ice table could be lower due to salts. Ice in and around the craters remains visibly bright for months to years, indicating that it is clean ice rather than ice-cemented regolith. Although some clean ice may be produced by the impact process, it is likely that the original ground ice was excess ice (exceeding dry soil pore space) in many cases. Observations of the craters suggest small-scale heterogeneities in this excess ice. The origin of such ice is uncertain. Ice lens formation by migration of thin films of liquid is most consistent with local heterogeneity in ice content and common surface boulders, but in some cases nearby thermokarst landforms suggest large amounts of excess ice that may be best explained by a degraded ice sheet.

Journal of Geophysical Research E: Planets↗

Stressing of the New Madrid seismic zone by a lower crust detachment fault

A new mechanical model for the cause of the New Madrid seismic zone in the central United States is analyzed. The model contains a subhorizontal detachment fault which is assumed to be near the domed top surface of locally thickened anomalous lower crust (“rift pillow”). Regional horizontal compression induces slip on the fault, and the slip creates a stress concentration in the upper crust above the rift pillow dome. In the coseismic stage of the model earthquake cycle, where the three largest magnitude 7–8 earthquakes in 1811–1812 are represented by a single model mainshock on a vertical northeast trending fault, the model mainshock has a moment equivalent to a magnitude 8 event. During the interseismic stage, corresponding to the present time, slip on the detachment fault exerts a right-lateral shear stress on the locked vertical fault whose failure produces the model mainshock. The sense of shear is generally consistent with the overall sense of slip of 1811–1812 and later earthquakes. Predicted rates of horizontal strain at the ground surface are about 10 −7 year −1 and are comparable to some observed rates. The model implies that rift pillow geometry is a significant influence on the maximum possible earthquake magnitude.

Journal of Geophysical Research B: Solid Earth↗

Inducing in situ, nonlinear soil response applying an active source

[1] It is well known that soil sites have a profound effect on ground motion during large earthquakes. The complex structure of soil deposits and the highly nonlinear constitutive behavior of soils largely control nonlinear site response at soil sites. Measurements of nonlinear soil response under natural conditions are critical to advancing our understanding of soil behavior during earthquakes. Many factors limit the use of earthquake observations to estimate nonlinear site response such that quantitative characterization of nonlinear behavior relies almost exclusively on laboratory experiments and modeling of wave propagation. Here we introduce a new method for in situ characterization of the nonlinear behavior of a natural soil formation using measurements obtained immediately adjacent to a large vibrator source. To our knowledge, we are the first group to propose and test such an approach. Employing a large, surface vibrator as a source, we measure the nonlinear behavior of the soil by incrementally increasing the source amplitude over a range of frequencies and monitoring changes in the output spectra. We apply a homodyne algorithm for measuring spectral amplitudes, which provides robust signal-to-noise ratios at the frequencies of interest. Spectral ratios are computed between the receivers and the source as well as receiver pairs located in an array adjacent to the source, providing the means to separate source and near-source nonlinearity from pervasive nonlinearity in the soil column. We find clear evidence of nonlinearity in significant decreases in the frequency of peak spectral ratios, corresponding to material softening with amplitude, observed across the array as the source amplitude is increased. The observed peak shifts are consistent with laboratory measurements of soil nonlinearity. Our results provide constraints for future numerical modeling studies of strong ground motion during earthquakes.

Journal of Geophysical Research B: Solid Earth↗

Overview of the Mars Pathfinder Mission: Launch through landing, surface operations, data sets, and science results

Mars Pathfinder successfully landed at Ares Vallis on July 4, 1997, deployed and navigated a small rover about 100 m clockwise around the lander, and collected data from three science instruments and ten technology experiments. The mission operated for three months and returned 2.3 Gbits of data, including over 16,500 lander and 550 rover images, 16 chemical analyses of rocks and soil, and 8.5 million individual temperature, pressure and wind measurements. Path‐finder is the best known location on Mars, having been clearly identified with respect to other features on the surface by correlating five prominent horizon features and two small craters in lander images with those in high‐resolution orbiter images and in inertial space from two‐way ranging and Doppler tracking. Tracking of the lander has fixed the spin pole of Mars, determined the precession rate since Viking 20 years ago, and indicates a polar moment of inertia, which constrains a central metallic core to be between 1300 and ∼2000 km in radius. Dark rocks appear to be high in silica and geochemically similar to anorogenic andesites; lighter rocks are richer in sulfur and lower in silica, consistent with being coated with various amounts of dust. Rover and lander images show rocks with a variety of morphologies, fabrics and textures, suggesting a variety of rock types are present. Rounded pebbles and cobbles on the surface as well as rounded bumps and pits on some rocks indicate these rocks may be conglomerates (although other explanations are also possible), which almost definitely require liquid water to form and a warmer and wetter past. Air‐borne dust is composed of composite silicate particles with a small fraction of a highly magnetic mineral, interpreted to be most likely maghemite; explanations suggest iron was dissolved from crustal materials during an active hydrologic cycle with maghemite freeze dried onto silicate dust grains. Remote sensing data at a scale of a kilometer or greater and an Earth analog correctly predicted a rocky plain safe for landing and roving with a variety of rocks deposited by catstrophic floods, which are relatively dust free. The surface appears to have changed little since it formed billions of years ago, with the exception that eolian activity may have deflated the surface by ∼3–7 cm, sculpted wind tails, collected sand into dunes, and eroded ventifacts (fluted and grooved rocks). Pathfinder found a dusty lower atmosphere, early morning water ice clouds, and morning near‐surface air temperatures that changed abruptly with time and height. Small scale vortices, interpreted to be dust devils, were observed repeatedly in the afternoon by the meteorology instruments and have been imaged.

Journal of Geophysical Research E: Planets↗

Revised groundwater-flow model of the glacial aquifer system north of Aberdeen, South Dakota, through water year 2015

The city of Aberdeen, in northeastern South Dakota, requires an expanded and sustainable supply of water to meet current and future demands. Conceptual and numerical models of the glacial aquifer system in the area north of Aberdeen were developed by the U.S. Geological Survey in cooperation with the City of Aberdeen in 2012. The U.S. Geological Survey, in cooperation with the City of Aberdeen, completed a study to revise the original numerical groundwater-flow model using data through water year (WY) 2015 to aid the City of Aberdeen in their development of plans and strategies for a sustainable water supply and to increase understanding of the glacial aquifer system and groundwater-flow system near Aberdeen. The original model was revised to improve the fit between model-simulated values and observed (measured or estimated) data, provide greater insight into surface-water interactions, and improve the usefulness of the model for water-supply planning. The revised groundwater-flow model (hereafter referred to as the “revised model”) presented in this report supersedes the original model. The purpose of this report is to describe a revised groundwater-flow model including data collection, model calibration, and model results for the glacial aquifer system including the Elm, Middle James, and Deep James aquifers north of Aberdeen, South Dakota, using updated hydrologic data through WY 2015. The original numerical model was revised in several ways. The model was modified by adding four new layers, which included a surficial layer, two intervening confining layers, and a shale bedrock layer. The revised model provides an improved understanding of the groundwater-flow system in comparison to the original model. The principal aquifers of the model area include portions of the Elm, Middle James, and Deep James aquifers. The lithologic information used to define and describe the aquifers in the model area was unaltered; however, aquifer properties and boundary conditions were reviewed and updated using geological information reported by the South Dakota Department of Environmental and Natural Resources and information obtained from geophysical investigations for this study. The horizontal extent of the Elm, Middle James, and Deep James aquifers was unaltered from the original model. The thickness of the Deep James aquifer was modified based on interpretations from the geophysical investigations. In general, groundwater in the Elm aquifer flowed from northwest to southeast and locally towards rivers and streams. Similarly, in the Middle James and Deep James aquifers, groundwater also typically flowed southeast. The revisions made to the original model include use of the following MODFLOW stress packages: Recharge, Evapotranspiration, Time-Variant Specified Head, Wells, Drains, and Stream Flow Routing, all of which were updated from the original model except for the Stream Flow Routing Package, which replaced the River Package used in the original model. Model calibration is the process of estimating model parameters to minimize the differences, or residuals, between observed data and simulated values; therefore, Parameter ESTimation (PEST) software was used to optimize model input parameters by matching model-simulated values to observed data. Calibration parameters included horizontal hydraulic conductivity, vertical hydraulic conductivity, specific yield, specific storage, and vertical streambed conductance for stream and drain cells. Multipliers were used to calibrate the recharge and evapotranspiration stresses. Evapotranspiration extinction depth also was adjusted during model calibration. Comparisons to the original model are described to highlight the changes made in the revised model. In general, the revised model adequately simulates the natural system and compares favorably with observed hydrologic data. Simulated water levels were evaluated by comparing them to single water-level observations at selected well locations. The selected wells were the same wells used in the original model. The coefficient of determination value between simulated and observed water levels for the revised model was 0.89 and included simulated and observed values from October 1, 1974 (WY 1975), through September 30, 2015 (WY 2015). The coefficient of determination value for the original model was 0.94 and included simulated and observed values from October 1, 1974, through September 30, 2009. The difference may indicate that the original model could have been overfit to hydraulic head observations because base flow was not simulated. The additional data used in the revised model included some climatically wetter, more extreme periods, such as 2011, in which annual precipitation was 30.9 inches. Average annual precipitation for the original model timeframe, which included data from WYs 1975–2009, was 20.26 inches. Additional precipitation data for WYs 2010–15, included in the revised model timeframe, resulted in an average annual precipitation for WYs 1975–2015 in the model area of 20.6 inches. The larger variability in climate data coupled with the additional water-level data could explain the lower coefficient of determination for water levels in the revised model. The revised model was used to calculate various groundwater-budget components for steady-state and transient conditions for WYs 1975–2015. The time-variant specified-head cells in the revised model had the largest change when compared to the original steady-state model for inflows and outflows. Comparing the transient budget components between the original and the revised models indicated that inflow from recharge and time-variant specified-head cells had the greatest effect on groundwater inflows, and outflow from storage had the greatest effect on groundwater outflows. The simulated potentiometric contours from the revised model were compared with (1) the observed (interpreted) potentiometric surface (layer 2) and the hydraulic head values (layers 4 and 6) and (2) the simulated contours from the original model. The simulated hydraulic gradients and general direction of groundwater flow in the Elm aquifer in the revised model generally matched the observed potentiometric contours, the simulated potentiometric contours from the original model, and general flow directions interpreted to be perpendicular to the contours. Minor discrepancies between simulated potentiometric contours from the revised model and the observed potentiometric contours may be due to the lack of observed data in the model area. The revised model was designed to reduce the limitations of the original model. The revisions were validated by comparing the results of the original model with the revised model. A primary benefit of the revised model is the inclusion of the surficial deposits and the confining units as explicit layers in the model. The addition of the surficial layer was beneficial for three primary reasons: (1) more accurate representation of recharge from precipitation, (2) more accurate representation of groundwater evapotranspiration, and (3) more accurate representation of groundwater and surface-water interactions. The groundwater model is a numeric approximation of a complex physical hydrologic system, and the revised model data were interpolated in regions with sparse data. Additionally, model discretization included averaged and interpolated values for water use, withdrawal rates, and hydraulic conductivity. The revised model provides a useful estimate for hydraulic gradients, groundwater-flow directions, and aquifer response to groundwater withdrawals.

South Dakota↗

Complex mean circulation over the inner shelf south of Martha's Vineyard revealed by observations and a high-resolution model

Inner-shelf circulation is governed by the interaction between tides, baroclinic forcing, winds, waves, and frictional losses; the mean circulation ultimately governs exchange between the coast and ocean. In some cases, oscillatory tidal currents interact with bathymetric features to generate a tidally rectified flow. Recent observational and modeling efforts in an overlapping domain centered on the Martha's Vineyard Coastal Observatory (MVCO) provided an opportunity to investigate the spatial and temporal complexity of circulation on the inner shelf. ADCP and surface radar observations revealed a mean circulation pattern that was highly variable in the alongshore and cross-shore directions. Nested modeling incrementally improved representation of the mean circulation as grid resolution increased and indicated tidal rectification as the generation mechanism of a counter-clockwise gyre near the MVCO. The loss of model skill with decreasing resolution is attributed to insufficient representation of the bathymetric gradients (Δh/h), which is important for representing nonlinear interactions between currents and bathymetry. The modeled momentum balance was characterized by large spatial variability of the pressure gradient and horizontal advection terms over short distances, suggesting that observed inner-shelf momentum balances may be confounded. Given the available observational and modeling data, this work defines the spatially variable mean circulation and its formation mechanism&mdash;tidal rectification&mdash;and illustrates the importance of model resolution for resolving circulation and constituent exchange near the coast. The results of this study have implications for future observational and modeling studies near the MVCO and other inner-shelf locations with alongshore bathymetric variability.

Massachusetts↗

Near real-time monitoring of volcanic surface deformation from GPS measurements at Long Valley Caldera, California

Long Valley Caldera in eastern California is an active volcanic area and has shown continued unrest in the last three decades. We have monitored surface deformation from Global Positioning System (GPS) data by using a projection method that we call Targeted Projection Operator (TPO). TPO projects residual time series with secular rates and periodic terms removed onto a predefined spatial pattern. We used the 2009–2010 slow deflation as a target spatial pattern. The resulting TPO time series shows a detailed deformation history including the 2007–2009 inflation, the 2009–2010 deflation, and a recent inflation that started in late-2011 and is continuing at the present time (November 2012). The recent inflation event is about four times faster than the previous 2007–2009 event. A Mogi source of the recent event is located beneath the resurgent dome at about 6.6 km depth at a rate of 0.009 km3/yr volume change. TPO is simple and fast and can provide a near real-time continuous monitoring tool without directly looking at all the data from many GPS sites in this potentially eruptive volcanic system.

California↗

ARRA-funded V S30 measurements using multi-technique approach at strong-motion stations in California and central-eastern United States

Funded by the 2009 American Recovery and Reinvestment Act (ARRA), we conducted geophysical site characterizations at 191 strong-motion stations: 187 in California and 4 in the Central-Eastern United States (CEUS). The geophysical methods used at each site included passive and active surface-wave and body-wave techniques. Multiple techniques were used at most sites, with the goal of robustly determining VS (shear-wave velocity) profiles and V S30 (the time-averaged shear-wave velocity in the upper 30 meters depth). These techniques included: horizontal-to-vertical spectral ratio (HVSR), two-dimensional (2-D) array microtremor (AM), refraction microtremor (ReMi™), spectral analysis of surface wave (SASW), multi-channel analysis of surface waves (Rayleigh wave: MAS R W; and Love wave: MAS L W), and compressional- and shear-wave refraction. Of the selected sites, 47 percent have crystalline, volcanic, or sedimentary rock at the surface or at relatively shallow depth, and 53 percent are of Quaternary sediments located in either rural or urban environments. Calculated values of VS30 span almost the full range of the National Earthquake Hazards Reduction Program (NEHRP) Site Classes, from D (stiff soils) to B (rock). The NEHRP Site Classes based on V S30 range from being consistent with the Class expected from analysis of surficial geology, to being one or two Site Classes below expected. In a few cases where differences between the observed and expected Site Class occurred, it was the consequence of inaccurate or coarse geologic mapping, as well as considerable degradation of the near-surface rock. Additionally, several sites mapped as rock have Site Class D (stiff soil) velocities, which is due to the extensive weathering of the surficial rock.

California↗

The earthquake prediction experiment at Parkfield, California

Since 1985, a focused earthquake prediction experiment has been in progress along the San Andreas fault near the town of Parkfield in central California. Parkfield has experienced six moderate earthquakes since 1857 at average intervals of 22 years, the most recent a magnitude 6 event in 1966. The probability of another moderate earthquake soon appears high, but studies assigning it a 95% chance of occurring before 1993 now appear to have been over-simplified. The identification of a Parkfield fault “segment” was initially based on geometric features in the surface trace of the San Andreas fault, but more recent microearthquake studies have demonstrated that those features do not extend to seismogenic depths. On the other hand, geodetic measurements are consistent with the existence of a “locked” patch on the fault beneath Parkfield that has presently accumulated a slip deficit equal to the slip in the 1966 earthquake. A magnitude 4.7 earthquake in October 1992 brought the Parkfield experiment to its highest level of alert, with a 72-hour public warning that there was a 37% chance of a magnitude 6 event. However, this warning proved to be a false alarm. Most data collected at Parkfield indicate that strain is accumulating at a constant rate on this part of the San Andreas fault, but some interesting departures from this behavior have been recorded. Here we outline the scientific arguments bearing on when the next Parkfield earthquake is likely to occur and summarize geophysical observations to date.

California↗

Application of ground-penetrating-radar methods in hydrogeologic studies

A ground-penetrating-radar system was used to study selected stratified-drift deposits in Connecticut. Ground-penetrating radar is a surface-geophysical method that depends on the emission, transmission, reflection, and reception of an electromagnetic pulse and can produce continuous high-resolution profiles of the subsurface rapidly and efficiently. Traverse locations on land included a well field in the town of Mansfield, a sand and gravel pit and a farm overlying a potential aquifer in the town of Coventry, and Haddam Meadows State Park in the town of Haddam. Traverse locations on water included the Willimantic River in Coventry and Mansfield Hollow Lake in Mansfield. The penetration depth of the radar signal ranged from about 20 feet in fine-grained glaciolacustrine sediments to about 70 feet in coarse sand and gravel. Some land records in coarse-grained sediments show a distinct, continuous reflection from the water table about 5 to 11 feet below land surface. Parallel reflectors on the records are interpreted as fine-grained sediments. Hummocky or chaotic reflectors are interpreted as cross-bedded or coarse-grained sediments. Other features observed on some of the radar records include the till and bedrock surface. Records collected on water had distinct water-bottom multiples (more than one reflection) and diffraction patterns from boulders. The interpretation of the radar records, which required little or no processing, was verified by using lithologic logs from test holes located along some of the land traverses and near the water traverses.

Connecticut↗