SEARCH · Search USGS

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 199 records · Page 11Linked to original sources

USGS2016

Variation of energy and carbon fluxes from a restored temperate freshwater wetland and implications for carbon market verification protocols

Temperate freshwater wetlands are among the most productive terrestrial ecosystems, stimulating interest in using restored wetlands as biological carbon sequestration projects for greenhouse gas reduction programs. In this study, we used the eddy covariance technique to measure surface energy carbon fluxes from a constructed, impounded freshwater wetland during two annual periods that were 8 years apart: 2002–2003 and 2010–2011. During 2010–2011, we measured methane (CH 4 ) fluxes to quantify the annual atmospheric carbon mass balance and its concomitant influence on global warming potential (GWP). Peak growing season fluxes of latent heat and carbon dioxide (CO 2 ) were greater in 2002–2003 compared to 2010–2011. In 2002, the daily net ecosystem exchange reached as low as −10.6 g C m −2  d −1 , which was greater than 3 times the magnitude observed in 2010 (−2.9 g C m −2  d −1 ). CH 4 fluxes during 2010–2011 were positive throughout the year and followed a strong seasonal pattern, ranging from 38.1 mg C m −2  d −1 in the winter to 375.9 mg C m −2  d −1 during the summer. The results of this study suggest that the wetland had reduced gross ecosystem productivity in 2010–2011, likely due to the increase in dead plant biomass (standing litter) that inhibited the generation of new vegetation growth. In 2010–2011, there was a net positive GWP (675.3 g C m −2  yr −1 ), and when these values are evaluated as a sustained flux, the wetland will not reach radiative balance even after 500 years.

California↗

USGS2015

Incorporating microbial dormancy dynamics into soil decomposition models to improve quantification of soil carbon dynamics of northern temperate forests

Soil carbon dynamics of terrestrial ecosystems play a significant role in the global carbon cycle. Microbial-based decomposition models have seen much growth recently for quantifying this role, yet dormancy as a common strategy used by microorganisms has not usually been represented and tested in these models against field observations. Here we developed an explicit microbial-enzyme decomposition model and examined model performance with and without representation of microbial dormancy at six temperate forest sites of different forest types. We then extrapolated the model to global temperate forest ecosystems to investigate biogeochemical controls on soil heterotrophic respiration and microbial dormancy dynamics at different temporal-spatial scales. The dormancy model consistently produced better match with field-observed heterotrophic soil CO 2 efflux ( R H ) than the no dormancy model. Our regional modeling results further indicated that models with dormancy were able to produce more realistic magnitude of microbial biomass (<2% of soil organic carbon) and soil R H (7.5 ± 2.4 Pg C yr −1 ). Spatial correlation analysis showed that soil organic carbon content was the dominating factor (correlation coefficient = 0.4–0.6) in the simulated spatial pattern of soil R H with both models. In contrast to strong temporal and local controls of soil temperature and moisture on microbial dormancy, our modeling results showed that soil carbon-to-nitrogen ratio (C:N) was a major regulating factor at regional scales (correlation coefficient = −0.43 to −0.58), indicating scale-dependent biogeochemical controls on microbial dynamics. Our findings suggest that incorporating microbial dormancy could improve the realism of microbial-based decomposition models and enhance the integration of soil experiments and mechanistically based modeling.

Journal of Geophysical Research G: Biogeosciences↗

USGS2015

Abundance, stable isotopic composition, and export fluxes of DOC, POC, and DIC from the Lower Mississippi River during 2006–2008

