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.

SourceResearch connectionFrom yearThrough year

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

At least 685 records · Page 38Linked to original sources

USGS2026

Re-examination of the 1958 Huslia earthquake sequence and regional tectonics of the northwestern Koyukuk Basin, Alaska

In 1958 a M w 6.8 (M S = 7.3) mainshock and two M w ∼ 5.9 aftershocks occurred near the village of Huslia in the remote Koyukuk Basin of northern interior Alaska, with no surface faulting recognized afterward. This region lies in a transitional zone between compressional tectonics associated with the collision of the Pacific and Yakutat plates with North America and the extrusion of western Alaska. Since 1969 several earthquake swarms with M w 4–5 events have occurred ∼50 km northwest of Huslia in the western Purcell Mountains, including a swarm in 2019–2022 with over 9500 events. Previous waveform modeling studies conducted 30 years ago suggest the 1958 mainshock occurred on a northeast-trending normal fault. This sense of rupture is inconsistent with focal mechanisms and trends of post-1973 seismicity in the 1958 epicentral region that indicate strike-slip motion along northwest or southeast trending faults. A recently mapped east-west striking fault scarp along the southern front of the Purcell Mountains is also inconsistent with northeast-striking normal faulting. To better characterize the 1958 sequence and its relation to the regional tectonics, we measured scarp profiles across the Purcell Mountains scarp, relocated pre-1974 seismicity, analyzed first motions of pre-1974 events, and modeled potential fields (gravity, magnetics) geophysical data. Our results suggest complex regional deformation occurs along east-west, north-south and northwest-southeast trending strike-slip faults and northwest-southeast trending normal faults. The Purcell Mountains swarm may be related to hydrothermal activity as several hot springs surround the region.

Alaska↗

USGS2001

Delineation of tidal scour through marine geophysical techniques at Sloop Channel and Goose Creek bridges, Jones Beach State Park, Long Island, New York

Inspection of the Goose Creek Bridge in southeastern Nassau County in April 1998 by the New York State Department of Transportation (NYSDOT) indicated a separation of bridge piers from the road bed as a result of pier instability due to apparent seabed scouring by tidal currents. This prompted a cooperative study by the U.S. Geological Survey with the NYSDOT to delineate the extent of tidal scour at this bridge and at the Sloop Channel Bridge, about 0.5 mile to the south, through several marine- geophysical techniques. These techniques included use of a narrow-beam, 200-kilohertz, research-grade fathometer, a global positioning system accurate to within 3 feet, a 3.5 to 7-kilohertz seismic-reflection profiler, and an acoustic Doppler current profiler (ADCP). The ADCP was used only at the Sloop Channel Bridge; the other techniques were used at both bridges. Results indicate extensive tidal scour at both bridges. The fathometer data indicate two major scour holes nearly parallel to the Sloop Channel Bridge—one along the east side, and one along the west side (bridge is oriented north-south). The scour-hole depths are as much as 47 feet below sea level and average more than 40 feet below sea level; these scour holes also appear to have begun to connect beneath the bridge. The deepest scour is at the north end of the bridge beneath the westernmost piers. The east-west symmetry of scour at Sloop Channel Bridge suggests that flood and ebb tides produce extensive scour. The thickness of sediment that has settled within scour holes could not be interpreted from fathometer data alone because fathometer frequencies cannot penetrate beneath the sea-floor surface. The lower frequencies used in seismic-reflection profiling can penetrate the sea floor and underlying sediments, and indicate the amount of infilling of scour holes, the extent of riprap under the bridge, and the assemblages of clay, sand, and silt beneath the sea floor. The seismic- reflection surveys detected 2 to 5 feet of sediment filling the scour holes at both bridges; this indicates that the fathometer surveys were undermeasuring the effective depth of bridge scour by 2 to 5 feet through their inability to penetrate the infilled sediment. Several clay layers with thicknesses of 3 to 5 feet were detected beneath the sea floor at both bridges. Most of the piers beneath Sloop Channel Bridge appear to be surrounded by riprap, but, in several areas the riprap appears to be slumping or sliding into adjacent scour holes. Similar slumping was indicated at the Goose Creek Bridge. Most of the sediment underlying the sea floor at both bridges is interpreted as a fine-grained, cross-bedded sand. ADCP data from Sloop Channel indicate that the constricted flow beneath the bridge increases the horizontal current velocities from 2 to 6 feet per second. Total measured discharge beneath Sloop Channel Bridge was 41,800 cubic feet per second at flood tide and 27,600 cubic feet per second at ebb tide.

Water-Resources Investigations Report↗

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↗

USGS1995

Water resources of the Grand Portage Indian Reservation, northeastern Minnesota

