Search USGS⌕ Search

SEARCH · Search USGS

Results for “International 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 487 records · Page 27Linked to original sources
USGS2001

Travel times of P and S from the global digital seismic networks: Implications for the relative variation of P and S velocity in the mantle

We present new data sets of P and S arrival times which have been handpicked from long-period vertical and transverse component recordings of the various global seismic networks. Using events which occurred from 1976 to 1994 results in ∼38,000 globally well-distributed measurements of teleseismic P and ∼41,000 measurements of S . These data are particularly useful for looking at the relative variation of S and P velocities in the lower mantle. We describe both the measurement techniques and the gross characteristics of the data sets. The size of our data sets allows us to exploit the internal consistency of the data to identify outliers using a summary ray analysis. Since the polarity of each arrival is also known, we can construct fault plane solutions and/or compare with polarities predicted by the Harvard centroid moment tensor solutions to further diagnose phase misidentification. This analysis results in ∼5% of the data being identified as outliers. An analysis of variance indicates that the S residual travel times are dominated by the effects of three-dimensional structure but the P data have comparable contributions from noise and source mislocation effects. The summary ray analysis reveals the basic character of lower mantle structure, and there are large-scale patterns in both the S and P data sets that correlate quite well with each other. This analysis suggests that on average, d ln v s / d ln v p is an increasing function of depth in the mantle going from a value of ∼1.7 at the top of the lower mantle to an apparent value of 4 near the base of the mantle. This latter extreme value of R seems to result mainly from data which sample one region in the lowermost mantle under the central Pacific, where large positive S residuals are associated with very small P residuals. Such an anomaly cannot be thermal in origin.

Journal of Geophysical Research B: Solid Earth↗
USGS2016

San Pedro River Aquifer Binational Report

The United States and Mexico share waters in a number of hydrological basins and aquifers that cross the international boundary. Both countries recognize that, in a region of scarce water resources and expanding populations, a greater scientific understanding of these aquifer systems would be beneficial. In light of this, the Mexican and U.S. Principal Engineers of the International Boundary and Water Commission (IBWC) signed the “Joint Report of the Principal Engineers Regarding the Joint Cooperative Process United States-Mexico for the Transboundary Aquifer Assessment Program" on August 19, 2009 (IBWC-CILA, 2009). This IBWC “Joint Report” serves as the framework for U.S.-Mexico coordination and dialogue to implement transboundary aquifer studies. The document clarifies several details about the program such as background, roles, responsibilities, funding, relevance of the international water treaties, and the use of information collected or compiled as part of the program. In the document, it was agreed by the parties involved, which included the IBWC, the Mexican National Water Commission (CONAGUA), the U.S. Geological Survey (USGS), and the Universities of Arizona and Sonora, to study two priority binational aquifers, one in the San Pedro River basin and the other in the Santa Cruz River basin. This report focuses on the Binational San Pedro Basin (BSPB). Reasons for the focus on and interest in this aquifer include the fact that it is shared by the two countries, that the San Pedro River has an elevated ecological value because of the riparian ecosystem that it sustains, and that water resources are needed to sustain the river, existing communities, and continued development. This study describes the aquifer’s characteristics in its binational context; however, most of the scientific work has been undertaken for many years by each country without full knowledge of the conditions on the other side of the border. The general objective of this study is to use new and existing research to define the general hydrologic framework of the Binational San Pedro Aquifer (BSPA), to gather hydrogeological and other relevant data in preparation for future work such as an updated groundwater conceptual model and budget and to establish the basis for a binational numerical model. The specific objectives are as follows: Understand the current state of knowledge with respect to climate, geology, soils, land cover, land use, and hydrology of the aquifer in its binational context; Compile and create a database of scientific information from both countries; Identify data gaps and identify what data would be necessary to update, in a subsequent phase, the hydrologic model of the aquifer system, including surface- and groundwater interactions on a binational level. The BSPB is one of the most studied basins in the region, and a database of publications has been compiled as part of this project. Previous studies include topics that range from geophysics and hydrogeology to biology and ecosystem services. The economic drivers on each side of the border are quite different. In the Arizona 4 portion of the basin military and tourism dominate while in the Sonoran portion, mining is the most important industry. Water management is also different in the two countries. In Mexico, primary authority for management of water resources devolves from the federal government. In the United States, primary authority rests with the states except in cases of interstate surface waters. Binational waters are not currently jointly managed by the two countries except in cases where treaties have been negotiated such as for the Rio Grande and Colorado Rivers. Thus, there is currently no binational coordination or treaty governing the management of groundwater.

Report↗
USGS1994

An image of the Columbia Plateau from inversion of high‐resolution seismic data