Sources, abundance, isotopic compositions, and export fluxes of dissolved inorganic carbon (DIC), dissolved and colloidal organic carbon (DOC and COC), and particulate organic carbon (POC), and their response to hydrologic regimes were examined through monthly sampling from the Lower Mississippi River during 2006–2008. DIC was the most abundant carbon species, followed by POC and DOC. Concentration and δ 13 C of DIC decreased with increasing river discharge, while those of DOC remained fairly stable. COC comprised 61 ± 3% of the bulk DOC with similar δ 13 C abundances but higher percentages of hydrophobic organic acids than DOC, suggesting its aromatic and diagenetically younger status. POC showed peak concentrations during medium flooding events and at the rising limb of large flooding events. While δ 13 C-POC increased, δ 15 N of particulate nitrogen decreased with increasing discharge. Overall, the differences in δ 13 C between DOC or DIC and POC show an inverse correlation with river discharge. The higher input of soil organic matter and respired CO 2 during wet seasons was likely the main driver for the convergence of δ 13 C between DIC and DOC or POC, whereas enhanced in situ primary production and respiration during dry seasons might be responsible for their isotopic divergence. Carbon export fluxes from the Mississippi River were estimated to be 13.6 Tg C yr −1 for DIC, 1.88 Tg C yr −1 for DOC, and 2.30 Tg C yr −1 for POC during 2006–2008. The discharge-normalized DIC yield decreased during wet seasons, while those of POC and DOC increased and remained constant, respectively, implying variable responses in carbon export to the increasing discharge.

Louisiana↗

USGS2016

Wetland tree transpiration modified by river-floodplain connectivity

Hydrologic connectivity provisions water and nutrient subsidies to floodplain wetlands and may be particularly important in floodplains with seasonal water deficits through its effects on soil moisture. In this study, we measured sapflow in 26 trees of two dominant floodplain forest species ( Celtis laevigata and Quercus lyrata ) at two hydrologically distinct sites in the lower White River floodplain in Arkansas, USA. Our objective was to investigate how connectivity-driven water table variations affected water use, an indicator of tree function. Meteorological variables (photosynthetically active radiation and vapor pressure deficit) were the dominant controls over water use at both sites; however, water table variations explained some site differences. At the wetter site, highest sapflow rates were during a late-season overbank flooding event, and no flood stress was apparent. At the drier site, sapflow decreased as the water table receded. The late-season flood pulse that resulted in flooding at the wetter site did not affect the water table at the drier site; accordingly, higher water use was not observed at the drier site. The species generally associated with wetter conditions ( Q. lyrata ) was more positively responsive to the flood pulse. Flood water subsidy lengthened the effective growing season, demonstrating ecological implications of hydrologic connectivity for alleviating water deficits that otherwise reduce function in this humid floodplain wetland.

Arkansas↗

USGS2016

Component greenhouse gas fluxes and radiative balance from two deltaic marshes in Louisiana: Pairing chamber techniques and eddy covariance

Coastal marshes take up atmospheric CO 2 while emitting CO 2 , CH 4 , and N 2 O. This ability to sequester carbon (C) is much greater for wetlands on a per-area basis than from most ecosystems, facilitating scientific, political, and economic interest in their value as greenhouse gas sinks. However, the greenhouse gas balance of Gulf of Mexico wetlands is particularly understudied. We describe the net ecosystem exchange (NEE c ) of CO 2 and CH 4 using eddy covariance (EC) in comparison with fluxes of CO 2 , CH 4 , and N 2 O using chambers from brackish and freshwater marshes in Louisiana, USA. From EC, we found that 182 g C m -2 y -1 was lost through NEE c from the brackish marsh. Of this, 11 g C m -2 y -1 resulted from net CH 4 emissions and the remaining 171 g C m -2 y -1 resulted from net CO 2 emissions. In contrast, -290 g C m 2 y -1 was taken up through NEE c by the freshwater marsh, with 47 g C m -2 y -1 emitted as CH 4 and -337 g C m -2 y -1 taken up as CO 2 . From chambers, we discovered that neither site had large fluxes of N 2 O. Sustained-flux greenhouse gas accounting metrics indicated that both marshes had a positive (warming) radiative balance, with the brackish marsh having a substantially greater warming effect than the freshwater marsh. That net respiratory emissions of CO 2 and CH 4 as estimated through chamber techniques were 2-4 times different from emissions estimated through EC requires additional understanding of the artifacts created by different spatial and temporal sampling footprints between techniques.

