USGS Science⌕ Search

SEARCH · USGS Science

Results for “International Journal of 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 253 records · Page 14Linked 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↗
USGS1980

Geophysical investigations in deep horizontal holes drilled ahead of tunnelling

Deep horizontal drill holes have been used since 1967 by the Defense Nuclear Agency as a primary exploration tool for siting nuclear events in tunnels at the Nevada Test Site. The U.S. Geological Survey had developed geophysical logging techniques for obtaining resistivity and velocity in these holes, and to date 33 horizontal drill holes in excess of 300 m in depth have been successfully logged. The deepest hole was drilled to a horizontal depth of 1125 m. The purposes of the logging measurements are to define clay zones, because of the unstable ground conditions such zones can present to tunnelling, and to define zones of partially saturated rock, because of the attenuating effects such zones have on the shock wave generated by the nuclear detonation. Excessive attenuation is undesirable because the shock wave is used as a tunnel closure mechanism to contain debris and other undesirable explosion products. Measurements are made by pumping resistivity, sonic and geophone probes down the drill string and out of the bit into the open hole. Clay zones are defined by the electrical resistivity technique based on empirical data relating the magnitude of the resistivity measurement to qualitative clay content. Rock exhibiting resistivity of less than 20 ??-m is considered potentially unstable, and resistivities less than 10 ??-m indicate appreciable amounts of clay are present in the rock. Partially saturated rock zones are defined by the measurement of the rock sound speed. Zones in the rock which exhibit velocities less than 2450 m/sec are considered of potential concern. ?? 1980.

International Journal of Rock Mechanics and Mining↗
USGS2001

Flow of variably fluidized granular masses across three-dimensional terrain I. Coulomb mixture theory

Rock avalanches, debris flows, and related phenomena consist of grain-fluid mixtures that move across three-dimensional terrain. In all these phenomena the same basic forces, govern motion, but differing mixture compositions, initial conditions, and boundary conditions yield varied dynamics and deposits. To predict motion of diverse grain-fluid masses from initiation to deposition, we develop a depth-averaged, threedimensional mathematical model that accounts explicitly for solid- and fluid-phase forces and interactions. Model input consists of initial conditions, path topography, basal and internal friction angles of solid grains, viscosity of pore fluid, mixture density, and a mixture diffusivity that controls pore pressure dissipation. Because these properties are constrained by independent measurements, the model requires little or no calibration and yields readily testable predictions. In the limit of vanishing Coulomb friction due to persistent high fluid pressure the model equations describe motion of viscous floods, and in the limit of vanishing fluid stress they describe one-phase granular avalanches. Analysis of intermediate phenomena such as debris flows and pyroclastic flows requires use of the full mixture equations, which can simulate interaction of high-friction surge fronts with more-fluid debris that follows. Special numerical methods (described in the companion paper) are necessary to solve the full equations, but exact analytical solutions of simplified equations provide critical insight. An analytical solution for translational motion of a Coulomb mixture accelerating from rest and descending a uniform slope demonstrates that steady flow can occur only asymptotically. A solution for the asymptotic limit of steady flow in a rectangular channel explains why shear may be concentrated in narrow marginal bands that border a plug of translating debris. Solutions for static equilibrium of source areas describe conditions of incipient slope instability, and other static solutions show that nonuniform distributions of pore fluid pressure produce bluntly tapered vertical profiles at the margins of deposits. Simplified equations and solutions may apply in additional situations identified by a scaling analysis. Assessment of dimensionless scaling parameters also reveals that miniature laboratory experiments poorly simulate the dynamics of full-scale flows in which fluid effects are significant. Therefore large geophysical flows can exhibit dynamics not evident at laboratory scales.

Journal of Geophysical Research B: Solid Earth↗
USGS2009

Snow and ice in a changing hydrological world

On the occasion of the 60th anniversary of the International Association of Hydrological Sciences, the 100th and 50th anniversaries of the First and Second International Polar Years, and the 25th anniversary of the International Geophysical Year, it seems appropriate to re-examine the world's water balance and the role of snow and ice in the global hydrological climatic system. Snow cover on land (especially in the Northern Hemisphere) and sea ice (especially in the Southern Hemisphere) vary seasonally, and this seasonal change has an important effect on the world climate because snow and sea ice reflect solar radiation efficiently and affect other heat flow processes between atmosphere and land or ocean. Glaciers, including ice sheets, store most of the fresh water on Earth, but change dimensions relatively slowly. There is no clear evidence that the glacier ice volume currently is declining, but more needs to be known about mountain glacier and ice sheet mass balances. The current rise in sea level poses an enigma: thermal expansion of the oceans may account for half of the present rise, but the other half is unexplained. Although major changes in the large ice sheets take place over time scales of 10 4 to 10 5 years, marine ice sheets may be subject to rapid disintegration due to grounding line instability, perhaps accompanied by surging. Ice cores may produce remarkably complete histories of air temperature, precipitation, fallout, and atmospheric composition. A recent core through the Greenland Ice Sheet shows an abrupt transition from glacial to modern climate just over 10000 years ago, suggesting that climate is an “almost intransitive” system. Because of the possibility of abrupt climate transitions and the uncertain stability of the West Antarctic Ice Sheet, future climatic variations are difficult to predict. The calculated heating of the atmosphere in the polar regions due to CO 2 increase is, therefore, of special interest.

