USGS Science⌕ Search

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

Ground surface deformation patterns, magma supply, and magma storage at Okmok volcano, Alaska, from InSAR analysis: 1. Intereruption deformation, 1997–2008

Starting soon after the 1997 eruption at Okmok volcano and continuing until the start of the 2008 eruption, magma accumulated in a storage zone centered ~3.5 km beneath the caldera floor at a rate that varied with time. A Mogi-type point pressure source or finite sphere with a radius of 1 km provides an adequate fit to the deformation field portrayed in time-sequential interferometric synthetic aperture radar images. From the end of the 1997 eruption through summer 2004, magma storage increased by 3.2–4.5 × 10 7 m 3 , which corresponds to 75–85% of the magma volume erupted in 1997. Thereafter, the average magma supply rate decreased such that by 10 July 2008, 2 days before the start of the 2008 eruption, magma storage had increased by 3.7–5.2 × 10 7 m 3 or 85–100% of the 1997 eruption volume. We propose that the supply rate decreased in response to the diminishing pressure gradient between the shallow storage zone and a deeper magma source region. Eventually the effects of continuing magma supply and vesiculation of stored magma caused a critical pressure threshold to be exceeded, triggering the 2008 eruption. A similar pattern of initially rapid inflation followed by oscillatory but generally slowing inflation was observed prior to the 1997 eruption. In both cases, withdrawal of magma during the eruptions depressurized the shallow storage zone, causing significant volcano-wide subsidence and initiating a new intereruption deformation cycle.

Alaska↗
USGS2010

Ground surface deformation patterns, magma supply, and magma storage at Okmok volcano, Alaska, from InSAR analysis: 2. Coeruptive deflation, July-August 2008

A hydrovolcanic eruption near Cone D on the floor of Okmok caldera, Alaska, began on 12 July 2008 and continued until late August 2008. The eruption was preceded by inflation of a magma reservoir located beneath the center of the caldera and ∼3 km below sea level (bsl), which began immediately after Okmok's previous eruption in 1997. In this paper we use data from several radar satellites and advanced interferometric synthetic aperture radar (InSAR) techniques to produce a suite of 2008 coeruption deformation maps. Most of the surface deformation that occurred during the eruption is explained by deflation of a Mogi-type source located beneath the center of the caldera and 2–3 km bsl, i.e., essentially the same source that inflated prior to the eruption. During the eruption the reservoir deflated at a rate that decreased exponentially with time with a 1/ e time constant of ∼13 days. We envision a sponge-like network of interconnected fractures and melt bodies that in aggregate constitute a complex magma storage zone beneath Okmok caldera. The rate at which the reservoir deflates during an eruption may be controlled by the diminishing pressure difference between the reservoir and surface. A similar mechanism might explain the tendency for reservoir inflation to slow as an eruption approaches until the pressure difference between a deep magma production zone and the reservoir is great enough to drive an intrusion or eruption along the caldera ring-fracture system.

Alaska↗
USGS2010

Anisotropy, repeating earthquakes, and seismicity associated with the 2008 eruption of Okmok Volcano, Alaska

We use shear wave splitting (SWS) analysis and double-difference relocation to examine temporal variations in seismic properties prior to and accompanying magmatic activity associated with the 2008 eruption of Okmok volcano, Alaska. Using bispectrum cross-correlation, a multiplet of 25 earthquakes is identified spanning five years leading up to the eruption, each event having first motions compatible with a normal fault striking NE–SW. Cross-correlation differential times are used to relocate earthquakes occurring between January 2003 and February 2009. The bulk of the seismicity prior to the onset of the eruption on 12 July 2008 occurred southwest of the caldera beneath a geothermal field. Earthquakes associated with the onset of the eruption occurred beneath the northern portion of the caldera and started as deep as 13 km. Subsequent earthquakes occurred predominantly at 3 km depth, coinciding with the depth at which the magma body has been modeled using geodetic data. Automated SWS analysis of the Okmok catalog reveals radial polarization outside the caldera and a northwest-southeast polarization within. We interpret these polarizations in terms of a magma reservoir near the center of the caldera, which we model with a Mogi point source. SWS analysis using the same input processing parameters for each event in the multiplet reveals no temporal changes in anisotropy over the duration of the multiplet, suggesting either a short-term or small increase in stress just before the eruption that was not detected by GPS, or eruption triggering by a mechanism other than a change of stress in the system.

