USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Geophysical Research: Earth Surface”

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 343 records · Page 19Linked to original sources
USGS2018

Combining multiphase groundwater flow and slope stability models to assess stratovolcano flank collapse in the Cascade Range

Hydrothermal alteration can create low‐permeability zones, potentially resulting in elevated pore‐fluid pressures, within a volcanic edifice. Strength reduction by rock alteration and high pore‐fluid pressures have been suggested as a mechanism for edifice flank instability. Here we combine numerical models of multiphase heat transport and groundwater flow with a slope‐stability code that incorporates three‐dimensional distributions of strength and pore‐water pressure to address the following questions: (1) What permeability distributions and contrasts produce elevated pore‐fluid pressures in a stratovolcano? (2) What are the effects of these elevated pressures on flank stability? (3) Finally, what are the effects of magma intrusion on potential flank failure in an edifice? Simulation results show that under a range of plausible parameters, water tables in a stratovolcano can be elevated or perched. These elevated water tables result in universally lower stability (lower factor of safety) compared with equivalent dry edifices, indicating a higher likelihood of flank collapse. Low‐permeability (<1 × 10 −17 m 2 ) layers such as altered pyroclastic deposits or breccias can result in locally saturated regions (perched water) and lower factors of safety near the ground surface but may actually reduce liquid water saturation and pore pressures in the core of the edifice and thus may favor small, shallow collapses over larger, deeper collapses. Magma intrusion into the base of the edifice increases pore‐fluid pressures and decreases the factor of safety. However, the shear strength of edifice rocks also exerts a significant control on stability, so both mechanical properties and pore‐fluid pressures are important for stability assessments.

Washington↗
USGS2001

Viscoelastic shear zone model of a strike-slip earthquake cycle

I examine the behavior of a two-dimensional (2-D) strike-slip fault system embedded in a 1-D elastic layer (schizosphere) overlying a uniform viscoelastic half-space (plastosphere) and within the boundaries of a finite width shear zone. The viscoelastic coupling model of Savage and Prescott [1978] considers the viscoelastic response of this system, in the absence of the shear zone boundaries, to an earthquake occurring within the upper elastic layer, steady slip beneath a prescribed depth, and the superposition of the responses of multiple earthquakes with characteristic slip occurring at regular intervals. So formulated, the viscoelastic coupling model predicts that sufficiently long after initiation of the system, (1) average fault-parallel velocity at any point is the average slip rate of that side of the fault and (2) far-field velocities equal the same constant rate. Because of the sensitivity to the mechanical properties of the schizosphere-plastosphere system (i.e., elastic layer thickness, plastosphere viscosity), this model has been used to infer such properties from measurements of interseismic velocity. Such inferences exploit the predicted behavior at a known time within the earthquake cycle. By modifying the viscoelastic coupling model to satisfy the additional constraint that the absolute velocity at prescribed shear zone boundaries is constant, I find that even though the time-averaged behavior remains the same, the spatiotemporal pattern of surface deformation (particularly its temporal variation within an earthquake cycle) is markedly different from that predicted by the conventional viscoelastic coupling model. These differences are magnified as plastosphere viscosity is reduced or as the recurrence interval of periodic earthquakes is lengthened. Application to the interseismic velocity field along the Mojave section of the San Andreas fault suggests that the region behaves mechanically like a ???600-km-wide shear zone accommodating 50 mm/yr fault-parallel motion distributed between the San Andreas fault system and Eastern California Shear Zone. Copyright 2001 by the American Geophysical Union.

California↗
USGS2018

Interseismic ground deformation and fault slip rates in the greater San Francisco Bay Area from two decades of space geodetic data

The detailed spatial variations of strain accumulation and creep on major faults in the northern San Francisco Bay Area (North Bay), which are important for seismic potential and evaluation of natural hazards, remain poorly understood. Here we combine interferometric synthetic aperture radar data from the ERS‐1/2 and Envisat satellites between 1992 and 2010 with continuous and campaign GPS data to obtain a high spatial and temporal coverage of ground deformation of the North Bay. The SAR data from both ascending and descending orbits are combined to separate horizontal and vertical components of the deformation. We jointly invert the horizontal component of the mean velocities derived from these data to infer the deep strike‐slip rates on major locked faults. We use the estimated deep rates to simulate the long‐wavelength deformation due to interseismic elastic strain accumulation along these locked faults. After removing the long‐wavelength signal from the InSAR horizontal mean velocity field, we estimate fault‐parallel surface creep rates of up to 2 mm/year along the central section of the Rodgers Creek fault and surface creep rates ranging between 2 and 4 mm/year along the Concord fault. No surface creep is geodetically resolved along the West Napa and Green Valley fault zones. We identified characteristically repeating earthquakes on the Rodgers Creek fault, the West Napa fault, the Green Valley fault, and the Concord fault. Nontectonic deformation in the Geysers geothermal field and in Late Cenozoic basins (Rohnert Park and Sonoma basins) are also observed, likely due to hydrological and sediment‐compaction processes, respectively.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

