USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Geophysics Research”

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 1,027 records · Page 57Linked to original sources
USGS2012

Slab1.0: A three-dimensional model of global subduction zone geometries

We describe and present a new model of global subduction zone geometries, called Slab1.0. An extension of previous efforts to constrain the two-dimensional non-planar geometry of subduction zones around the focus of large earthquakes, Slab1.0 describes the detailed, non-planar, three-dimensional geometry of approximately 85% of subduction zones worldwide. While the model focuses on the detailed form of each slab from their trenches through the seismogenic zone, where it combines data sets from active source and passive seismology, it also continues to the limits of their seismic extent in the upper-mid mantle, providing a uniform approach to the definition of the entire seismically active slab geometry. Examples are shown for two well-constrained global locations; models for many other regions are available and can be freely downloaded in several formats from our new Slab1.0 website, http://on.doi.gov/d9ARbS . We describe improvements in our two-dimensional geometry constraint inversion, including the use of ‘average’ active source seismic data profiles in the shallow trench regions where data are otherwise lacking, derived from the interpolation between other active source seismic data along-strike in the same subduction zone. We include several analyses of the uncertainty and robustness of our three-dimensional interpolation methods. In addition, we use the filtered, subduction-related earthquake data sets compiled to build Slab1.0 in a reassessment of previous analyses of the deep limit of the thrust interface seismogenic zone for all subduction zones included in our global model thus far, concluding that the width of these seismogenic zones is on average 30% larger than previous studies have suggested.

Journal of Geophysical Research B: Solid Earth↗
USGS2011

Shallow conduit system at Kilauea Volcano, Hawaii, revealed by seismic signals associated with degassing bursts

Eruptive activity at the summit of Kilauea Volcano, Hawaii, beginning in March, 2008 and continuing to the present time is characterized by episodic explosive bursts of gas and ash from a vent within Halemaumau Pit Crater. These bursts are accompanied by seismic signals that are well recorded by a broadband network deployed in the summit caldera. We investigate in detail the dimensions and oscillation modes of the source of a representative burst in the 1−10 s band. An extended source is realized by a set of point sources distributed on a grid surrounding the source centroid, where the centroid position and source geometry are fixed from previous modeling of very-long-period (VLP) data in the 10–50 s band. The source time histories of all point sources are obtained simultaneously through waveform inversion carried out in the frequency domain. Short-scale noisy fluctuations of the source time histories between adjacent sources are suppressed with a smoothing constraint, whose strength is determined through a minimization of the Akaike Bayesian Information Criterion (ABIC). Waveform inversions carried out for homogeneous and heterogeneous velocity structures both image a dominant source component in the form of an east trending dike with dimensions of 2.9 × 2.9 km. The dike extends ∼2 km west and ∼0.9 km east of the VLP centroid and spans the depth range 0.2–3.1 km. The source model for a homogeneous velocity structure suggests the dike is hinged at the source centroid where it bends from a strike E 27°N with northern dip of 85° west of the centroid, to a strike E 7°N with northern dip of 80° east of the centroid. The oscillating behavior of the dike is dominated by simple harmonic modes with frequencies ∼0.2 Hz and ∼0.5 Hz, representing the fundamental mode ν 11 and first degenerate mode ν 12 = ν 21 of the dike. Although not strongly supported by data in the 1–10 s band, a north striking dike segment is required for enhanced compatibility with the model elaborated in the 10–50 s band. This dike provides connectivity between the east trending dike and the new vent within Halemaumau Pit Crater. Waveform inversions with a dual-dike model suggest dimensions of 0.7 × 0.7 km to 2.6 × 2.6 km for this segment. Further elaboration of the complex dike system under Halemaumau does not appear to be feasible with presently available data.

Hawai'i↗
USGS2012

Geometry and subsidence history of the Dead Sea basin: A case for fluid-induced mid-crustal shear zone?