Journal of Geophysical Research: Biogeosciences↗

USGS2016

The differing biogeochemical and microbial signatures of glaciers and rock glaciers

Glaciers and rock glaciers supply water and bioavailable nutrients to headwater mountain lakes and streams across all regions of the American West. Here we present a comparative study of the metal, nutrient, and microbial characteristics of glacial and rock glacial influence on headwater ecosystems in three mountain ranges of the contiguous U.S.: The Cascade Mountains, Rocky Mountains, and Sierra Nevada. Several meltwater characteristics (water temperature, conductivity, pH, heavy metals, nutrients, complexity of dissolved organic matter (DOM), and bacterial richness and diversity) differed significantly between glacier and rock glacier meltwaters, while other characteristics (Ca 2+ , Fe 3+ , SiO 2 concentrations, reactive nitrogen, and microbial processing of DOM) showed distinct trends between mountain ranges regardless of meltwater source. Some characteristics were affected both by glacier type and mountain range (e.g. temperature, ammonium (NH 4 + ) and nitrate (NO 3 - ) concentrations, bacterial diversity). Due to the ubiquity of rock glaciers and the accelerating loss of the low latitude glaciers our results point to the important and changing influence that these frozen features place on headwater ecosystems.

California, Colorado, Oregon, Washington, Wyoming↗

USGS2016

Climate change and dissolved organic carbon export to the Gulf of Maine

Ongoing climate change is affecting the concentration, export (flux), and timing of dissolved organic carbon (DOC) exported to the Gulf of Maine (GoM) through changes in hydrologic regime. DOC export was calculated for water years 1950 through 2013 for 20 rivers and for water years 1930 through 2013 for 14 rivers draining to the GoM. DOC export was also estimated for the 21st century based on climate and hydrologic modeling in a previously published study. DOC export was calculated by using the regression model LOADEST to fit seasonally adjusted concentration discharge (C-Q) relations. Our results are an analysis of the sensitivity of DOC export to changes in hydrologic conditions over time since land cover and vegetation were held constant over time. Despite large interannual variability, all rivers had increasing DOC export during winter and these trends were significant ( p < 0.05) in 10 out of 20 rivers for 1950 to 2013 and in 13 out of 14 rivers for 1930 to 2013. All rivers also had increasing annual export of DOC although fewer trends were statistically significant than for winter export. Projections for DOC export during the 21st century were variable depending on the climate model and greenhouse gas emission scenario that affected future river discharge through effects on precipitation and evapotranspiration. The most consistent result was a significant increase in DOC export in winter in all model-by-emission scenarios. DOC export was projected to decrease during the summer in all model-by-emission scenarios, with statistically significant decreases in half of the scenarios.

Maine↗

USGS2016

Induced earthquake magnitudes are as large as (statistically) expected

A major question for the hazard posed by injection-induced seismicity is how large induced earthquakes can be. Are their maximum magnitudes determined by injection parameters or by tectonics? Deterministic limits on induced earthquake magnitudes have been proposed based on the size of the reservoir or the volume of fluid injected. However, if induced earthquakes occur on tectonic faults oriented favorably with respect to the tectonic stress field, then they may be limited only by the regional tectonics and connectivity of the fault network. In this study, we show that the largest magnitudes observed at fluid injection sites are consistent with the sampling statistics of the Gutenberg-Richter distribution for tectonic earthquakes, assuming no upper magnitude bound. The data pass three specific tests: (1) the largest observed earthquake at each site scales with the log of the total number of induced earthquakes, (2) the order of occurrence of the largest event is random within the induced sequence, and (3) the injected volume controls the total number of earthquakes rather than the total seismic moment. All three tests point to an injection control on earthquake nucleation but a tectonic control on earthquake magnitude. Given that the largest observed earthquakes are exactly as large as expected from the sampling statistics, we should not conclude that these are the largest earthquakes possible. Instead, the results imply that induced earthquake magnitudes should be treated with the same maximum magnitude bound that is currently used to treat seismic hazard from tectonic earthquakes.