The 1979 Homestead Valley earthquake sequence, California: Control of aftershocks and postseismic deformation

The coseismic slip and geometry of the March 15, 1979, Homestead Valley, California, earthquake sequence are well constrained by precise horizontal and vertical geodetic observations and by data from a dense local seismic network. These observations indicate 0.52±0.10 m of right-lateral slip and 0.17±0.04 m of reverse slip on a buried vertical 6-km-long and 5-km-deep fault and yield a mean static stress drop of 7.2±1.3 MPa. The largest shock had M S = 5.6. Observations of the ground rupture revealed up to 0.1 m of right-lateral slip on two mapped faults that are subparallel to the modeled seismic slip plane. In the 1.9 years since the earthquakes, geodetic network displacements indicate that an additional 60±10 mm of postseismic creep took place. The rate of postseismic shear strain (0.53±0.13 μrad/yr) measured within a 30×30-km network centered on the principal events was anomalously high compared to its preearthquake value and the postseismic rate in the adjacent network. This transient cannot be explained by postseismic slip on the seismic fault but rather indicates that broadscale release of strain followed the earthquake sequence. We have calculated the postearthquake stress field caused by the modeled coseismic slip. We assume that failure is promoted when the sum of the shear stress plus 0.75 times the faultopening stress increases. Most aftershocks concentrate at points where the stresses are enhanced by 0.3 MPa (3 bars) or more; aftershocks are nearly absent where postearthquake stresses decrease by 0.3–0.5 MPa. Isolated off-fault clusters of aftershocks that locate at one fault length from the rupture plane are explainable by this hypothesis. We find that ground rupture and postseismic creep take place where near-surface stresses are calculated to increase within the preexisting fault zones. Two patches that extend 4 km from both ends of the seismic fault exhibited neither aftershocks nor measurable postseismic creep. The sensitivity of aftershocks and ground rupture to changes in stress that are less than 5% of the earthquake stress drop demonstrates that the region around the earthquakes was within a few percent of its failure threshold before the main shocks. The preearthquake stress field and the stress required for failure must also have been nearly uniform.

Journal of Geophysical Research Solid Earth↗
USGS2007

Earthquake-by-earthquake fold growth above the Puente Hills blind thrust fault, Los Angeles, California: Implications for fold kinematics and seismic hazard

Boreholes and high-resolution seismic reflection data collected across the forelimb growth triangle above the central segment of the Puente Hills thrust fault (PHT) beneath Los Angeles, California, provide a detailed record of incremental fold growth during large earthquakes on this major blind thrust fault. These data document fold growth within a discrete kink band that narrows upward from ???460 m at the base of the Quaternary section (200-250 m depth) to <150 m at 2.5 m depth, with most growth during the most recent folding event occurring within a zone only ???60 m wide. These observations, coupled with evidence from petroleum industry seismic reflection data, demonstrate that most (>82% at 250 m depth) folding and uplift occur within discrete kink bands, thereby enabling us to develop a paleoseismic history of the underlying blind thrust fault. The borehole data reveal that the youngest part of the growth triangle in the uppermost 20 m comprises three stratigraphically discrete growth intervals marked by southward thickening sedimentary strata that are separated by intervals in which sediments do not change thickness across the site. We interpret the intervals of growth as occurring after the formation of now-buried paleofold scarps during three large PHT earthquakes in the past 8 kyr. The intervening intervals of no growth record periods of structural quiescence and deposition at the regional, near-horizontal stream gradient at the study site. Minimum uplift in each of the scarp-forming events, which occurred at 0.2-2.2 ka (event Y), 3.0-6.3 ka (event X), and 6.6-8.1 ka (event W), ranged from ???1.1 to ???1.6 m, indicating minimum thrust displacements of ???2.5 to 4.5 m. Such large displacements are consistent with the occurrence of large-magnitude earthquakes (Mw > 7). Cumulative, minimum uplift in the past three events was 3.3 to 4.7 m, suggesting cumulative thrust displacement of ???7 to 10.5 m. These values yield a minimum Holocene slip rate for the PHT of ???0.9 to 1.6 mm/yr. The borehole and seismic reflection data demonstrate that dip within the kink band is acquired incrementally, such that older strata that have been deformed by more earthquakes dip more steeply than younger strata. Specifically, strata dip 0.4?? at 4 m depth, 0.7?? at 20 m depth, 8?? at 90 m, 16?? at 110 m, and 17?? at 200 m. Moreover, structural restorations of the borehole data show that the locus of active folding (the anticlinal active axial surface) does not extend to the surface in exactly the same location from earthquake to earthquake. Rather, that the axial surfaces migrate from earthquake to earthquake, reflecting a component of fold growth by kink band migration. The incremental acquisition of bed dip in the growth triangle may reflect some combination of fold growth by limb rotation in addition to kink band migration, possibly through a component of trishear or shear fault bend folding. Alternatively, the component of limb rotation may result from curved hinge fault bend folding, and/or the mechanical response of loosely consolidated granular sediments in the shallow subsurface to folding at depth. Copyright 2007 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Analysis of time series of glacier speed: Columbia Glacier, Alaska