Alaska↗
USGS2010

Lessons from (triggered) tremor

I test a “clock-advance” model that implies triggered tremor is ambient tremor that occurs at a sped-up rate as a result of loading from passing seismic waves. This proposed model predicts that triggering probability is proportional to the product of the ambient tremor rate and a function describing the efficacy of the triggering wave to initiate a tremor event. Using data mostly from Cascadia, I have compared qualitatively a suite of teleseismic waves that did and did not trigger tremor with ambient tremor rates. Many of the observations are consistent with the model if the efficacy of the triggering wave depends on wave amplitude. One triggered tremor observation clearly violates the clock-advance model. The model prediction that larger triggering waves result in larger triggered tremor signals also appears inconsistent with the measurements. I conclude that the tremor source process is a more complex system than that described by the clock-advance model predictions tested. Results of this and previous studies also demonstrate that (1) conditions suitable for tremor generation exist in many tectonic environments, but, within each, only occur at particular spots whose locations change with time; (2) any fluid flow must be restricted to less than a meter; (3) the degree to which delayed failure and secondary triggering occurs is likely insignificant; and 4) both shear and dilatational deformations may trigger tremor. Triggered and ambient tremor rates correlate more strongly with stress than stressing rate, suggesting tremor sources result from time-dependent weakening processes rather than simple Coulomb failure.

Cascadia↗
USGS2010

Volcano monitoring using GPS: Developing data analysis strategies based on the June 2007 Kīlauea Volcano intrusion and eruption

The global positioning system (GPS) is one of the most common techniques, and the current state of the art, used to monitor volcano deformation. In addition to slow (several centimeters per year) displacement rates, GPS can be used to study eruptions and intrusions that result in much larger (tens of centimeters over hours-days) displacements. It is challenging to resolve precise positions using GPS at subdaily time intervals because of error sources such as multipath and atmospheric refraction. In this paper, the impact of errors due to multipath and atmospheric refraction at subdaily periods is examined using data from the GPS network on Kīlauea Volcano, Hawai'i. Methods for filtering position estimates to enhance precision are both simulated and tested on data collected during the June 2007 intrusion and eruption. Comparisons with tiltmeter records show that GPS instruments can precisely recover the timing of the activity.

Hawai'i↗
USGS2010

Cyclic ground tilt associated with the 2004–2008 eruption of Mount St. Helens

The long‐term behavior of the 2004–2008 effusive eruption of Mount St. Helens was characterized by a gradual decline in the rates of seismicity, dome growth, and broad‐scale ground deformation, but shallow near‐periodic “drumbeat” earthquakes over timescales of minutes indicated episodic short‐term behavior. In part to better characterize this behavior and any associated ground deformation, a network of tiltmeters was installed and recorded thousands of cyclic tilt events within the crater. The duration of these events, from minutes to hours, was too long for them to be directly caused by the drumbeat seismicity. Tilt events were characterized by a recoverable, asymmetric pattern of rapid tilt away from the vent followed by a more gradual reversal, were highly correlated between different tiltmeters in the crater, sometimes occurred in association with volcanic seismicity, and ceased at the end of the eruption. Tilt vectors converged on a point just south of the center of the preexisting 1980s lava dome, and the absence of detectable tilt outside the crater suggests a shallow source (<1 km). We examine several models, including cycles of conduit pressurization and plug slip or gas loss and stick‐slip behavior on the interface between the ascending plug and the 1980s lava dome. The small number of stations within the crater prevents a unique determination of source type or geometry, but results are consistent with a mechanism involving extrusion of the semisolid dacite plug and/or cycles of conduit pressurization.

Washington↗
USGS2010

Coulomb stress interactions among M&ge;5.9 earthquakes in the Gorda deformation zone and on the Mendocino Fracture Zone, Cascadia megathrust, and northern San Andreas fault