Pull‐apart basins are narrow zones of crustal extension bounded by strike‐slip faults that can serve as analogs to the early stages of crustal rifting. We use seismic tomography, 2‐D ray tracing, gravity modeling, and subsidence analysis to study crustal extension of the Dead Sea basin (DSB), a large and long‐lived pull‐apart basin along the Dead Sea transform (DST). The basin gradually shallows southward for 50 km from the only significant transverse normal fault. Stratigraphic relationships there indicate basin elongation with time. The basin is deepest (8–8.5 km) and widest (∼15 km) under the Lisan about 40 km north of the transverse fault. Farther north, basin depth is ambiguous, but is 3 km deep immediately north of the lake. The underlying pre‐basin sedimentary layer thickens gradually from 2 to 3 km under the southern edge of the DSB to 3–4 km under the northern end of the lake and 5–6 km farther north. Crystalline basement is ∼11 km deep under the deepest part of the basin. The upper crust under the basin has lower P wave velocity than in the surrounding regions, which is interpreted to reflect elevated pore fluids there. Within data resolution, the lower crust below ∼18 km and the Moho are not affected by basin development. The subsidence rate was several hundreds of m/m.y. since the development of the DST ∼17 Ma, similar to other basins along the DST, but subsidence rate has accelerated by an order of magnitude during the Pleistocene, which allowed the accumulation of 4 km of sediment. We propose that the rapid subsidence and perhaps elongation of the DSB are due to the development of inter‐connected mid‐crustal ductile shear zones caused by alteration of feldspar to muscovite in the presence of pore fluids. This alteration resulted in a significant strength decrease and viscous creep. We propose a similar cause to the enigmatic rapid subsidence of the North Sea at the onset the North Atlantic mantle plume. Thus, we propose that aqueous fluid flux into a slowly extending continental crust can cause rapid basin subsidence that may be erroneously interpreted as an increased rate of tectonic activity.

Dead Sea shoreline↗
USGS2012

Fixed recurrence and slip models better predict earthquake behavior than the time- and slip-predictable models 1: repeating earthquakes

The behavior of individual events in repeating earthquake sequences in California, Taiwan and Japan is better predicted by a model with fixed inter-event time or fixed slip than it is by the time- and slip-predictable models for earthquake occurrence. Given that repeating earthquakes are highly regular in both inter-event time and seismic moment, the time- and slip-predictable models seem ideally suited to explain their behavior. Taken together with evidence from the companion manuscript that shows similar results for laboratory experiments we conclude that the short-term predictions of the time- and slip-predictable models should be rejected in favor of earthquake models that assume either fixed slip or fixed recurrence interval. This implies that the elastic rebound model underlying the time- and slip-predictable models offers no additional value in describing earthquake behavior in an event-to-event sense, but its value in a long-term sense cannot be determined. These models likely fail because they rely on assumptions that oversimplify the earthquake cycle. We note that the time and slip of these events is predicted quite well by fixed slip and fixed recurrence models, so in some sense they are time- and slip-predictable. While fixed recurrence and slip models better predict repeating earthquake behavior than the time- and slip-predictable models, we observe a correlation between slip and the preceding recurrence time for many repeating earthquake sequences in Parkfield, California. This correlation is not found in other regions, and the sequences with the correlative slip-predictable behavior are not distinguishable from nearby earthquake sequences that do not exhibit this behavior.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Holocene faulting in the Bellingham forearc basin: Upper-plate deformation at the northern end of the Cascadia subduction zone

The northern Cascadia forearc takes up most of the strain transmitted northward via the Oregon Coast block from the northward-migrating Sierra Nevada block. The north-south contractional strain in the forearc manifests in upper-plate faults active during the Holocene, the northern-most components of which are faults within the Bellingham Basin. The Bellingham Basin is the northern of four basins of the actively deforming northern Cascadia forearc. A set of Holocene faults, Drayton Harbor, Birch Bay, and Sandy Point faults, occur within the Bellingham Basin and can be traced from onshore to offshore using a combination of aeromagnetic lineaments, paleoseismic investigations and scarps identified using LiDAR imagery. With the recognition of such Holocene faults, the northernmost margin of the actively deforming Cascadia forearc extends 60 km north of the previously recognized limit of Holocene forearc deformation. Although to date no Holocene faults are recognized at the northern boundary of the Bellingham Basin, which is 15 km north of the international border, there is no compelling tectonic reason to expect that Holocene faults are limited to south of the international border.

Washington↗
USGS2012

Impact of reduced near-field entrainment of overpressured volcanic jets on plume development

