USGS Science⌕ Search

SEARCH · USGS Science

Results for “Special Publication”

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

Shallow-conduit dynamics at Stromboli Volcano, Italy, imaged from waveform inversions

Modelling of Very-Long-Period (VLP) seismic data recorded during explosive activity at Stromboli in 1997 provides an image of the uppermost 1 km of its volcanic plumbing system. Two distinct dyke-like conduit structures are identified, each representative of explosive eruptions from two different vents located near the northern and southern perimeters of the summit crater. Observed volumetric changes in the dykes are viewed as the result of a piston-like action of the magma associated with the disruption of a gas slug transiting through discontinuities in the dyke apertures. Accompanying these volumetric source components are single vertical forces resulting from an exchange of linear momentum between the source and the Earth. In the dyke system underlying the northern vent, a primary disruption site is observed at an elevation near 440 m where a bifurcation in the conduit occurs. At a depth of 80 m below sea level, a sharp corner in the conduit marks another location where the elastic response of the solid to the action of the upper source induces pressure and momentum changes in the magma. In the conduit underlying the southern vent, the junction of two inclined dykes with a sub-vertical dyke at 520 m elevation is a primary site of gas slug disruption, and another conduit corner 280 m below sea level represents a coupling location between the elastic response of the solid and fluid motion.

Stromboli↗
USGS2009

Implications of estimated magmatic additions and recycling losses at the subduction zones of accretionary (non-collisional) and collisional (suturing) orogens

Arc magmatism at subduction zones (SZs) most voluminously supplies juvenile igneous material to build rafts of continental and intra-oceanic or island arc (CIA) crust. Return or recycling of accumulated CIA material to the mantle is also most vigorous at SZs. Recycling is effected by the processes of sediment subduction, subduction erosion, and detachment and sinking of deeply underthrust sectors of CIA crust. Long-term (>10-20 Ma) rates of additions and losses can be estimated from observational data gathered where oceanic crust underruns modern, long-running (Cenozoic to mid-Mesozoic) ocean-margin subduction zones (OMSZs, e.g. Aleutian and South America SZs). Long-term rates can also be observationally assessed at Mesozoic and older crust-suturing subduction zone (CSSZs) where thick bodies of CIA crust collided in tectonic contact (e.g. Wopmay and Appalachian orogens, India and SE Asia). At modern OMSZs arc magmatic additions at intra-oceanic arcs and at continental margins are globally estimated at c. 1.5 AU and c. 1.0 AU, respectively (1 AU, or Armstrong Unit,= 1 km 3 a -1 of solid material). During collisional suturing at fossil CSSZs, global arc magmatic addition is estimated at 0.2 AU. This assessment presumes that in the past the global length of crustal collision zones averaged c. 6000 km, which is one-half that under way since the early Tertiary. The average long-term rate of arc magmatic additions extracted from modern OMSZs and older CSSZs is thus evaluated at 2.7 AU. Crustal recycling at Mesozoic and younger OMSZs is assessed at c. 60 km 3 Ma -1 km -1 (c. 60% by subduction erosion). The corresponding global recycling rate is c. 2.5 AU. At CSSZs of Mesozoic, Palaeozoic and Proterozoic age, the combined upper and lower plate losses of CIA crust via subduction erosion, sediment subduction, and lower plate crustal detachment and sinking are assessed far less securely at c. 115 km 3 Ma -1 km -1 . At a global length of 6000 km, recycling at CSSZs is accordingly c. 0.7 AU. The collective loss of CIA crust estimated for modern OMSZs and for older CSSZs is thus estimated at c. 3.2 AU. SZ additions (2.7 AU) and subtractions (23.2 AU) are similar. Because many uncertainties and assumptions are involved in assessing and applying them to the deep past, the net growth of CIA crust during at least Phanerozoic time is viewed as effectively nil. With increasing uncertainty, the long-term balance can be applied to the Proterozoic, but not before the initiation of the present style of subduction at c. 3 Ga. Allowing that since this time a rounded-down rate of recycling of 3 AU is applicable, a startlingly high volume of CIA crust equal to that existing now has been recycled to the mantle. Although the recycled volume (c. 9 ?? 10 9 km 3 ) is small (c. 1%) compared with that of the mantle, it is large enough to impart to the mantle the signature of recycled CIA crust. Because subduction zones are not spatially fixed, and their average global lengths have episodically been less or greater than at present, recycling must have contributed significantly to creating recognized heterogeneities in mantle geochemistry. ?? The Geological Society of London 2009.

