USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Geophysical Research - Atmospheres”

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

Monte Carlo method for determining earthquake recurrence parameters from short paleoseismic catalogs: Example calculations for California

Paleoearthquake observations often lack enough events at a given site to directly define a probability density function (PDF) for earthquake recurrence. Sites with fewer than 10-15 intervals do not provide enough information to reliably determine the shape of the PDF using standard maximum-likelihood techniques (e.g., Ellsworth et al., 1999). In this paper I present a method that attempts to fit wide ranges of distribution parameters to short paleoseismic series. From repeated Monte Carlo draws, it becomes possible to quantitatively estimate most likely recurrence PDF parameters, and a ranked distribution of parameters is returned that can be used to assess uncertainties in hazard calculations. In tests on short synthetic earthquake series, the method gives results that cluster around the mean of the input distribution, whereas maximum likelihood methods return the sample means (e.g., NIST/SEMATECH, 2006). For short series (fewer than 10 intervals), sample means tend to reflect the median of an asymmetric recurrence distribution, possibly leading to an overestimate of the hazard should they be used in probability calculations. Therefore a Monte Carlo approach may be useful for assessing recurrence from limited paleoearthquake records. Further, the degree of functional dependence among parameters like mean recurrence interval and coefficient of variation can be established. The method is described for use with time-independent and time-dependent PDFs, and results from 19 paleoseismic sequences on strike-slip faults throughout the state of California are given.

Journal of Geophysical Research B: Solid Earth↗
USGS2008

Finding concealed active faults: Extending the southern Whidbey Island fault across the Puget Lowland, Washington

The southern Whidbey Island fault zone (SWIF), as previously mapped using borehole data, potential field anomalies, and marine seismic reflection surveys, consists of three subparallel, northwest trending strands extending ∼100 km from near Vancouver Island to the northern Puget Lowland. East of Puget Sound, the SWIF makes landfall between the cities of Seattle and Everett but is concealed beneath a thick mantle of young glacial deposits and vegetation. A ∼20‐km‐wide, northwest trending swath of subparallel, low‐amplitude aeromagnetic anomalies crosses this region of the Puget Lowland and is on strike with the SWIF. The most prominent aeromagnetic anomaly, the Cottage Lake lineament, extends at least 18 km and lies approximately on strike with the SWIF on Whidbey Island. Subtle scarps and topographic lineaments on Pleistocene surfaces, visible on high‐resolution lidar topography at a number of locations along the SWIF, lie on or near these magnetic anomalies. In the field, scarps exhibit northeast‐side‐up and vertical relief of 1 to 5 m. Excavations across several lidar scarps lying on or near magnetic anomalies show evidence for multiple folding and faulting events since deglaciation, most likely above buried reverse/oblique faults. Excavations in areas away from magnetic anomalies do not show evidence of tectonic deformation. In total, paleoseismological evidence suggests that the SWIF produced at least four earthquakes since deglaciation about 16,400 years ago, the most recent less than 2700 years ago.

Washington↗
USGS2008

Variability of passive gas emissions, seismicity, and deformation during crater lake growth at White Island Volcano, New Zealand, 2002-2006

We report on 4 years of airborne measurements of CO2, SO2, and H2S emission rates during a quiescent period at White Island volcano, New Zealand, beginning in 2003. During this time a significant crater lake emerged, allowing scrubbig processes to be investigated. CO2 emissions varied from a baseline of 250 to >2000 t d-1 and demonstrated clear annual cycling that was consistent with numbers of earthquake detections and annual changes in sea level. The annual variability was found to be most likely related to increases in the strain on the volcano during sea level highs, temporarily causing fractures to reduce in size in the upper conduit. SO2 emissions varied from 0 to >400 t d-1 and were clearly affected by scrubbing processes within the first year of take development. Scrubbing caused increases of SO42- and Cl- in lake waters, and the ratio of carbon to total sulphur suggested that elemental sulphur deposition was also significant in the lake during the first year. Careful measurements of the lake level and chemistry allowed estimates of the rate of H2O(g) and HCl(g) input into the lake and suggested that the molar abundances of major gas species (H2O, CO2, SO2, and HCl) during this quiescent phase were similar to fumarolic ratios observed between earlier eruptive periods. The volume of magma estimated from CO2 emissions (0.0 15-0.04 km3) was validated by Cl- increases in the lake, suggesting that the gas and magma are transported from deep to shallow depths as a closed system and likely become open in the upper conduit region. The absence of surface deformation further leads to a necessity of magma convection to supply and remove magma from the degassing depths. Two models of convection configurations are discussed. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2008