Volcanic plumes are often studied using one-dimensional analytical models, which use an empirical entrainment ratio to close the equations. Although this ratio is typically treated as constant, its value near the vent is significantly reduced due to flow development and overpressured conditions. To improve the accuracy of these models, a series of experiments was performed using particle image velocimetry, a high-accuracy, full-field velocity measurement technique. Experiments considered a high-speed jet with Reynolds numbers up to 467,000 and exit pressures up to 2.93 times atmospheric. Exit gas densities were also varied from 0.18 to 1.4 times that of air. The measured velocity was integrated to determine entrainment directly. For jets with exit pressures near atmospheric, entrainment was approximately 30% less than the fully developed level at 20 diameters from the exit. At pressures nearly three times that of the atmosphere, entrainment was 60% less. These results were introduced into Plumeria, a one-dimensional plume model, to examine the impact of reduced entrainment. The maximum column height was only slightly modified, but the critical radius for collapse was significantly reduced, decreasing by nearly a factor of two at moderate eruptive pressures.

Journal of Geophysical Research↗
USGS2012

Toward an understanding of disequilibrium dihedral angles in mafic rocks

The median dihedral angle at clinopyroxene-plagioclase-plagioclase junctions in mafic rocks, Θcpp, is generally lower than equilibrium (109˚ {plus minus} 2˚). Observation of a wide range of mafic bodies demonstrates that previous work on systematic variations of Θcpp is incorrect in several important respects. Firstly, the spatial distribution of plagioclase compositional zoning demonstrates that the final geometry of three-grain junctions, and hence Θcpp, is formed during solidification (the igneous process): sub-solidus textural modification in most dolerites and gabbros, previously thought to be the dominant control on Θcpp, is insignificant. Θcpp is governed by mass transport constraints, the inhibiting effects of small pore size on crystallization, and variation in relative growth rates of pyroxene and plagioclase. During rapid cooling, pyroxene preferentially fills wider pores while the narrower pores remain melt-filled, resulting in an initial value of Θcpp of 78˚, rather than 60˚ which would be expected if all melt-filled pores were filled with pyroxene. Lower cooling rates create a higher initial Θcpp due to changes in relative growth rates of the two minerals at the nascent three-grain junction. Low Θcpp (associated with cuspate clinopyroxene grains at triple junctions) can also be diagnostic of infiltration of previously melt-free rocks by late-stage evolved liquids (the metasomatic process). Modification of Θcpp by sub-solidus textural equilibration (the metamorphic process) is only important for fine-grained mafic rocks such as chilled margins and intra-plutonic chill zones. In coarse-grained gabbros from shallow crustal intrusions the metamorphic process occurs only in the centres of oikocrysts, associated with rounding of chadacrysts.

Journal of Geophysical Research↗
USGS2012

Why the 2002 Denali fault rupture propagated onto the Totschunda fault: implications for fault branching and seismic hazards

The propagation of the rupture of the M w 7.9 Denali fault earthquake from the central Denali fault onto the Totschunda fault has provided a basis for dynamic models of fault branching in which the angle of the regional or local prestress relative to the orientation of the main fault and branch plays a principal role in determining which fault branch is taken. GeoEarthScope LiDAR and paleoseismic data allow us to map the structure of the Denali-Totschunda fault intersection and evaluate controls of fault branching from a geological perspective. LiDAR data reveal the Denali-Totschunda fault intersection is structurally simple with the two faults directly connected. At the branch point, 227.2 km east of the 2002 epicenter, the 2002 rupture diverges southeast to become the Totschunda fault. We use paleoseismic data to propose that differences in the accumulated strain on each fault segment, which express differences in the elapsed time since the most recent event, was one important control of the branching direction. We suggest that data on event history, slip rate, paleo offsets, fault geometry and structure, and connectivity, especially on high slip rate-short recurrence interval faults, can be used to assess the likelihood of branching and its direction. Analysis of the Denali-Totschunda fault intersection has implications for evaluating the potential for a rupture to propagate across other types of fault intersections and for characterizing sources of future large earthquakes.

Alaska↗
USGS2012

Ash3d: A finite-volume, conservative numerical model for ash transport and tephra deposition

We develop a transient, 3-D Eulerian model (Ash3d) to predict airborne volcanic ash concentration and tephra deposition during volcanic eruptions. This model simulates downwind advection, turbulent diffusion, and settling of ash injected into the atmosphere by a volcanic eruption column. Ash advection is calculated using time-varying pre-existing wind data and a robust, high-order, finite-volume method. Our routine is mass-conservative and uses the coordinate system of the wind data, either a Cartesian system local to the volcano or a global spherical system for the Earth. Volcanic ash is specified with an arbitrary number of grain sizes, which affects the fall velocity, distribution and duration of transport. Above the source volcano, the vertical mass distribution with elevation is calculated using a Suzuki distribution for a given plume height, eruptive volume, and eruption duration. Multiple eruptions separated in time may be included in a single simulation. We test the model using analytical solutions for transport. Comparisons of the predicted and observed ash distributions for the 18 August 1992 eruption of Mt. Spurr in Alaska demonstrate to the efficacy and efficiency of the routine.

