Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Journal of Geophysical Research - Atmospheres”

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 649 records · Page 36Linked to original sources
USGS2008

A prelanding assessment of the ice table depth and ground ice characteristics in Martian permafrost at the Phoenix landing site

We review multiple estimates of the ice table depth at potential Phoenix landing sites and consider the possible state and distribution of subsurface ice. A two-layer model of ice-rich material overlain by ice-free material is consistent with both the observational and theoretical lines of evidence. Results indicate ground ice to be shallow and ubiquitous, 2-6 cm below the surface. Undulations in the ice table depth are expected because of the thermodynamic effects of rocks, slopes, and soil variations on the scale of the Phoenix Lander and within the digging area, which can be advantageous for analysis of both dry surficial soils and buried ice-rich materials. The ground ice at the ice table to be sampled by the Phoenix Lander is expected to be geologically young because of recent climate oscillations. However, estimates of the ratio of soil to ice in the ice-rich subsurface layer suggest that that the ice content exceeds the available pore space, which is difficult to reconcile with existing ground ice stability and dynamics models. These high concentrations of ice may be the result of either the burial of surface snow during times of higher obliquity, initially high-porosity soils, or the migration of water along thin films. Measurement of the D/H ratio within the ice at the ice table and of the soil-to-ice ratio, as well as imaging ice-soil textures, will help determine if the ice is indeed young and if the models of the effects of climate change on the ground ice are reasonable.

Journal of Geophysical Research E: Planets↗
USGS2008

Compositional trends in aeolian dust along a transect across the southwestern United States

Aeolian dust strongly influences ecology and landscape geochemistry over large areas that span several desert ecosystems of the southwestern United States. This study evaluates spatial and temporal variations and trends of the physical and chemical properties of dust in the southwestern United States by examining dust deposited in natural depressions on high isolated surfaces along a transect from the Mojave Desert to the central Colorado Plateau. Aeolian dust is recognized in these depressions on the basis of textural, chemical, isotopic, and mineralogical characteristics and comparisons of those characteristics to the underlying bedrock units. Spatial and temporal trends suggest that although local dust sources are important to the accumulated material in these depressions, Mojave Desert dust sources may also contribute. Depth trends in the depressions suggest that Mojave sources may have contributed more dust to the Colorado Plateau recently than in the past. These interpretations point to the important roles of far-traveled aeolian dust for landscape geochemistry and imply future changes to soil geochemistry under changing conditions in far-distant dust source areas. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research F: Earth Surface↗
USGS2008

Magnetic fabric of sheared till: A strain indicator for evaluating the bed deformation model of glacier flow

Wet-based portions of ice sheets may move primarily by shearing their till beds, resting in high sediment fluxes and the development of subglacial landforms. This model of glacier movement, which requires high bed shear strains, can be tested using till microstructural characteristics that evolve during till deformation. Here we examine the development of magnetic fabric using a ring shear device to defom two Wisconsin-age basal tills to shear strains as high as 70. Hysteresis experiments and the dependence of magnetic susceptibility of these tills on temperature demonstrate that anisotropy of magnetic susceptibility (AMS) develops during shear due to the rotation of primarily magnetite particles that are silt sized or smaller. At moderate shear strains (???6-25), principal axes of maximum magnetic susceptibility develop a strong fabric (S1 eignevalues of 0.83-0.96), without further strengthening at higher strains, During deformation, directions of maximum susceptibility cluster strongly in the direction of shear and plunge 'up-glacier,' consistent with the behavior of pebbles and sand particles studied in earlier experiments. In contrast, the magnitude of AMS does not vary systematically with strain and is small relative to its variability among samples; this is because most magnetite grains are contained as inclusions in larger particles and hence do not align during shear. Although processes other than pervasive bed deformation may result in strong flow parallel fabrics, AMS fabrics provide a rapid and objective means of identifying basal tills that have not been sheared sufficiently to be compatible with the bed deformation model. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research F: Earth Surface↗
USGS2009

Transport slopes, sediment cover, and bedrock channel incision in the Henry Mountains, Utah