During the summers of 1984 and 1985, laser measurements were made of the distance from a reference location to markers on the surface of the lower reach of Columbia Glacier, Alaska. The data have numerous gaps, mostly because of inclement weather. The laser measurements were corrected for variations in atmospheric temperature and pressure and for the trajectory of the marker. The marker speed was calculated by using a combination of cubic splines and digital filters to arrive at a data set with 1-hour time intervals. The speed varies from 7 to 15 m/d and has three noteworthy components: (1) a low-frequency perturbation in speed with a time scale of days related to increased precipitation, (2) semidiurnal and diurnal variations related to sea tides, and (3) diurnal variations related to glacier surface melt. The low frequency and tidal period variations are separated by low-pass filtering and subtracting the filtered data from the original data. The variations from melt runoff and tides have the same frequency range and are examined using harmonic analysis and multiple regression with the various forcing functions. The high-frequency portion of the ice speed signal is dominated by tidal influences and to a lesser extent by meltwater influences. The low-frequency portion of the signal is dominated by the effects of precipitation.

Journal of Geophysical Research Solid Earth↗
USGS2012

Venus cartography

The entire surface of the planet Venus is being mapped at global and regional scales (1:50,000,000 through 1:1,500,000) with synthetic aperture radar (SAR), radar altimeter, and radiometer measurements of physical properties from the Magellan spacecraft. The mapping includes SAR image mosaics, shaded relief maps, and topographic contour overlays made from altimetry data and by radargrammetric methods. Methods used include new techniques of radar image processing that became operational as a result of the Magellan mission. The area being mapped (at a resolution of 75 m/pixel) is roughly equivalent to that of Earth, including seafloors. The mapping is designed to support geologic and geophysical investigations.

Journal of Geophysical Research E: Planets↗
USGS2000

Viscoelastic-coupling model for the earthquake cycle driven from below

In a linear system the earthquake cycle can be represented as the sum of a solution which reproduces the earthquake cycle itself (viscoelastic-coupling model) and a solution that provides the driving force. We consider two cases, one in which the earthquake cycle is driven by stresses transmitted along the schizosphere and a second in which the cycle is driven from below by stresses transmitted along the upper mantle (i.e., the schizosphere and upper mantle, respectively, act as stress guides in the lithosphere). In both cases the driving stress is attributed to steady motion of the stress guide, and the upper crust is assumed to be elastic. The surface deformation that accumulates during the interseismic interval depends solely upon the earthquake-cycle solution (viscoelastic-coupling model) not upon the driving source solution. Thus geodetic observations of interseismic deformation are insensitive to the source of the driving forces in a linear system. In particular, the suggestion of Bourne et al . [1998] that the deformation that accumulates across a transform fault system in the interseismic interval is a replica of the deformation that accumulates in the upper mantle during the same interval does not appear to be correct for linear systems.

Journal of Geophysical Research B: Solid Earth↗
USGS2017

Integrating puffing and explosions in a general scheme for Strombolian-style activity