The Gorda deformation zone, a 50,000 km 2 area of diffuse shear and rotation offshore northernmost California, has been the site of 20 M ≥ 5.9 earthquakes on four different fault orientations since 1976, including four M ≥ 7 shocks. This is the highest rate of large earthquakes in the contiguous United States. We calculate that the source faults of six recent M ≥ 5.9 earthquakes had experienced ≥0.6 bar Coulomb stress increases imparted by earthquakes that struck less than 9 months beforehand. Control tests indicate that ≥0.6 bar Coulomb stress interactions between M ≥ 5.9 earthquakes separated by <9 months are unlikely to occur by random chance, suggesting that the multiple short-term stress interactions observed among the recent Gorda zone earthquakes are not an apparent effect. In all well-constrained ≥0.2 bar Coulomb stress interactions between earthquakes that occurred within 4 years of each other, the second earthquake is promoted. On longer timescales, calculated stress changes imparted by the 1980 M w = 7.3 Trinidad earthquake are consistent with the locations of M ≥ 5.9 earthquakes in the Gorda zone until at least 1995, as well as earthquakes on the Mendocino Fault Zone in 1994 and 2000. Coulomb stress changes imparted by the 1980 earthquake are also consistent with its distinct elbow-shaped aftershock pattern. From these observations, we derive generalized static stress interactions among right-lateral, left-lateral and thrust faults near triple junctions.

California↗
USGS2010

Origins of large-volume, compositionally zoned volcanic eruptions: New constraints from U-series isotopes and numerical thermal modeling for the 1912 Katmai-Novarupta eruption

We present the results of a combined U-series isotope and numerical modeling study of the 1912 Katmai-Novarupta eruption in Alaska. A stratigraphically constrained set of samples have compositions that range from basalt through basaltic andesite, andesite, dacite, and rhyolite. The major and trace element range can be modeled by 80–90% closed-system crystal fractionation over a temperature interval from 1279°C to 719°C at 100 MPa, with an implied volume of parental basalt of ∼65 km 3 . Numerical models suggest, for wall rock temperatures appropriate to this depth, that 90% of this volume of magma would cool and crystallize over this temperature interval within a few tens of kiloyears. However, the range in 87 Sr/ 86 Sr, ( 230 Th/ 238 U), and ( 226 Ra/ 230 Th) requires open-system processes. Assimilation of the host sediments can replicate the range of Sr isotopes. The variation of ( 226 Ra/ 230 Th) ratios in the basalt to andesite compositional range requires that these were generated less than several thousand years before eruption. Residence times for dacites are close to 8000 years, whereas the rhyolites appear to be 50–200 kyr old. Thus, the magmas that erupted within only 60 h had a wide range of crustal residence times. Nevertheless, they were emplaced in the same thermal regime and evolved along similar liquid lines of descent from parental magmas with similar compositions. The system was built progressively with multiple inputs providing both mass and heat, some of which led to thawing of older silicic material that provided much of the rhyolite.

Alaska↗
USGS2010

Variability and trends in dry day frequency and dry event length in the southwestern United States

Daily precipitation from 22 National Weather Service first-order weather stations in the southwestern United States for water years 1951 through 2006 are used to examine variability and trends in the frequency of dry days and dry event length. Dry events with minimum thresholds of 10 and 20 consecutive days of precipitation with less than 2.54 mm are analyzed. For water years and cool seasons (October through March), most sites indicate negative trends in dry event length (i.e., dry event durations are becoming shorter). For the warm season (April through September), most sites also indicate negative trends; however, more sites indicate positive trends in dry event length for the warm season than for water years or cool seasons. The larger number of sites indicating positive trends in dry event length during the warm season is due to a series of dry warm seasons near the end of the 20th century and the beginning of the 21st century. Overall, a large portion of the variability in dry event length is attributable to variability of the El Niño–Southern Oscillation, especially for water years and cool seasons. Our results are consistent with analyses of trends in discharge for sites in the southwestern United States, an increased frequency in El Niño events, and positive trends in precipitation in the southwestern United States.

Journal of Geophysical Research D: Atmospheres↗
USGS2009

Identification of hydrated silicate minerals on Mars using MRO-CRISM: Geologic context near Nili Fossae and implications for aqueous alteration