We use a method of traveltime inversion of high‐resolution seismic data to provide the first reliable images of internal details of the Columbia River Basalt Group (CRBG), the subsurface basalt/sediment interface, and the deeper sediment/basement interface. Velocity structure within the basalts, delineated on the order of 1 km horizontally and 0.2 km vertically, is constrained to within ±0.1 km/s for most of the seismic profile. Over 5000 observed traveltimes fit our model with an rms error of 0.018 s. The maximum depth of penetration of the basalt diving waves (truncated by underlying low‐velocity sediments) provides a reliable estimate of the depth to the base of the basalt, which agrees with well‐log measurements to within 0.05 km (165 ft). We use image blurring, calculated from the resolution matrix, to estimate the aspect ratio of imaged velocity anomaly widths to true widths for velocity features within the basalt. From our calculations of image blurring, we interpret low velocity zones (LVZ) within the basalts at Boylston Mountain and the Whiskey Dick anticline to have widths of 4.5 and 3 km, respectively, within the upper 1.5 km of the model. At greater depth, the widths of these imaged LVZs thin to approximately 2 km or less. We interpret these linear, subparallel, low‐velocity zones imaged adjacent to anticlines of the Yakima Fold Belt to be brecciated fault zones. These fault zones dip to the south at angles between 15 to 45 degrees.

Geophysics↗
USGS2014

Differential response of carbon fluxes to climate in three peatland ecosystems that vary in the presence and stability of permafrost

Changes in vegetation and soil properties following permafrost degradation and thermokarst development in peatlands may cause changes in net carbon storage. To better understand these dynamics, we established three sites in Alaska that vary in permafrost regime, including a black spruce peat plateau forest with stable permafrost, an internal collapse scar bog formed as a result of thermokarst, and a rich fen without permafrost. Measurements include year-round eddy covariance estimates of carbon dioxide (CO 2 ), water, and energy fluxes, associated environmental variables, and methane (CH 4 ) fluxes at the collapse scar bog. The ecosystems all acted as net sinks of CO 2 in 2011 and 2012, when air temperature and precipitation remained near long-term means. In 2013, under a late snowmelt and late leaf out followed by a hot, dry summer, the permafrost forest and collapse scar bog were sources of CO 2 . In this same year, CO 2 uptake in the fen increased, largely because summer inundation from groundwater inputs suppressed ecosystem respiration. CO 2 exchange in the permafrost forest and collapse scar bog was sensitive to warm air temperatures, with 0.5 g C m −2 lost each day when maximum air temperature was very warm (≥29°C). The bog lost 4981 ± 300 mg CH 4 m −2 between April and September 2013, indicating that this ecosystem acted as a significant source of both CO 2 and CH 4 to the atmosphere in 2013. These results suggest that boreal peatland responses to warming and drying, both of which are expected to occur in a changing climate, will depend on permafrost regime.

Alaska↗
USGS2000

Geologic mapping of Europa

Galileo data enable the major geological units, structures, and surface features to be identified on Europa. These include five primary units (plains, chaos, band, ridge, and crater materials) and their subunits, along with various tectonic structures such as faults. Plains units are the most widespread. Ridged plains material spans a wide range of geological ages, including the oldest recognizable features on Europa, and appears to represent a style of tectonic resurfacing, rather than cryovolcanism. Smooth plains material typically embays other terrains and units, possibly as a type of fluid emplacement, and is among the youngest material units observed. At global scales, plains are typically mapped as undifferentiated plains material, although in some areas differences can be discerned in the near infrared which might be related to differences in ice grain size. Chaos material is composed of plains and other preexisting materials that have been severely disrupted by inferred internal activity; chaos is characterized by blocks of icy material set in a hummocky matrix. Band material is arrayed in linear, curvilinear, wedge-shaped, or cuspate zones with contrasting albedo and surface textures with respect to the surrounding terrain. Bilateral symmetry observed in some bands and the relationships with the surrounding units suggest that band material forms by the lithosphere fracturing, spreading apart, and infilling with material derived from the subsurface. Ridge material is mapped as a unit on local and some regional maps but shown with symbols at global scales. Ridge material includes single ridges, doublet ridges, and ridge complexes. Ridge materials are considered to represent tectonic processes, possibly accompanied by the extrusion or intrusion of subsurface materials, such as diapirs. The tectonic processes might be related to tidal flexing of the icy lithosphere on diurnal or longer timescales. Crater materials include various interior (smooth central, rough inner, and annular massif) and exterior (continuous ejecta) subunits. Structural features and landforms are shown with conventional symbols. Type localities for the units are identified, along with suggestions for portraying the features on geological maps, including colors and letter abbreviations for material units. Implementing these suggestions by the planetary mapping community would facilitate comparisons of maps for different parts of Europa and contribute to an eventual global synthesis of its complex geology. On the basis of initial mapping results, a stratigraphic sequence is suggested in which ridged plains form the oldest unit on Europa, followed by development of band material and individual ridges. Band materials tend to be somewhat older than ridges, but in many areas the two units formed simultaneously. Similarly, the formation of most chaos follows the development of ridged plains; although chaos is among the youngest materials on Europa, some chaos units might have formed contemporaneously with ridged plains. Smooth plains generally embay all other units and are late-stage in the evolution of the surface. C1 craters are superposed on ridged plains but are crosscut by other materials, including bands and ridges. Most c2 craters postdate all other units, but a few c2 craters are cut by ridge material. C3 craters constitute the youngest recognizable material on Europa. Copyright 2000 by the American Geophysical Union.