Strombolian eruptions are among the most common subaerial styles of explosive volcanism worldwide. Distinctive features of each volcano lead to a correspondingly wide range of variations of magnitude and erupted products, but most papers focus on a single type of event at a single volcano. Here, in order to emphasize the common features underlying this diversity of styles, we scrutinize a database from 35 different erupting vents, including 21 thermal infrared videos from Stromboli (Italy), Etna (Italy), Yasur (Vanuatu), and Batu Tara (Indonesia), from puffing, through rapid explosions to normal explosions, with variable ejection parameters and relative abundance of gas, ash, and bombs. Using field observations and high-speed thermal infrared videos processed by a new algorithm, we identify the distinguishing characteristics of each type of activity and how they may relate and interact. In particular, we record that ash-poor normal explosions may be preceded and followed by the onset or the increase of the puffing activity, while ash-rich explosions are emergent, i.e., with inflation of the free surface followed directly by emission of increasingly large gas pockets. Overall, we see that all Strombolian activities form a continuum arising from a common mechanism and are modulated by the combination of two well-established controls: (1) the length of the bursting gas pocket with respect to the vent diameter and (2) the presence and thickness of a high-viscosity layer in the uppermost part of the volcanic conduit.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Three-dimensional gravity modeling of the geologic structure of Long Valley caldera

A 48-mGal gravity low coincides with Long Valley caldera and is mainly attributed to low-density caldera fill. Gravity measurements by Unocal Geothermal have been integrated with U.S. Geological Survey data, vastly improving gravity station coverage throughout the caldera. A strong regional gravity trend is mainly attributed to isostasy. A “best fitting” (based on regional control of basement densities) Airy-Heiskanen isostatic model was used for the regional correction. A three-dimensional, multiple-unit gravity modeling program with iterative capabilities was developed to model the residual gravity. The density structure of Long Valley caldera and vicinity was modeled with 22 discrete density units, most of which were based on geologic units. Information from drill hole lithologies, surface geology, and structural geology interpretations constrain the model. Some important points revealed by the three-dimensional gravity modeling are that (1) the volume of ejected magma associated with the Bishop Tuff eruption is greater than previously thought, (2) the caldera structure is strongly influenced by precaldera topography and the extensions of major, active faults, (3) the main west ring fracture is coincident with the Inyo Domes-Mono Craters fracture system, (4) a relatively low-density region probably underlies the caldera, and (5) a silicic magma chamber may underlie Devils Postpile.

Journal of Geophysical Research Solid Earth↗
USGS2012

Variations in fault slip and strain accumulation at Parkfield, California: Initial results using two-color geodimeter measurements 1984-1988

Repeated length measurements of several geodetic baselines near Parkfield, California, have revealed significant variations in the local rates of shallow fault slip and strain. This network of baselines, surveyed several times each week, straddles the San Andreas fault in the transition zone between the creeping section to the northwest and the locked section to the southeast. The length measurements, characterized by a precision approaching 0.1 ppm, reveal large fluctuations in the rates of baseline extension. Principal mode analysis of the length change data indicates that the two largest components of the signal are (1) secular extensions and contractions consistent with surface slip on the main strand of the San Andreas fault, and (2) a large seasonal oscillation with no obvious spatial coherence. On most of these baselines, the second component appears to be in phase with seasonal rainfall. When data from the baselines with the largest amplitude of the seasonal signal are excluded, the remaining data can be modeled in terms of both spatial and temporal variations in surface slip, variations in the components of the changes in uniform strain, and the possible displacement of the central monument in this radial network. In parameterizing this model, the spatial variation of slip beneath the near surface is reflected by changes in shear strain. Although the computed secular shear is highly dependent upon the specified parameterization of surface slip, the data are consistent with the hypothesis that slip at intermediate depths lags behind the surface slip rate. However, the range in models that fit the data does not necessarily imply that there is a deficit in slip at depth relative to the surface. Comparison of the inferred values of surface slip from the model with the observed fault slip measured by very short baseline creep meters indicates close agreement in secular rates, but the short-term variations observed with the creep meters are either highly attenuated or nonexistent in the modeled slip since the modeled slip is a spatial average which smooths out possible short-wavelength variations in the surface slip for which the creep instruments are most sensitive. An interesting conclusion from the two-color data is that surface slip on the San Andreas fault appears to be spread over a 2-km-wide zone on the south flank of Middle Mountain but is confined to a very narrow zone to the south as the fault passes through the center of the network. This conclusion is dependent upon the assumption that a few critical monuments are stable and track tectonic displacements in the long term. Finally, the largest observed strain change is an extensional strain coincident with the Kettleman Hills earthquake M 5.5 in August 1985.

Journal of Geophysical Research Solid Earth↗
USGS2012

Lithospheric loading by the 1896 Riku-u earthquake, northern Japan: Implications for plate flexure and asthenospheric rheology