The Grand Portage Indian Reservation Tribal Council needs information about the availability and quality of the ground water in the Reservation to develop, protect, and manage this resource for future use. The U.S. Geological Survey, in cooperation with the Grand Portage Indian Reservation Tribal Council, did a three-year study of the ground water in the Reservation to provide this needed information. This report presents the results of that study. Presently, ground water from bedrock is the principal source of supply for municipal, commercial, and residential water use. The bedrock aquifers are the (1) North Shore Volcanic Group basalt, (2) Keweenawan Volcanic and intrusive rocks, which are gabbro and diabase, and (3) Rove Formation argillite, slate, and graywacke. Sand and gravel aquifers are a small source of ground water. The storage coefficient of the bedrock aquifers was estimated to be 1 x 10 -4 , which is a small value typical of confined, fractured rock aquifers. The median estimate of transmissivity determined from data for 17 wells completed in bedrock was 20 feet squared per day; the range was from 3 to 500 feet squared per day. Reported yield of 19 wells completed in bedrock had a range of 1 to 100 gallons per minute and a median of 7 gallons per minute. The median yield of 11 wells completed in the North Shore Volcanic Group was 16 gallons per minute; the median yield of 8 wells completed in the Keweenawan Volcanic and intrusive rocks and Rove Formation was 4 gallons per minute. Geophysical logs and televiewer images of two wells completed in bedrock indicated the boreholes penetrated many fractures. Hydrofracturing of the two wells increased their yield from about 0.05 and 0.25 gallons per minute to about 1.5 and 1.2 gallons per minute, respectively. Although the estimated yield from the two wells was increased by 30 and by nearly 5 times after hydrofracturing, the well yield after hydrofracturing was still small. Water types determined from analyses of water from nine wells completed in bedrock were sodium-chloride, calcium-chloride, sodium-bicarbonate, and calcium-bicarbonate. Water from three wells had concentrations of dissolved solids (800 to 3,110 milligrams per liter) and dissolved chloride (410 to 1,600 milligrams per liter) that were higher than their respective Secondary Maximum Contaminant Levels of 500 and 250 milligrams per liter established by the U.S. Environmental Protection Agency. Water from two wells had concentrations of dissolved iron (1,600 and 1,300 micrograms per liter) that were higher than the Secondary Maximum Contaminant Level of 300 micrograms per liter. Water from an observation well located about 200 feet downgradient from an abandoned landfill and screened from 79 to 84 feet below land surface in a gravel aquifer had a trace amount of toluene (0.2 micrograms per liter). The presence of toluene suggested possible contamination.

Minnesota↗

USGS1999

Areal studies aid protection of ground-water quality in Illinois, Indiana, and Wisconsin

In 1991, the U.S. Geological Survey, in cooperation with the U.S. Environmental Protection Agency, initiated studies designed to characterize the ground-water quality and hydrogeology in northern Illinois, and southern and eastern Wisconsin (with a focus on the north-central Illinois cities of Belvidere and Rockford, and the Calumet region of northeastern Illinois and northwestern Indiana). These areas are considered especially susceptible to ground-water contamination because of the high density of industrial and waste-disposal sites and the shallow depth to the unconsolidated sand and gravel aquifers and the fractured, carbonate bedrock aquifers that underlie the areas. The data and conceptual models of ground-water flow and contaminant distribution and movement developed as part of the studies have allowed Federal, State, and local agencies to better manage, protect, and restore the water supplies of the areas. Water-quality, hydrologic, geologic, and geophysical data collected as part of these areal studies indicate that industrial contaminants are present locally in the aquifers underlying the areas. Most of the contaminants, particularly those at concentrations that exceeded regulatory water-quality levels, were detected in the sand and gravel aquifers near industrial or waste-disposal sites. In water from water-supply wells, the contaminants that were present generally were at concentrations below regulatory levels. The organic compounds detected most frequently at concentrations near or above regulatory levels varied by area. Trichloroethene, tetrachloroethene, and 1,1,1-trichloroethane (volatile chlorinated compounds) were most prevalent in north-central Illinois; benzene (a petroleum-related compound) was most prevalent in the Calumet region. Differences in the type of organic compounds that were detected in each area likely reflect differences in the types of industrial sites that predominate in the areas. Nickel and aluminum were the trace metals detected most frequently at concentrations above regulatory levels in both areas. Contaminants in the shallow sand and gravel aquifers and carbonate aquifers appear to have moved with ground water discharging to local lakes, streams, and wetlands. Ground-water flow and possibly contaminant movement is concentrated in the weathered surface zones and in deeper fractures of the carbonate aquifers underlying both areas.

Illinois, Indiana, Wisconsin↗

USGS2001

