USGS ScienceSearch

SEARCH · USGS Science

Results for “Journal of Geophysical Research - Oceans”

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 271 records · Page 15Linked to original sources

Functional group, biomass, and climate change effects on ecological drought in semiarid grasslands

Water relations in plant communities are influenced both by contrasting functional groups (grasses, shrubs) and by climate change via complex effects on interception, uptake and transpiration. We modelled the effects of functional group replacement and biomass increase, both of which can be outcomes of invasion and vegetation management, and climate change on ecological drought (soil water potential below which photosynthesis stops) in 340 semiarid grassland sites over 30‐year periods. Relative to control vegetation (climate and site‐determined mixes of functional groups), the frequency and duration of drought were increased by shrubs and decreased by annual grasses. The rankings of shrubs, control vegetation, and annual grasses in terms of drought effects were generally consistent in current and future climates, suggesting that current differences among functional groups on drought effects predict future differences. Climate change accompanied by experimentally‐increased biomass (i.e. the effects of invasions that increase community biomass, or management that increases productivity through fertilization or respite from grazing) increased drought frequency and duration, and advanced drought onset. Our results suggest that the replacement of perennial temperate semiarid grasslands by shrubs, or increased biomass, can increase ecological drought both in current and future climates.

Journal of Geophysical Research G: Biogeosciences

Nutrient dynamics in partially drained arctic thaw lakes

Thaw lakes are ubiquitous on arctic coastal plains (ACPs). While many thaw lakes have steep banks, stable water levels, and static surface areas, others only partially fill their basins and vary in area over the summer. These partially drained lakes (PDLs) are hydrologically connected to the wetlands immediately surrounding them. Heat and nutrient availability limit aquatic productivity on ACPs, and we hypothesized that shallow shorelines and greater hydrologic connectivity with the landscape should result in greater nutrient concentrations and biogeochemical cycling in PDLs. We tested this by monitoring water chemistry in lakes with varying levels of seasonal drainage in sandy and silty peaty lowland sites on the ACP of Alaska. One highly drained lake (N1) was significantly warmer than minimally drained lakes (minDLs) related to earlier ice off, reaching temperatures as high as 16 °C in June when minDLs still contained ice. Ammonia, total dissolved phosphorus, and dissolved organic carbon and nitrogen concentrations were higher in lakes with greater drainage, and concentrations in N1 rivaled those in the small, biologically productive ponds. Many PDLs displayed a midsummer decrease in nutrients consistent with assimilation by the aquatic ecosystem, and a late‐summer increase most likely related to runoff from drained lake margins following precipitation. N1 exported kilograms of ammonium and total dissolved phosphorus to the stream network over the summer. Given increased warming and drying in the arctic, the proportion of PDLs may be changing, which in turn may affect nutrient and organic matter availability in arctic lakes and export to downstream environments.

Journal of Geophysical Research: Biogeosciences

Quantifying postfire aeolian sediment transport using rare earth element tracers

Grasslands, which provide fundamental ecosystem services in many arid and semiarid regions of the world, are undergoing rapid increases in fire activity and are highly susceptible to postfire-accelerated soil erosion by wind. A quantitative assessment of physical processes that integrates fire-wind erosion feedbacks is therefore needed relative to vegetation change, soil biogeochemical cycling, air quality, and landscape evolution. We investigated the applicability of a novel tracer technique—the use of multiple rare earth elements (REE)—to quantify soil transport by wind and to identify sources and sinks of wind-blown sediments in both burned and unburned shrub-grass transition zone in the Chihuahuan Desert, NM, USA. Results indicate that the horizontal mass flux of wind-borne sediment increased approximately threefold following the fire. The REE tracer analysis of wind-borne sediments shows that the source of the horizontal mass flux in the unburned site was derived from bare microsites (88.5%), while in the burned site it was primarily sourced from shrub (42.3%) and bare (39.1%) microsites. Vegetated microsites which were predominantly sinks of aeolian sediments in the unburned areas became sediment sources following the fire. The burned areas showed a spatial homogenization of sediment tracers, highlighting a potential negative feedback on landscape heterogeneity induced by shrub encroachment into grasslands. Though fires are known to increase aeolian sediment transport, accompanying changes in the sources and sinks of wind-borne sediments may influence biogeochemical cycling and land degradation dynamics. Furthermore, our experiment demonstrated that REEs can be used as reliable tracers for field-scale aeolian studies.