[1] The Noachian terrain west of the Isidis basin hosts a diverse collection of alteration minerals in rocks comprising varied geomorphic units within a 100,000 km 2 region in and near the Nili Fossae. Prior investigations in this region by the Observatoire pour l'Minéralogie, l'Eau, les Glaces, et l'Activité (OMEGA) instrument on Mars Express revealed large exposures of both mafic minerals and iron magnesium phyllosilicates in stratigraphic context. Expanding on the discoveries of OMEGA, the Compact Reconnaissance Imaging Spectrometer for Mars (CRISM) aboard the Mars Reconnaissance Orbiter (MRO) has found more spatially widespread and mineralogically diverse alteration minerals than previously realized, which represent multiple aqueous environments. Using CRISM near-infrared spectral data, we detail the basis for identification of iron and magnesium smectites (including both nontronite and more Mg-rich varieties), chlorite, prehnite, serpentine, kaolinite, potassium mica (illite or muscovite), hydrated (opaline) silica, the sodium zeolite analcime, and magnesium carbonate. The detection of serpentine and analcime on Mars is reported here for the first time. We detail the geomorphic context of these minerals using data from high-resolution imagers onboard MRO in conjunction with CRISM. We find that the distribution of alteration minerals is not homogeneous; rather, they occur in provinces with distinctive assemblages of alteration minerals. Key findings are (1) a distinctive stratigraphy, in and around the Nili Fossae, of kaolinite and magnesium carbonate in bedrock units always overlying Fe/Mg smectites and (2) evidence for mineral phases and assemblages indicative of low-grade metamorphic or hydrothermal aqueous alteration in cratered terrains. The alteration minerals around the Nili Fossae are more typical of those resulting from neutral to alkaline conditions rather than acidic conditions, which appear to have dominated much of Mars. Moreover, the mineralogic diversity and geologic context of alteration minerals found in the region around the Nili Fossae indicates several episodes of aqueous activity in multiple distinct environments.

Journal of Geophysical Research E: Planets↗
USGS2010

Seasonal H2O and CO2 ice cycles at the Mars Phoenix landing site: 1. Prelanding CRISM and HiRISE observations

The condensation, evolution, and sublimation of seasonal water and carbon dioxide ices were characterized at the Mars Phoenix landing site from Martian northern midsummer to midspring (L s ∼ 142° – L s ∼ 60°) for the year prior to the Phoenix landing on 25 May 2008. Ice relative abundances and grain sizes were estimated using data from the Compact Reconnaissance Imaging Spectrometer for Mars and High Resolution Imaging Science Experiment aboard Mars Reconnaissance Orbiter and a nonlinear mixing model. Water ice first appeared at the Phoenix landing site during the afternoon in late summer (L s ∼ 167°) as an optically thin layer on top of soil. CO 2 ice appeared after the fall equinox. By late winter (L s ∼ 344°), the site was covered by relatively pure CO 2 ice (∼30 cm thick), with a small amount of ∼100 μ m diameter water ice and soil. As spring progressed, CO 2 ice grain sizes gradually decreased, a change interpreted to result from granulation during sublimation losses. The combined effect of CO 2 sublimation and decreasing H 2 O ice grain sizes allowed H 2 O ice to dominate spectra during the spring and significantly brightened the surface. CO 2 ice disappeared by early spring (L s ∼ 34°) and H 2 O ice by midspring (L s ∼ 59°). Spring defrosting was not uniform and occurred more rapidly over the centers of polygons and geomorphic units with relatively higher thermal inertia values.

Journal of Geophysical Research E: Planets↗
USGS2010

Silica in a Mars analog environment: Ka u Desert, Kilauea Volcano, Hawaii

Airborne Visible/Near-Infrared Imaging Spectrometer (AVIRIS) data acquired over the Ka u Desert are atmospherically corrected to ground reflectance and used to identify the mineralogic components of relatively young basaltic materials, including 250-700 and 200-400 year old lava flows, 1971 and 1974 flows, ash deposits, and solfatara incrustations. To provide context, a geologic surface units map is constructed, verified with field observations, and supported by laboratory analyses. AVIRIS spectral endmembers are identified in the visible (0.4 to 1.2 ??m) and short wave infrared (2.0 to 2.5 ??m) wavelength ranges. Nearly all the spectral variability is controlled by the presence of ferrous and ferric iron in such minerals as pyroxene, olivine, hematite, goethite, and poorly crystalline iron oxides or glass. A broad, nearly ubiquitous absorption feature centered at 2.25 ??m is attributed to opaline (amorphous, hydrated) silica and is found to correlate spatially with mapped geologic surface units. Laboratory analyses show the silica to be consistently present as a deposited phase, including incrustations downwind from solfatara vents, cementing agent for ash duricrusts, and thin coatings on the youngest lava flow surfaces. A second, Ti-rich upper coating on young flows also influences spectral behavior. This study demonstrates that secondary silica is mobile in the Ka u Desert on a variety of time scales and spatial domains. The investigation from remote, field, and laboratory perspectives also mimics exploration of Mars using orbital and landed missions, with important implications for spectral characterization of coated basalts and formation of opaline silica in arid, acidic alteration environments. Copyright 2010 by the American Geophysical Union.