Hydrology of Yucca Mountain, Nevada

Yucca Mountain, located in southern Nevada in the Mojave Desert, is being considered as a geologic repository for high-level radioactive waste. Although the site is arid, previous studies indicate net infiltration rates of 5-10 mm yr-1 under current climate conditions. Unsaturated flow of water through the mountain generally is vertical and rapid through the fractures of the welded tuffs and slow through the matrix of the nonwelded tuffs. The vitric-zeolitic boundary of the nonwelded tuffs below the potential repository, where it exists, causes perching and substantial lateral flow that eventually flows through faults near the eastern edge of the potential repository and recharges the underlying groundwater system. Fast pathways are located where water flows relatively quickly through the unsaturated zone to the water table. For the bulk of the water a large part of the travel time from land surface to the potential repository horizon (~300 m below land surface) is through the interlayered, low fracture density, nonwelded tuff where flow is predominately through the matrix. The unsaturated zone at Yucca Mountain is being modeled using a three-dimensional, dual-continuum numerical model to predict the results of measurements and observations in new boreholes and excavations. The interaction between experimentalists and modelers is providing confidence in the conceptual model and the numerical model and is providing researchers with the ability to plan further testing and to evaluate the usefulness or necessity of further data collection.

Reviews of Geophysics↗

USGS2000

Geophysical constraints on the Virgin River Depression, Nevada, Utah, and Arizona

Gravity and aeromagnetic data provide insights into the subsurface lithology and structure of the Virgin River Depression (VRD) of Nevada, Utah, and Arizona. The gravity data indicate that the Quaternary and Tertiary sedimentary deposits hide a complex pre-Cenozoic surface. A north-northwest-trending basement ridge separates the Mesquite and Mormon basins, as evidenced by seismic-reflection, gravity, and aeromagnetic data. The Mesquite basin is very deep, reaching depths of 8?10 km. The Mormon basin reaches thicknesses of 5 km. Its northern margin is very steep and may be characterized by right steps, although this interpretation could change with additional gravity stations. Most of the young (Quaternary), small-displacement faults trend within 10? of due north and occur within the deeper parts of the Mesquite basin north of the Virgin River. South of the Virgin River, only a few, young, small-displacement faults are mapped; the trend of these faults is more northeasterly and parallels the basement topography and is distinct from that of the faults to the north. The Virgin River appears to follow the margin of the basin as it emerges from the plateau. The high-resolution aeromagnetic data outline the extent of shallow volcanic rocks in the Mesquite basin. The north-northwest alignment of volcanic rocks east of Toquop Wash appear to be structurally controlled because of faults imaged on seismic-reflection profiles and because the alignment is nearly perpendicular to the direction of Cenozoic extension. More buried volcanics likely exist to the north and east of the high-resolution aeromagnetic survey. Broader aeromagnetic anomalies beneath pre-Cenozoic basement in the Mormon Mountains and Tule Springs Hills reflect either Precambrian basement or Tertiary intrusions. These rocks are probably barriers to groundwater flow, except where fractured.

Arizona, Nevada, Utah↗

USGS2002

Seismic velocity structure of the continental lithosphere from controlled source data

The purpose of this chapter is to provide a summary of the seismic velocity structure of the continental lithosphere, i.e., the crust and uppermost mantle. We define the crust as the outer layer of the Earth that is separated from the underlying mantle by the Mohorovičić discontinuity (Moho). We adopted the usual convention of defining the seismic Moho as the level in the Earth where the seismic compressional-wave ( P -wave) velocity increases rapidly or gradually to a value greater than or equal to 7.6 km sec −1 (Steinhart, 1967), defined in the data by the so-called “ Pn ” phase ( P-n ormal). Here we use the term uppermost mantle to refer to the 50–200+ km thick lithospheric mantle that forms the root of the continents and that is attached to the crust (i.e., moves with the continental plates). This summary has been preceded by 90 y of intense scientific activity. Mohorovičić (1910) was the first to publish an estimate for crustal thickness (54 km near Zagreb, Croatia), and to describe the seismically defined boundary between the crust and mantle that now bears his name (often shortened to “the Moho”). The fact that the oceanic crust is significantly thinner than continental crust (5 km versus about 40 km) was documented 40 y later (e.g., Hersey et al., 1952). Numerous later studies demonstrated that the continental crust varies in thickness from about 15 km to greater than 70 km beneath the Tibet plateau. Jarchow and Thompson (1989) provide a useful summary of early crustal studies, and Table 1 provides additional references. We emphasize results from active-source seismic refraction profiles that provide detailed P -wave velocity information. The shear-wave ( S -wave) structure of the crust and uppermost mantle, as determined by surface waves, is discussed, for example, by Ekström et al. (1997) and Ritzwoller et al. (1998).