Journal of Geophysical Research: Biogeosciences

On factors influencing air-water gas exchange in emergent wetlands

Knowledge of gas exchange in wetlands is important in order to determine fluxes of climatically and biogeochemically important trace gases and to conduct mass balances for metabolism studies. Very few studies have been conducted to quantify gas transfer velocities in wetlands, and many wind speed/gas exchange parameterizations used in oceanographic or limnological settings are inappropriate under conditions found in wetlands. Here six measurements of gas transfer velocities are made with SF 6 tracer release experiments in three different years in the Everglades, a subtropical peatland with surface water flowing through emergent vegetation. The experiments were conducted under different flow conditions and with different amounts of emergent vegetation to determine the influence of wind, rain, water flow, waterside thermal convection, and vegetation on air-water gas exchange in wetlands. Measured gas transfer velocities under the different conditions ranged from 1.1 cm h −1 during baseline conditions to 3.2 cm h −1 when rain and water flow rates were high. Commonly used wind speed/gas exchange relationships would overestimate the gas transfer velocity by a factor of 1.2 to 6.8. Gas exchange due to thermal convection was relatively constant and accounted for 14 to 51% of the total measured gas exchange. Differences in rain and water flow among the different years were responsible for the variability in gas exchange, with flow accounting for 37 to 77% of the gas exchange, and rain responsible for up to 40%.

Florida

Seasonal, tidal, and geomorphic controls on sediment export to Amazon River tidal floodplains

Mainstem–floodplain material exchange in the tidal freshwater reach of ma jor rivers may lead to significant sequestration of riverine sediment, but this zone remains understudied compared to adjacent fluvial and marine environ ments. This knowledge gap prompts investigation of floodplain-incising tidalchannels found along the banks of tidal rivers and their role in facilitating water and suspended-sediment fluxes between mainstem and floodplain. To evaluate this role, and how it evolves along the tidal river and with time, we measured water level, flow velocity, temperature, and suspended-sediment concentration (SSC) in four tidal channels along the tidal Amazon River, Brazil. Eleven deployments were made during low, rising, high, and falling seasonal Amazon discharge. Generally, channels export high-SSC water from the mainstem to the tidal floodplain on flood tides and transfer low-SSC water back to the mainstem on ebbs. Along the length of the tidal river, the interaction between tidal and seasonal water-level variations and channel–floodplain morphology is a primary control on tidal-channel sediment dynamics. Close to the river mouth, where tides are large, this interaction produces transient flow features and current induced sediment resuspension, but the importance of these processes decreases with distance upstream. Although the magnitude of the exchange of water and sediment between mainstem and floodplain via tidal channels is a small percentage of the total mainstem discharge in this large tidal-river system, tidal channels are important conduits for material flux between these two environments. This flux is critical to resisting floodplain submergence during times of rising sea level.

Amazon River

Missing great earthquakes

The occurrence of three earthquakes with moment magnitude (Mw) greater than 8.8 and six earthquakes larger than Mw 8.5, since 2004, has raised interest in the long-term global rate of great earthquakes. Past studies have focused on the analysis of earthquakes since 1900, which roughly marks the start of the instrumental era in seismology. Before this time, the catalog is less complete and magnitude estimates are more uncertain. Yet substantial information is available for earthquakes before 1900, and the catalog of historical events is being used increasingly to improve hazard assessment. Here I consider the catalog of historical earthquakes and show that approximately half of all Mw ≥ 8.5 earthquakes are likely missing or underestimated in the 19th century. I further present a reconsideration of the felt effects of the 8 February 1843, Lesser Antilles earthquake, including a first thorough assessment of felt reports from the United States, and show it is an example of a known historical earthquake that was significantly larger than initially estimated. The results suggest that incorporation of best available catalogs of historical earthquakes will likely lead to a significant underestimation of seismic hazard and/or the maximum possible magnitude in many regions, including parts of the Caribbean.