Journal of Geophysical Research B: Solid Earth↗

USGS2016

Subsidence rates at the southern Salton Sea consistent with reservoir depletion

Space geodetic measurements from the Envisat satellite between 2003 and 2010 show that subsidence rates near the southeastern shoreline of the Salton Sea in Southern California are up to 52mmyr −1 greater than the far-field background rate. By comparing these measurements with model predictions, we find that this subsidence appears to be dominated by poroelastic contraction associated with ongoing geothermal fluid production, rather than the purely fault-related subsidence proposed previously. Using a simple point source model, we suggest that the source of this proposed volumetric strain is at depths between 1.0 km and 2.4 km (95% confidence interval), comparable to generalized boundaries of the Salton Sea geothermal reservoir. We find that fault slip on two previously imaged tectonic structures, which are part of a larger system of faults in the Brawley Seismic Zone, is not an adequate predictor of surface velocity fields because the magnitudes of the best fitting slip rates are often greater than the full plate boundary rate and at least 2 times greater than characteristic sedimentation rates in this region. Large-scale residual velocity anomalies indicate that spatial patterns predicted by fault slip are incompatible with the observations.

California↗

USGS2016

A large refined catalog of earthquake relocations and focal mechanisms for the Island of Hawai'i and its seismotectonic implications

We present high-quality focal mechanisms based on a refined earthquake location catalog for the Island of Hawai'i, focusing on Mauna Loa and Kīlauea volcanoes. The relocation catalog is based on first-arrival times and waveform data of both compressional and shear waves for about 180,000 events on and near the Island of Hawai'i between 1986 and 2009 recorded by the seismic stations at the Hawaiian Volcano Observatory. We relocate all the earthquakes by applying ray tracing through an existing three-dimensional velocity model, similar event cluster analysis, and a differential-time relocation method. The resulting location catalog represents an expansion of previous relocation studies, covering a longer time period and consisting of more events with well-constrained absolute locations. The focal mechanisms are obtained based on the compressional-wave first-motion polarities and compressional-to-shear wave amplitude ratios by applying the HASH program to the waveform cross correlation relocated earthquakes. Overall, the good-quality (defined by the HASH parameters) focal solutions are dominated by normal faulting in our study area, especially in the active Ka'ōiki and Hīlea seismic zones. Kīlauea caldera is characterized by a mixture of approximately equal numbers of normal, strike-slip, and reverse faults, whereas its south flank has slightly fewer strike-slip events. Our relocation and focal mechanism results will be useful for mapping the seismic stress and strain fields and for understanding the seismic-volcanic-tectonic relationships within the magmatic systems.

Hawai'i↗

USGS2016

Volcano dome dynamics at Mount St. Helens: Deformation and intermittent subsidence monitored by seismicity and camera imagery pixel offsets

The surface deformation field measured at volcanic domes provides insights into the effects of magmatic processes, gravity- and gas-driven processes, and the development and distribution of internal dome structures. Here we study short-term dome deformation associated with earthquakes at Mount St. Helens, recorded by a permanent optical camera and seismic monitoring network. We use Digital Image Correlation (DIC) to compute the displacement field between successive images and compare the results to the occurrence and characteristics of seismic events during a 6 week period of dome growth in 2006. The results reveal that dome growth at Mount St. Helens was repeatedly interrupted by short-term meter-scale downward displacements at the dome surface, which were associated in time with low-frequency, large-magnitude seismic events followed by a tremor-like signal. The tremor was only recorded by the seismic stations closest to the dome. We find a correlation between the magnitudes of the camera-derived displacements and the spectral amplitudes of the associated tremor. We use the DIC results from two cameras and a high-resolution topographic model to derive full 3-D displacement maps, which reveals internal dome structures and the effect of the seismic activity on daily surface velocities. We postulate that the tremor is recording the gravity-driven response of the upper dome due to mechanical collapse or depressurization and fault-controlled slumping. Our results highlight the different scales and structural expressions during growth and disintegration of lava domes and the relationships between seismic and deformation signals.

