Search USGS⌕ Search

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.

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

At least 955 records · Page 53Linked to original sources
USGS2019

Interseismic quiescence and triggered slip of active normal faults of Kīlauea Volcano’s south flank during 2001-2018

The mobile south flank of Kīlauea Volcano hosts two normal fault systems, the Koa'e fault system (KFS) and the Hilina fault system (HFS). In historical time, at least three M>6.5 earthquakes have occurred on the basal detachment of the Kīlauea Volcano's south flank, with the most recent being the 4 May 2018 M6.9 earthquake. Here we analyze kinematic Global Positioning System data collected from 2001 to 2017 and interferometric synthetic aperture radar data before, during, and after the 2018 M6.9 earthquake to determine the crustal motion across the HFS and KFS faults. Our results indicate that the HFS faults did not significantly slip during the interseismic period from 2007 to 2011. Despite its substantial magnitude, interferometric synthetic aperture radar (InSAR) data show that the 2018 M6.9 earthquake triggered subcentimeter level slip along sections of the previously mapped HFS branches. Up to 20 cm of offset occurred on what appears to be a newly formed (or previously unknown) fault near the eastern end of the HFS. During the 3 months following the M6.9 earthquake, up to ~30 cm of slip occurred along the KFS, which helps accommodate rapid large‐scale subsidence of Kīlauea's summit region as large volumes of summit reservoir magma fed the lower East Rift Zone eruption. The HFS appears to activate only in concert with large earthquakes on the basal detachment. The KFS, on the other hand, moves both seismically during small local earthquakes and aseismically in response to nearby earthquakes and caldera subsidence.

Hawaii↗
USGS2019

Insights into pāhoehoe lava emplacement using visible and thermal structure-from-motion photogrammetry

We present the evolution over 3 months of a 2016–2017 pāhoehoe flow at Kīlauea as it changed from a narrow sheet flow into a compound lava field fed by a stable system of tubes. The portion of the flow located on Kīlauea's coastal plain was characterized using helicopter‐based visible and thermal structure‐from‐motion photogrammetry to construct a series of georeferenced digital surface models and thermal maps on eight different days. Results reveal key influences on the emplacement and evolution of such long‐lived pāhoehoe flows. This region of the flow grew by ~12 × 10 6 m 3 with a near‐constant time‐average discharge rate of 1.2–2.7 m 3 /s. The development of two tube systems is captured and shows an initial nascent tube enhanced by a narrow topographic confinement, which later inflated and created a topographic inversion that modulated the emplacement of a second flow lobe with its own tube system. The analysis of breakouts at various stages of the field's life suggests that the evolution of the thermal and morphological properties of the flow surface reflect its maturity. Thermal properties of breakouts were used to expand the empirical relationship of breakout cooling to longer timescales. This study contributes to the long‐term development and validation of more accurate predictive models for pāhoehoe, required during the management of long‐lasting lava flow crises in Hawai'i and elsewhere.

Hawaii↗
USGS2019

Slow-growing and extended-duration seismicity swarms: Reactivating joints or foliations in the Cahuilla Valley Pluton, Central Peninsular Ranges, Southern California

Three prolific earthquake swarms and numerous smaller ones have occurred since 1980 in the Mesozoic igneous plutonic rocks of the Perris block of the Peninsular Ranges, Southern California. The major swarms occurred in 1980–1981, 1983–1984, and 2016–2018, with the latest swarm still ongoing. These swarms have no clear mainshock, with the largest events of M L 3.6, M L 3.7, and M w 4.4. Each successive swarm had larger cumulative seismic moment release with about 314 and 411 events of M ≥ 1.5, while the third swarm has produced about 451 events of M ≥ 1.5 (as of September 2018). The concurrent strike‐slip faulting occurred on north to northwest striking planes but with no orthogonal northeast trending seismicity alignments. These shallow swarms are probably driven by intrablock Pacific‐North America plate boundary stress loading of the two bounding major late Quaternary strike‐slip faults, the Elsinore and San Jacinto faults. The state of stress within the Cahuilla Valley pluton has a ~40° angle between the maximum principal stress and the average trend of the swarms, suggesting that migrating pore fluid pressures aid in the formation and growth of zones of weakness. These swarms, which last more than 600 days each, exhibit clear bilateral spatial migration for distances of up to ~7–8 km and reach their full length in about 20 months. The slow spatial‐temporal development of the swarms corresponds to a fluid diffusivity of 0.006 to 0.01 m 2 /s, consistent with very low permeability rocks as expected for this block. There is no geodetic or other evidence for a slow slip event driving the swarms.