Under favorable circumstances the time-dependent aseismic deformation resulting from the loading of the lithosphere by the stress drop of large dip slip earthquakes can be used to determine both the effective elastic plate thickness and the asthenospheric viscosity. The deformation has several similarities with the deflection of the lithosphere by surface loads and with movements due to postglacial rebound. Level changes obtained in the 80 years since the M = 7.5, 1896 Riku-u earthquake, an intraplate thrust event in northern Honshu, provide convincing evidence that asthenospheric readjustments are responsible for the observed movements. Leveling surveys crossing the zone of surface faulting have been repeated five times since 1900 and delineate a localized depression that has subsided at a continually decreasing rate. The depression is centered close to the 1896 faulting, and its shape and width, about 75 km, are matched by our model using a plate thickness of 30 km. The decaying subsidence rate constrains the viscosity of the uppermost asthenosphere to be 1×10 20 P. A linear viscous rheology matches the observed decay quite well, although measurements are sparse during the several decades following the earthquake.

Journal of Geophysical Research Solid Earth↗
USGS2001

Permeability reduction in granite under hydrothermal conditions

The formation of impermeable fault seals between earthquake events is a feature of many models of earthquake generation, suggesting that earthquake recurrence may depend in part on the rate of permeability reduction of fault zone materials under hydrothermal conditions. In this study, permeability measurements were conducted on intact, fractured, and gouge-bearing Westerly granite at an effective pressure of 50 MPa and at temperatures from 150° to 500°C, simulating conditions in the earthquake-generating portions of fault zones. Pore fluids were cycled back and forth under a 2 MPa pressure differential for periods of up to 40 days. Permeability of the granite decreased with time t , following the exponential relation k = c (10 − rt ). For intact samples run between 250° and 500°C the time constant for permeability decrease r was proportional to temperature and ranged between 0.001 and 0.1 days −1 (i.e., between 0.4 and 40 decades year −1 loss of permeability). Values of r for the lower-temperature experiments differed little from the 250°C runs. In contrast, prefractured samples showed higher rates of permeability decrease at a given temperature. The surfaces of the fractured samples showed evidence of dissolution and mineral growth that increased in abundance with both temperature and time. The experimentally grown mineral assemblages varied with temperature and were consistent with a rock-dominated hydrothermal system. As such mineral deposits progressively seal the fractured samples, their rates of permeability decrease approach the rates for intact rocks at the same temperature. These results place constraints on models of precipitation sealing and suggest that fault rocks may seal at a rate consistent with earthquake recurrence intervals of typical fault zones.

California↗
USGS2013

Annual modulation of non-volcanic tremor in northern Cascadia

Two catalogs of episodic tremor events in northern Cascadia, one from 2006 to 2012 and the other from 1997 to 2011, reveal two systematic patterns of tremor occurrence in southern Vancouver Island: (1) most individual events tend to occur in the third quarter of the year; (2) the number of events in prolonged episodes (i.e., episodic tremor and slip events), which generally propagate to Vancouver Island from elsewhere along the Cascadia subduction zone, is inversely correlated with the amount of precipitation that occurred in the preceding 2 months. We rationalize these patterns as the product of hydrologic loading of the crust of southern Vancouver Island and the surrounding continental region, superimposed with annual variations from oceanic tidal loading. Loading of the Vancouver Island crust in the winter (when the land surface receives ample precipitation) and unloading in the summer tends to inhibit and enhance downdip shear stress, respectively. Quantitatively, for an annually variable surface load, the predicted stress perturbation depends on mantle viscoelastic rheology. A mechanical model of downdip shear stress on the transition zone beneath Vancouver Island—driven predominantly by the annual hydrologic cycle—is consistent with the 1997–2012 tremor observations, with peak-to-peak downdip shear stress of about 0.4 kPa. This seasonal dependence of tremor occurrence appears to be restricted to southern Vancouver Island because of its unique situation as an elongated narrow-width land mass surrounded by ocean, which permits seasonal perturbations in shear stress at depth.

Journal of Geophysical Research B: Solid Earth↗
USGS2020

Peak ground velocity spatial variability revealed by dense seismic array in southern California