Journal of Geophysical Research B: Solid Earth

The utility of atmospheric analyses for the mitigation of artifacts in InSAR

The numerical weather models (NWMs) developed by the meteorological community are able to provide accurate analyses of the current state of the atmosphere in addition to the predictions of the future state. To date, most attempts to apply the NWMs to estimate the refractivity of the atmosphere at the time of satellite synthetic aperture radar (SAR) data acquisitions have relied on predictive models. We test the hypothesis that performing a final assimilative routine, ingesting all available meteorological observations for the times of SAR acquisitions, and generating customized analyses of the atmosphere at those times will better mitigate atmospheric artifacts in differential interferograms. We find that, for our study area around Mount St. Helens (Amboy, Washington, USA), this approach is unable to model the refractive changes and provides no mean benefit for interferogram analysis. The performance is improved slightly by ingesting atmospheric delay estimates derived from the limited local GPS network; however, the addition of water vapor products from the GOES satellites reduces the quality of the corrections. We interpret our results to indicate that, even with this advanced approach, NWMs are not a reliable mitigation technique for regions such as Mount St. Helens with highly variable moisture fields and complex topography and atmospheric dynamics. It is possible, however, that the addition of more spatially dense meteorological data to constrain the analyses might significantly improve the performance of weather modeling of atmospheric artifacts in satellite radar interferograms.

Journal of Geophysical Research B: Solid Earth

Heterogeneous rupture in the great Cascadia earthquake of 1700 inferred from coastal subsidence estimates

Past earthquake rupture models used to explain paleoseismic estimates of coastal subsidence during the great A.D. 1700 Cascadia earthquake have assumed a uniform slip distribution along the megathrust. Here we infer heterogeneous slip for the Cascadia margin in A.D. 1700 that is analogous to slip distributions during instrumentally recorded great subduction earthquakes worldwide. The assumption of uniform distribution in previous rupture models was due partly to the large uncertainties of then available paleoseismic data used to constrain the models. In this work, we use more precise estimates of subsidence in 1700 from detailed tidal microfossil studies. We develop a 3-D elastic dislocation model that allows the slip to vary both along strike and in the dip direction. Despite uncertainties in the updip and downdip slip extensions, the more precise subsidence estimates are best explained by a model with along-strike slip heterogeneity, with multiple patches of high-moment release separated by areas of low-moment release. For example, in A.D. 1700, there was very little slip near Alsea Bay, Oregon (~44.4°N), an area that coincides with a segment boundary previously suggested on the basis of gravity anomalies. A probable subducting seamount in this area may be responsible for impeding rupture during great earthquakes. Our results highlight the need for more precise, high-quality estimates of subsidence or uplift during prehistoric earthquakes from the coasts of southern British Columbia, northern Washington (north of 47°N), southernmost Oregon, and northern California (south of 43°N), where slip distributions of prehistoric earthquakes are poorly constrained.

British Columbia;Washington;Oregon;California