Washington↗

USGS2016

Insights into shallow magmatic processes at Kīlauea Volcano, Hawaiʻi, from a multiyear continuous gravity time series

Continuous gravity data collected near the summit eruptive vent at Kīlauea Volcano, Hawaiʻi, during 2011–2015 show a strong correlation with summit-area surface deformation and the level of the lava lake within the vent over periods of days to weeks, suggesting that changes in gravity reflect variations in volcanic activity. Joint analysis of gravity and lava level time series data indicates that over the entire time period studied, the average density of the lava within the upper tens to hundreds of meters of the summit eruptive vent remained low—approximately 1000–1500 kg/m 3 . The ratio of gravity change (adjusted for Earth tides and instrumental drift) to lava level change measured over 15 day windows rose gradually over the course of 2011–2015, probably reflecting either (1) a small increase in the density of lava within the eruptive vent or (2) an increase in the volume of lava within the vent due to gradual vent enlargement. Superimposed on the overall time series were transient spikes of mass change associated with inflation and deflation of Kīlauea's summit and coincident changes in lava level. The unexpectedly strong mass variations during these episodes suggest magma flux to and from the shallow magmatic system without commensurate deformation, perhaps indicating magma accumulation within, and withdrawal from, void space—a process that might not otherwise be apparent from lava level and deformation data alone. Continuous gravity data thus provide unique insights into magmatic processes, arguing for continued application of the method at other frequently active volcanoes.

Hawaii↗

USGS2016

Coseismic slip and early afterslip of the 2015 Illapel, Chile, earthquake: Implications for frictional heterogeneity and coastal uplift

Great subduction earthquakes are thought to rupture portions of the megathrust, where interseismic coupling is high and velocity-weakening frictional behavior is dominant, releasing elastic deformation accrued over a seismic cycle. Conversely, postseismic afterslip is assumed to occur primarily in regions of velocity-strengthening frictional characteristics that may correlate with lower interseismic coupling. However, it remains unclear if fixed frictional properties of the subduction interface, coseismic or aftershock-induced stress redistribution, or other factors control the spatial distribution of afterslip. Here we use interferometric synthetic aperture radar and Global Position System observations to map the distribution of coseismic slip of the 2015 M w 8.3 Illapel, Chile, earthquake and afterslip within the first 38 days following the earthquake. We find that afterslip overlaps the coseismic slip area and propagates along-strike into regions of both high and moderate interseismic coupling. The significance of these observations, however, is tempered by the limited resolution of geodetic inversions for both slip and coupling. Additional afterslip imaged deeper on the fault surface bounds a discrete region of deep coseismic slip, and both contribute to net uplift of the Chilean Coastal Cordillera. A simple partitioning of the subduction interface into regions of fixed frictional properties cannot reconcile our geodetic observations. Instead, stress heterogeneities, either preexisting or induced by the earthquake, likely provide the primary control on the afterslip distribution for this subduction zone earthquake. We also explore the occurrence of coseismic and postseismic coastal uplift in this sequence and its implications for recent hypotheses concerning the source of permanent coastal uplift along subduction zones.

Journal of Geophysical Research B: Solid Earth↗

USGS2016

Continental rupture and the creation of new crust in the Salton Trough rift, southern California and northern Mexico: Results from the Salton Seismic Imaging Project