Effects of supershear rupture speed on the high-frequency content of S waves investigated using spontaneous dynamic rupture models and isochrone theory

In this paper we achieve three goals: (1) We demonstrate that crack tips governed by friction laws, including slip weakening, rate- and state-dependent laws, and thermal pressurization of pore fluids, propagating at supershear speed have slip velocity functions with reduced high-frequency content compared to crack tips traveling at subshear speeds. This is demonstrated using a fully dynamic, spontaneous, three-dimensional earthquake model, in which we calculate fault slip velocity at nine points (locations) distributed along a quarter circle on the fault where the rupture is traveling at supershear speed in the in-plane direction and subshear speed in the antiplane direction. This holds for a fault governed by the linear slip-weakening constitutive equation, by slip weakening with thermal pressurization of pore fluid, and by rate- and state-dependent laws with thermal pressurization. The same is also true even assuming a highly heterogeneous initial shear stress field on the fault. (2) Using isochrone theory, we derive a general expression for the spectral characteristics and geometric spreading of two pulses arising from supershear rupture, the well-known Mach wave, and a second lesser known pulse caused by rupture acceleration. (3) We demonstrate that the Mach cone amplification of high frequencies overwhelms the de-amplification of high-frequency content in the slip velocity functions in supershear ruptures. Consequently, when earthquake ruptures travel at supershear speed, a net enhancement of high-frequency radiation is expected, and the alleged “low” peak accelerations observed for the 2002 Denali and other large earthquakes are probably not caused by diminished high-frequency content in the slip velocity function, as has been speculated.

Journal of Geophysical Research Solid Earth↗
USGS2008

Dislocation models of interseismic deformation in the western United States

The GPS-derived crustal velocity field of the western United States is used to construct dislocation models in a viscoelastic medium of interseismic crustal deformation. The interseismic velocity field is constrained by 1052 GPS velocity vectors spanning the ???2500-km-long plate boundary zone adjacent to the San Andreas fault and Cascadia subduction zone and extending ???1000 km into the plate interior. The GPS data set is compiled from U.S. Geological Survey campaign data, Plate Boundary Observatory data, and the Western U.S. Cordillera velocity field of Bennett et al. (1999). In the context of viscoelastic cycle models of postearthquake deformation, the interseismic velocity field is modeled with a combination of earthquake sources on ???100 known faults plus broadly distributed sources. Models that best explain the observed interseismic velocity field include the contributions of viscoelastic relaxation from faulting near the major plate margins, viscoelastic relaxation from distributed faulting in the plate interior, as well as lateral variations in depth-averaged rigidity in the elastic lithosphere. Resulting rigidity variations are consistent with reduced effective elastic plate thickness in a zone a few tens of kilometers wide surrounding the San Andreas fault (SAF) system. Primary deformation characteristics are captured along the entire SAF system, Eastern California Shear Zone, Walker Lane, the Mendocino triple junction, the Cascadia margin, and the plate interior up to ???1000 km from the major plate boundaries.

Journal of Geophysical Research B: Solid Earth↗
USGS2008

In situ Raman spectroscopic investigation of the structure of subduction-zone fluids

In situ Raman spectra of synthetic subduction-zone fluids (KAlSi3O8-H2O system) were measured to 900?? and 2.3 GPa using a hydrothermal diamond-anvil cell. The structures of aqueous fluid and hydrous melt become closer when conditions approach the second critical endpoint. Almost no three-dimensional network was observed in the supercritical fluid above 2 GPa although a large amount of silicate component is dissolved, suggesting that the physical and chemical properties of these phases change drastically at around the second critical endpoint. Our experimental results indicate that the fluids released from a subducting slab change from aqueous fluid to supercritical fluid with increasing depth under the volcanic arcs. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2008

Probabilistic seismic hazard in the San Francisco Bay area based on a simplified viscoelastic cycle model of fault interactions