Journal of Geophysical Research↗
USGS2012

Tidal triggering of low frequency earthquakes near Parkfield, California: Implications for fault mechanics within the brittle-ductile transition

Studies of nonvolcanic tremor (NVT) have established the significant impact of small stress perturbations on NVT generation. Here we analyze the influence of the solid earth and ocean tides on a catalog of ∼550,000 low frequency earthquakes (LFEs) distributed along a 150 km section of the San Andreas Fault centered at Parkfield. LFE families are identified in the NVT data on the basis of waveform similarity and are thought to represent small, effectively co-located earthquakes occurring on brittle asperities on an otherwise aseismic fault at depths of 16 to 30 km. We calculate the sensitivity of each of these 88 LFE families to the tidally induced right-lateral shear stress (RLSS), fault-normal stress (FNS), and their time derivatives and use the hypocentral locations of each family to map the spatial variability of this sensitivity. LFE occurrence is most strongly modulated by fluctuations in shear stress, with the majority of families demonstrating a correlation with RLSS at the 99% confidence level or above. Producing the observed LFE rate modulation in response to shear stress perturbations requires low effective stress in the LFE source region. There are substantial lateral and vertical variations in tidal shear stress sensitivity, which we interpret to reflect spatial variation in source region properties, such as friction and pore fluid pressure. Additionally, we find that highly episodic, shallow LFE families are generally less correlated with tidal stresses than their deeper, continuously active counterparts. The majority of families have weaker or insignificant correlation with positive (tensile) FNS. Two groups of families demonstrate a stronger correlation with fault-normal tension to the north and with compression to the south of Parkfield. The families that correlate with fault-normal clamping coincide with a releasing right bend in the surface fault trace and the LFE locations, suggesting that the San Andreas remains localized and contiguous down to near the base of the crust. The deep families that have high sensitivity to both shear and tensile normal stress perturbations may be indicative of an increase in effective fault contact area with depth. Synthesizing our observations with those of other LFE-hosting localities will help to develop a comprehensive understanding of transient fault slip below the “seismogenic zone” by providing constraints on parameters in physical models of slow slip and LFEs.

California↗
USGS2011

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

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

Massachusetts↗
USGS2012

Utilizing multichannel electrical resistivity methods to examine the dynamics of the fresh water–seawater interface in two Hawaiian groundwater systems

Multichannel electrical resistivity (ER) measurements were conducted at two contrasting coastal sites in Hawaii to obtain new information on the spatial scales and dynamics of the fresh water–seawater interface and rates of coastal groundwater exchange. At Kiholo Bay (located on the dry, Kona side of the Big Island) and at a site in Maunalua Bay (Oahu), there is an evidence for abundant submarine groundwater discharge (SGD). However, the hydrologic and geologic controls on coastal groundwater discharge are likely to be different at these two sites. While at Kiholo Bay SGD is predominantly through lava tubes, at the Maunalua Bay site exchange occurs mostly through nearshore submarine springs. In order to calculate SGD fluxes, it is important to understand the spatial and temporal scales of coastal groundwater exchange. From ER time series data, subsurface salinity distributions were calculated using site-specific formation factors. A salinity mass balance box model was then used to calculate rates of point source (i.e., spatially discreet) and total fresh water discharge. From these data, mean SGD rates were calculated for Kiholo Bay (∼9,200 m 3 /d) and for the Maunalua Bay site (∼5,900 m 3 /d). While such results are on the same order of magnitude to geochemical tracer-derived SGD rates, the ER SGD rates provide enhanced details of coastal groundwater exchange that can enable a more cohesive whole watershed perspective.

Hawai'i↗
USGS2012

Use of vertical temperature gradients for prediction of tidal flat sediment characteristics