A refraction and wide-angle reflection seismic profile along the axis of the Salton Trough, California and Mexico, was analyzed to constrain crustal and upper mantle seismic velocity structure during active continental rifting. From the northern Salton Sea to the southern Imperial Valley, the crust is 17-18 km thick and approximately one-dimensional. The transition at depth from Colorado River sediment to underlying crystalline rock is gradual and is not a depositional surface. The crystalline rock from ~3 to ~8 km depth is interpreted as sediment metamorphosed by high heat flow. Deeper felsic crystalline rock could be stretched pre-existing crust or higher grade metamorphosed sediment. The lower crust below ~12 km depth is interpreted to be gabbro emplaced by rift-related magmatic intrusion by underplating. Low upper-mantle velocity indicates high temperature and partial melting. Under the Coachella Valley, sediment thins to the north and the underlying crystalline rock is interpreted as granitic basement. Mafic rock does not exist at 12-18 depth as it does to the south, and a weak reflection suggests Moho at ~28 km depth. Structure in adjacent Mexico has slower mid-crustal velocity and rocks with mantle velocity must be much deeper than in the Imperial Valley. Slower velocity and thicker crust in the Coachella and Mexicali valleys define the rift zone between them to be >100 km wide in the direction of plate motion. North American lithosphere in the central Salton Trough has been rifted apart and is being replaced by new crust created by magmatism, sedimentation, and metamorphism.

Journal of Geophysical Research↗

USGS2016

Origin of the pulse-like signature of shallow long-period volcano seismicity

Short-duration, pulse-like long-period (LP) events are a characteristic type of seismicity accompanying eruptive activity at Mount Etna in Italy in 2004 and 2008 and at Turrialba Volcano in Costa Rica and Ubinas Volcano in Peru in 2009. We use the discrete wave number method to compute the free surface response in the near field of a rectangular tensile crack embedded in a homogeneous elastic half space and to gain insights into the origin of the LP pulses. Two source models are considered, including (1) a vertical fluid-driven crack and (2) a unilateral tensile rupture growing at a fixed sub-Rayleigh velocity with constant opening on a vertical crack. We apply cross correlation to the synthetics and data to demonstrate that a fluid-driven crack provides a natural explanation for these data with realistic source sizes and fluid properties. Our modeling points to shallow sources (<1 km depth), whose signatures are representative of the Rayleigh pulse sampled at epicentral distances >∼1 km. While a slow-rupture failure provides another potential model for these events, the synthetics and resulting fits to the data are not optimal in this model compared to a fluid-driven source. We infer that pulse-like LP signatures are parts of the continuum of responses produced by shallow fluid-driven sources in volcanoes.

Journal of Geophysical Research B: Solid Earth↗

USGS2017

Low stress drops observed for aftershocks of the 2011 Mw 5.7 Prague, Oklahoma, earthquake

In November 2011, three M w ≥ 4.8 earthquakes and thousands of aftershocks occurred along the structurally complex Wilzetta fault system near Prague, Oklahoma. Previous studies suggest that wastewater injection induced a M w 4.8 foreshock, which subsequently triggered a M w 5.7 mainshock. We examine source properties of aftershocks with a standard Brune-type spectral model and jointly solve for seismic moment ( M 0 ), corner frequency ( f 0 ), and kappa ( κ ) with an iterative Gauss-Newton global downhill optimization method. We examine 934 earthquakes with initial moment magnitudes ( M w ) between 0.33 and 4.99 based on the pseudospectral acceleration and recover reasonable M 0 , f 0 , and κ for 87 earthquakes with M w 1.83–3.51 determined by spectral fit. We use M 0 and f 0 to estimate the Brune-type stress drop, assuming a circular fault and shear-wave velocity at the hypocentral depth of the event. Our observations suggest that stress drops range between 0.005 and 4.8 MPa with a median of 0.2 MPa (0.03–26.4 MPa with a median of 1.1 MPa for Madariaga-type), which is significantly lower than typical eastern United States intraplate events (>10 MPa). We find that stress drops correlate weakly with hypocentral depth and magnitude. Additionally, we find the stress drops increase with time after the mainshock, although temporal variation in stress drop is difficult to separate from spatial heterogeneity and changing event locations. The overall low median stress drop suggests that the fault segments may have been primed to fail as a result of high pore fluid pressures, likely related to nearby wastewater injection.

Oklahoma↗

USGS2016

The role of crystallization-driven exsolution on the sulfur mass balance in volcanic arc magmas

