Search USGS⌕ Search

SEARCH · Search USGS

Results for “Journal of Geophysical 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 829 records · Page 46Linked to original sources
USGS2011

InSAR observations of aseismic slip associated with an earthquake swarm in the Columbia River flood basalts

In 2009 a swarm of small shallow earthquakes occurred within the basalt flows of the Columbia River Basalt Group (CRBG). The swarm occurred within a dense seismic network in the U.S. Department of Energys Hanford Site. Data from the seismic network along with interferometric synthetic aperture radar (InSAR) data from the European Space Agencys (ESA) ENVISAT satellite provide insight into the nature of the swarm. By modeling the InSAR deformation data we constructed a model that consists of a shallow thrust fault and a near horizontal fault. We suggest that the near horizontal lying fault is a bedding-plane fault located between basalt flows. The geodetic moment of the modeled fault system is about eight times the cumulative seismic moment of the swarm. Precise location estimates of the swarm earthquakes indicate that the area of highest slip on the thrust fault, ???70mm of slip less than ???0.5km depth, was not located within the swarm cluster. Most of the slip on the faults appears to have progressed aseismically and we suggest that interbed sediments play a central role in the slip process. Copyright 2011 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2011

Estimates of stress drop and crustal tectonic stress from the 27 February 2010 Maule, Chile, earthquake: Implications for fault strength

The great 27 February 2010 M w 8.8 earthquake off the coast of southern Chile ruptured a ∼600 km length of subduction zone. In this paper, we make two independent estimates of shear stress in the crust in the region of the Chile earthquake. First, we use a coseismic slip model constrained by geodetic observations from interferometric synthetic aperture radar (InSAR) and GPS to derive a spatially variable estimate of the change in static shear stress along the ruptured fault. Second, we use a static force balance model to constrain the crustal shear stress required to simultaneously support observed fore‐arc topography and the stress orientation indicated by the earthquake focal mechanism. This includes the derivation of a semianalytic solution for the stress field exerted by surface and Moho topography loading the crust. We find that the deviatoric stress exerted by topography is minimized in the limit when the crust is considered an incompressible elastic solid, with a Poisson ratio of 0.5, and is independent of Young's modulus. This places a strict lower bound on the critical stress state maintained by the crust supporting plastically deformed accretionary wedge topography. We estimate the coseismic shear stress change from the Maule event ranged from −6 MPa (stress increase) to 17 MPa (stress drop), with a maximum depth‐averaged crustal shear‐stress drop of 4 MPa. We separately estimate that the plate‐driving forces acting in the region, regardless of their exact mechanism, must contribute at least 27 MPa trench‐perpendicular compression and 15 MPa trench‐parallel compression. This corresponds to a depth‐averaged shear stress of at least 7 MPa. The comparable magnitude of these two independent shear stress estimates is consistent with the interpretation that the section of the megathrust fault ruptured in the Maule earthquake is weak, with the seismic cycle relieving much of the total sustained shear stress in the crust.

Maule↗
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↗
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

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↗
USGS2011

Wave-current interaction in Willapa Bay

This paper describes the importance of wave-current interaction in an inlet-estuary system. The three-dimensional, fully coupled, Coupled Ocean-Atmosphere-Wave-Sediment Transport (COAWST) modeling system was applied in Willapa Bay (Washington State) from 22 to 29 October 1998 that included a large storm event. To represent the interaction between waves and currents, the vortex-force method was used. Model results were compared with water elevations, currents, and wave measurements obtained by the U.S. Army Corp of Engineers. In general, a good agreement between field data and computed results was achieved, although some discrepancies were also observed in regard to wave peak directions in the most upstream station. Several numerical experiments that considered different forcing terms were run in order to identify the effects of each wind, tide, and wave-current interaction process. Comparison of the horizontal momentum balances results identified that wave-breaking-induced acceleration is one of the leading terms in the inlet area. The enhancement of the apparent bed roughness caused by waves also affected the values and distribution of the bottom shear stress. The pressure gradient showed significant changes with respect to the pure tidal case. During storm conditions the momentum balance in the inlet shares the characteristics of tidal-dominated and wave-dominated surf zone environments. The changes in the momentum balance caused by waves were manifested both in water level and current variations. The most relevant effect on hydrodynamics was a wave-induced setup in the inner part of the estuary.

Washington↗
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↗
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↗