International Geophysics↗

USGS2003

A comprehensive approach to monitoring volcano deformation as a window on the eruption cycle

Since the 1980 eruption of Mount St. Helens, volcanologists have made considerable progress toward predicting eruptions on the basis of precursors that typically start a few days to several months in advance. Although accurate eruption prediction is by no means routine, it may now be possible in some cases to extend the effective warning period by anticipating the onset of short-term precursors. Three promising indicators of deep magmatic processes are (1) deep, long-period earthquakes and tremor that indicate the ascent of magma through the crust, (2) magmatic CO 2 emission rate as a proxy for magma supply rate, and (3) relatively broad, generally aseismic surface uplift caused by magmatic intrusions. In the latter case it is essential to sample the deformation field thoroughly in both time and space to adequately constrain source models. Until recently, this has been nearly impossible because high-precision sensors could not be deployed in sufficient numbers, nor could extensive geodetic surveys be conducted often enough. Advances in instrumentation, interferometric synthetic aperture radar (InSAR), and telecommunications are helping to overcome these limitations. As a result, comprehensive geodetic monitoring of selected volcanoes is now feasible. A combination of InSAR, large-aperture GPS surveys, microgravity surveys, and dense arrays of continuous GPS stations, strain meters, and tiltmeters can reveal both spatial and temporal patterns of ground deformation throughout the eruption cycle. Improved geodetic monitoring of many of the world's volcanoes would be a major stride toward better understanding of magmatic processes and longer-term eruption forecasts.

Washington↗

USGS2009

THEMIS high-resolution digital terrain: Topographic and thermophysical mapping of Gusev Crater, Mars

We discuss a new technique to generate high-resolution digital terrain models (DTMs) and to quantitatively derive and map slope-corrected thermophysical properties such as albedo, thermal inertia, and surface temperatures. This investigation is a continuation of work started by Kirk et al. (2005), who empirically deconvolved Thermal Emission Imaging System (THEMIS) visible and thermal infrared data of this area, isolating topographic information that produced an accurate DTM. Surface temperatures change as a function of many variables such as slope, albedo, thermal inertia, time, season, and atmospheric opacity. We constrain each of these variables to construct a DTM and maps of slope-corrected albedo, slope- and albedo-corrected thermal inertia, and surface temperatures across the scene for any time of day or year and at any atmospheric opacity. DTMs greatly facilitate analyses of the Martian surface, and the MOLA global data set is not finely scaled enough (128 pixels per degree, ???0.5 km per pixel near the equator) to be combined with newer data sets (e.g., High Resolution Imaging Science Experiment, Context Camera, and Compact Reconnaissance Imaging Spectrometer for Mars at ???0.25, ???6, and ???20 m per pixel, respectively), so new techniques to derive high-resolution DTMs are always being explored. This paper discusses our technique of combining a set of THEMIS visible and thermal infrared observations such that albedo and thermal inertia variations within the scene are eliminated and only topographic variations remain. This enables us to produce a high-resolution DTM via photoclinometry techniques that are largely free of albedo-induced errors. With this DTM, THEMIS observations, and a subsurface thermal diffusion model, we generate slope-corrected maps of albedo, thermal inertia, and surface temperatures. In addition to greater accuracy, these products allow thermophysical properties to be directly compared with topography.

Journal of Geophysical Research E: Planets↗

USGS2012

Interpretation of gravity data in a complex volcano-tectonic setting, southwestern Nevada

This regional gravity study, based on an irregular 2-km data grid, was conducted during the past few years at Yucca Mountain, southern Nye County, Nevada, as part of a program to locate a suitable repository for high-level nuclear waste. About 100 surface rock samples, three borehole gamma-gamma logs, and one borehole gravity study provide excellent density control. A nearly linear increase in density of 0.26 g/cm 3 per kilometer of depth is indicated in the thick tuff sequences that underlie the mountain. Isostatic and 2.0-g/cm 3 Bouguer corrections were applied to the observed gravity values to remove regional gradients and topographic effects, respectively. The Bare Mountain gravity high, with an isostatic anomaly maximum of 48 mGal, is connected with a greater gravity high over the Funeral Mountains, to the southwest; together, these highs result from a continuous block of dense, metamorphosed Precambrian and Paleozoic rocks that stretches across much of the Walker Lane from the east edge of Death Valley to Bare Mountain. The Calico Hills gravity high appears more likely to originate from a northeast trending buried ridge of Paleozoic rocks that extends southwestward beneath Busted Butte, 5 km southeast of the proposed repository, where two- and three-dimensional modeling indicates that the pre-Cenozoic rocks lie less than 1000 m beneath the surface. Tuff, at least 4000 m thick, fills a large steep-sided depression in the pretuff rocks beneath Yucca Mountain and Crater Flat. The gravity low and the thick tuff section lie within a large collapse area that includes the Crater Flat-Timber Mountain-Silent Canyon caldera complexes. Gravity lows in Crater Flat itself are interpreted to coincide with the source areas of the Prow Pass Member, the Bullfrog Member, and the Tram Member of the Crater Flat Tuff; these source areas add nearly 350 km 2 to the previously recognized extent of the local caldera complexes. Southward extension of the broad gravity low associated with Crater Flat into the Amargosa Desert is evidence for sector graben-type collapse segments related to the formation of the Timber Mountain caldera and superimposed on the other volcanic and extensional structures within Crater Flat.