Journal of Geophysical Research E: Planets↗
USGS1999

Formation of natural gas hydrates in marine sediments 1. Conceptual model of gas hydrate growth conditioned by host sediment properties

The stability of submarine gas hydrates is largely dictated by pressure and temperature, gas composition, and pore water salinity. However, the physical properties and surface chemistry of deep marine sediments may also affect the thermodynamic state, growth kinetics, spatial distributions, and growth forms of clathrates. Our conceptual model presumes that gas hydrate behaves in a way analogous to ice in a freezing soil. Hydrate growth is inhibited within fine-grained sediments by a combination of reduced pore water activity in the vicinity of hydrophilic mineral surfaces, and the excess internal energy of small crystals confined in pores. The excess energy can be thought of as a "capillary pressure" in the hydrate crystal, related to the pore size distribution and the state of stress in the sediment framework. The base of gas hydrate stability in a sequence of fine sediments is predicted by our model to occur at a lower temperature (nearer to the seabed) than would be calculated from bulk thermodynamic equilibrium. Capillary effects or a build up of salt in the system can expand the phase boundary between hydrate and free gas into a divariant field extending over a finite depth range dictated by total methane content and pore-size distribution. Hysteresis between the temperatures of crystallization and dissociation of the clathrate is also predicted. Growth forms commonly observed in hydrate samples recovered from marine sediments (nodules, and lenses in muds; cements in sands) can largely be explained by capillary effects, but kinetics of nucleation and growth are also important. The formation of concentrated gas hydrates in a partially closed system with respect to material transport, or where gas can flush through the system, may lead to water depletion in the host sediment. This "freeze-drying" may be detectable through physical changes to the sediment (low water content and overconsolidation) and/or chemical anomalies in the pore waters and metastable presence of free gas within the normal zone of hydrate stability.

Journal of Geophysical Research B: Solid Earth↗
USGS2010

After a century: Revised Paleogene coal stratigraphy, correlation, and deposition, Powder River Basin, Wyoming and Montana

The stratigraphy, correlation, mapping, and depositional history of coal-bearing strata in the Paleogene Fort Union and Wasatch Formations in the Powder River Basin were mainly based on measurement and description of outcrops during the early 20th century. Subsequently, the quality and quantity of data improved with (1) exploration and development of oil, gas, and coal during the middle 20th century and (2) the onset of coalbed methane (CBM) development during the late 20th and early 21st centuries that resulted in the drilling of more than 26,000 closely spaced wells with accompanying geophysical logs. The closeness of the data control points, which average 0.5 mi (805 m) apart, made for better accuracy in the subsurface delineation and correlation of coal beds that greatly facilitated the construction of regional stratigraphic cross sections and the assessment of resources. The drillhole data show that coal beds previously mapped as merged coal zones, such as the Wyodak coal zone in the Wyoming part of the Powder River Basin, gradually thinned into several discontinuous beds and sequentially split into as many as 7 hierarchical orders westward and northward. The thinning and splitting of coal beds in these directions were accompanied by as much as a ten-fold increase in the thicknesses of sandstone-dominated intervals within the Wyodak coal zone. This probably resulted from thrust loading by the eastern front of the Bighorn uplift accompanied by vertical displacement along lineaments that caused subsidence of the western axial part of the Powder River Basin during Laramide deformation in Late Cretaceous and early Tertiary time. Accommodation space was thereby created for synsedimentary alluvial infilling that controlled thickening, thinning, splitting, pinching out, and areal distribution of coal beds. Equally important was differential subsidence between this main accommodation space and adjoining areas, which influenced the overlapping, for example, of the Dietz coal zone in Montana, over the Wyodak coal zone in Wyoming. Correlation in a circular track of the Wyodak coal zone in the southern part of the basin also demonstrates overlapping with lower coal zones. Recognition of this stratigraphic relationship has led to revision of the correlations and nomenclature of coal beds because of inconsistency within these zones as well as those below and above them, which have long been subjects of controversy. Also, it significantly changes the traditional coal bed-to-bed correlations, and estimates of coal and coalbed methane resources of these coal zones due to thinning and pinching out of beds. More notably, thickness isopach, orientation, and distribution of the merged Wyodak coal bodies in the south-southeast part of the basin suggest that differential movement of lineament zones active during the Cretaceous was not a major influence on coal accumulation during the Paleocene. Improved knowledge of alluvial depositional environments as influenced by external and internal paleotectonic conditions within the Powder River Basin permits more accurate correlation, mapping, and resource estimation of the Fort Union and Wasatch coal beds. The result is a better understanding of the sedimentology of the basin infill deposits in relation to peat bog accumulation.