Evolution of dike opening during the March 2011 Kamoamoa fissure eruption, Kīlauea Volcano, Hawai`i

The 5–9 March 2011 Kamoamoa fissure eruption along the east rift zone of Kīlauea Volcano, Hawai`i, followed months of pronounced inflation at Kīlauea summit. We examine dike opening during and after the eruption using a comprehensive interferometric synthetic aperture radar (InSAR) data set in combination with continuous GPS data. We solve for distributed dike displacements using a whole Kīlauea model with dilating rift zones and possibly a deep décollement. Modeled surface dike opening increased from nearly 1.5 m to over 2.8 m from the first day to the end of the eruption, in agreement with field observations of surface fracturing. Surface dike opening ceased following the eruption, but subsurface opening in the dike continued into May 2011. Dike volumes increased from 15, to 16, to 21 million cubic meters (MCM) after the first day, eruption end, and 2 months following, respectively. Dike shape is distinctive, with a main limb plunging from the surface to 2–3 km depth in the up‐rift direction toward Kīlauea's summit, and a lesser projection extending in the down‐rift direction toward Pu`u `Ō`ō at 2 km depth. Volume losses beneath Kīlauea summit (1.7 MCM) and Pu`u `Ō`ō (5.6 MCM) crater, relative to dike plus erupted volume (18.3 MCM), yield a dike to source volume ratio of 2.5 that is in the range expected for compressible magma without requiring additional sources. Inflation of Kīlauea's summit in the months before the March 2011 eruption suggests that the Kamoamoa eruption resulted from overpressure of the volcano's magmatic system.

Hawai'i

The 2011 Mw 7.1 Van (Eastern Turkey) earthquake

We use interferometric synthetic aperture radar (InSAR), body wave seismology, satellite imagery, and field observations to constrain the fault parameters of the M w 7.1 2011 Van (Eastern Turkey) reverse-slip earthquake, in the Turkish-Iranian plateau. Distributed slip models from elastic dislocation modeling of the InSAR surface displacements from ENVISAT and COSMO-SkyMed interferograms indicate up to 9 m of reverse and oblique slip on a pair of en echelon NW 40 °–54 ° dipping fault planes which have surface extensions projecting to just 10 km north of the city of Van. The slip remained buried and is relatively deep, with a centroid depth of 14 km, and the rupture reaching only within 8–9 km of the surface, consistent with the lack of significant ground rupture. The up-dip extension of this modeled WSW striking fault plane coincides with field observations of weak ground deformation seen on the western of the two fault segments and has a dip consistent with that seen at the surface in fault gouge exposed in Quaternary sediments. No significant coseismic slip is found in the upper 8 km of the crust above the main slip patches, except for a small region on the eastern segment potentially resulting from the M w 5.9 aftershock on the same day. We perform extensive resolution tests on the data to confirm the robustness of the observed slip deficit in the shallow crust. We resolve a steep gradient in displacement at the point where the planes of the two fault segments ends are inferred to abut at depth, possibly exerting some structural control on rupture extent.

Journal of Geophysical Research

InSAR Evidence for an active shallow thrust fault beneath the city of Spokane Washington, USA

In 2001, a nearly five month long sequence of shallow, mostly small magnitude earthquakes occurred beneath the city of Spokane, a city with a population of about 200,000, in the state of Washington. During most of the sequence, the earthquakes were not well located because seismic instrumentation was sparse. Despite poor-quality locations, the earthquake hypocenters were likely very shallow, because residents near the city center both heard and felt many of the earthquakes. The combination of poor earthquake locations and a lack of known surface faults with recent movement make assessing the seismic hazards related to the earthquake swarm difficult. However, the potential for destruction from a shallow moderate-sized earthquake is high, for example Christchurch New Zealand in 2011, so assessing the hazard potential of a seismic structure involved in the Spokane earthquake sequence is important. Using interferometric synthetic aperture radar (InSAR) data from the European Space Agency ERS2 and ENVISAT satellites and the Canadian Space Agency RADARSAT-1, satellite we are able to show that slip on a shallow previously unknown thrust fault, which we name the Spokane Fault, is the source of the earthquake sequence. The part of the Spokane Fault that slipped during the 2001 earthquake sequence underlies the north part of the city, and slip on the fault was concentrated between ~0.3 and 2 km depth. Projecting the buried fault plane to the surface gives a possible surface trace for the Spokane Fault that strikes northeast from the city center into north Spokane.