Journal of Geophysical Research Solid Earth↗

USGS1995

Results of geophysical surveys of glacial deposits near a former waste-disposal site, Nashua, New Hampshire

Geophysical investigations were done near a former waste-disposal site in Nashua, New Hampshire to determine the thickness and infer hydraulic characteristics of the glacial sediments that underlie the area. Approximately 5 miles of ground- penetrating radar (GPR) data were collected in the study area by use of dual-80 Megahertz antennas. Three distinct radar-reflection signatures were evident from the data and are interpreted to represent (1) glacial lake-bottom sediments, (2) coarse sand and gravel and (or) sandy glacial till, and (3) bedrock. The GPR signal penetrated as much as 70 feet of sediment in coarse-grained areas, but penetration depth was generally less than 40 feet in extensive areas of fine-grained deposits. Geologic features were evident in many of the profiles. Glacial-lake-bottom sediments were the most common features identified. Other features include deltas deposited in glacial Lake Nashua and lobate fans of sediment deposited subaqueously at the distal end of deltaic sediments. Cross-bedded sands were often identifiable in the deltaic sediments. Seismic-refraction data were also collected at five of the GPR data sites. In most cases, depths to the water table and to the till and (or) bedrock surface indicated by the seismic-refraction data compared favorably with depths calculated from the GPR data. Test holes were drilled at three locations to determine the true depths to radar reflectors and to determine the types of geologic material represented by the various reflectors.

Open-File Report↗

USGS2010

Complexities in barrier island response to sea level rise: Insights from numerical model experiments, North Carolina Outer Banks

Using a morphological-behavior model to conduct sensitivity experiments, we investigate the sea level rise response of a complex coastal environment to changes in a variety of factors. Experiments reveal that substrate composition, followed in rank order by substrate slope, sea level rise rate, and sediment supply rate, are the most important factors in determining barrier island response to sea level rise. We find that geomorphic threshold crossing, defined as a change in state (e.g., from landward migrating to drowning) that is irreversible over decadal to millennial time scales, is most likely to occur in muddy coastal systems where the combination of substrate composition, depth-dependent limitations on shoreface response rates, and substrate erodibility may prevent sand from being liberated rapidly enough, or in sufficient quantity, to maintain a subaerial barrier. Analyses indicate that factors affecting sediment availability such as low substrate sand proportions and high sediment loss rates cause a barrier to migrate landward along a trajectory having a lower slope than average barrier island slope, thereby defining an “effective” barrier island slope. Other factors being equal, such barriers will tend to be smaller and associated with a more deeply incised shoreface, thereby requiring less migration per sea level rise increment to liberate sufficient sand to maintain subaerial exposure than larger, less incised barriers. As a result, the evolution of larger/less incised barriers is more likely to be limited by shoreface erosion rates or substrate erodibility making them more prone to disintegration related to increasing sea level rise rates than smaller/more incised barriers. Thus, the small/deeply incised North Carolina barriers are likely to persist in the near term (although their long-term fate is less certain because of the low substrate slopes that will soon be encountered). In aggregate, results point to the importance of system history (e.g., previous slopes, sediment budgets, etc.) in determining migration trajectories and therefore how a barrier island will respond to sea level rise. Although simple analytical calculations may predict barrier response in simplified coastal environments (e.g., constant slope, constant sea level rise rate, etc.), our model experiments demonstrate that morphological-behavior modeling is necessary to provide critical insights regarding changes that may occur in environments having complex geometries, especially when multiple parameters change simultaneously.

North Carolina↗

USGS2006

Frictional strength heterogeneity and surface heat flow: Implications for the strength of the creeping San Andreas fault