Montana, Wyoming↗
USGS2019

Protracted multipulse emplacement of a post-resurgent pluton: The case of Platoro caldera complex (Southern Rocky Mountain volcanic field, Colorado)

Many eroded calderas expose associated postcollapse plutons, but detailed fieldwork‐supported studies have rarely focused on the internal structure that can contribute to understanding of emplacement dynamics. The Alamosa River monzonite pluton is a postcollapse intrusion at the Platoro caldera complex that erupted six large ignimbrites between 30.2 and 28.8 Ma in the Southern Rocky Mountains volcanic field. Magnetic fabrics in this intrusion indicate the pulsed emplacement of a vertically extensive pluton. The magmatic pulses are documented by three concentric domains of magnetic foliations elongated in ~NE‐SW direction, corresponding to structural trends at the Platoro caldera complex and preexisting regional structures. As no evidence for deformation of wall rocks and the adjacent resurgent block has been identified, we interpret the Alamosa River pluton as a postresurgent intrusion. The space‐opening process involved magmatic stoping and small‐scale magma wedging. New SHRIMP‐RG U/Pb zircon dates (28.98 ±0.18, 27.42 ± 0.35, and 27.32 ± 0.38 Ma) suggest a magmatic lifespan of ~1.7 My for the Alamosa River pluton. Our results indicate that postcaldera magmatism includes pulsed and protracted activity from large intracaldera resurgent plutons to smaller postresurgent stocks and sheeted complexes. As demonstrated by the Alamosa River pluton, some intrusions are emplaced shortly after collapse and resurgence, but postcaldera volcano‐plutonic systems may remain active for several million years or more. We also suggest that subvolcanic magma bodies may be assembled incrementally and that the record of early composite magma lenses preserved as magma wedges are later obliterated by convective flowage and crystallization.

Colorado↗
USGS1993

Relationship of geological and geothermal field properties: Midcontinent area, USA, an example

Quantitative approaches to data analysis in the last decade have become important in basin modeling and mineral-resource estimation. The interrelation of geological, geophysical, geochemical, and geohydrological variables is important in adjusting a model to a real-world situation. Revealing the interdependences of variables can contribute in understanding the processes interacting in sedimentary basins. It is reasonably simple to compare spatial data of the same type but more difficult if different properties are involved. Statistical techniques, such as cluster analysis or principal components analysis, or some algebraic approaches can be used to ascertain the relations of standardized spatial data. In this example, structural configuration on five different stratigraphic horizons, one total sediment thickness map, and four maps of geothermal data were copared. As expected, the structural maps are highly related because all had undergone about the same deformation with differing degrees of intensity. The temperature gradients derived (1) from shallow borehole logging measurements under equilibrium conditions with the surrounding rock, and (2) from non-equilibrium bottom-hole temperatures (BHT) from deeper depths are mainly independent of each other. This was expected and confirmed also for the two temperature maps at 1000 ft which were constructed using both types of gradient values. Thus, it is evident that the use of a 2-point (BHT and surface temperature) straightline calculation of a mean temperature gradient gives different information about the geothermal regime than using gradients from temperatures logged under equilibrium conditions. Nevertheless, it is useful to determine to what a degree the larger dataset of nonequilibrium temperatures could reflect quantitative relationships to geologic conditions. Comparing all maps of geothermal information vs. the structural and the sediment thickness maps, it was determined that all correlations are moderately negative or slightly positive. These results are clearly shown by the cluster analysis and the principal components. Considering a close relationship between temperature and thermal conductivity of the sediments as observed for most of the Midcontinent area and relatively homogeneous heat-flow density conditions for the study area these results support the following assumptions: (1) undifferentiated geothermal gradients, computed from temperatures of different depth intervals and differing sediment properties, cannot contribute to an improved understanding of the temperature structure and its controls within the sedimentary cover, and (2) the quantitative approach of revealing such relations needs refined datasets of temperature information valid for the different depth levels or stratigraphic units. ?? 1993 International Association for Mathematical Geology.

Mathematical Geology↗
USGS2012

Internal structure of the Sierra Nevada batholith based on specific gravity and gravity measurements