California↗
USGS2020

Near-fault velocity spectra from laboratory failures and their relation to natural ground motion

We compared near-fault velocity spectra recorded during laboratory experiments to that of natural earthquakes. We fractured crystalline rock samples at room temperature and intermediate confining pressure (50 MPa). Subsequent slip events were generated on the fracture surfaces under higher confinement (300 MPa). Velocity spectra from rock fracture resemble the inverse frequency (1/ f ) decay of natural earthquake velocity. This spectrum can be attributed to fault creation via seismic fracturing over a wide range of spatial scales. In contrast, subsequent slips on the rough fracture surfaces are depleted in high frequency energy and falloff approximately as 1/ f 2 . The 1/ f 2 spectrum is more consistent with a slider-block model obeying static-kinetic friction than a natural earthquake. The depleted high frequency content precludes the rough fault experiments from being directly analogous to natural sources. The suppression of high frequencies may have resulted from two possible factors: (1) the presence of a well-developed shear zone and coseismic damping of the fault motion by dissipation within it or, in our favored interpretation, (2) a smaller amount of energy dissipated by shearing relative to the total energy release at elevated confining pressure. In context of the latter explanation, a unifying concept that applies to these experiments, earthquakes, ground motion, and models of complex radiated motion is that high frequency radiated energy is relatively enhanced when total energy release is nearly balanced within the source region by dissipative processes. This near-critical energy release condition can be accessed at low normal stress in laboratory experiments.

Journal of Geophysical Research↗
USGS2020

How often can Earthquake Early Warning systems alert sites with high intensity ground motion?

Although numerous Earthquake Early Warning (EEW) algorithms have been developed we still lack a detailed understanding of how often and under what circumstances useful ground motion alerts can be provided to end-users. Here we analyze the alerting performance of the PLUM, EPIC and FinDer algorithms by running them retrospectively on the seismic strong motion data of the 219 earthquakes in Japan since 1996 that exceeded Modified Mercalli Intensity (MMI) of 4.5 on at least 10 sites (Mw 4.5-9.1). Our analysis suggests that, irrespective of the algorithm, EEW end-users should be prepared that EEW can often but not always provide useful ground motion alerts. A majority of sites with moderate-strong ground motion (MMI 5-6) can generally get at least a few seconds of warning time from all algorithms. If such shaking is caused by a shallow crustal event, around 50% of such sites receive alerts with warning times >5 s. Many sites with severe-extreme ground motion (MMI >=8) can be alerted successfully in the case of very large offshore earthquakes, but less than 20% can be alerted ahead of time if such shaking is caused by a shallow crustal event. Our results provide detailed quantitative insight into the expected alerting performance for EEW algorithms under realistic conditions. The main caveat is that the largest shallow crustal event in our data set has Mw7.0, i.e. the data set does not contain very large strike slip events.

Journal of Geophysical Research↗
USGS2019

Crustal magmatism and anisotropy beneath the Arabian Shield - A cautionary tale

Volcanism in Saudi Arabia includes a historic eruption close to the holy city of Al Madinah. As part of a volcanic hazard assessment of this area, magnetotelluric (MT) data were collected to investigate the structural setting, the distribution of melt within the crust, and the mantle source of volcanism. Interpretation of a new 3‐D resistivity model includes a shallow graben beneath thin lava fields (Harrats), a melt‐free upper crust, and decompression melting in the asthenosphere below thin lithosphere. Within the lower crust the model images elongate conductivity anomalies, one of which was attributed in a previous MT study to melt. The regional MT data, combined with perspective from geology and geophysical modeling, suggest the lower crust is anisotropic with no interconnected melt zones. These divergent interpretations have distinct hazard implications and highlight the importance of large survey aperture and anisotropic modeling to MT studies of volcanic regions. Lower‐crustal anisotropy extends beyond the Harrat, with the most conductive direction oriented N10°E and a factor of 3–5, determined from 2‐D anisotropic inversion, between the most and least conductive directions. The enhanced conductivity is likely due to interconnected grain boundary graphite, while the anisotropy direction reflects either frozen‐in fabric from Neoproterozoic stabilization of the Arabian Shield or modern ductile deformation driven by channelized asthenospheric flow coupled through a thin rigid mantle lid. Asthenospheric melt is interpreted to transect the crust primarily through diking, with limited melt storage and short residence times in the crust.