Sediment characteristics largely govern tidal flat morphologic evolution; however, conventional methods of investigating spatial variability in lithology on tidal flats are difficult to employ in these highly dynamic regions. In response, a series of laboratory experiments was designed to investigate the use of temperature diffusion toward sediment characterization. A vertical thermistor array was used to quantify temperature gradients in simulated tidal flat sediments of varying compositions. Thermal conductivity estimates derived from these arrays were similar to measurements from a standard heated needle probe, which substantiates the thermistor methodology. While the thermal diffusivities of dry homogeneous sediments were similar, diffusivities for saturated homogeneous sediments ranged approximately one order of magnitude. The thermal diffusivity of saturated sand was five times the thermal diffusivity of saturated kaolin and more than eight times the thermal diffusivity of saturated bentonite. This suggests that vertical temperature gradients can be used for distinguishing homogeneous saturated sands from homogeneous saturated clays and perhaps even between homogeneous saturated clay types. However, experiments with more realistic tidal flat mixtures were less discriminating. Relationships between thermal diffusivity and percent fines for saturated mixtures varied depending upon clay composition, indicating that clay hydration and/or water content controls thermal gradients. Furthermore, existing models for the bulk conductivity of sediment mixtures were improved only through the use of calibrated estimates of homogeneous end-member conductivity and water content values. Our findings suggest that remotely sensed observations of water content and thermal diffusivity could only be used to qualitatively estimate tidal flat sediment characteristics.

Journal of Geophysical Research C: Oceans↗
USGS2011

The Moon Mineralogy Mapper (M3) imaging spectrometer for lunar science: Instrument description, calibration, on‐orbit measurements, science data calibration and on‐orbit validation

The NASA Discovery Moon Mineralogy Mapper imaging spectrometer was selected to pursue a wide range of science objectives requiring measurement of composition at fine spatial scales over the full lunar surface. To pursue these objectives, a broad spectral range imaging spectrometer with high uniformity and high signal-to-noise ratio capable of measuring compositionally diagnostic spectral absorption features from a wide variety of known and possible lunar materials was required. For this purpose the Moon Mineralogy Mapper imaging spectrometer was designed and developed that measures the spectral range from 430 to 3000 nm with 10 nm spectral sampling through a 24 degree field of view with 0.7 milliradian spatial sampling. The instrument has a signal-to-noise ratio of greater than 400 for the specified equatorial reference radiance and greater than 100 for the polar reference radiance. The spectral cross-track uniformity is >90% and spectral instantaneous field-of-view uniformity is >90%. The Moon Mineralogy Mapper was launched on Chandrayaan-1 on the 22nd of October. On the 18th of November 2008 the Moon Mineralogy Mapper was turned on and collected a first light data set within 24 h. During this early checkout period and throughout the mission the spacecraft thermal environment and orbital parameters varied more than expected and placed operational and data quality constraints on the measurements. On the 29th of August 2009, spacecraft communication was lost. Over the course of the flight mission 1542 downlinked data sets were acquired that provide coverage of more than 95% of the lunar surface. An end-to-end science data calibration system was developed and all measurements have been passed through this system and delivered to the Planetary Data System (PDS.NASA.GOV). An extensive effort has been undertaken by the science team to validate the Moon Mineralogy Mapper science measurements in the context of the mission objectives. A focused spectral, radiometric, spatial, and uniformity validation effort has been pursued with selected data sets including an Earth-view data set. With this effort an initial validation of the on-orbit performance of the imaging spectrometer has been achieved, including validation of the cross-track spectral uniformity and spectral instantaneous field of view uniformity. The Moon Mineralogy Mapper is the first imaging spectrometer to measure a data set of this kind at the Moon. These calibrated science measurements are being used to address the full set of science goals and objectives for this mission.

Journal of Geophysical Research E: Planets↗
USGS2011

Stability of Mg-sulfates at-10C and the rates of dehydration/rehydration processes under conditions relevant to Mars

We report the results of low temperature (−10°C) experiments on the stability fields and phase transition pathways of five hydrous Mg-sulfates. A low temperature form of MgSO 4 ·7H 2 O (LT-7w) was found to have a wide stability field that extends to low relative humidity (∼13% RH at −10°C). Using information on the timing of phase transitions, we extracted information on the reaction rates of five important dehydration and rehydration processes. We found that the temperature dependencies of rate constants for dehydration processes differ from those of rehydration, which reflect differences in reaction mechanisms. By extrapolating these rate constants versus T correlations into the T range relevant to Mars, we can evaluate the possibility of occurrence of specific processes and the presence of common Mg-sulfate species present on Mars in different periods and locations. We anticipate in a moderate obliquity period, starkeyite and LH-MgSO 4 ·H 2 O should be two common Mg-sulfates at the surface, another polymorph MH-MgSO 4 ·H 2 O can exist at the locations where hydrothermal processes may have occurred. In polar regions or within the subsurface of other regions, meridianiite (coexisting with water ice, near 100% RH) and LT-7w (over a large RH range) are the stable phases. During a high obliquity period, meridianiite and LT-7w should exhibit widespread occurrence. The correlations of reaction rates versus temperature found in this study imply that dehydration and rehydration of hydrous Mg-sulfates would always be slower than the sublimation and crystallization of water ice, which would be supported by mission observations from Odyssey and by Mars Exploration Rovers.