Journal of Geophysical Research E: Planets↗
USGS2010

Mineralogy and stratigraphy of phyllosilicate-bearing and dark mantling units in the greater Mawrth Vallis/west Arabia Terra area: Constraints on geological origin

Analyses of MRO/CRISM images of the greater Mawrth Vallis region of Mars affirm the presence of two primary phyllosilicate assemblages throughout a region ∼1000 × 1000 km. These two units consist of an Fe/Mg-phyllosilicate assemblage overlain by an Al-phyllosilicate and hydrated silica assemblage. The lower unit contains Fe/Mg-smectites, sometimes combined with one or more of these other Fe/Mg-phyllosilicates: serpentine, chlorite, biotite, and/or vermiculite. It is more than 100 m thick and finely layered at meter scales. The upper unit includes Al-smectite, kaolin group minerals, and hydrated silica. It is tens of meters thick and finely layered as well. A common phyllosilicate stratigraphy and morphology is observed throughout the greater region wherever erosional windows are present. This suggests that the geologic processes forming these units must have occurred on at least a regional scale. Sinuous ridges (interpreted to be inverted channels) and narrow channels cut into the upper clay-bearing unit suggesting that aqueous processes were prevalent after, and possibly during, the deposition of the layered units. We propose that layered units may have been deposited at Mawrth Vallis and then subsequently altered to form the hydrated units. The Fe/Mg-phyllosilicate assemblage is consistent with hydrothermal alteration or pedogenesis of mafic to ultramafic rocks. The Al-phyllosilicate/hydrated silica unit may have formed through alteration of felsic material or via leaching of basaltic material through pedogenic alteration or a mildly acidic environment. These phyllosilicate-bearing units are overlain by a darker, relatively unaltered, and indurated material that has probably experienced a complex geological history.

Journal of Geophysical Research↗
USGS2009

Global distribution of bedrock exposures on Mars using THEMIS high-resolution thermal inertia

We investigate high thermal inertia surfaces using the Mars Odyssey Thermal Emission Imaging System (THEMIS) nighttime temperature images (100 m/pixel spatial sampling). For this study, we interpret any pixel in a THEMIS image with a thermal inertia over 1200 J m -2 K -1 s -1/2 as "bedrock" which represents either in situ rock exposures or rock-dominated surfaces. Three distinct morphologies, ranked from most to least common, are associated with these high thermal inertia surfaces: (1) valley and crater walls associated with mass wasting and high surface slope angles; (2) floors of craters with diameters >25 km and containing melt or volcanics associated with larger, high-energy impacts; and (3) intercrater surfaces with compositions significantly more mafic than the surrounding regolith. In general, bedrock instances on Mars occur as small exposures (less than several square kilometers) situated in lower-albedo (<0.18), moderate to high thermal inertia (>350 J m -2 K -1 s -1/2 ), and relatively dust-free (dust cover index <0.95) regions; however, there are instances that do not follow these generalizations. Most instances are concentrated in the southern highlands, with very few located at high latitudes (poleward of 45 o N and 58 o S), suggesting enhanced mechanical breakdown probably associated with permafrost. Overall, Mars has very little exposed bedrock with only 960 instances identified from 75 o S to 75 o N with likely <3500 km 2 exposed, representing???1% of the total surface area. These data indicate that Mars has likely undergone large-scale surface processing and reworking, both chemically and mechanically, either destroying or masking a majority of the bedrock exposures on the planet. Copyright 2009 by the American Geophysical Union.

Journal of Geophysical Research E: Planets↗
USGS2010

Detection and mapping of hydrocarbon deposits on Titan