Arabian Shield↗
USGS2020

Submarine landslide kinematics derived from high-resolution imaging in Port Valdez, Alaska

Submarine landslides caused by strong ground shaking during the M9.2 1964 Great Alaska earthquake generated a tsunami that destroyed much of the old town of Valdez, Alaska, and was responsible for 32 deaths at that location. We explore structural details of the 1964 landslide deposit, as well as landslide deposits from earlier events, in order to characterize kinematics of the landslide process. We present a new high‐resolution seismic reflection data set that images the 1964 landslide deposit and six pre‐1964 deposits with great detail. These deposits are represented by thick packages (~7–23 m) of debris within >500 m of fjord sedimentation above basement. Internal slide structures are associated with distinctive landslide failure mechanisms, including detailed erosional and depositional features and structures resolved within both landslide blocks and distal debris flow layers. Based on comparisons of deposit volume from subbottom structure and differenced bathymetry, we refine prior interpretations of the source of failed material. New data show evidence for basal erosion and reworking of fjord‐floor sedimentation. Additionally, material comprising the 1964 landslide appears to have been translated and deformed by lateral thrusting, rather than having been sourced entirely from upslope evacuation zones. Taking into account these complexities in depositional patterns, we show variations in slide size through Holocene time and relate the history of landslides to the paleoseismic record. Collectively, these new observations demonstrate that Port Valdez has a repeated history of large submarine landslides, which are likely associated with large megathrust earthquakes.

Alaska↗
USGS2019

Introduction to special issue on gas hydrate in porous media: Linking laboratory and field‐scale phenomena

The proliferation of drilling expeditions focused on characterizing natural gas hydrate as a potential energy resource has spawned widespread interest in gas hydrate reservoir properties and associated porous media phenomena. Between 2017 and 2019, a Special Section of this journal compiled contributed papers elucidating interactions between gas hydrate and sediment based on laboratory, numerical modeling, and field studies. Motivated mostly by field observations in the northern Gulf of Mexico and offshore Japan, several papers focus on the mechanisms for gas hydrate formation and accumulation, particularly with vapor phase gas, not dissolved gas, as the precursor to hydrate. These studies rely on numerical modeling or laboratory experiments using sediment packs or benchtop micromodels. A second focus of the Special Section is the role of fines in inhibiting production of gas from methane hydrate, controlling the distribution of hydrate at a pore scale, and influencing the bulk behavior of seafloor sediments. Other papers fill knowledge gaps related to the physical properties of hydrate-bearing sediments and advance new approaches in coupled thermal-mechanical modeling of these sediments during hydrate dissociation. Finally, one study addresses the long-standing question about the fate of methane hydrate at the molecular level when CO2 is injected into natural reservoirs under hydrate-forming conditions.

Journal of Geophysical Research B: Solid Earth↗
USGS2020

Magma intrusion and volatile ascent beneath Norris Geyser Basin, Yellowstone National Park

Recent activity has provided new insights into the causes of surface deformation in and around the Yellowstone Caldera, a topic that has been debated since the discovery of caldera floor uplift more than four decades ago. An episode of unusually rapid uplift (>15 cm/yr) centered near Norris Geyser Basin along the north caldera rim began in late 2013 and continued until a M w 4.9 earthquake on 30 March 2014; thereafter, uplift abruptly switched to subsidence. Uplift at rates of several centimeters per year resumed in 2016 and continued at least through the end of 2018. Modeling of Global Positioning System and interferometric synthetic aperture radar data suggests an evolving process of deep magma intrusion during 1996–2001 followed by volatile ascent and accumulation at shallow levels, perhaps as shallow as a few hundred meters depth. The depth of shallow volatile accumulation appears to have shallowed from the 2014 to the 2016 deformation episode, and frequent eruptions of Steamboat Geyser since March 2018 are likely a surface manifestation of this ongoing process. Hydrothermal explosion features are prominent in the Norris Geyser Basin area, and the apparent shallow nature of the volatile accumulation implies an increased risk of hydrothermal explosions.