Heat flow measurements along much of the San Andreas fault (SAF) constrain the apparent coefficient of friction (??app) of the fault to <0.2, much lower than laboratory-derived friction values for most geologic materials. However, heat flow data are sparse near the creeping section of the SAF, a frictional "asperity" where the fault slips almost exclusively by aseismic creep. We test the hypothesis that the creeping section has a substantially higher or lower ?? app than adjacent sections of the SAF. We use numerical models to explore the effects of faults with spatially and temporally heterogeneous frictional strength on the spatial distribution of surface heat flow. Heat flow from finite length asperities is uniformly lower than predicted by assuming an infinitely long fault. Over geologic time, lateral offset from strike-slip faulting produces heat flow patterns that are asymmetric across the fault and along strike. We explore a range of asperity sizes, slip rates, and displacement histories for comparing predicted spatial patterns of heat flow with existing measurements. Models with ??app ??? 0.1 fit the data best. For most scenarios, heat flow anomalies from a frictional asperity with ??app > 0.2 should be detectable even with the sparse existing observations, implying that ??app for the creeping section is as low as the surrounding SAF. Because the creeping section does not slip in large earthquakes, the mechanism controlling its weakness is not related to dynamic processes resulting from high slip rate earthquake ruptures. Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

USGS2009

The 2006-2007 Kuril Islands great earthquake sequence

[1] The southwestern half of a ∼500 km long seismic gap in the central Kuril Island arc subduction zone experienced two great earthquakes with extensive preshock and aftershock sequences in late 2006 to early 2007. The nature of seismic coupling in the gap had been uncertain due to the limited historical record of prior large events and the presence of distinctive upper plate, trench and outer rise structures relative to adjacent regions along the arc that have experienced repeated great interplate earthquakes in the last few centuries. The intraplate region seaward of the seismic gap had several shallow compressional events during the preceding decades (notably an M S 7.2 event on 16 March 1963), leading to speculation that the interplate fault was seismically coupled. This issue was partly resolved by failure of the shallow portion of the interplate megathrust in an M W = 8.3 thrust event on 15 November 2006. This event ruptured ∼250 km along the seismic gap, just northeast of the great 1963 Kuril Island (M w = 8.5) earthquake rupture zone. Within minutes of the thrust event, intense earthquake activity commenced beneath the outer wall of the trench seaward of the interplate rupture, with the larger events having normal-faulting mechanisms. An unusual double band of interplate and intraplate aftershocks developed. On 13 January 2007, an M W = 8.1 extensional earthquake ruptured within the Pacific plate beneath the seaward edge of the Kuril trench. This event is the third largest normal-faulting earthquake seaward of a subduction zone on record, and its rupture zone extended to at least 33 km depth and paralleled most of the length of the 2006 rupture. The 13 January 2007 event produced stronger shaking in Japan than the larger thrust event, as a consequence of higher short-period energy radiation from the source. The great event aftershock sequences were dominated by the expected faulting geometries; thrust faulting for the 2006 rupture zone, and normal faulting for the 2007 rupture zone. A large intraplate compressional event occurred on 15 January 2009 (M w = 7.4) near 45 km depth, below the rupture zone of the 2007 event and in the vicinity of the 16 March 1963 compressional event. The fault geometry, rupture process and slip distributions of the two great events are estimated using very broadband teleseismic body and surface wave observations. The occurrence of the thrust event in the shallowest portion of the interplate fault in a region with a paucity of large thrust events at greater depths suggests that the event removed most of the slip deficit on this portion of the interplate fault. This great earthquake doublet demonstrates the heightened seismic hazard posed by induced intraplate faulting following large interplate thrust events. Future seismic failure of the remainder of the seismic gap appears viable, with the northeastern region that has also experienced compressional activity seaward of the megathrust warranting particular attention.

Journal of Geophysical Research B: Solid Earth↗

USGS2012

The glass transition in basalt