Hydrological Sciences Journal↗
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↗
USGS2003

GCIP water and energy budget synthesis (WEBS)

As part of the World Climate Research Program's (WCRPs) Global Energy and Water-Cycle Experiment (GEWEX) Continental-scale International Project (GCIP), a preliminary water and energy budget synthesis (WEBS) was developed for the period 1996-1999 fromthe "best available" observations and models. Besides this summary paper, a companion CD-ROM with more extensive discussion, figures, tables, and raw data is available to the interested researcher from the GEWEX project office, the GAPP project office, or the first author. An updated online version of the CD-ROM is also available at http://ecpc.ucsd.edu/gcip/webs.htm/. Observations cannot adequately characterize or "close" budgets since too many fundamental processes are missing. Models that properly represent the many complicated atmospheric and near-surface interactions are also required. This preliminary synthesis therefore included a representative global general circulation model, regional climate model, and a macroscale hydrologic model as well as a global reanalysis and a regional analysis. By the qualitative agreement among the models and available observations, it did appear that we now qualitatively understand water and energy budgets of the Mississippi River Basin. However, there is still much quantitative uncertainty. In that regard, there did appear to be a clear advantage to using a regional analysis over a global analysis or a regional simulation over a global simulation to describe the Mississippi River Basin water and energy budgets. There also appeared to be some advantage to using a macroscale hydrologic model for at least the surface water budgets. Copyright 2003 by the American Geophysical Union.

Journal of Geophysical Research D: Atmospheres↗
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↗
USGS2020

SPEAR: The next generation GFDL modeling system for seasonal to multidecadal prediction and projection

We document the development and simulation characteristics of the next generation modeling system for seasonal to decadal prediction and projection at the Geophysical Fluid Dynamics Laboratory (GFDL). SPEAR ( S eamless System for P rediction and EA rth System R esearch) is built from component models recently developed at GFDL—the AM4 atmosphere model, MOM6 ocean code, LM4 land model, and SIS2 sea ice model. The SPEAR models are specifically designed with attributes needed for a prediction model for seasonal to decadal time scales, including the ability to run large ensembles of simulations with available computational resources. For computational speed SPEAR uses a coarse ocean resolution of approximately 1.0° (with tropical refinement). SPEAR can use differing atmospheric horizontal resolutions ranging from 1° to 0.25°. The higher atmospheric resolution facilitates improved simulation of regional climate and extremes. SPEAR is built from the same components as the GFDL CM4 and ESM4 models but with design choices geared toward seasonal to multidecadal physical climate prediction and projection. We document simulation characteristics for the time mean climate, aspects of internal variability, and the response to both idealized and realistic radiative forcing change. We describe in greater detail one focus of the model development process that was motivated by the importance of the Southern Ocean to the global climate system. We present sensitivity tests that document the influence of the Antarctic surface heat budget on Southern Ocean ventilation and deep global ocean circulation. These findings were also useful in the development processes for the GFDL CM4 and ESM4 models.

Journal of Advances in Modeling Earth Systems↗
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↗
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↗
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↗
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↗
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↗
USGS2014

Methane hydrates in nature - Current knowledge and challenges