Wyoming↗
USGS2019

The effect of brine on the electrical properties of methane hydrate

Gas hydrates possess lower electrical conductivity (inverse of resistivity) than either seawater or ice, but higher than clastic silts and sands, such that electromagnetic methods can be employed to help identify their natural formation in marine and permafrost environments. Controlled laboratory studies offer a means to isolate and quantify the effects of changing individual components within gas‐hydrate‐bearing systems, in turn yielding insight into the behavior of natural systems. Here we investigate the electrical properties of polycrystalline methane hydrate with ≥25% gas‐filled porosity and in mixture with brine. Initially, pure methane hydrate was synthesized from H 2 O ice and CH 4 gas while undergoing electrical impedance measurement, then partially dissociated to assess the effects of pure pore water accumulation on electrical conductivity. Methane hydrate + brine mixtures were then formed by either adding NaCl (0.25–2.5 wt %) to high‐purity ice or by using frozen seawater as a reactant. Conductivity was obtained from impedance measurements made in situ throughout synthesis while temperature cycled between +15 °C and −25 °C. Several possible conduction mechanisms were subsequently determined using equivalent circuit modeling. Samples with low NaCl concentration show a doping/impurity effect and a log linear conductivity response as a function of temperature. For higher salt content samples, conductivity increases exponentially with temperature and the log linear relationship no longer holds; instead, we observe phase changes within the samples that follow NaCl–H 2 O–CH 4 phase equilibrium predictions. Final samples were quenched in liquid nitrogen and imaged by cryogenic scanning electron microscopy (cryo‐SEM) to assess grain‐scale characteristics.

Journal of Geophysical Research↗
USGS2019

Coseismic slip and early afterslip of the M6.0 August 24, 2014 South Napa, California, earthquake

We employ strong motion seismograms and static offsets from the Global Positioning System, Interferometric Synthetic Aperture Radar, and other measurements in order to derive a coseismic slip and afterslip model of the M6.0 24 August 2014 South Napa earthquake. This earthquake ruptured an ∼13‐km‐long portion of the West Napa fault with predominantly right‐lateral strike slip. In the kinematic seismic slip inversions, we couple the coseismic slip and afterslip distributions by requiring both distributions to involve right‐lateral strike slip with positive amplitude, with the net static slip being the sum of the two. We consider several candidate fault geometries: a first involving two steeply east dipping fault planes that reach Earth's surface at the western surface trace (STW), where most surface rupture was observed, a second involving a steeply west dipping plane that also reaches Earth's surface at the STW, and a third involving a combination of two variably west dipping planes constrained to pass through the locus of postseismic seismicity located ∼1 km west of the STW. The data are best fit using the model of two east dipping fault planes, with coseismic slip up to ∼1.2 m on a dominant shallow asperity about 10 km north of the hypocenter and on deeper asperities on the southern part of the rupture. Afterslip up to 1 m is concentrated along the southern part of the rupture at depths 5 km, consistent with surface observations of afterslip. Seismic moments associated with coseismic slip and afterslip are 1.13×10 18 N m (Mw 6.00) and 3.64×10 17 N m, respectively.

California↗
USGS2020

Semiautomated estimates of directivity and related source properties of small to moderate southern California earthquakes using second seismic moments