Geological Society Special Publication↗
USGS2009

Gas hydrate drilling transect across northern Cascadia margin - IODP Expedition 311

A transect of four sites (U1325, U1326, U1327 and U1329) across the northern Cascadia margin was established during Integrated Ocean Drilling Program Expedition 311 to study the occurrence and formation of gas hydrate in accretionary complexes. In addition to the transect sites, a fifth site (U1328) was established at a cold vent with active fluid flow. The four transect sites represent different typical geological environments of gas hydrate occurrence across the northern Cascadia margin from the earliest occurrence on the westernmost first accreted ridge (Site U1326) to the eastward limit of the gas hydrate occurrence in shallower water (Site U1329). Expedition 311 complements previous gas hydrate studies along the Cascadia accretionary complex, especially ODP Leg 146 and Leg 204 by extending the aperture of the transect sampled and introducing new tools to systematically quantify the gas hydrate content of the sediments. Among the most significant findings of the expedition was the occurrence of up to 20 m thick sand-rich turbidite intervals with gas hydrate concentrations locally exceeding 50% of the pore space at Sites U1326 and U1327. Moreover, these anomalous gas hydrate intervals occur at unexpectedly shallow depths of 50-120 metres below seafloor, which is the opposite of what was expected from previous models of gas hydrate formation in accretionary complexes, where gas hydrate was predicted to be more concentrated near the base of the gas hydrate stability zone just above the bottom-simulating reflector. Gas hydrate appears to be mainly concentrated in turbidite sand layers. During Expedition 311, the visual correlation of gas hydrate with sand layers was clearly and repeatedly documented, strongly supporting the importance of grain size in controlling gas hydrate occurrence. The results from the transect sites provide evidence for a structurally complex, lithology-controlled gas hydrate environment on the northern Cascadia margin. Local shallow occurrences of high gas hydrate concentrations contradict the previous model of gas hydrate formation at an accretionary prism. However, long-lived fluid flow (part of the old model) is still required to explain the shallow high gas hydrate concentrations, although it is most likely not pervasive throughout the entire accretionary prism, but rather localized and focused by the tectonic processes. Differences in the fluid flow regime across all of the transect drill sites indicate site-specific and probably disconnected (compartmented) deeper fluid sources in the various parts of the accretionary prism. The data and future analyses will yield a better understanding of the geologic controls, evolution and ultimate fate of gas hydrate in an accretionary prism as an important contribution to the role of gas hydrate methane gas in slope stability and possibly in climate change. ?? The Geological Society of London 2009.

Geological Society Special Publication↗
USGS2011

Cold-climate slope deposits and landscape modifications of the Mid-Atlantic Coastal Plain, Eastern USA

The effects of Pleistocene cold-climate geomorphology are distributed across the weathered and eroded Mid-Atlantic Coastal Plain uplands from the Wisconsinan terminal moraine south to Tidewater Virginia. Cold-climate deposits and landscape modifications are superimposed on antecedent landscapes of old, weathered Neogene upland gravels and Pleistocene marine terraces that had been built during warm periods and sea-level highstands. In New Jersey, sequences of surficial deposits define a long history of repeating climate change events. To the south across the Delmarva Peninsula and southern Maryland, most antecedent topography has been obscured by Late Pleistocene surficial deposits. These are spatially variable and are collectively described as a cold-climate alloformation. The cold-climate alloformation includes time-transgressive details of climate deterioration from at least marine isotope stage (MIS) 4 through the end of MIS 2. Some deposits and landforms within the alloformation may be as young as the Younger Dryas. Southwards along the trend of the Potomac River, these deposits and their climatic affinities become diffused. In Virginia, a continuum of erosion and surficial deposits appears to be the product of ‘normal’ temperate, climate-forced processes. The cold-climate alloformation and more temperate deposits in Virginia are being partly covered by Holocene alluvium and bay mud.

Geological Society Special Publication↗
USGS2011

The indication of Martian gully formation processes by slope-area analysis