We construct a viscoelastic cycle model of plate boundary deformation that includes the effect of time-dependent interseismic strain accumulation, coseismic strain release, and viscoelastic relaxation of the substrate beneath the seismogenic crust. For a given fault system, time-averaged stress changes at any point (not on a fault) are constrained to zero; that is, kinematic consistency is enforced for the fault system. The dates of last rupture, mean recurrence times, and the slip distributions of the (assumed) repeating ruptures are key inputs into the viscoelastic cycle model. This simple formulation allows construction of stress evolution at all points in the plate boundary zone for purposes of probabilistic seismic hazard analysis (PSHA). Stress evolution is combined with a Coulomb failure stress threshold at representative points on the fault segments to estimate the times of their respective future ruptures. In our PSHA we consider uncertainties in a four-dimensional parameter space: the rupture peridocities, slip distributions, time of last earthquake (for prehistoric ruptures) and Coulomb failure stress thresholds. We apply this methodology to the San Francisco Bay region using a recently determined fault chronology of area faults. Assuming single-segment rupture scenarios, we find that fature rupture probabilities of area faults in the coming decades are the highest for the southern Hayward, Rodgers Creek, and northern Calaveras faults. This conclusion is qualitatively similar to that of Working Group on California Earthquake Probabilities, but the probabilities derived here are significantly higher. Given that fault rupture probabilities are highly model-dependent, no single model should be used to assess to time-dependent rupture probabilities. We suggest that several models, including the present one, be used in a comprehensive PSHA methodology, as was done by Working Group on California Earthquake Probabilities.

Journal of Geophysical Research B: Solid Earth↗
USGS2008

Noise in GPS displacement measurements from southern California and southern Nevada

Time series of position changes estimated from data from 236 continuously recording GPS receivers operating in Southern California and Southern Nevada are evaluated for noise models that characterize their temporal correlations. The lengths of the time series range between 3.5 and 10 years. After adjusting these data for postseismic deformation, offsets, and annual periodicities, I find that about one-half of the time series have temporal correlations that are categorized as either flicker or random-walk noise. The remaining time series can be best categorized as either a combination of flicker and random-walk; power law noise; first-order Gauss-Markov plus random-walk noise; or power law plus broadband, seasonal noise. A variety of geodetic monuments are used in Southern California and Nevada, including deeply braced designs, cement piers, pins drilled in outcrop, and buildings. When I evaluate the noise for each time series in terms of an estimate of the standard error in velocity, I find that the sites with the smallest errors are those located in Nevada using deeply braced monuments. Sites that are installed within regions of active pumping, both for groundwater and oil, had the largest standard errors in velocity. Comparison of monument stability, as measured by standard error in rate, with average, annual rainfall nearby indicates a marginally significant correlation. In addition, even though regional filtering removed much of the common-mode signals in these time series, there still remains a common-mode seasonal signal which can and should be removed.

California, Nevada↗
USGS2008

Eruption dynamics at Mount St. Helens imaged from broadband seismic waveforms: Interaction of the shallow magmatic and hydrothermal systems

The current eruption at Mount St. Helens is characterized by dome building and shallow, repetitive, long-period (LP) earthquakes. Waveform cross-correlation reveals remarkable similarity for a majority of the earthquakes over periods of several weeks. Stacked spectra of these events display multiple peaks between 0.5 and 2 Hz that are common to most stations. Lower-amplitude very-long-period (VLP) events commonly accompany the LP events. We model the source mechanisms of LP and VLP events in the 0.5-4 s and 8-40 s bands, respectively, using data recorded in July 2005 with a 19-station temporary broadband network. The source mechanism of the LP events includes: 1) a volumetric component modeled as resonance of a gently NNW-dipping, steam-filled crack located directly beneath the actively extruding part of the new dome and within 100 m of the crater floor and 2) a vertical single force attributed to movement of the overlying dome. The VLP source, which also includes volumetric and single-force components, is 250 m deeper and NNW of the LP source, at the SW edge of the 1980s lava dome. The volumetric component points to the compression and expansion of a shallow, magma-filled sill, which is subparallel to the hydrothermal crack imaged at the LP source, coupled with a smaller component of expansion and compression of a dike. The single-force components are due to mass advection in the magma conduit. The location, geometry and timing of the sources suggest the VLP and LP events are caused by perturbations of a common crack system.

Journal of Geophysical Research B: Solid Earth↗
USGS2008

A global search for stress shadows