We report the identification of compounds on Titan's surface by spatially resolved imaging spectroscopy methods through Titan's atmosphere, and set upper limits to other organic compounds. We present evidence for surface deposits of solid benzene (C 6 H 6 ), solid and/or liquid ethane (C 2 H 6 ), or methane (CH 4 ), and clouds of hydrogen cyanide (HCN) aerosols using diagnostic spectral features in data from the Cassini Visual and Infrared Mapping Spectrometer (VIMS). Cyanoacetylene (2-propynenitrile, IUPAC nomenclature, HC 3 N) is indicated in spectra of some bright regions, but the spectral resolution of VIMS is insufficient to make a unique identification although it is a closer match to the feature previously attributed to CO 2 . We identify benzene, an aromatic hydrocarbon, in larger abundances than expected by some models. Acetylene (C 2 H 2 ), expected to be more abundant on Titan according to some models than benzene, is not detected. Solid acetonitrile (CH 3 CN) or other nitriles might be candidates for matching other spectral features in some Titan spectra. An as yet unidentified absorption at 5.01- μ m indicates that yet another compound exists on Titan's surface. We place upper limits for liquid methane and ethane in some locations on Titan and find local areas consistent with millimeter path lengths. Except for potential lakes in the southern and northern polar regions, most of Titan appears “dry.” Finally, we find there is little evidence for exposed water ice on the surface. Water ice, if present, must be covered with organic compounds to the depth probed by 1–5- μ m photons: a few millimeters to centimeters.

Journal of Geophysical Research - Planets↗
USGS2009

Aeolian dunes as ground truth for atmospheric modeling on Mars

Martian aeolian dunes preserve a record of atmosphere/surface interaction on a variety of scales, serving as ground truth for both Global Climate Models (GCMs) and mesoscale climate models, such as the Mars Regional Atmospheric Modeling System (MRAMS). We hypothesize that the location of dune fields, expressed globally by geographic distribution and locally by dune centroid azimuth (DCA), may record the long-term integration of atmospheric activity across a broad area, preserving GCM-scale atmospheric trends. In contrast, individual dune morphology, as expressed in slipface orientation (SF), may be more sensitive to localized variations in circulation, preserving topographically controlled mesoscale trends. We test this hypothesis by comparing the geographic distribution, DCA, and SF of dunes with output from the Ames Mars GCM and, at a local study site, with output from MRAMS. When compared to the GCM: 1) dunes generally lie adjacent to areas with strongest winds, 2) DCA agrees fairly well with GCM modeled wind directions in smooth-floored craters, and 3) SF does not agree well with GCM modeled wind directions. When compared to MRAMS modeled winds at our study site: 1) DCA generally coincides with the part of the crater where modeled mean winds are weak, and 2) SFs are consistent with some weak, topographically influenced modeled winds. We conclude that: 1) geographic distribution may be valuable as ground truth for GCMs, 2) DCA may be useful as ground truth for both GCM and mesoscale models, and 3) SF may be useful as ground truth for mesoscale models. Copyright 2009 by the American Geophysical Union.

Journal of Geophysical Research E: Planets↗
USGS2010

A comparison of Martian north and south polar cold spots and the long‐term effects of the 2001 global dust storm

In the 1970s, Viking and Mariner observed areas in the polar regions of Mars with winter brightness temperatures below the expected kinetic temperatures for CO 2 ice sublimation. These areas have since been termed “cold spots” and have been identified as surface deposits of CO 2 atmospheric condensates and, occasionally, active CO 2 storms. Three Mars years of data from the Mars Global Surveyor Thermal Emission Spectrometer were used to observe autumn and winter cold spot activity. In this study, cold spots that occur near and on the southern perennial cap were compared to those found near or on the northern perennial cap. On the southern perennial cap, cold spots associated with topographic features (induced by orographic lifting) were less common than cold spots independent of topography, similar to the north. However, the cold spots in the south lasted longer than those observed in the north. There is also evidence that cold spot formation in the south was affected by the global dust storm of 2001, even though the dust storm occurred during the southern spring and summer seasons. Prior to the dust storm, the amount of overall cold spot activity closer to the perennial cap increased and the average CO 2 grain size for most of the cold spots increased as well. Following the dust storm, the majority of cold spots in the south increased in size and duration but they did not form north of 62°S latitude, whereas, in other years, cold spots formed as far north as 48°S.

Journal of Geophysical Research E: Planets↗