[1] Field data from channels in the Henry Mountains of Utah demonstrate that abundant coarse sediment can inhibit fluvial incision into bedrock by armoring channel beds (the cover effect). We compare several small channels that share tributary junctions and have incised into the same sedimentary bedrock unit (Navajo Sandstone) but contain differing amounts of coarse diorite clasts owing to the spatial distribution of localized sediment sources. Bedrock channels that contain abundant clasts (diorite-rich) have steeper longitudinal slopes than tributaries of these channels with smaller drainage areas and less sediment (diorite-poor). The diorite-poor tributaries have incised more deeply to lower average slopes and have more reach-scale slope variability, which may reflect bedrock properties, longitudinal sediment sorting, and incision at lower sediment supply. Diorite-rich channels have less bedrock exposed and smoother longitudinal profiles than diorite-poor channels. We find that (1) coarse sediment can mantle bedrock channel beds and reduce the efficiency of incision, validating the hypothesized cover effect in fluvial incision models; (2) the channel slope needed to transport the sediment load can be larger than that needed to erode bedrock, suggesting that the slope of incising bedrock channels can become adjusted to the sediment load; (3) when abundant sediment is available, transport capacity rather than thresholds of motion can be dominant in setting bedrock channel slope; and (4) cover effects can be important even when moderate amounts of bedrock are exposed in channel beds.

Utah↗
USGS2008

Modeling the evolution of channel shape: Balancing computational efficiency with hydraulic fidelity

The cross-sectional shape of a natural river channel controls the capacity of the system to carry water off a landscape, to convey sediment derived from hillslopes, and to erode its bed and banks. Numerical models that describe the response of a landscape to changes in climate or tectonics therefore require formulations that can accommodate evolution of channel cross-sectional geometry. However, fully two-dimensional (2-D) flow models are too computationally expensive to implement in large-scale landscape evolution models, while available simple empirical relationships between width and discharge do not adequately capture the dynamics of channel adjustment. We have developed a simplified 2-D numerical model of channel evolution in a cohesive, detachment-limited substrate subject to steady, unidirectional flow. Erosion is assumed to be proportional to boundary shear stress, which is calculated using an approximation of the flow field in which log-velocity profiles are assumed to apply along vectors that are perpendicular to the local channel bed. Model predictions of the velocity structure, peak boundary shear stress, and equilibrium channel shape compare well with predictions of a more sophisticated but more computationally demanding ray-isovel model. For example, the mean velocities computed by the two models are consistent to within ???3%, and the predicted peak shear stress is consistent to within ???7%. Furthermore, the shear stress distributions predicted by our model compare favorably with available laboratory measurements for prescribed channel shapes. A modification to our simplified code in which the flow includes a high-velocity core allows the model to be extended to estimate shear stress distributions in channels with large width-to-depth ratios. Our model is efficient enough to incorporate into large-scale landscape evolution codes and can be used to examine how channels adjust both cross-sectional shape and slope in response to tectonic and climatic forcing. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research F: Earth Surface↗
USGS2007

Hydrogeologic controls on nitrate transport in a small agricultural catchment, Iowa