About 6,000 specific‐gravity (SG) measurements of samples collected from nearly 200 granitic plutons comprising the central Sierra Nevada batholith yield a SG contour map across the batholith from 36.25° to 38° north latitude. With notable exceptions, SG decreases from values generally greater than 2.7 in the west to less than 2.6 over a few small areas of high‐silica, high‐potassium granites near the east edge. A good correlation between measured SG and analyzed weight percent SiO 2 enables estimation of average silica variations across the batholith. The average SG is 2.69 corresponding to an average of 68 wt. % SiO 2 for the 18,000 km² central part of the batholith. A 1‐km gridded version of the SG measurements has been used to generate a series of synthetic gravity maps, assuming that the SG of rocks at the surface extends unchanged to various depths. The synthetic map computed for a depth of 10 km shows the best correspondence with the isostatic residual gravity map indicating that variations in the observed gravity residuals are largely caused by SG variations of the plutonic rocks exposed at the surface that apparently extend downward to an average depth of about 10 km. Although the 10‐km synthetic gravity map gives the best overall fit to the observed gravity data, comparison of individual anomalies indicates that the bottoms of the plutons as defined by SG variations at the surface are generally shallower along the west edge of the Sierra Nevada (7±2 km) and deeper in the younger and more felsic eastern part (12±3 km). These depths do not necessarily represent a distinct base of the Sierra Nevada batholith. They may indicate the depth below which density homogenization occurs, either by igneous, or possibly, structural processes.

California↗
USGS2016

Lithospheric flexure under the Hawaiian volcanic load: Internal stresses and a broken plate revealed by earthquakes