[1] Debate continues regarding the relative proportion of earthquakes triggered by passing seismic waves versus static stress changes from a main shock. Static stress changes are expected to have long‐term effects on earthquake probabilities, whereas dynamic stress changes due to the passing of seismic waves should not. Both mechanisms are expected to raise seismicity rates in some areas, but only static stress change calculations predict rate decrease shadows. Thus, identification of post‐main‐shock earthquake suppression is diagnostic of a static stress change process. We note that in principle, static stress change theory predicts suppression of particular earthquake mechanisms in a shadow zone rather than an overall rate reduction. A stress shadow can therefore be characterized by a change in the average earthquake focal mechanism before and after a main shock that results from suppression of a given mechanism type. We examined average mechanisms from ±2° radii and 5‐year periods before and after 119 M s ≥ 7 main shock earthquakes drawn from the Harvard Centroid Moment Tensor (CMT) catalog. Significant average mechanism changes caused by earthquake suppression were found in only two cases. However, by stacking the data, we were able to resolve statistically significant suppression of particular post‐main‐shock focal mechanisms. This indicates that, while static stress shadows are subtle, they are indeed present in the global catalog.

Journal of Geophysical Research Solid Earth↗
USGS2008

Exponential decline of aftershocks of the M7.9 1868 great Kau earthquake, Hawaii, through the 20th century

The remarkable catalog of Hawaiian earthquakes going back to the 1820s is based on missionary diaries, newspaper accounts, and instrumental records and spans the great M 7.9 Kau earthquake of April 1868 and its aftershock sequence. The earthquake record since 1868 defines a smooth curve complete to M 5.2 of the declining rate into the 21st century, after five short volcanic swarms are removed. A single aftershock curve fits the earthquake record, even with numerous M 6 and 7 main shocks and eruptions. The timing of some moderate earthquakes may be controlled by magmatic stresses, but their overall long‐term rate reflects one of aftershocks of the Kau earthquake. The 1868 earthquake is, therefore, the largest and most controlling stress event in the 19th and 20th centuries. We fit both the modified Omori (power law) and stretched exponential (SE) functions to the earthquakes. We found that the modified Omori law is a good fit to the M ≥ 5.2 earthquake rate for the first 10 years or so and the more rapidly declining SE function fits better thereafter, as supported by three statistical tests. The switch to exponential decay suggests that a possible change in aftershock physics may occur from rate and state fault friction, with no change in the stress rate, to viscoelastic stress relaxation. The 61‐year exponential decay constant is at the upper end of the range of geodetic relaxation times seen after other global earthquakes. Modeling deformation in Hawaii is beyond the scope of this paper, but a simple interpretation of the decay suggests an effective viscosity of 10 19 to 10 20 Pa s pertains in the volcanic spreading of Hawaii's flanks. The rapid decline in earthquake rate poses questions for seismic hazard estimates in an area that is cited as one of the most hazardous in the United States.

Hawaii↗
USGS2008

Diffuse flow hydrothermal manganese mineralization along the active Mariana and southern Izu-Bonin arc system, western Pacific