Effects of subsurface deposits on nitrate loss in stream riparian zones are recognized, but little attention has been focused on similar processes occurring in upland agricultural settings. In this paper, we evaluated hydrogeologic controls on nitrate transport processes occurring in a small 7.6 ha Iowa catchment. Subsurface deposits in the catchment consisted of upland areas of loess overlying weathered pre-Illinoian till, drained by two ephemeral drainageways that consisted of Holocene-age silty and organic rich alluvium. Water tables in upland areas fluctuated more than 4 m per year compared to less than 0.3 m in the drainageway. Water quality patterns showed a distinct spatial pattern, with groundwater in the drainageways having lower nitrate concentrations (<0.5 mg L-1 compared to upland areas (>10 mg L-1) as wells as lower pH, dissolved oxygen and redox, and higher ammonium and dissolved organic carbon levels. Several lines of evidence suggested that conditions are conducive for denitrification of groundwater flowing from uplands through the drainageways. Field-measured nitrate decay rates in the drainageways (???0.02 day-1) were consistent with other laboratory studies and regional patterns. Results from MODFLOW and MT3DMS simulations indicated that soils in the ephemeral drainageways could process all upland groundwater nitrate flowing through them. However, model-simulated tile drainage increased both water flux and nitrate loss from the upland catchment. Study results suggest that ephemeral drainageways can provide a natural nitrate treatment system in our upland glaciated catchments, offering management opportunities to reduce nitrate delivery to streams. Copyright 2007 by the American Geophysical Union.

Journal of Geophysical Research G: Biogeosciences↗
USGS2007

Influence of disturbance on carbon exchange in a permafrost collapse and adjacent burned forest

We measured CO 2 and CH 4 exchange from the center of a Sphagnum ‐dominated permafrost collapse, through an aquatic moat, and into a recently burned black spruce forest on the Tanana River floodplain in interior Alaska. In the anomalously dry growing season of 2004, both the collapse and the surrounding burned area were net sinks for CO 2 , with a mean daytime net ecosystem exchange of −1.4 μ mol CO 2 m −2 s −1 , while the moat was a CH 4 source with a mean flux of 0.013 μ mol CH 4 m −2 s −1 . Regression analyses identified temperature as the dominant factor affecting intragrowing season variation in CO 2 exchange and soil moisture as the primary control influencing CH 4 emissions. CH 4 emissions during the wettest portion of the growing season were four times higher than during the driest periods. If temperatures continue to warm, peatland vegetation will likely expand with permafrost degradation, resulting in greater carbon accumulation and methane emissions for the landscape as a whole.

Alaska↗
USGS2008

Oceanic loading of wildfire-derived organic compounds from a small mountainous river

Small mountainous rivers (SMRs) export substantial amounts of sediment into the world's oceans. The concomitant yield of organic carbon (OC) associated with this class of rivers has also been shown to be significant and compositionally unique. We report here excessively high loadings of polycyclic aromatic hydrocarbons (PAHs), lignin, and levoglucosan, discharged from the Santa Clara River into the Santa Barbara Channel. The abundance of PAHs, levoglucosan, and lignin in Santa Barbara Channel sediments ranged from 201.7 to 1232.3 ng gdw-1, 1.3 to 6.9 ??g gdw-1, and 0.3 to 2.2 mg per 100 mg of the sedimentary OC, respectively. Assuming a constant rate of sediment accumulation, the annual fluxes of PAHs, levoglucosan, and lignin, to the Santa Barbara Channel were respectively, 885.5 ?? 170.2 ng cm-2 a-1, 3.5 ?? 1.9 ??g cm-2 a-1 and 1.4 ?? 0.3 mg per 100 mg OC cm-2 a-1, over ???30 years. The close agreement between PAHs, levoglucosan, and lignin abundance suggests that the depositional flux of these compounds is largely biomass combustion-derived. To that end, use of the Santa Clara River as a model for SMRs suggests this class of rivers may be one of the largest contributors of pyrolyzed carbon to coastal systems and the open ocean. Wildfire associated carbon discharged from other high yield fluvial systems, when considered collectively, may be a significant source of lignin, pyrolytic PAHs, and other pyrogenic compounds to the ocean. Extrapolating these methods over geologic time may offer useful historical information about carbon sequestration and burial in coastal sediments and affect coastal carbon budgets. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research G: Biogeosciences↗
USGS2008

Changes in the soil C cycle at the arid-hyperarid transition in the Atacama Desert

We examined soil organic C (OC) turnover and transport across the rainfall transition from a biotic, arid site to a largely abiotic, hyperarid site. With this transition, OC concentrations decrease, and C cycling slows precipitously, both in surface horizons and below ground. The concentration and isotopic character of soil OC across this transition reflect decreasing rates of inputs, decomposition, and downward transport. OC concentrations in the arid soil increase slightly with depth in the upper meter, but are generally low and variable (???0.05%; total inventory of 1.82 kg m-2); OC-??14C values decrease from modern (+7???) to very 14C-depleted (-966???) with depth; and OC-??13C values are variable (-23.7??? to -14.1???). Using a transport model, we show that these trends reflect relatively rapid cycling in the upper few centimeters, and spatially variable preservation of belowground OC from root inputs, possibly during a previous, wetter climate supporting higher soil OC concentrations. In the driest soil, the OC inventory is the lowest among the sites (0.19 kg m-2), and radiocarbon values are 14C-depleted (-365??? to -696???) but show no trend with depth, indicating belowground OC inputs and long OC residence times throughout the upper meter (104 y). A distinct depth trend in ??13C values and OC/ON values within the upper 40 cm at the driest site may reflect photochemical alteration of organic matter at the soil surface, combined with limited subsurface decomposition and downward transport. We argue that while root inputs are preserved at the wetter sites, C cycling in the most hyperarid soil occurs through infrequent, rapid dissolved transport of highly photodegraded organic matter during rare rain events, each followed by a pulse of decomposition and subsequent prolonged drought. These belowground inputs are likely a primary control on the character, activity, and depth distribution of small microbial populations. While the lack of water is the dominant control on C cycling, very low C/N ratios of organic matter suggest that when rainfall occurs, hyperarid soils are effectively C limited. The preservation of fossil root fragments in the sediment beneath the driest soil indicates that wetter climate conditions preceded formation of this soil, and that vadose zone microbial activity has been extremely limited for the past 2 My. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research G: Biogeosciences↗
USGS2008

Short-term response of methane fluxes and methanogen activity to water table and soil warming manipulations in an Alaskan peatland

Growing season CH 4 fluxes were monitored over a two year period following the start of ecosystem-scale manipulations of water table position and surface soil temperatures in a moderate rich fen in interior Alaska. The largest CH 4 fluxes occurred in plots that received both flooding (raised water table position) and soil warming, while the lowest fluxes occurred in unwarmed plots in the lowered water table treatment. A combination of treatment and soil hydroclimate variables explained more than 70% of the variation in ln-transformed CH 4 fluxes, with mean daily water table position representing the strongest predictor. We used quantitative PCR of the α -subunit of mcr operon to explore the influence of soil climate manipulations on methanogen abundances. Methanogen abundances were greatest in warmed plots, and showed a positive relationship with mean daily CH 4 fluxes. Our results show that water table manipulations that led to soil inundation (flooding) had a stronger effect on CH 4 fluxes than water table drawdown. Seasonal CH 4 fluxes increased by 80–300% under the combined wetter and warmer soil climate treatments. Thus, while warming is expected to increase CH 4 emissions from Alaskan wetlands, higher water table positions caused by increases in precipitation or disturbances such as permafrost thaw that lead to thermokarst and flooding in wetlands will stimulate CH 4 emissions beyond the effects of soil warming alone. Consequently, we argue that modeling the effects of climate change on Alaskan wetland CH 4 emissions needs to consider the interactive effects of soil warming and water table position on CH 4 production and transport.

Alaska↗
USGS2007

Sulfur species behavior in soil organic matter during decomposition

Soil organic matter (SOM) is a primary re??servoir of terrestrial sulfur (S), but its role in the global S cycle remains poorly understood. We examine S speciation by X-ray absorption near-edge structure (XANES) spectroscopy to describe S species behavior during SOM decomposition. Sulfur species in SOM were best represented by organic sulfide, sulfoxide, sulfonate, and sulfate. The highest fraction of S in litter was organic sulfide, but as decomposition progressed, relative fractions of sulfonate and sulfate generally increased. Over 6-month laboratory incubations, organic sulfide was most reactive, suggesting that a fraction of this species was associated with a highly labile pool of SOM. During humification, relative concentrations of sulfoxide consistently decreased, demonstrating the importance of sulfoxide as a reactive S phase in soil. Sulfonate fractional abundance increased during humification irrespective of litter type, illustrating its relative stability in soils. The proportion of S species did not differ systematically by litter type, but organic sulfide became less abundant in conifer SOM during decomposition, while sulfate fractional abundance increased. Conversely, deciduous SOM exhibited lesser or nonexistent shifts in organic sulfide and sulfate fractions during decomposition, possibly suggesting that S reactivity in deciduous litter is coupled to rapid C mineralization and independent of S speciation. All trends were consistent in soils across study sites. We conclude that S reactivity is related to spqciation in SOM, particularly in conifer forests, and S species fractions in SOM change, during decomposition. Our data highlight the importance of intermediate valence species (sulfoxide and sulfonate) in the pedochemical cycling of organic bound S. Copyright 2007 by the American Geophysical Union.

Journal of Geophysical Research G: Biogeosciences↗
USGS2009

Three‐dimensional model of Hellenic Arc deformation and origin of the Cretan uplift

[1] The Hellenic Arc of Greece is the most seismically active part of Europe, but little is know about its mechanics. We modeled deformation along the arc using a finite element model. The model was intended to capture large‐scale 3‐D structure of Nubian plate subduction beneath the Aegean block and its deformational consequences. The shape of the interface was developed using mapped traces at the surface and earthquake hypocenters at depth. Model block motions were constrained by recent compilations of GPS velocity vectors. We simulated a 10 ka period of convergence between Nubia and the Aegean and calculated the strain field in the overriding plate as well as the spatial distribution and orientation of differential stress (∣ σ 1 − σ 3 ∣). From these calculations we derived testable quantities such as the expected seismic moment rate on the interplate contact, uplift pattern, and distribution of strain modes. Our relatively simple model broadly reproduced observed uplift patterns, earthquake activity, and loci of extension and contraction. The model showed a localization of uplift near the island of Crete, where the fastest Aegean uplift rates are well documented. Comparison of calculated expected seismic moment and observed earthquake catalogs implies a nearly fully coupled interplate contact. On the basis of our modeling results, we suggest that south Aegean deformation is driven primarily by the fast moving (∼33 mm a −1 ) Aegean upper plate overriding a nearly stalled (∼5 mm a −1 ) Nubian lower plate. This tectonic setting thus more closely resembles a continental thrust than it does a typical oceanic subduction zone.

Aegean Sea, Hellenic Arc↗
USGS2008

Postearthquake relaxation after the 2004 M6 Parkfield, California, earthquake and rate-and-state friction

An unusually complete set of measurements (including rapid rate GPS over the first 10 days) of postseismic deformation is available at 12 continuous GPS stations located close to the epicenter of the 2004 M6.0 Parkfield earthquake. The principal component modes for the relaxation of the ensemble of those 12 GPS stations were determined. The first mode alone furnishes an adequate approximation to the data. Thus, the relaxation at all stations can be represented by the product of a common temporal function and distinct amplitudes for each component (north or east) of relaxation at each station. The distribution in space of the amplitudes indicates that the relaxation is dominantly strike slip. The temporal function, which spans times from about 5 min to 900 days postearthquake, can be fit by a superposition of three creep terms, each of the form ??l loge(1 + t/??l), with characteristic times ??, = 4.06, 0.11, and 0.0001 days. It seems likely that what is actually involved is a broad spectrum of characteristic times, the individual components of which arise from afterslip on different fault patches. Perfettini and Avouac (2004) have shown that an individual creep term can be explained by the spring-slider model with rate-dependent (no state variable) friction. The observed temporal function can also be explained using a single spring-slider model (i.e., single fault patch) that includes rate-and-state-dependent friction, a single-state variable, and either of the two commonly used (aging and slip) state evolution laws. In the latter fits, the rate-and-state friction parameter b is negative.

Journal of Geophysical Research B: Solid Earth↗
USGS2009

Constraining earthquake source inversions with GPS data: 2. A two-step approach to combine seismic and geodetic data sets

We present a new method to combine static and wavefield data to image earthquake ruptures. Our combined inversion is a two-step procedure, following the work of Hernandez et al. (1999), and takes into account the differences between the resolutions of the two data sets. The first step consists of an inversion of the static field, which yields a map of slip amplitude. This inversion exploits a special irregular grid that takes into account the resolution of the static data. The second step is an inversion of the radiated wavefield; it results in the determination of the time evolution of slip on the fault. In the second step, the slip amplitude is constrained to resemble the static slip amplitude map inferred from the GPS inversion. Using this combined inversion, we study the source process of the 2004 M6 Parkfield, California, earthquake. We conclude that slip occurred in two main regions of the fault, each of which displayed distinct rupture behaviors. Slip initiated at the hypocenter with a very strong bilateral burst of energy. Here, slip was localized in a narrow area approximately 10 km long, the rupture velocity was very fast (???3.5 km/s), and slip only lasted a short period of time (<1 s). Then the rupture proceeded to a wider region 12-20 km northwest of the hypocenter. Here, the earthquake developed in a more moderated way: the rupture velocity slowed to ???3.0 km/s and slip lasted longer (1-2 s). The maximum slip amplitude was 0.45 m. Copyright 2009 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2009

Premonitory acoustic emissions and stick-slip in natural and smooth-faulted Westerly granite

A stick-slip event was induced in a cylindrical sample of Westerly granite containing a preexisting natural fault by loading at constant confining pressure of 150 MPa. Continuously recorded acoustic emission (AE) data and computer tomography (CT)-generated images of the fault plane were combined to provide a detailed examination of microscale processes operating on the fault. The dynamic stick-slip event, considered to be a laboratory analog of an earthquake, generated an ultrasonic signal that was recorded as a large-amplitude AE event. First arrivals of this event were inverted to determine the nucleation site of slip, which is associated with a geometric asperity on the fault surface. CT images and AE locations suggest that a variety of asperities existed in the sample because of the intersection of branch or splay faults with the main fault. This experiment is compared with a stick-slip experiment on a sample prepared with a smooth, artificial saw-cut fault surface. Nearly a thousand times more AE were observed for the natural fault, which has a higher friction coefficient (0.78 compared to 0.53) and larger shear stress drop (140 compared to 68 MPa). However at the measured resolution, the ultrasonic signal emitted during slip initiation does not vary significantly between the two experiments, suggesting a similar dynamic rupture process. We propose that the natural faulted sample under triaxial compression provides a good laboratory analogue for a field-scale fault system in terms of the presence of asperities, fault surface heterogeneity, and interaction of branching faults. ?? 2009.

Journal of Geophysical Research B: Solid Earth↗
USGS2009

Regional three-dimensional seismic velocity model of the crust and uppermost mantle of northern California

We present a three-dimensional (3D) tomographic model of the P wave velocity (Vp) structure of northern California. We employed a regional-scale double-difference tomography algorithm that incorporates a finite-difference travel time calculator and spatial smoothing constraints. Arrival times from earthquakes and travel times from controlled-source explosions, recorded at network and/or temporary stations, were inverted for Vp on a 3D grid with horizontal node spacing of 10 to 20 km and vertical node spacing of 3 to 8 km. Our model provides an unprecedented, comprehensive view of the regional-scale structure of northern California, putting many previously identified features into a broader regional context and improving the resolution of a number of them and revealing a number of new features, especially in the middle and lower crust, that have never before been reported. Examples of the former include the complex subducting Gorda slab, a steep, deeply penetrating fault beneath the Sacramento River Delta, crustal low-velocity zones beneath Geysers-Clear Lake and Long Valley, and the high-velocity ophiolite body underlying the Great Valley. Examples of the latter include mid-crustal low-velocity zones beneath Mount Shasta and north of Lake Tahoe. Copyright 2009 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2009

Inducing in situ, nonlinear soil response applying an active source

[1] It is well known that soil sites have a profound effect on ground motion during large earthquakes. The complex structure of soil deposits and the highly nonlinear constitutive behavior of soils largely control nonlinear site response at soil sites. Measurements of nonlinear soil response under natural conditions are critical to advancing our understanding of soil behavior during earthquakes. Many factors limit the use of earthquake observations to estimate nonlinear site response such that quantitative characterization of nonlinear behavior relies almost exclusively on laboratory experiments and modeling of wave propagation. Here we introduce a new method for in situ characterization of the nonlinear behavior of a natural soil formation using measurements obtained immediately adjacent to a large vibrator source. To our knowledge, we are the first group to propose and test such an approach. Employing a large, surface vibrator as a source, we measure the nonlinear behavior of the soil by incrementally increasing the source amplitude over a range of frequencies and monitoring changes in the output spectra. We apply a homodyne algorithm for measuring spectral amplitudes, which provides robust signal-to-noise ratios at the frequencies of interest. Spectral ratios are computed between the receivers and the source as well as receiver pairs located in an array adjacent to the source, providing the means to separate source and near-source nonlinearity from pervasive nonlinearity in the soil column. We find clear evidence of nonlinearity in significant decreases in the frequency of peak spectral ratios, corresponding to material softening with amplitude, observed across the array as the source amplitude is increased. The observed peak shifts are consistent with laboratory measurements of soil nonlinearity. Our results provide constraints for future numerical modeling studies of strong ground motion during earthquakes.

Journal of Geophysical Research B: Solid Earth↗
USGS2009

Tracking magma volume recovery at Okmok Volcano using GPS and an unscented kalman filter

Changes beneath a volcano can be observed through position changes in a GPS network, but distinguishing the source of site motion is not always straightforward. The records of continuous GPS sites provide a favorable data set for tracking magma migration. Dense campaign observations usually provide a better spatial picture of the overall deformation field, at the expense of an episodic temporal record. Combining these observations provides the best of both worlds. A Kalman filter provides a means for integrating discrete and continuous measurements and for interpreting subtle signals. The unscented Kalman filter (UKF) is a nonlinear method for time‐dependent observations. We demonstrate the application of this technique to deformation data by applying it to GPS data collected at Okmok volcano. Seven years of GPS observations at Okmok are analyzed using a Mogi source model and the UKF. The deformation source at Okmok is relatively stable at 2.5 km depth below sea level, located beneath the center of the caldera, which means the surface deformation is caused by changes in the strength of the source. During the 7 years of GPS observations more than 0.5 m of uplift has occurred, a majority of that during the time period January 2003 to July 2004. The total volume recovery at Okmok since the last eruption in 1997 is ∼60–80%. The UKF allows us to solve simultaneously for the time‐dependence of the source strength and for the location without a priori information about the source. 2009 by the American Geophysical Union.

Alaska↗