The formation process of recent gullies on Mars is currently under debate. This study aims to discriminate between the proposed formation processes - pure water flow, debris flow and dry mass wasting - through the application of geomorphological indices commonly used in terrestrial geomorphology. High-resolution digital elevation models (DEMs) of Earth and Mars were used to evaluate the drainage characteristics of small slope sections. Data from Earth were used to validate the hillslope, debris-flow and alluvial process domains previously found for large fluvial catchments on Earth, and these domains were applied to gullied and ungullied slopes on Mars. In accordance with other studies, our results indicate that debris flow is one of the main processes forming the Martian gullies that were being examined. The source of the water is predominantly distributed surface melting, not an underground aquifer. Evidence is also presented indicating that other processes may have shaped Martian crater slopes, such as ice-assisted creep and solifluction, in agreement with the proposed recent Martian glacial and periglacial climate. Our results suggest that, within impact craters, different processes are acting on differently oriented slopes, but further work is needed to investigate the potential link between these observations and changes in Martian climate.

Geological Society Special Publication↗
USGS2011

Polyphase Neoproterozoic orogenesis within the east Africa- Antarctica orogenic belt in central and northern Madagascar

Our recent geological survey of the basement of central and northern Madagascar allowed us to re-evaluate the evolution of this part of the East Africa-Antarctica Orogen (EAAO). Five crustal domains are recognized, characterized by distinctive lithologies and histories of sedimentation, magmatism, deformation and metamorphism, and separated by tectonic and/or unconformable contacts. Four consist largely of Archaean metamorphic rocks (Antongil, Masora and Antananarivo Cratons, Tsaratanana Complex). The fifth (Bemarivo Belt) comprises Proterozoic meta-igneous rocks. The older rocks were intruded by plutonic suites at c. 1000 Ma, 820-760 Ma, 630-595 Ma and 560-520 Ma. The evolution of the four Archaean domains and their boundaries remains contentious, with two end-member interpretations evaluated: (1) all five crustal domains are separate tectonic elements, juxtaposed along Neoproterozoic sutures and (2) the four Archaean domains are segments of an older Archaean craton, which was sutured against the Bemarivo Belt in the Neoproterozoic. Rodinia fragmented during the early Neoproterozoic with intracratonic rifts that sometimes developed into oceanic basins. Subsequent Mid- Neoproterozoic collision of smaller cratonic blocks was followed by renewed extension and magmatism. The global 'Terminal Pan-African' event (560-490 Ma) finally stitched together the Mid-Neoproterozoic cratons to form Gondwana. ?? The Geological Society of London 2011.

Geological Society Special Publication↗
USGS2011

Deep rock damage in the San Andreas Fault revealed by P- and S-type fault-zone-guided waves

Damage to fault-zone rocks during fault slip results in the formation of a channel of low seismic-wave velocities. Within such channels guided seismic waves, denoted by Fg, can propagate. Here we show with core samples, well logs and Fg-waves that such a channel is crossed by the SAFOD (San Andreas Fault Observatory at Depth) borehole at a depth of 2.7 km near Parkfield, California, USA. This laterally extensive channel extends downwards to at least half way through the seismogenic crust, more than about 7 km. The channel supports not only the previously recognized Love-type- (FL) and Rayleigh-type- (FR) guided waves, but also a new fault-guided wave, which we name FF. As recorded 2.7 km underground, FF is normally dispersed, ends in an Airy phase, and arrives between the P- and S-waves. Modelling shows that FF travels as a leaky mode within the core of the fault zone. Combined with the drill core samples, well logs and the two other types of guided waves, FF at SAFOD reveals a zone of profound, deep, rock damage. Originating from damage accumulated over the recent history of fault movement, we suggest it is maintained either by fracturing near the slip surface of earthquakes, such as the 1857 Fort Tejon M 7.9, or is an unexplained part of the fault-creep process known to be active at this site.

California↗
USGS2012

Advances in carbonate exploration and reservoir analysis