Abundant ferromanganese oxides were collected along 1200 km of the active Izu-Bonin-Mariana arc system. Chemical compositions and mineralogy show that samples were collected from two deposit types: Fe-Mn crusts of mixed hydrogenetic/hydrothermal origin and hydrothermal Mn oxide deposits; this paper addresses only the second type. Mn oxides cement volcaniclastic and biogenic sandstone and breccia layers (Mn sandstone) and form discrete dense stratabound layers along bedding planes and within beds (stratabound Mn). The Mn oxide was deposited within coarse-grained sediments from diffuse flow systems where precipitation occurred below the seafloor. Deposits were exposed at the seabed by faulting, mass wasting, and erosion. Scanning electron microscopy and microprobe analyses indicate the presence of both amorphous and crystalline 10 ?? and 7 ?? manganate minerals, the fundamental chemical difference being high water contents in the amorphous Mn oxides. Alternation of amorphous and crystalline laminae occurs in many samples, which likely resulted from initial rapid precipitation of amorphous Mn oxides from waxing pulses of hydrothermal fluids followed by precipitation of slow forming crystallites during waning stages. The chemical composition is characteristic of a hydrothermal origin including strong fractionation between Fe (mean 0.9 wt %) and Mn (mean 48 wt %) for the stratabound Mn, generally low trace metal contents, and very low rare earth element and platinum group element contents. However, Mo, Cd, Zn, Cu, Ni, and Co occur in high concentrations in some samples and may be good indicator elements for proximity to the heat source or to massive sulfide deposits. For the Mn sandstones, Fe (mean-8.4%) and Mn (12.4%) are not significantly fractionated because of high Fe contents in the volcaniclastic material. However, the proportion of hydrothermal Fe (nondetrital Fe) to total Fe is remarkably constant (49-58%) for all the sample groups, regardless of the degree of Mn mineralization. Factor analyses indicate various mixtures of two dominant components: hydrothermal Mn oxide for the stratabound Mn and detrital aluminosilicate for the Mn-cemented sandstone; and two minor components, hydrothermal Fe oxyhydroxide and biocarbonate/biosilica. Our conceptual model shows that Mn mineralization was produced by hydrothermal convection cells within arc volcanoes and sedimentary prisms that occur along, the flanks and within calderas. The main source of hydrothermal fluid was seawater that penetrated through fractures, faults, and permeable volcanic edifices. The fluids were heated by magma, enriched in metals by leaching of basement rocks and sediments, and mixed with magmatic fluids and gases. Dikes and sills may have been another source of heat that drove small-scale circulation within sedimentary prisms. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2008

Observations and interpretation of fundamental mode Rayleigh wavefields recorded by the Transportable Array (USArray)

Broadband recordings of the dense Transportable Array (TA) in the western United States provide unparalleled detailed images of long-period seismic surface wavefields. With 400 stations spanning most of the western United States, wavefronts of fundamental mode Rayleigh waves may be visualized coherently across the array at periods ???40 s. In order to constrain the Rayleigh wave phase velocity structure in the western United States, I assemble a data set of vertical component seismograms from 53 teleseismic events recorded by the TA from April 2006 to October 2007. Complex amplitude spectra from these recordings at peni ods 27-100 s are interpreted using the multiplane wave tomographic method of Friederich and Wielandt (1995) and Pollitz (1999). This analysis yields detailed surface wave phase velocity and three-dimensional shear wave velocity patterns across the North American plate boundary zone, elucidating the active processes in the highly heterogeneous western U.S. upper mantle.

Journal of Geophysical Research B: Solid Earth↗
USGS2007

Tidal asymmetry and residual circulation over linear sandbanks and their implication on sediment transport: a process-oriented numerical study

A series of process-oriented numerical simulations is carried out in order to evaluate the relative role of locally generated residual flow and overtides on net sediment transport over linear sandbanks. The idealized bathymetry and forcing are similar to those present in the Norfolk Sandbanks, North Sea. The importance of bottom drag parameterization and bank orientation with respect to the ambient flow is examined in terms of residual flow and overtide generation, and subsequent sediment transport implications are discussed. The results show that although the magnitudes of residual flow and overtides are sensitive to bottom roughness parameterization and bank orientation, the magnitude of the generated residual flow is always larger than that of the locally generated overtides. Also, net sediment transport is always dominated by the nonlinear interaction of the residual flow and the semidiurnal tidal currents, although cross-bank sediment transport can occur even in the absence of a cross-shore residual flow. On the other hand, net sediment divergence/convergence increases as the bottom drag decreases and as bank orientation increases. The sediment erosion/deposition is not symmetric about the crest of the bank, suggesting that originally symmetric banks would have the tendency to become asymmetric.

Journal of Geophysical Research C: Oceans↗
USGS2007

Role of aquifer heterogeneity in fresh groundwater discharge and seawater recycling: An example from the Carmel coast, Israel