Journal of Geophysical Research E: Planets↗
USGS2012

Characterization of previously unidentified lunar pyroclastic deposits using Lunar Reconnaissance Orbiter Camera (LROC) data

We used a Lunar Reconnaissance Orbiter Camera (LROC) global monochrome Wide-angle Camera (WAC) mosaic to conduct a survey of the Moon to search for previously unidentified pyroclastic deposits. Promising locations were examined in detail using LROC multispectral WAC mosaics, high-resolution LROC Narrow Angle Camera (NAC) images, and Clementine multispectral (ultraviolet-visible or UVVIS) data. Out of 47 potential deposits chosen for closer examination, 12 were selected as probable newly identified pyroclastic deposits. Potential pyroclastic deposits were generally found in settings similar to previously identified deposits, including areas within or near mare deposits adjacent to highlands, within floor-fractured craters, and along fissures in mare deposits. However, a significant new finding is the discovery of localized pyroclastic deposits within floor-fractured craters Anderson E and F on the lunar farside, isolated from other known similar deposits. Our search confirms that most major regional and localized low-albedo pyroclastic deposits have been identified on the Moon down to ~100 m/pix resolution, and that additional newly identified deposits are likely to be either isolated small deposits or additional portions of discontinuous, patchy deposits.

Journal of Geophysical Research E: Planets↗
USGS2011

The Sculptured Hills of the Taurus Highlands: Implications for the relative age of Serenitatis, basin chronologies and the cratering history of the Moon

New images from the Lunar Reconnaissance Orbiter Camera show the distribution and geological relations of the Sculptured Hills, a geological unit widespread in the highlands between the Serenitatis and Crisium basins. The Sculptured Hills shows knobby, undulating, radially textured, and plains-like morphologies and in many places is indistinguishable from the similarly knobby Alpes Formation, a facies of ejecta from the Imbrium basin. The new LROC image data show that the Sculptured Hills in the Taurus highlands is Imbrium ejecta and not directly related to the formation of the Serenitatis basin. This occurrence and the geological relations of this unit suggests that the Apollo 17 impact melts may not be not samples of the Serenitatis basin-forming impact, leaving their provenance undetermined and origin unexplained. If the Apollo 17 melt rocks are Serenitatis impact melt, up to half of the basin and large crater population of the Moon was created within a 30 Ma interval around 3.8 Ga in a global impact "cataclysm." Either interpretation significantly changes our view of the impact process and history of the Earth-Moon system. Copyright 2011 by the American Geophysical Union.

Journal of Geophysical Research E: Planets↗
USGS2012

Advances in the simulation and automated measurement of well-sorted granular material: 1. Simulation

1. In this, the first of a pair of papers which address the simulation and automated measurement of well-sorted natural granular material, a method is presented for simulation of two-phase (solid, void) assemblages of discrete non-cohesive particles. The purpose is to have a flexible, yet computationally and theoretically simple, suite of tools with well constrained and well known statistical properties, in order to simulate realistic granular material as a discrete element model with realistic size and shape distributions, for a variety of purposes. The stochastic modeling framework is based on three-dimensional tessellations with variable degrees of order in particle-packing arrangement. Examples of sediments with a variety of particle size distributions and spatial variability in grain size are presented. The relationship between particle shape and porosity conforms to published data. The immediate application is testing new algorithms for automated measurements of particle properties (mean and standard deviation of particle sizes, and apparent porosity) from images of natural sediment, as detailed in the second of this pair of papers. The model could also prove useful for simulating specific depositional structures found in natural sediments, the result of physical alterations to packing and grain fabric, using discrete particle flow models. While the principal focus here is on naturally occurring sediment and sedimentary rock, the methods presented might also be useful for simulations of similar granular or cellular material encountered in engineering, industrial and life sciences.

Journal of Geophysical Research F: Earth Surface↗