The glass transition has been experimentally detected in basalt as (1) an increase in the aggregate linear thermal expansion coefficient α L , (2) an abrupt change in the temperature dependence of Young's modulus dE / dT , and (3) a change in stress relaxation behavior that effectively separates the T > T G and T < T G creep regimes. Transition temperatures determined by the respective experimental methods are (1) 730° ± 15°C, by expansion dilatometry, (2) 725° ± 10°C, by acoustic spectroscopy, and (3) 725° ± 25°C, by stress relaxation. For olivine tholeiite from Kilauea Iki lava lake, Hawaii, thermal expansivity rises through the glass transition interval, attaining α υ = 39.0 × 10 −6 °C −1 at 850°C. For olivine tholeiites with moderate glass contents, the temperature dependence of the Young's modulus changes on either side of the inferred T G of 725°C, and dE / dT = −0.071 kbar °C −1 for T < T G , while dE / dT = − 0.21 kbar °C −1 for T > T G . Collectively, the mechanical results suggest that for Hawaiian olivine tholeiite at 1-atm pressure, the principal material responses are (1) elastic ( T ≤ 600°C), (2) reduced creep (600 < T < 725°C), (3) glass transition ( T ≃ 725°C), (4) enhanced creep (725 < T < 980°C), and (5) partial melt ( T > 980°C). Fracture surface morphologies developed during solidification suggest that the presence of the supercooled melt grain boundary phase may participate in the regulation of the thermal cracking process. Well-preserved fracture surfaces formed by incremental crack growth are found to be covered with striations that correspond to the inferred sequential stopping positions of the advancing fracture front. These striae may be rationalized in terms of experimental analogues that have been produced in viscoelastic polymers by cyclic tension-tension loading. The fracture surface morphology, mechanical data, and the controlled crack growth analogues suggest that thermal fracture in solidifying basalt is an incremental and cyclic process, involving three steps: (1) the accumulation of elastic strain energy in cooling rock at temperatures below that required for stress relaxation due to viscous flow in the intercrystalline liquid phase, (2) fracture at a Δ T determined primarily by the aggregate thermal expansion coefficient α υ and Young's modulus E , (3) the penetration by the advancing crack tip, of the thermal horizon capable of relaxing stress due to the creep of intercrystalline supercooled melt, producing the rough surface texture associated with the termination of a striation. Further crack growth must now await the migration of the solidus. The cycle then repeats. Striations measured in deep Hawaiian lava lakes have been compared with the crack advance increments expected in the vicinity of the glass transition, based on two tests: (1) thermal gradients measured in Kilauea Iki, combined with the mechanical properties of olivine tholeiite evaluated near T G , and (2) the crack advance required to match the recorded seismic stopping phases for prexisting cracks of the dimensions expected for Kilauea Iki. The observed versus predicted comparisons are (1) 31 versus 36 cm; and (2) 31 versus 30 cm. We envision this incremental crack growth process as contributing to the control on the downward movement of the thermal cracking front—and its associated hydrothermal circulation zone—in the upper portions of solidifying subaerial and submarine ponded basalt.

Journal of Geophysical Research Solid Earth↗

USGS2002

An integrated surface-geophysical investigation of the University of Connecticut landfill, Storrs, Connecticut: 2000

A surface-geophysical investigation to characterize the hydrogeology and contaminant distribution of the former landfill area at the University of Connecticut in Storrs, Connecticut, was conducted in 2000 to supplement the preliminary hydrogeologic assessment of the contamination of soil, surface water, and ground water at the site. A geophysical-toolbox approach was used to characterize the hydrogeology and contaminant distribution of the former landfill. Two-dimensional direct-current resistivity, inductive terrain-conductivity, and seismic-refraction surface-geophysical data were collected and interpreted in an iterative manner with exploratory drilling, borehole geophysics, and hydraulic testing. In this investigation, a geophysical-toolbox approach was used to 1) further define previously identified conductive anomalies and leachate plumes; 2) identify additional leachate plumes, possible fracture zones, and (or) conductive lithologic layers in the bedrock; and 3) delineate bedrock-surface topography in the drainage valleys north and south of the landfill. Resistivity and terrain-conductivity surveys were used to further delineate previously identified geophysical anomalies to the north and southwest of the landfill. A conductive anomaly identified in the terrain-conductivity survey to the north of the landfill in 2000 had a similar location and magnitude as an anomaly identified in terrain-conductivity surveys conducted in 1998 and 1999. Collectively, these surveys indicated that the magnitude of the conductive anomaly decreased with depth and with distance from the landfill. These anomalies indicated landfill leachate in the overburden and shallow bedrock. Results of previous surface-geophysical investigations southwest of the landfill indicated a shallow conductive anomaly in the overburden that extended into the fractured-bedrock aquifer. This conductive anomaly had a sheet-like geometry that had a north-south strike, dipped to the west, and terminated abruptly about 450 feet southwest of the landfill. The sheet-like conductive anomaly was interpreted as a fractured, conductive lithologic feature filled with conductive fluids. To further delineate this anomaly, two two-dimensional resistivity profiles were collected west of the sheet-like conductive anomaly to assess the possibility that the sheet-like conductive anomaly continued to the west in its down-dip direction. Each of the north-south oriented resistivity profiles showed bullet-shaped rather than linear-shaped anomalies, with a relatively smaller magnitude of conductivity than the sheet-like conductive anomaly to the east. If these bullet-like features are spatially connected, they may represent a linear, or pipe-like, conductive anomaly in the bedrock with a trend of N290°E and a plunge of 12°. Additional surveys were conducted to assess the apparent southern termination of the sheet-like conductive feature. Terrain-conductivity surveys indicated the sheet-like feature was not continuous to the south. A two-dimensional resistivity line and a coincident terrain-conductivity profile indicated the presence of a steep, eastward dipping, low magnitude, electrically conductive anomaly on the eastern end of the profile. Although the sheet-like conductive anomaly apparently did not continue to the south, the survey conducted in 2000 identified an isolated, weak conductive anomaly south of the previously identified anomaly. Inductive terrain-conductivity surveys performed north of the sheet-like conductive anomaly and west of the landfill indicated the anomaly did not extend to the north into the area of the former chemical-waste disposal pits. No conductive plumes or conductive features were observed in the subsurface bedrock west of the landfill. A conductive anomaly was identified in the southern section of the new terrain-conductivity grid. The magnitude and distribution of the apparent conductivity of this anomaly was identified as a nearly vertical sheet-like conductive feature. The anomaly extended north-south for 150 feet, and the depth and (or) magnitude of the conductive anomaly decreased towards the south. Based on the location, orientation, and dip, this feature was interpreted as a separate anomaly rather than as a continuation of the sheet-like, conductive bedrock feature that was previously identified southwest of the landfill. No other new conductive anomalies were identified south or west of the landfill. Seismic-refraction surveys were used to delineate the depth to the water table and the depth and topography of the bedrock surface. The seismic-refraction surveys at the northern and southern ends of the landfill confirm the presence of shallow bedrock at 25 feet below land surface. Seismic-refraction surveys conducted southwest of the landfill in a minor topographic valley indicate the bedrock is about 10 feet below land surface.