Several lines of earthquake evidence indicate that the lithospheric plate is broken under the load of the island of Hawai`i, where the geometry of the lithosphere is circular with a central depression. The plate bends concave downward surrounding a stress-free hole, rather than bending concave upward as with past assumptions. Earthquake focal mechanisms show that the center of load stress and the weak hole is between the summits of Mauna Loa and Mauna Kea where the load is greatest. The earthquake gap at 21 km depth coincides with the predicted neutral plane of flexure where horizontal stress changes sign. Focal mechanism P axes below the neutral plane display a striking radial pattern pointing to the stress center. Earthquakes above the neutral plane in the north part of the island have opposite stress patterns; T axes tend to be radial. The M 6.2 Honomu and M 6.7 Kiholo main shocks (both at 39 km depth) are below the neutral plane and show radial compression, and the M 6.0 Kiholo aftershock above the neutral plane has tangential compression. Earthquakes deeper than 20 km define a donut of seismicity around the stress center where flexural bending is a maximum. The hole is interpreted as the soft center where the lithospheric plate is broken. Kilauea's deep conduit is seismically active because it is in the ring of maximum bending. A simplified two-dimensional stress model for a bending slab with a load at one end yields stress orientations that agree with earthquake stress axes and radial P axes below the neutral plane. A previous inversion of deep Hawaiian focal mechanisms found a circular solution around the stress center that agrees with the model. For horizontal faults, the shear stress within the bending slab matches the slip in the deep Kilauea seismic zone and enhances outward slip of active flanks.

Hawaii↗
USGS1998

Closing the gap between regional and global travel time tomography

Recent global travel time tomography studies by Zhou [1996] and van der Hilst et al . [1997] have been performed with cell parameterizations of the order of those frequently used in regional tomography studies (i.e., with cell sizes of 1°–2°). These new global models constitute a considerable improvement over previous results that were obtained with rather coarse parameterizations (5° cells). The inferred structures are, however, of larger scale than is usually obtained in regional models, and it is not clear where and if individual cells are actually resolved. This study aims at resolving lateral heterogeneity on scales as small as 0.6° in the upper mantle and 1.2°–3° in the lower mantle. This allows for the adequate mapping of expected small-scale structures induced by, for example, lithosphere subduction, deep mantle upwellings, and mid-ocean ridges. There are three major contributions that allow for this advancement. First, we employ an irregular grid of nonoverlapping cells adapted to the heterogeneous sampling of the Earth's mantle by seismic waves [ Spakman and Bijwaard , 1998]. Second, we exploit the global data set of Engdahl et al . [1998], which is a reprocessed version of the global data set of the International Seismological Centre. Their reprocessing included hypocenter redetermination and phase reidentification. Finally, we combine all data used ( P , pP , and pwP phases) into nearly 5 million ray bundles with a limited spatial extent such that averaging over large mantle volumes is prevented while the signal-to-noise ratio is improved. In the approximate solution of the huge inverse problem we obtain a variance reduction of 57.1%. Synthetic sensitivity tests indicate horizontal resolution on the scale of the smallest cells (0.6° or 1.2°) in the shallow parts of subduction zones decreasing to approximately 2°–3° resolution in well-sampled regions in the lower mantle. Vertical resolution can be worse (up to several hundreds of kilometers) in subduction zones with rays predominantly pointing along dip. Important features of the solution are as follows: 100–200 km thick high-velocity slabs beneath all major subduction zones, sometimes flattening in the transition zone and sometimes directly penetrating into the lower mantle; large high-velocity anomalies in the lower mantle that have been attributed to subduction of the Tethys ocean and the Farallon plate; and low-velocity anomalies continuing across the 660 km discontinuity to hotspots at the surface under Iceland, east Africa, the Canary Islands, Yellowstone, and the Society Islands. Our findings corroborate that the 660 km boundary may resist but not prevent (present day) large-scale mass transfer from upper to lower mantle or vice versa. This observation confirms the results of previous, global mantle studies that employed coarser parameterizations.

Journal of Geophysical Research B: Solid Earth↗
USGS2004

Paleomagnetism and 40Ar/39Ar ages from volcanics extruded during the Matuyama and Brunhes Chrons near McMurdo Sound, Antarctica

Maps of virtual geomagnetic poles derived from international geomagnetic reference field models show large lobes with significant departures from the spin axis. These lobes persist in field models for the last few millenia. The anomalous lobes are associated with observation sites at extreme southerly latitudes. To determine whether these features persist for millions of years, paleomagnetic vector data from the continent of Antarctica are essential. We present here new paleomagnetic vector data and 40 Ar/ 39 Ar ages from lava flows spanning the Brunhes and Matuyama Chrons from the vicinity of McMurdo Sound, Antarctica. Oriented paleomagnetic samples were collected from 50 lava flows by E. Mankinen and A. Cox in the 1965–1966 austral summer season. Preliminary data based largely on the natural remanent magnetization (NRM) directions were published by Mankinen and Cox [1988]. We have performed detailed paleomagnetic investigations of 37 sites with multiple fully oriented core samples to investigate the reliability of results from this unique sample collection. Of these, only one site fails to meet our acceptance criteria for directional data. Seven sites are reversely magnetized. The mean normal and reverse directions are antipodal. The combined mean direction has = 12, = −86, α = 4, κ = 37 and is indistinguishable from that expected from a GAD field. We obtained reproducible absolute paleointensity estimates from 15 lava flows with a mean dipole moment of 49 ZAm 2 and a standard deviation of 28 ZAm 2 . 40 Ar/ 39 Ar age determinations were successfully carried out on samples from 18 of the flows. Our new isotopic ages and paleomagnetic polarities are consistent with the currently accepted geomagnetic reversal timescales.

Geochemistry, Geophysics, Geosystems↗
USGS2021

Evaluating the state-of-the-art in remote volcanic eruption characterization Part II: Ulawun volcano, Papua New Guinea

Retrospective eruption characterization is valuable for advancing our understanding of volcanic systems and evaluating our observational capabilities, especially with remote technologies (defined here as a space-borne system or non-local, ground-based instrumentation which include regional and remote infrasound sensors). In June 2019, the open-system Ulawun volcano, Papua New Guinea, produced a VEI 4 eruption. We combined data from satellites (including Sentinel-2, TROPOMI, MODIS, Himawari-8), the International Monitoring System infrasound network, and GLD360 globally detected lightning with information from the local authorities and social media to characterize the pre-, syn - and post-eruptive behaviour. The Rabaul Volcano Observatory recorded ~24 h of seismicity and detected SO 2 emissions ~16 h before the visually-documented start of the Plinian phase on 26 June at 04:20 UTC. Infrasound and SO 2 detections suggest the eruption started during the night on 24 June 2019 at 10:39 UTC ~38 h before ash detections with a gas-dominated jetting phase. Local reports and infrasound detections show that the second phase of the eruption started on 25 June 19:28 UTC with ~6 h of jetting. The first detected lightning occurred on 26 June 00:14 UTC, and ash emissions were first detected by Himawari-8 at 01:00 UTC. Post-eruptive satellite imagery indicates new flow deposits to the south and north of the edifice and ash fall to the west and southwest. In particular, regional infrasound data provided novel insight into eruption onset and syn -eruptive changes in intensity. We conclude that, while remote observations are sufficient for detection and tracking of syn-eruptive changes, key challenges in data latency, acquisition, and synthesis must be addressed to improve future near-real-time characterization of eruptions at minimally-monitored or unmonitored volcanoes.

Bismark volcanic arc, New Britain Island, Ulawun v↗
USGS2012

Variations in the reflectivity of the moho transition zone beneath the Midcontinent Rift System of North America: results from true amplitude analysis of GLIMPCE data

True amplitude processing of The Great Lakes International Multidisciplinary Program on Crustal Evolution seismic reflection data from the Midcontinent Rift System of North America shows large differences in the reflectivity of the Moho transition zone beneath the axial rift, beneath the rift flanks, and outside of the rift. The Moho reflection from the axial rift has a discontinuous, diffractive character marginally stronger (several decibels) than an otherwise transparent lower crust and upper mantle. Beneath the axial rift, Moho is interpreted to be a synrift igneous feature. Beneath the rift flanks, the reflectivity of the Moho transition is generally well developed with two identifiable boundaries, although in places it is weakly reflective to nonreflective, similar to Moho outside the rift. The two boundaries are interpreted as the base of essentially intact, although stretched, prerift Archean crust (upper boundary) and new synrift Moho 1-2 s (6-7 km) deeper (lower boundary). Beneath the rift flanks, the layered reflection Moho transition results from the preexisting crustal composition and fabric modified by synrift igneous processes and extensional tectonic/metamorphic processes. The geologic evidence for extensive basaltic magmatism in the rift is the basis for interpreting the Moho signature as a Keweenawan structure that has been preserved for 1.1 b.y. Extension and magmatism appear to enhance reflectivity in the lower crust and Moho transition zone only where stretching factors are moderate (rift flanks) and not where they are extreme (axial rift). This leads to the prediction that the reflectivity across analogous volcanic passive continental margins should be greatest beneath the moderately stretched continental shelves and should decrease towards the ocean-continent boundary.

Journal of Geophysical Research B: Solid Earth↗
USGS2006

Metallogeny of the Great Basin: Crustal evolution, fluid flow, and ore deposits

The Great Basin physiographic province in the Western United States contains a diverse assortment of world-class ore deposits. It currently (2006) is the world’s second leading producer of gold, contains large silver and base metal (Cu, Zn, Pb, Mo, W) deposits, a variety of other important metallic (Fe, Ni, Be, REE’s, Hg, PGE) and industrial mineral (diatomite, barite, perlite, kaolinite, gallium) resources, as well as petroleum and geothermal energy resources. Ore deposits are most numerous and largest in size in linear mineral belts with complex geology. U.S. Geological Survey (USGS) scientists are in the final year of a research project initiated in the fall of 2001 to increase understanding of relations between crustal evolution, fluid flow, and ore deposits in the Great Basin. Because of its substantial past and current mineral production, this region has been the focus of numerous investigations over the past century and is the site of ongoing research by industry, academia, and state agencies. A variety of geoinformatic tools was used to organize, reinterpret, and display, in space and time, the large amounts of geologic, geophysical, geochemical, and hydrologic information deemed pertinent to this problem. This information, in combination with concentrated research on (1) critical aspects of the geologic history, (2) an area in northern Nevada that encompasses the major mineral belts, and (3) important mining districts and deposits, is producing new insights about the interplay between key tectonic events, hydrothermal fluid flow, and ore genesis in mineral belts. The results suggest that the Archean to Holocene history of the Great Basin was punctuated by several tectonic events that caused fluids of different origins (sea water, basinal brine, meteoric water, metamorphic water, magmatic water) to move through the crust. Basement faults reactivated during these events localized deformation, sedimentation, magmatism, and hydrothermal fluid flow in overlying rocks to form mineral belts that contain ore deposits of different types and ages that are locally superimposed (demonstrating inheritance). Fluid flow in these systems also was influenced by the distribution of permeable lithologies and paleotopographic highs and lows. Hydrothermal fluids evolved from their initial chemistries towards compositions that reflect the ƒ O 2 and ƒ S 2 buffering capacity of, and the ligands and metals present in, the rocks (±older mineralization) through which they moved. In northern Nevada, where gold deposits are relatively common, carbonaceous, pyritic strata buffered fluids of diverse origins to H 2 S-rich compositions so they could transport gold repeatedly over Paleozoic-Cenozoic time (convergent evolution). Ore formed where metal-laden fluids encountered effective physicochemical traps. Maps of Neogene basin fill and erosion surfaces identify areas where preexisting ore deposits have been progressively exposed or concealed. Comparisons with analogous terrains and deposit types in other parts of the world provide global context. The initial findings and some of the databases, geologic maps, sections, reconstructions, hydrogeologic models, topical syntheses, regional overviews, short courses, field guides, and deposit comparisons produced by project staff and associated managers, contractors, and collaborators have been presented in numerous abstracts, symposia, USGS publications, and professional journals over the last 5 years (see the extensive bibliography). Notable among these was the 2005 Geological Society of Nevada symposium in Reno, Nevada, and the 2005 Geological Society of America annual meeting in Salt Lake City, Utah, where project results were presented to audiences from around the nation and world. The final results of the project will be submitted for publication in 2007 to appropriate USGS and professional journals. A special issue of GEOSPHERE, scheduled for publication in 2007, will be devoted to the results of this project and related work. This special issue will reach an international audience and be available worldwide on the internet. Much of the research for this project has concentrated on areas that will receive the focused attention of the mining industry in the future. As such, the data and interpretations generated by this project have direct use for land-use managers in Federal, State, and local agencies. Improved hydrogeologic models developed by this project will considerably enhance ongoing and future water resource investigations in the region. The increased understanding of when, where, and how hydrothermal systems produce significant economic deposits has direct uses for mineral exploration and for future USGS mineral resource assessments in the Great Basin and other parts of the world.

Great Basin↗
USGS1999

Debris-flow origin for the Simud/Tiu deposit on Mars

A late Hesperian smooth plains deposit on Mars interpreted as a debris flow extends more than 2000 km from Hydraotes Chaos, through Simud and Tiu Valles, and into Chryse Planitia. The Simud/Tiu deposit widens out to >1000 km and embays streamlined landforms and knobs made up of sedimentary and perhaps volcanic deposits that were carved by earlier channeling activity. Morphologic features of the Simud/Tiu deposit observed in Viking and Pathfinder images are generally consistent with a debris-flow origin, but some of the deposit's salient features are not readily explained by catastrophic flooding or ice flow. Internal depressions appear to be bounded by linear scarps along flow margins where differential shearing may have occurred and in areas where flow spreading may have produced zones of extensional breakup and thinning within the flow. Possible flow lobes within the deposit may have formed by successive flow surges within the flow unit. The Pathfinder landing site is on the Simud/Tiu deposit, and the observations there are consistent with debris flow. The low, longitudinal ridges at the site may have formed by clast interactions as the flow ground to a halt. Imbricated, planar rocks on the ridges, such as in the Rock Garden, also may have been emplaced by debris or ice flow. However, stream energy calculations at Ares Vallis and channel geology indicate that flooding probably was incapable of emplacing the meter-size boulders observed at the Pathfinder site. Dewatering of pressurized zones in the debris flow or underlying material may be responsible for mud eruptions that formed a couple of patches of low pancakelike shields up to 5 km in diameter and for probable water flows that formed two small rille channels a few kilometers long. Local irregular grooves may be cracks that resulted from later desiccation and contraction of the flow material. The debris-flow unit apparently coalesced from outflows of water-fluidized debris originating from beneath chaotic and hummocky terrains within and along the margins of Simud and Tiu Valles. The deposit is onlapped from the north by another flow deposit originating from Acidalia Planitia. If the Simud/Tiu debris flow had entered a standing body of water, a turbidity current may have arisen from the debris flow and then backflowed over the debris flow to account for the Acidalia deposit.

Journal of Geophysical Research E: Planets↗
USGS2012

Lithology and evolution of the crust-mantle boundary region in the southwestern Basin and Range Province

Mantle and crustal xenoliths from volcanic rocks in the southwestern Basin and Range province and Colorado Plateau Transition Zone reveal histories of episodic magmatism and deformation that have profoundly influenced the crustal structure of this region. Seismic transects in this area show a strongly reflective Moho of generally low relief, which, in the area of modern transects, consists of a thin zone (<2 km thick) of short reflectors. The upper mantle is transparent and has a P n of 7.8–8.0 km/s similar to much of the western United States. A lower crustal zone, 2–13 km thick, has variable internal reflectivity and a relatively low velocity of 6.6–6.8 km/s. Upper mantle peridotite xenoliths show both ductile and brittle deformational features and have structures and compositions affected by magmatic intrusion; intrusions form complex dike systems and extensive zones of grain boundary infiltration in peridotite xenoliths. Whereas melt infiltration preceded and followed ductile deformation, brittle deformation, represented by closely spaced joint systems and faults, followed ductile deformation and is related to the youngest magmatic episodes. These structural characteristics and high uppermost mantle temperature (∼1000°C) may combine to explain the relatively low P n . Alternating layers of ductily deformed and undeformed peridotites, with or without igneous intrusions, may contribute to the reflectivity of the Moho. Lower crustal xenoliths are dominantly igneous-textured pyroxenites and mafic to intermediate gabbros identical to the dikes in peridotite xenoliths. The crustal xenoliths also commonly are jointed, and in addition many show partial melting and have abundant cavities that probably were filled with CO 2 -rich fluids. These rocks are interpreted as products of underplated magmas that were fed through the mantle dike systems and may represent the lowest crustal unit identified in the seismic records. The mafic compositions and high densities of the crustal xenoliths indicate that the low velocity of the lower crust may be caused in part by fracture systems, partial melts, and high temperatures. Garnet granulite xenoliths from a locality with no mantle peridotite xenoliths probably represent crust of the region before late Miocene extension. Felsic granulite xenoliths from two localities have velocities like those of the two lower crustal units identified seismically and could be present in the modern crust as unequilibrated remnants of old crust. The preferred model for the evolution of the lower lithosphere is one in which extension affects the upper mantle as well as the crust and is overlapped in time by multiple magmatic episodes. The earliest magmatic events preceded extension, and later events accompanied and followed extension.

Arizona, California, Nevada↗