We develop a semiautomated method for estimating with second seismic moments the directivity, rupture area, duration, and centroid velocity of earthquakes. The method is applied to 41 southern California earthquakes with magnitude in the range 3.5–5.2 and provides stable results for 28 events. Apparent source time functions (ASTFs) of P and S phases are derived using deconvolution with three stacked empirical Green's functions (seGf). The use of seGf suppresses nongeneric source effects, improves the focal mechanism correspondence to the analyzed earthquakes, and typically allows inclusion of 5 to 15 more ASTFs compared with analysis using a single eGf. Most analyzed earthquakes in the Trifurcation area of the San Jacinto Fault have directivities toward the northwest, while events around Cajon Pass and San Gabriel Mountain tend to propagate toward the southeast. These results are generally consistent with predictions for dynamic rupture on bimaterial interfaces associated with the imaged velocity contrasts in the area. The second moment inversions also provide constraints on the upper and lower bounds of rupture areas in our data set. Stress drops and uncertainties are estimated for elliptical ruptures using the derived characteristic rupture length and width. The semiautomated second moment method with seGfs can be used for routine application to moderate earthquakes in locations with good station coverage.

California↗
USGS2019

Offset channels may not accurately record strike-slip fault displacement: Evidence from landscape evolution models

Slip distribution, slip rate, and slip per event for strike‐slip faults are commonly determined by correlating offset stream channels—under the assumption that they record seismic slip—but offset channels are formed by the interplay of tectonic and geomorphic processes. To constrain offset channel development under known tectonic and geomorphic conditions, we use numerical landscape evolution simulations along a theoretical strike‐slip fault with uniform and steady uplift, erosion, and diffusion. We investigate the influence of four tectonic parameters (fault zone width, earthquake recurrence interval, variance of the recurrence interval, and total slip relative to channel spacing) on offset channel development through multiple earthquake cycles. Analysis of >3,000 automatically measured offsets from >135 simulations suggests ~30% variability in individual measurements, but modeled displacement is recovered by averaging multiple measurements. However, the average of multiple offset measurements systematically underestimates modeled slip except when the fault zone is less than ~5 m wide, total slip is less than channel spacing, and offsets are measured shortly after an earthquake. In these simulations, postearthquake landscape evolution widens the geomorphic expression of the fault zone and modifies apparent channel offsets. We distinguish this “geomorphic fault zone” from the tectonic fault zone (zone of coseismic distributed deformation). This study highlights the capability of landscape evolution models to explore a range of conditions not easily defined in natural examples and the importance of averaging multiple measurements. Our results verify that paleoseismic studies must consider how geomorphic change has modified offset markers and use caution interpreting slip histories with multiple earthquakes.

Journal of Geophysical Research Solid Earth↗
USGS2020

Are the stress drops of small earthquakes good predictors of the stress drops of moderate-to-large earthquakes?

The stress drops of small earthquakes often exhibit spatial patterns of variability. If moderate and large earthquakes follow the same spatial patterns, the stress drops of possible future damaging earthquakes could be better predicted by considering the stress drops of nearby small events. Better stress drop predictability could reduce ground-motion uncertainty in Probabilistic Seismic Hazard Assessment (PSHA) and Earthquake Early Warning (EEW). I find that for an internally consistent stress drop catalog of M1.8-3.1 events in southern California, the stress drops of the bigger earthquakes are predictable from the nearby smaller events. However, this catalog only weakly spatially correlates with another catalog of M3.0-5.8 earthquakes, and is spatially uncorrelated with five other stress drop catalogs of M≥3.4 earthquakes. For southern California events M5.5-7.5, stress drops compiled from the literature are weakly spatially correlated with the stress drops of the M1.8-3.1 events, although the correlations are not statistically significant. The lack of strong spatial correlation may be due to actual differences in the controlling factors of stress drop, for example dynamic weakening in moderate-to-large earthquakes. Alternatively, a stronger spatial correlation may exist that is obscured by methodological heterogeneity and large errors in the stress drop estimates. Either way, the stress drops of small earthquakes do not appear to be good predictors of the stress drops of nearby moderate-to-large earthquakes, at least for current techniques of stress drop estimation. If these results are representative, small-earthquake stress drops are not currently useful for substantially reducing uncertainty in PSHA and EEW.

California↗
USGS2020

Geodetic measurements of slow slip events southeast of Parkfield, CA