A case study is shown in which the pattern of submarine groundwater discharge and of seawater recycling is controlled by local hydrogeological variability. The coastal aquifer in Dor Bay is composed of two units: a partly confined calcaranitic sandstone (Kurkar) and an overlying loose sand. Groundwater in the Kurkar has elevated activities of 222 Rn (∼390 dpm/L) and relatively low 224 Ra/ 223 Ra activity ratios (3–4), while the sand groundwater is significantly less radiogenic (6–90 dpm/L) and shows higher 224 Ra/ 223 Ra ratios. Groundwater discharging from sand-covered areas of the bay has salinities of 16–31 and an average 222 Rn activity of 168 dpm/L, which lies on a mixing line between Rn-rich Kurkar fresh water and Rn-poor seawater. Another key observation is that seawater infiltrates to some extent into onshore sand groundwater, while the fresh water within the submarine Kurkar can be traced up to 40 m offshore. This implies that while fresh water mainly discharges from the Kurkar unit, seawater recycling is limited to the loose sand, and that the discharge from sand-covered areas is a mixture of Kurkar water with recycled seawater. Advection rates from the bay floor were calculated from Rn time series and found to vary between 0 and 36 cm/d, correlating negatively with bay water depth. The average flux was 8.1 cm/d, and it did not seem to change much during March, May, and July 2006. The average amount of fresh water discharging to the bay was 5.0 m 3 /d per meter of shoreline. Radon activity in the sand groundwater also fluctuates due to influx of Kurkar-type groundwater.

Carmel Coast, Dor Bay↗
USGS2008

Variational analysis of drifter positions and model outputs for the reconstruction of surface currents in the central Adriatic during fall 2002

In this paper we present an application of a variational method for the reconstruction of the velocity field in a coastal flow in the central Adriatic Sea, using in situ data from surface drifters and outputs from the ROMS circulation model. The variational approach, previously developed and tested for mesoscale open ocean flows, has been improved and adapted to account for inhomogeneities on boundary current dynamics over complex bathymetry and coastline and for weak Lagrangian persistence in coastal flows. The velocity reconstruction is performed using nine drifter trajectories over 45 d, and a hierarchy of indirect tests is introduced to evaluate the results as the real ocean state is not known. For internal consistency and impact of the analysis, three diagnostics characterizing the particle prediction and transport, in terms of residence times in various zones and export rates from the boundary current toward the interior, show that the reconstruction is quite effective. A qualitative comparison with sea color data from the MODIS satellite images shows that the reconstruction significantly improves the description of the boundary current with respect to the ROMS model first guess, capturing its main features and its exchanges with the interior when sampled by the drifters. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research C: Oceans↗
USGS2007

Generation and propagation of nonlinear internal waves in Massachusetts Bay

During the summer, nonlinear internal waves (NLIWs) are commonly observed propagating in Massachusetts Bay. The topography of the area is unique in the sense that the generation area (over Stellwagen Bank) is only 25 km away from the shoaling area, and thus it represents an excellent natural laboratory to study the life cycle of NLIWs. To assist in the interpretation of the data collected during the 1998 Massachusetts Bay Internal Wave Experiment (MBIWE98), a fully nonlinear and nonhydrostatic model covering the generation/shoaling region was developed, to investigate the response of the system to the range of background and driving conditions observed. Simplified models were also used to elucidate the role of nonlinearity and dispersion in shaping the NLIW field. This paper concentrates on the generation process and the subsequent evolution in the basin. The model was found to reproduce well the range of propagation characteristics observed (arrival time, propagation speed, amplitude), and provided a coherent framework to interpret the observations. Comparison with a fully nonlinear hydrostatic model shows that during the generation and initial evolution of the waves as they move away from Stellwagen Bank, dispersive effects play a negligible role. Thus the problem can be well understood considering the geometry of the characteristics along which the Riemann invariants of the hydrostatic problem propagate. Dispersion plays a role only during the evolution of the undular bore in the middle of Stellwagen Basin. The consequences for modeling NLIWs within hydrostatic models are briefly discussed at the end.

Massachusetts↗
USGS2008

Geomorphic evidence for former lobate debris aprons at low latitudes on Mars: Indicators of the Martian paleoclimate

Circumferential depressions enclosing mesas and plateaus in the northern Kasei Valles and in the Tartarus Colles regions of Mars are interpreted as indicators of the former extent of lobate debris aprons, thought to be mixtures of ice and elastic particles. These former lobate debris aprons existed about 1 Ga ago and were embayed by lavas or other flow deposits. After the lobate debris aprons had been removed by sublimation and deflation, topographic depressions with a depth of 50 m and a width of several kilometers were left behind between the mesa or plateau scarp and the solidified flow materials. These depressions or moats are located equatorward of ??30?? at significantly lower latitudes than generally observed for occurrences of modem, intact lobate debris aprons. This observation provides evidence that the paleoclimate at that time was different than today, probably due to a higher averaged obliquity of the planet's rotational axis. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research E: Planets↗