Understanding and modeling variability of ground motion is essential for building accurate and precise ground motion prediction equations, which can net site‐specific characterization and reduced hazard levels. Here, we explore the spatial variability in peak ground velocity (PGV) at Sage Brush Flats along the San Jacinto Fault in southern California. We use data from a dense array (0.6 x 0.6 km 2 , 1,108 geophones, station spacings 10‐30 m) deployed in 2014 for ~1‐month. These data offer an opportunity to study small scale variability in this region. We examine 38 2≤M L ≤4.2 earthquakes within 200 km of the array. Fault strands and a small basin impact the ground motions, producing PGV variations up to 22% of the mean and a 40% reduction in P and S wave near‐surface velocities. We find along‐fault rupture directivity, source, and path effects can increase PGVs by 167%. Surface PGV measurements exceed the co‐located borehole station (depth 148 m) PGV by factors of 3‐10, confirming the impact on PGV from near surface fault structures, basins, topography, and amplifications from soft sediments. Consistently we find high PGVs within the basin structure. A pair of ~co‐located M L 2.6 events produce repeatable PGV values with similar spatial patterns. The average corner frequencies of these two events are 11‐16 Hz and viable measurements of stress drop can differ by 6.45MPa. Within this small array, the PGV values are variable implying spatial extrapolation of PGV to regions of known faults and basins, even across a small area, should be done with caution.

California↗
USGS2022

Toward next-generation lava flow forecasting: Development of a fast, physics-based lava propagation model

During effusive volcanic crises, the eruption and propagation of lava flows pose a significant hazard to nearby populations, homes, and infrastructure. Consequently, timely lava flow forecasts are a critical need for volcano observatory and emergency management operations. Previous lava flow modeling tools are typically either too slow to produce timely forecasts, or are fast, but lack critical aspects of lava physics or important forecasting outputs. In particular, the strong thermal stratification present in laminar, high-Prandtl number flows has generally been neglected. Bulk rheological changes have previously been computed from cell-averaged temperatures, assuming that the flow is thermally mixed. Here, we detail the development and initial testing of Lava2d, a new two-dimensional depth-averaged finite volume model of lava flow propagation over natural terrain which accounts for bulk rheological changes due to thermorheological stratification. We use a novel approach to energy conservation based on tracking cooling and solidifying at the flow base and at the moving flow surface, allowing for the estimation of more realistic vertical thermorheological profiles, while maintaining computational efficiency, producing very timely model runs. We validate our approach with three examples: comparison with theoretical propagation of crust-dominated lava flows, comparison with a large-scale molten basalt experiment from the Syracuse University Lava Project, and efficiency testing and comparison with the initial phase of the 1984 Mauna Loa lava flows. Our model is shown to produce rapid, realistic forecasts, making it a good candidate for operationalization in active volcanic regions such as in Hawai'i.

Journal of Geophysical Research - Solid Earth↗
USGS2011

A wavelength-dependent visible and infrared spectrophotometric function for the Moon based on ROLO data

The USGS's Robotic Lunar Observatory (ROLO) dedicated ground‐based lunar calibration project obtained photometric observations of the Moon over the spectral range attainable from Earth (0.347–2.39 μ m) and over solar phase angles of 1.55°–97°. From these observations, we derived empirical lunar surface solar phase functions for both the highlands and maria that can be used for a wide range of applications. The functions can be used to correct for the effects of viewing geometry to produce lunar mosaics, spectra, and quick‐look products for future lunar missions and ground‐based observations. Our methodology can be used for a wide range of objects for which multiply scattered radiation is not significant, including all but the very brightest asteroids and moons.

Journal of Geophysical Research E: Planets↗
USGS2012

Gravity model studies of Newberry Volcano, Oregon

Newberry Volcano, a large Quaternary volcano located about 60 km east of the axis of the High Cascades volcanoes in central Oregon, has a coincident positive residual gravity anomaly of about 12 mGals. Model calculations of the gravity anomaly field suggest that the volcano is underlain by an intrusive complex of mafic composition of about 20-km diameter and 2-km thickness, at depths above 4 km below sea level. However, uplifted basement in a northwest trending ridge may form part of the underlying excess mass, thus reducing the volume of the subvolcanic intrusive. A ring dike of mafic composition is inferred to intrude to near-surface levels along the caldera ring fractures, and low-density fill of the caldera floor probably has a thickness of 0.7–0.9 km. The gravity anomaly attributable to the volcano is reduced to the east across a north-northwest trending gravity anomaly gradient through Newberry caldera and suggests that normal, perhaps extensional, faulting has occurred subsequent to caldera formation and may have controlled the location of some late-stage basaltic and rhyolitic eruptions. Significant amounts of felsic intrusive material may exist above the mafic intrusive zone but cannot be resolved by the gravity data.

Journal of Geophysical Research Solid Earth↗