Tremor and low-frequency earthquakes are presumed to be indicative of surrounding slow, aseismic slip that is often below geodetic detection thresholds. This study uses data from borehole seismometers and long-baseline laser strainmeters to observe both the seismic and geodetic signatures of episodic tremor and slip on the Parkfield region of the San Andreas Fault near Cholame, CA. The observed occurrence rates of both the tremors and co-located families of low-frequency earthquakes are not steady but instead exhibit quasiperiodic bursts of increased activity. We show that these periods of elevated seismic activity correlate with statistically significant stacked strain signals consisting of 44 slow-slip events. Modeled individual slow-slip events and their total summed moment, which are constrained by seismic signals and stacked strain, respectively, indicate that the individual moment magnitudes of these events range from 4.6–5.2. We find that the measured geodetic signal likely precedes the seismic signal by several hours, consistent with the aseismic slip preceding and driving the observed seismic tremor activity. We confirm that strike-slip faults, in addition to subduction zones, are capable of producing episodic tremor and slip.

California↗
USGS2020

Methods for rapidly estimating velocity precision from GNSS time series in the presence of temporal correlation: A new method and comparison of existing methods

Time series of position estimates from Global Navigational Satellite System (GNSS) are used to measure the velocities of points on the surface of the Earth. Along with the velocity estimates, a measure of the precision is needed to assess the quality of the velocity measurement. Here, I evaluate rate uncertainties provided by four different methods that have been applied to geodetic time series. The most rigorous approach uses a data covariance that incorporates a variety of noise processes relevant to geodetic time series but is computationally demanding. Two other approaches are efficient algorithms and are used widely, but both can provide less rigorous estimates of the rate uncertainty. I propose and evaluate a fourth method, which provides estimates of rate uncertainty closer to the rigorous approach but is significantly less computationally demanding. I have evaluated all three methods against the more rigorous method using both simulations and time series from 190 GNSS sites. For data best characterized as having a flicker type noise process, one of the widely used methods overestimates the uncertainty by up to a factor of 2, while the other widely used method underestimates the uncertainty by less than a factor of 2. For a random-walk process, both methods underestimate the rate uncertainty by a factor of 3 to 5.

Journal of Geophysical Research↗
USGS2020

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

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

California↗
USGS2019

Assessing the seasonal dynamics of nitrate and sulfate aerosols at the South Pole utilizing stable isotopes

Atmospheric nitrate (NO 3 − = particulate NO 3 − + gas‐phase nitric acid [HNO 3 ]) and sulfate (SO 4 2− ) are key molecules that play important roles in numerous atmospheric processes. Here, the seasonal cycles of NO 3 − and total suspended particulate sulfate (SO 4 2− (TSP) ) were evaluated at the South Pole from aerosol samples collected weekly for approximately 10 months (26 January to 25 October) in 2002 and analyzed for their concentration and isotopic compositions. Aerosol NO 3 − was largely affected by snowpack emissions in which [NO 3 − ] and δ 15 N(NO 3 − ) were highest (49.3 ± 21.4 ng/m 3 , n = 8) and lowest (−47.0 ± 11.7‰, n = 5), respectively, during periods of sunlight in the interior of Antarctica. The seasonal cycle of Δ 17 O(NO 3 − ) reflected tropospheric chemistry year‐round with lower values observed during sunlight periods and higher values observed during dark periods, reflecting shifts from HO x ‐ to O 3 ‐dominated oxidation chemistry. SO 4 2− (TSP) concentrations were highest during austral summer and fall (86.7 ± 73.7 ng/m 3 , n = 18) and are indicated to be derived from dimethyl sulfide (DMS) emissions, as δ 34 S(SO 4 2− ) (TSP) values (18.5 ± 1.0‰, n = 10) were similar to literature δ 34 S(DMS) values. The seasonal cycle of Δ 17 O(SO 4 2− ) (TSP) exhibited minima during austral summer (0.9 ± 0.1‰, n = 5) and maxima during austral fall (1.3 ± 0.3‰, n = 6) and austral spring (1.6 ± 0.1‰, n = 5), indicating a shift from HO x ‐ to O 3 ‐dominated chemistry in the atmospheric derived SO 4 2− component. Overall, the budgets of NO 3 − and SO 4 2− (TSP) at the South Pole were complex functions of transport, localized chemistry, biological activity, and meteorological conditions, and these results will be important for interpretations of oxyanions in ice core records in the interior of Antarctica.

Journal of Geophysical Research D: Atmospheres↗