Washington

Chemical controls on fault behavior: weakening of serpentinite sheared against quartz-bearing rocks and its significance for fault creep in the San Andreas system

The serpentinized ultramafic rocks found in many plate-tectonic settings commonly are juxtaposed against crustal rocks along faults, and the chemical contrast between the rock types potentially could influence the mechanical behavior of such faults. To investigate this possibility, we conducted triaxial experiments under hydrothermal conditions (200-350°C), shearing serpentinite gouge between forcing blocks of granite or quartzite. In an ultramafic chemical environment, the coefficient of friction, µ, of lizardite and antigorite serpentinite is 0.5-0.6, and µ increases with increasing temperature over the tested range. However, when either lizardite or antigorite serpentinite is sheared against granite or quartzite, strength is reduced to µ ~ 0.3, with the greatest strength reductions at the highest temperatures (temperature weakening) and slowest shearing rates (velocity strengthening). The weakening is attributed to a solution-transfer process that is promoted by the enhanced solubility of serpentine in pore fluids whose chemistry has been modified by interaction with the quartzose wall rocks. The operation of this process will promote aseismic slip (creep) along serpentinite-bearing crustal faults at otherwise seismogenic depths. During short-term experiments serpentine minerals reprecipitate in low-stress areas, whereas in longer experiments new Mg-rich phyllosilicates crystallize in response to metasomatic exchanges across the serpentinite-crustal rock contact. Long-term shear of serpentinite against crustal rocks will cause the metasomatic mineral assemblages, which may include extremely weak minerals such as saponite or talc, to play an increasingly important role in the mechanical behavior of the fault. Our results may explain the distribution of creep on faults in the San Andreas system.

California

Bayesian inversion of data from effusive volcanic eruptions using physics-based models: Application to Mount St. Helens 2004--2008

Physics-based models of volcanic eruptions can directly link magmatic processes with diverse, time-varying geophysical observations, and when used in an inverse procedure make it possible to bring all available information to bear on estimating properties of the volcanic system. We develop a technique for inverting geodetic, extrusive flux, and other types of data using a physics-based model of an effusive silicic volcanic eruption to estimate the geometry, pressure, depth, and volatile content of a magma chamber, and properties of the conduit linking the chamber to the surface. A Bayesian inverse formulation makes it possible to easily incorporate independent information into the inversion, such as petrologic estimates of melt water content, and yields probabilistic estimates for model parameters and other properties of the volcano. Probability distributions are sampled using a Markov-Chain Monte Carlo algorithm. We apply the technique using GPS and extrusion data from the 2004–2008 eruption of Mount St. Helens. In contrast to more traditional inversions such as those involving geodetic data alone in combination with kinematic forward models, this technique is able to provide constraint on properties of the magma, including its volatile content, and on the absolute volume and pressure of the magma chamber. Results suggest a large chamber of >40 km 3 with a centroid depth of 11–18 km and a dissolved water content at the top of the chamber of 2.6–4.9 wt%.

Washington

Comment on “Historical perspective on seismic hazard to Hispaniola and the northeast Caribbean region” by U. ten Brink et al.

The analysis of historical earthquakes in the northeastern Caribbean by ten Brink et al. [2011, hereafter TB11] addresses the occurrence of large and destructive historical earthquakes associated with the North American-Caribbean plate boundary. One conclusion presented in TB11 is that the recurrence interval for large earthquakes on the left-lateral, strike-slip Septentrional Fault (SF) (Figure 1a) is approximately 300 years. Their Figure 7 shows rupture of the SF across the entire island of Hispaniola in CE 1200, 1542, and 1842. Our comment challenges this model for SF earthquake recurrence because it is inconsistent with our published paleoseismic data that show no large historical earthquake is associated with surface rupture along the SF east of Santiago (Figure 1a)[Prentice et al., 1993; Mann et al., 1998; Prentice et al., 2003].