The development of innovative techniques and concepts, and the emergence of new plays in carbonate rocks are creating a resurgence of oil and gas discoveries worldwide. The maturity of a basin and the application of exploration concepts have a fundamental influence on exploration strategies. Exploration success often occurs in underexplored basins by applying existing established geological concepts. This approach is commonly undertaken when new basins ‘open up’ owing to previous political upheavals. The strategy of using new techniques in a proven mature area is particularly appropriate when dealing with unconventional resources (heavy oil, bitumen, stranded gas), while the application of new play concepts (such as lacustrine carbonates) to new areas (i.e. ultra-deep South Atlantic basins) epitomizes frontier exploration. Many low-matrix-porosity hydrocarbon reservoirs are productive because permeability is controlled by fractures and faults. Understanding basic fracture properties is critical in reducing geological risk and therefore reducing well costs and increasing well recovery. The advent of resource plays in carbonate rocks, and the long-standing recognition of naturally fractured carbonate reservoirs means that new fracture and fault analysis and prediction techniques and concepts are essential.

Geological Society Special Publication↗
USGS2012

Late Pleistocene and Holocene uplift history of Cyprus: implications for active tectonics along the southern margin of the Anatolian microplate

The nature of the southern margin of the Anatolian microplate during the Neogene is complex, controversial and fundamental in understanding active plate-margin tectonics and natural hazards in the Eastern Mediterranean region. Our investigation provides new insights into the Late Pleistocene uplift history of Cyprus and the Troodos Ophiolite. We provide isotopic ( 14 C) and radiogenic (luminescence) dates of outcropping marine sediments in eastern Cyprus that identify periods of deposition during marine isotope stages (MIS) 3, 4, 5 and 6. Past sea-levels indicated by these deposits are c. 95±25 m higher in elevation than estimates of worldwide eustatic sea-level. An uplift rate of c. 1.8 mm/year and possibly as much as c. 4.1 mm/year in the past c. 26–40 ka is indicated. Holocene marine deposits also occur at elevations higher than those expected for past SL and suggest uplift rates of c. 1.2–2.1 mm/year. MIS-3 marine deposits that crop out in southern and western Cyprus indicate uniform island-wide uplift. We propose a model of tectonic wedging at a plate-bounding restraining bend as a mechanism for Late Pleistocene to Holocene uplift of Cyprus; uplift is accommodated by deformation and seismicity along the margins of the Troodos Ophiolite and re-activation of its low-angle, basal shear zone.

Geological Society, London, Special Publications: ↗
USGS2012

Magnetostratigraphy susceptibility for the Guadalupian Series GSSPs (Middle Permian) in Guadalupe Mountains National Park and adjacent areas in West Texas

Here we establish a magnetostratigraphy susceptibility zonation for the three Middle Permian Global boundary Stratotype Sections and Points (GSSPs) that have recently been defined, located in Guadalupe Mountains National Park, West Texas, USA. These GSSPs, all within the Middle Permian Guadalupian Series, define (1) the base of the Roadian Stage (base of the Guadalupian Series), (2) the base of the Wordian Stage and (3) the base of the Capitanian Stage. Data from two additional stratigraphic successions in the region, equivalent in age to the Kungurian–Roadian and Wordian–Capitanian boundary intervals, are also reported. Based on low-field, mass specific magnetic susceptibility (χ) measurements of 706 closely spaced samples from these stratigraphic sections and time-series analysis of one of these sections, we (1) define the magnetostratigraphy susceptibility zonation for the three Guadalupian Series Global boundary Stratotype Sections and Points; (2) demonstrate that χ datasets provide a proxy for climate cyclicity; (3) give quantitative estimates of the time it took for some of these sediments to accumulate; (4) give the rates at which sediments were accumulated; (5) allow more precise correlation to equivalent sections in the region; (6) identify anomalous stratigraphic horizons; and (7) give estimates for timing and duration of geological events within sections.

Texas↗
USGS2016

A stochastic bioenergetics model based approach to translating large river flow and temperature in to fish population responses: The pallid sturgeon example

In managing fish populations, especially at-risk species, realistic mathematical models are needed to help predict population response to potential management actions in the context of environmental conditions and changing climate while effectively incorporating the stochastic nature of real world conditions. We provide a key component of such a model for the endangered pallid sturgeon ( Scaphirhynchus albus ) in the form of an individual-based bioenergetics model influenced not only by temperature but also by flow. This component is based on modification of a known individual-based bioenergetics model through incorporation of: the observed ontogenetic shift in pallid sturgeon diet from marcroinvertebrates to fish; the energetic costs of swimming under flowing-water conditions; and stochasticity. We provide an assessment of how differences in environmental conditions could potentially alter pallid sturgeon growth estimates, using observed temperature and velocity from channelized portions of the Lower Missouri River mainstem. We do this using separate relationships between the proportion of maximum consumption and fork length and swimming cost standard error estimates for fish captured above and below the Kansas River in the Lower Missouri River. Critical to our matching observed growth in the field with predicted growth based on observed environmental conditions was a two-step shift in diet from macroinvertebrates to fish.