Connecticut↗

USGS2012

The M2 tide on the Amazon Shelf

As part of A Multidisciplinary Amazon Shelf Sediment Study (AMASSEDS), moored and shipboard current measurements made over the Amazon shelf during 1990–1991 have been analyzed to determine the dominant semidiurnal tidal constituent, the M 2 . These results have been combined with coastal sea level data from within the Amazon and Para Rivers, the adjacent shelf, and with satellite‐derived tidal elevation data from off the shelf to provide a more complete description of the M 2 tide in this complex river/shelf system. Near the Amazon River mouth the M 2 tide propagates across the shelf and through the mouth as a damped progressive wave, with its amplitude decreasing and phase increasing upriver. Over the adjacent shelf north of Cabo Norte, the M 2 tide approaches a damped standing wave, with large amplitudes (greater than 1.5 m) near the coast due to near resonance within the coastal embayment formed by the Cabo Norte shoal to the south and Cabo Cassipore to the north. The observed M 2 tidal currents are nearly rectilinear and oriented primarily across the local isobaths. Comparisons between tidal observations in both the North Channel and the Cabo Norte‐Cabo Cassipore embayment and a simple variable‐width channel tidal model indicate that (1) most of the M 2 tidal energy dissipation occurs over the mid‐ and inner shelf (in water depths less than 20 m) and (2) fluid muds found there cause a significant reduction (of order 50%) in the effective bottom friction felt by the M 2 tide. The approximate resonant period of the Cabo Norte‐Cabo Cassipore embayment is 11.9 hours, and at resonance the average energy dissipation per forcing period is roughly 2.2 times the average mechanical energy in the embayment. This damping rate is large enough that the tidal amplification is rather insensitive to forcing frequency, so that the response of the embayment to forcing over the semidiurnal band should be essentially the same. The vertical structure of the M 2 tidal current is examined at one outer shelf site located in 65‐m water depth. The observed semimajor axis increases logarithmically with height above bottom within the lowest 1–2 m and reaches a maximum in excess of 0.5 m/s at approximately 11 m above bottom. The mean ellipticity is small (less than 0.1) and positive, indicating clockwise rotation of a nearly rectilinear current, and the semimajor axis is oriented within 10° of the local cross‐isobath direction. The M 2 phase increases with height above bottom, with flood at the bottom leading flood at the surface by about 1 hour. A simple, local homogeneous tidal model with time‐ and space‐dependent eddy viscosity simulates the observed near‐bottom velocity reasonably well, however, the model suggests that stratification above the lowest few meters may significantly affect the tidal boundary layer structure at this site. The M 2 energy flux onto the Amazon shelf and into the Amazon and Para Rivers has been estimated using current and surface elevation data and the best fit variable‐width channel model results. The net M 2 energy flux into the mouths of the Amazon and Para Rivers is 0.47×10 10 W and 0.19×10 10 W, respectively. A net M 2 energy flux of about 3.3×10 10 W occurs onto the shelf between the North Channel of the Amazon River and Cabo Cassipore. This stretch of the Amazon shelf accounts for about 1.3% of the global dissipation of the M 2 tide.

Journal of Geophysical Research C: Oceans↗