Hispaniola

The role of viscous magma mush spreading in volcanic flank motion at Kīlauea Volcano, Hawai‘i

Multiple mechanisms have been suggested to explain seaward motion of the south flank of Kīlauea Volcano, Hawai‘i. The consistency of flank motion during both waxing and waning magmatic activity at Kīlauea suggests that a continuously acting force, like gravity body force, plays a substantial role. Using finite element models, we test whether gravity is the principal driver of long-term motion of Kīlauea's flank. We compare our model results to geodetic data from Global Positioning System and interferometric synthetic aperture radar during a time period with few magmatic and tectonic events (2000-2003), when deformation of Kīlauea was dominated by summit subsidence and seaward motion of the south flank. We find that gravity-only models can reproduce the horizontal surface velocities if we incorporate a regional décollement fault and a deep, low-viscosity magma mush zone. To obtain quasi steady state horizontal surface velocities that explain the long-term seaward motion of the flank, we find that an additional weak zone is needed, which is an extensional rift zone above the magma mush. The spreading rate in our model is mainly controlled by the magma mush viscosity, while its density plays a less significant role. We find that a viscosity of 2.5 × 1017–2.5 × 1019 Pa s for the magma mush provides an acceptable fit to the observed horizontal surface deformation. Using high magma mush viscosities, such as 2.5 × 1019 Pa s, the deformation rates remain more steady state over longer time scales. These models explain a significant amount of the observed subsidence at Kīlauea's summit. Some of the remaining subsidence is probably a result of magma withdrawal from subsurface reservoirs.

Hawai'i

A domain decomposition approach to implementing fault slip in finite-element models of quasi-static and dynamic crustal deformation

We employ a domain decomposition approach with Lagrange multipliers to implement fault slip in a finite-element code, PyLith, for use in both quasi-static and dynamic crustal deformation applications. This integrated approach to solving both quasi-static and dynamic simulations leverages common finite-element data structures and implementations of various boundary conditions, discretization schemes, and bulk and fault rheologies. We have developed a custom preconditioner for the Lagrange multiplier portion of the system of equations that provides excellent scalability with problem size compared to conventional additive Schwarz methods. We demonstrate application of this approach using benchmarks for both quasi-static viscoelastic deformation and dynamic spontaneous rupture propagation that verify the numerical implementation in PyLith.

Journal of Geophysical Research B: Solid Earth

Foreshocks during the nucleation of stick-slip instability

We report on laboratory experiments which investigate interactions between aseismic slip, stress changes, and seismicity on a critically stressed fault during the nucleation of stick-slip instability. We monitor quasi-static and dynamic changes in local shear stress and fault slip with arrays of gages deployed along a simulated strike-slip fault (2 m long and 0.4 m deep) in a saw cut sample of Sierra White granite. With 14 piezoelectric sensors, we simultaneously monitor seismic signals produced during the nucleation phase and subsequent dynamic rupture. We observe localized aseismic fault slip in an approximately meter-sized zone in the center of the fault, while the ends of the fault remain locked. Clusters of high-frequency foreshocks (M w ~ −6.5 to −5.0) can occur in this slowly slipping zone 5–50 ms prior to the initiation of dynamic rupture; their occurrence appears to be dependent on the rate at which local shear stress is applied to the fault. The meter-sized nucleation zone is generally consistent with theoretical estimates, but source radii of the foreshocks (2 to 70 mm) are 1 to 2 orders of magnitude smaller than the theoretical minimum length scale over which earthquake nucleation can occur. We propose that frictional stability and the transition between seismic and aseismic slip are modulated by local stressing rate and that fault sections, which would typically slip aseismically, may radiate seismic waves if they are rapidly stressed. Fault behavior of this type may provide physical insight into the mechanics of foreshocks, tremor, repeating earthquake sequences, and a minimum earthquake source dimension.

Journal of Geophysical Research B: Solid Earth