Missouri River↗
USGS2015

Hierarchical stochastic modeling of large river ecosystems and fish growth across spatio-temporal scales and climate models: the Missouri River endangered pallid sturgeon example

We present a hierarchical series of spatially decreasing and temporally increasing models to evaluate the uncertainty in the atmosphere – ocean global climate model (AOGCM) and the regional climate model (RCM) relative to the uncertainty in the somatic growth of the endangered pallid sturgeon (Scaphirhynchus albus). For effects on fish populations of riverine ecosystems, cli- mate output simulated by coarse-resolution AOGCMs and RCMs must be downscaled to basins to river hydrology to population response. One needs to transfer the information from these climate simulations down to the individual scale in a way that minimizes extrapolation and can account for spatio-temporal variability in the intervening stages. The goal is a framework to determine whether, given uncertainties in the climate models and the biological response, meaningful inference can still be made. The non-linear downscaling of climate information to the river scale requires that one realistically account for spatial and temporal variability across scale. Our down- scaling procedure includes the use of fixed/calibrated hydrological flow and temperature models coupled with a stochastically parameterized sturgeon bioenergetics model. We show that, although there is a large amount of uncertainty associated with both the climate model output and the fish growth process, one can establish significant differences in fish growth distributions between models, and between future and current climates for a given model.

Missouri River basin↗
USGS2015

Integrated Environmental Modelling: Human decisions, human challenges

Integrated Environmental Modelling (IEM) is an invaluable tool for understanding the complex, dynamic ecosystems that house our natural resources and control our environments. Human behaviour affects the ways in which the science of IEM is assembled and used for meaningful societal applications. In particular, human biases and heuristics reflect adaptation and experiential learning to issues with frequent, sharply distinguished, feedbacks. Unfortunately, human behaviour is not adapted to the more diffusely experienced problems that IEM typically seeks to address. Twelve biases are identified that affect IEM (and science in general). These biases are supported by personal observations and by the findings of behavioural scientists. A process for critical analysis is proposed that addresses some human challenges of IEM and solicits explicit description of (1) represented processes and information, (2) unrepresented processes and information, and (3) accounting for, and cognizance of, potential human biases. Several other suggestions are also made that generally complement maintaining attitudes of watchful humility, open-mindedness, honesty and transparent accountability. These suggestions include (1) creating a new area of study in the behavioural biogeosciences, (2) using structured processes for engaging the modelling and stakeholder communities in IEM, and (3) using ‘red teams’ to increase resilience of IEM constructs and use.

Geological Society of London Special Publications↗
USGS2015

Operational thermal remote sensing and lava flow monitoring at the Hawaiian Volcano Observatory

Hawaiian volcanoes are highly accessible and well monitored by ground instruments. Nevertheless, observational gaps remain and thermal satellite imagery has proven useful in Hawai‘i for providing synoptic views of activity during intervals between field visits. Here we describe the beginning of a thermal remote sensing programme at the US Geological Survey Hawaiian Volcano Observatory (HVO). Whereas expensive receiving stations have been traditionally required to achieve rapid downloading of satellite data, we exploit free, low-latency data sources on the internet for timely access to GOES, MODIS, ASTER and EO-1 ALI imagery. Automated scripts at the observatory download these data and provide a basic display of the images. Satellite data have been extremely useful for monitoring the ongoing lava flow activity on Kīlauea's East Rift Zone at Pu‘u ‘Ō‘ō over the past few years. A recent lava flow, named Kahauale‘a 2, was upslope from residential subdivisions for over a year. Satellite data helped track the slow advance of the flow and contributed to hazard assessments. Ongoing improvement to thermal remote sensing at HVO incorporates automated hotspot detection, effusion rate estimation and lava flow forecasting, as has been done in Italy. These improvements should be useful for monitoring future activity on Mauna Loa.

Hawaii↗
USGS2016

New insights into the Kawah Ijen hydrothermal system from geophysical data