Recognizing the importance of methane hydrate research and the need for a coordinated effort, the United States Congress enacted the Methane Hydrate Research and Development Act of 2000. At the same time, the Ministry of International Trade and Industry in Japan launched a research program to develop plans for a methane hydrate exploratory drilling project in the Nankai Trough. India, China, the Republic of Korea, and other nations also have established large methane hydrate research and development programs. Government-funded scientific research drilling expeditions and production test studies have provided a wealth of information on the occurrence of methane hydrates in nature. Numerous studies have shown that the amount of gas stored as methane hydrates in the world may exceed the volume of known organic carbon sources. However, methane hydrates represent both a scientific and technical challenge, and much remains to be learned about their characteristics and occurrence in nature. Methane hydrate research in recent years has mostly focused on: (1) documenting the geologic parameters that control the occurrence and stability of methane hydrates in nature, (2) assessing the volume of natural gas stored within various methane hydrate accumulations, (3) analyzing the production response and characteristics of methane hydrates, (4) identifying and predicting natural and induced environmental and climate impacts of natural methane hydrates, (5) analyzing the methane hydrate role as a geohazard, (6) establishing the means to detect and characterize methane hydrate accumulations using geologic and geophysical data, and (7) establishing the thermodynamic phase equilibrium properties of methane hydrates as a function of temperature, pressure, and gas composition. The U.S. Department of Energy (DOE) and the Consortium for Ocean Leadership (COL) combined their efforts in 2012 to assess the contributions that scientific drilling has made and could continue to make to advance our understanding of methane hydrates in nature. COL assembled a Methane Hydrate Project Science Team with members from academia, industry, and government. This Science Team worked with COL and DOE to develop and host the Methane Hydrate Community Workshop, which surveyed a substantial cross section of the methane hydrate research community for input on the most important research developments in our understanding of methane hydrates in nature and their potential role as an energy resource, a geohazard, and/or as an agent of global climate change. Our understanding of how methane hydrates occur in nature is still growing and evolving, and it is known with certainty that field, laboratory, and modeling studies have contributed greatly to our understanding of hydrates in nature and will continue to be a critical source of the information needed to advance our understanding of methane hydrates.

Journal of Chemical and Engineering Data↗
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↗
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↗
USGS2012

The mechanics and three-dimensional internal structure of active magmatic systems: Kilauea volcano, Hawaii

Interpretation of abundant seismic data suggests that Kilauea's primary conduit within the upper mantle is concentrically zoned to about 34-km depth. This zoned structure is inferred to contain a central core region of relatively higher permeability, surrounded by numerous dikes that are in intermittent hydraulic communication with each other and with the central core. During periods of relatively high magma transport, the entire cross section of the conduit is utilized. During periods of relatively low to moderate transport, however, only the central core is active. As the conduit penetrates the oceanic crust and enters the volcanic shield, it simultaneously supplies the deeper sections of the rift zones (6-to 10-km depth) and the roots of the summit reservoir with picritic magma. The rift zones at depth are inferred to be almost wholly molten and to possess a high degree of fluid continuity from Heiheiahulu in the East Rift Zone, 45 km westward through the roots of the summit magma reservoir, and well into the Southwest Rift Zone. Higher in the shield, the subcaldera magma reservoir and the shallow rift zones occupy the 2-to 4-km depth interval. Summit-differentiated olivine tholeiite (ρ ≈ 2.62 g cm −3 ) is periodically injected laterally along a horizon of neutral buoyancy within the rift zones, where the density of the magma is just balanced by the in situ density of the shield (Ryan, 1987 a , b ). Deep rift zone intrusions push seaward the deep tectonic blocks of the volcano's south flank. Shallow rift intrusions build a sheeted dike complex, inferred to be in isostatic equilibrium with the higher-density deep rift cores below. General finite element analyses are presented for the deformation and stress fields surrounding such dikes in the horizontal and vertical planes. The dike tip in two and three dimensions is surrounded by a tubular core of tensile (σ 1 , σ 2 ) and shear stress (τ max ). The displacement field is characterized by counterrotating cells on either side of the dike tip which, in vertical orientation, produce the characteristic subsidence above the dike complex, with uplift on either side, forming a ridge-trough-ridge structure. A finite element model of Kilauea's shield computes the displacement fields and principal stress (σ 1 ) distributions resulting from intrusive activity on each or both of the rift zones. Within the summit region, tensile stress lobes produced by the three-dimensional upward extension of the intrusions superpose constructively to produce calderawide regimes of tensile stress, conducive to caldera development. Parametric studies of (1) intrusion in the East Rift Zone only, (2) intrusion in the Southwest Rift Zone only, and (3) intrusion in both rift zones demonstrate their unique kinematic contributions. For case 1, the caldera undergoes a counterclockwise rotation (torque up state) conducive to the development of rightstepping en echelon eruptive fissures, as exemplified by the August 14, 1971, eruption. For case 2, the caldera undergoes a clockwise rotation (torque down state) conducive to the development of left-stepping eruptive fissures, as occurred during the December 31, 1974, eruption. For case 3, the caldera substructure is driven due southward, producing the southward migration of the upper portions of the summit magma reservoir.

Journal of Geophysical Research Solid Earth↗