The release of large amounts of sulfur to the stratosphere during explosive eruptions affects the radiative balance in the atmosphere and consequentially impacts climate for up to several years after the event. Quantitative estimations of the processes that control the mass balance of sulfur between melt, crystals, and vapor bubbles is needed to better understand the potential sulfur yield of individual eruption events and the conditions that favor large sulfur outputs to the atmosphere. The processes that control sulfur partitioning in magmas are (1) exsolution of volatiles (dominantly H 2 O) during decompression (first boiling) and during isobaric crystallization (second boiling), (2) the crystallization and breakdown of sulfide or sulfate phases in the magma, and (3) the transport of sulfur-rich vapor (gas influx) from deeper unerupted regions of the magma reservoir. Vapor exsolution and the formation/breakdown of sulfur-rich phases can all be considered as closed-system processes where mass balance arguments are generally easier to constrain, whereas the contribution of sulfur by vapor transport (open system process) is more difficult to quantify. The ubiquitous “excess sulfur” problem, which refers to the much higher sulfur mass released during eruptions than what can be accounted for by amount of sulfur originally dissolved in erupted melt, as estimated from melt inclusion sulfur concentrations (the “petrologic estimate”), reflects the challenges in closing the sulfur mass balance between crystals, melt, and vapor before and during a volcanic eruption. In this work, we try to quantify the relative importance of closed- and open-system processes for silicic arc volcanoes using kinetic models of sulfur partitioning during exsolution. Our calculations show that crystallization-induced exsolution (second boiling) can generate a significant fraction of the excess sulfur observed in crystal-rich arc magmas. This result does not negate the important role of vapor migration in sulfur mass balance but rather points out that second boiling (in situ exsolution) can provide the necessary yield to drive the excess sulfur to the levels observed for crystal-rich systems. In contrast, in crystal-poor systems, magma recharge that releases sulfur-rich bubbles is necessary and most likely the primary contributor to sulfur mass balance. Finally, we apply our model to account for the effect of sulfur partitioning during second boiling and its impact on sulfur released during the Cerro Galan supereruption in Argentina (2.08 Ma) and show the potential importance of second boiling in releasing a large amount of sulfur to the atmosphere during the eruption of large crystal-rich ignimbrites.

Journal of Geophysical Research B: Solid Earth↗

USGS2017

Reexamination of the magnitudes for the 1906 and 1922 Chilean earthquakes using Japanese tsunami amplitudes: Implications for source depth constraints

Far-field tsunami records from the Japanese tide gauge network allow the reexamination of the moment magnitudes ( M w ) for the 1906 and 1922 Chilean earthquakes, which to date rely on limited information mainly from seismological observations alone. Tide gauges along the Japanese coast provide extensive records of tsunamis triggered by six great ( M w >8) Chilean earthquakes with instrumentally determined moment magnitudes. These tsunami records are used to explore the dependence of tsunami amplitudes in Japan on the parent earthquake magnitude of Chilean origin. Using the resulting regression parameters together with tide gauge amplitudes measured in Japan we estimate apparent moment magnitudes of M w 8.0–8.2 and M w 8.5–8.6 for the 1906 central and 1922 north-central Chile earthquakes. The large discrepancy of the 1906 magnitude estimated from the tsunami observed in Japan as compared with those previously determined from seismic waves ( M s 8.4) suggests a deeper than average source with reduced tsunami excitation. A deep dislocation along the Chilean megathrust would favor uplift of the coast rather than beneath the sea, giving rise to a smaller tsunami and producing effects consistent with those observed in 1906. The 1922 magnitude inferred from far-field tsunami amplitudes appear to better explain the large extent of damage and the destructive tsunami that were locally observed following the earthquake than the lower seismic magnitudes ( M s 8.3) that were likely affected by the well-known saturation effects. Thus, a repeat of the large 1922 earthquake poses seismic and tsunami hazards in a region identified as a mature seismic gap.

Journal of Geophysical Research B: Solid Earth↗