Volcanoes with crater lakes and/or extensive hydrothermal systems pose significant challenges with respect to monitoring and forecasting eruptions, but they also provide new opportunities to enhance our understanding of magmatic–hydrothermal processes. Their lakes and hydrothermal systems serve as reservoirs for magmatic heat and fluid emissions, filtering and delaying the surface expressions of magmatic unrest and eruption, yet they also enable sampling and monitoring of geochemical tracers. Here, we describe the outcomes of a highly focused international experimental campaign and workshop carried out at Kawah Ijen volcano, Indonesia, in September 2014, designed to answer fundamental questions about how to improve monitoring and eruption forecasting at wet volcanoes.

Kawah Ijen↗
USGS2017

A global probabilistic tsunami hazard assessment from earthquake sources

Large tsunamis occur infrequently but have the capacity to cause enormous numbers of casualties, damage to the built environment and critical infrastructure, and economic losses. A sound understanding of tsunami hazard is required to underpin management of these risks, and while tsunami hazard assessments are typically conducted at regional or local scales, globally consistent assessments are required to support international disaster risk reduction efforts, and can serve as a reference for local and regional studies. This study presents a global-scale probabilistic tsunami hazard assessment (PTHA), extending previous global-scale assessments based largely on scenario analysis. Only earthquake sources are considered, as they represent about 80% of the recorded damaging tsunami events. Globally extensive estimates of tsunami run-up height are derived at various exceedance rates, and the associated uncertainties are quantified. Epistemic uncertainties in the exceedance rates of large earthquakes often lead to large uncertainties in tsunami run-up. Deviations between modelled tsunami run-up and event observations are quantified, and found to be larger than suggested in previous studies. Accounting for these deviations in PTHA is important, as it leads to a pronounced increase in predicted tsunami run-up for a given exceedance rate.

Geological Society, London, Special Publications↗
USGS2017

The application of microtextural and heavy mineral analysis to discriminate between storm and tsunami deposits

Recent work has applied microtextural and heavy mineral analyses to sandy storm and tsunami deposits from Portugal, Scotland, Indonesia and the USA. We looked at the interpretation of microtextural imagery (scanning electron microscopy) of quartz grains and heavy mineral compositions. We consider inundation events of different chronologies and sources (the AD 1755 Lisbon and 2004 Indian Ocean tsunamis, the Great Storm of 11 January 2005 in Scotland, and Hurricane Sandy in 2012) that affected contrasting coastal and hinterland settings with different regional oceanographic conditions. Storm and tsunami deposits were examined along with potential source sediments (alluvial, beach, dune and nearshore sediments) to determine provenance. Results suggest that tsunami deposits typically exhibit a significant spatial variation in grain sizes, microtextures and heavy minerals. Storm deposits show less variability, especially in vertical profiles. Tsunami and storm quartz grains had more percussion marks and fresh surfaces compared to potential source material. Moreover, in the studied cases, tsunami samples had fewer fresh surfaces than storm deposits. Heavy mineral assemblages are typically site-specific. The concentration of heavy minerals decreases upwards in tsunamigenic units, whereas storm sediments show cyclic concentrations of heavy minerals, reflected in the laminations observed macroscopically in the deposits.

Geological Society, London, Special Publications↗
USGS2017

The formation of gullies on Mars today

A decade of high-resolution monitoring has revealed extensive activity in fresh Martian gullies. Flows within the gullies are diverse: they can be relatively light, neutral or dark, colourful or bland, and range from superficial deposits to 10 m-scale topographic changes. We observed erosion and transport of material within gullies, new terraces, freshly eroded channel segments, migrating sinuous curves, channel abandonment, and lobate deposits. We also observed early stages of gully initiation, demonstrating that these processes are not merely modifying pre-existing landforms. The timing of activity closely correlates with the presence of seasonal CO 2 frost, so the current changes must be part of ongoing gully formation that is driven largely by its presence. We suggest that the cumulative effect of many flows erodes alcoves and channels, and builds lobate aprons, with no involvement of liquid water. Instead, flows may be fluidized by sublimation of entrained CO 2 ice or other mechanisms. The frequent activity is likely to have erased any features dating from high-obliquity periods, so fresh gully geomorphology at middle and high latitudes is not evidence for past liquid water. CO 2 ice-driven processes may have been important throughout Martian geological history and their deposits could exist in the rock record, perhaps resembling debris-flow sediments.

Geological Society, London, Special Publications↗