SEARCH · USGS Science

Results for “Geophysical Journal”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

SourceResearch connectionFrom yearThrough year

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

At least 955 records · Page 53Linked to original sources

USGS2004

A physical model for strain accumulation in the San Francisco Bay region: Stress evolution since 1838

Understanding of the behavior of plate boundary zones has progressed to the point where reasonably comprehensive physical models can predict their evolution. The San Andreas fault system in the San Francisco Bay region (SFBR) is dominated by a few major faults whose behavior over about one earthquake cycle is fairly well understood. By combining the past history of large ruptures on SFBR faults with a recently proposed physical model of strain accumulation in the SFBR, we derive the evolution of regional stress from 1838 until the present. This effort depends on (1) an existing compilation of the source properties of historic and contemporary SFBR earthquakes based on documented shaking, geodetic data, and seismic data (Bakun, 1999) and (2) a few key parameters of a simple regional viscoelastic coupling model constrained by recent GPS data (Pollitz and Nyst, 2004). Although uncertainties abound in the location, magnitude, and fault geometries of historic ruptures and the physical model relies on gross simplifications, the resulting stress evolution model is sufficiently detailed to provide a useful window into the past stress history. In the framework of Coulomb failure stress, we find that virtually all M ??? 5.8 earthquakes prior to 1906 and M ??? 5.5 earthquakes after 1906 are consistent with stress triggering from previous earthquakes. These events systematically lie in zones of predicted stress concentration elevated 5-10 bars above the regional average. The SFBR is predicted to have emerged from the 1906 "shadow" in about 1980, consistent with the acceleration in regional seismicity at that time. The stress evolution model may be a reliable indicator of the most likely areas to experience M ??? 5.5 shocks in the future.

Journal of Geophysical Research B: Solid Earth↗

USGS2005

Interferometric synthetic aperture radar study of Okmok volcano, Alaska, 1992-2003: Magma supply dynamics and postemplacement lava flow deformation

Okmok volcano, located in the central Aleutian arc, Alaska, is a dominantly basaltic complex topped with a 10-km-wide caldera that formed circa 2.05 ka. Okmok erupted several times during the 20th century, most recently in 1997; eruptions in 1945, 1958, and 1997 produced lava flows within the caldera. We used 80 interferometric synthetic aperture radar (InSAR) images (interferograms) to study transient deformation of the volcano before, during, and after the 1997 eruption. Point source models suggest that a magma reservoir at a depth of 3.2 km below sea level, located beneath the center of the caldera and about 5 km northeast of the 1997 vent, is responsible for observed volcano-wide deformation. The preeruption uplift rate decreased from about 10 cm yr −1 during 1992–1993 to 2 ∼ 3 cm yr −1 during 1993–1995 and then to about −1 ∼ −2 cm yr −1 during 1995–1996. The posteruption inflation rate generally decreased with time during 1997–2001, but increased significantly during 2001–2003. By the summer of 2003, 30 ∼ 60% of the magma volume lost from the reservoir in the 1997 eruption had been replenished. Interferograms for periods before the 1997 eruption indicate consistent subsidence of the surface of the 1958 lava flows, most likely due to thermal contraction. Interferograms for periods after the eruption suggest at least four distinct deformation processes: (1) volcano-wide inflation due to replenishment of the shallow magma reservoir, (2) subsidence of the 1997 lava flows, most likely due to thermal contraction, (3) deformation of the 1958 lava flows due to loading by the 1997 flows, and (4) continuing subsidence of 1958 lava flows buried beneath 1997 flows. Our results provide insights into the postemplacement behavior of lava flows and have cautionary implications for the interpretation of inflation patterns at active volcanoes.

Alaska↗

USGS2005

Significance of stress transfer in time-dependent earthquake probability calculations

A sudden change in stress is seen to modify earthquake rates, but should it also revise earthquake probability? Data used to derive input parameters permits an array of forecasts; so how large a static stress change is require to cause a statistically significant earthquake probability change? To answer that question, effects of parameter and philosophical choices are examined through all phases of sample calculations, Drawing at random from distributions of recurrence-aperiodicity pairs identifies many that recreate long paleoseismic and historic earthquake catalogs. Probability density funtions built from the recurrence-aperiodicity pairs give the range of possible earthquake forecasts under a point process renewal model. Consequences of choices made in stress transfer calculations, such as different slip models, fault rake, dip, and friction are, tracked. For interactions among large faults, calculated peak stress changes may be localized, with most of the receiving fault area changed less than the mean. Thus, to avoid overstating probability change on segments, stress change values should be drawn from a distribution reflecting the spatial pattern rather than using the segment mean. Disparity resulting from interaction probability methodology is also examined. For a fault with a well-understood earthquake history, a minimum stress change to stressing rate ratio of 10:1 to 20:1 is required to significantly skew probabilities with >80-85% confidence. That ratio must be closer to 50:1 to exceed 90-95% confidence levels. Thus revision to earthquake probability is achievable when a perturbing event is very close to the fault in question or the tectonic stressing rate is low.

Journal of Geophysical Research B: Solid Earth↗

USGS2005

Rupture dynamics with energy loss outside the slip zone

Energy loss in a fault damage zone, outside the slip zone, contributes to the fracture energy that determines rupture velocity of an earthquake. A nonelastic two-dimensional dynamic calculation is done in which the slip zone is modeled as a fault plane and material off the fault is subject to a Coulomb yield condition. In a mode 2 crack-like solution in which an abrupt uniform drop of shear traction on the fault spreads from a point, Coulomb yielding occurs on the extensional side of the fault. Plastic strain is distributed with uniform magnitude along the fault, and it has a thickness normal to the fault proportional to propagation distance. Energy loss off the fault is also proportional to propagation distance, and it can become much larger than energy loss on the fault specified by the fault constitutive relation. The slip velocity function could be produced in an equivalent elastic problem by a slip-weakening friction law with breakdown slip Dc increasing with distance. Fracture energy G and equivalent Dc will be different in ruptures with different initiation points and stress drops, so they are not constitutive properties; they are determined by the dynamic solution that arrives at a particular point. Peak slip velocity is, however, a property of a fault location. Nonelastic response can be mimicked by imposing a limit on slip velocity on a fault in an elastic medium.

Journal of Geophysical Research B: Solid Earth↗

USGS2005

New constraints on mechanisms of remotely triggered seismicity at Long Valley Caldera

Regional-scale triggering of local earthquakes in the crust by seismic waves from distant main shocks has now been robustly documented for over a decade. Some of the most thoroughly recorded examples of repeated triggering of a single site from multiple, large earthquakes are measured in geothermal fields of the western United States like Long Valley Caldera. As one of the few natural cases where the causality of an earthquake sequence is apparent, triggering provides fundamental constraints on the failure processes in earthquakes. We show here that the observed triggering by seismic waves is inconsistent with any mechanism that depends on cumulative shaking as measured by integrated energy density. We also present evidence for a frequency-dependent triggering threshold. On the basis of the seismic records of 12 regional and teleseismic events recorded at Long Valley Caldera, long-period waves (>30 s) are more effective at generating local seismicity than short-period waves of comparable amplitude. If the properties of the system are stationary over time, the failure threshold for long-period waves is ~0.05 cm/s vertical shaking. Assuming a phase velocity of 3.5 km/s and an elastic modulus of 3.5 x 10 10 Pa, the threshold in terms of stress is 5 kPa. The frequency dependence is due in part to the attenuation of the surface waves with depth. Fluid flow through a porous medium can produce the rest of the observed frequency dependence of the threshold. If the threshold is not stationary with time, pore pressures that are >99.5% of lithostatic and vary over time by a factor of 4 could explain the observations with no frequency dependence of the triggering threshold.

California↗

USGS2005

Evolution of melt-vapor surface tension in silicic volcanic systems: Experiments with hydrous melts

We evaluate the melt‐vapor surface tension (σ) of natural, water‐saturated dacite melt at 200 MPa, 950–1055°C, and 4.8–5.7 wt % H 2 O. We experimentally determine the critical supersaturation pressure for bubble nucleation as a function of dissolved water and then solve for σ at those conditions using classical nucleation theory. The solutions obtained give dacite melt‐vapor surface tensions that vary inversely with dissolved water from 0.042 (±0.003) J m −2 at 5.7 wt % H 2 O to 0.060 (±0.007) J m −2 at 5.2 wt % H 2 O to 0.073 (±0.003) J m −2 at 4.8 wt % H 2 O. Combining our dacite results with data from published hydrous haplogranite and high‐silica rhyolite experiments reveals that melt‐vapor surface tension also varies inversely with the concentration of mafic melt components (e.g., CaO, FeO total , MgO). We develop a thermodynamic context for these observations in which melt‐vapor surface tension is represented by a balance of work terms controlled by melt structure. Overall, our results suggest that cooling, crystallization, and vapor exsolution cause systematic changes in σ that should be considered in dynamic modeling of magmatic processes.

Journal of Geophysical Research B: Solid Earth↗

USGS2004

Stress orientations at intermediate angles to the San Andreas Fault, California

There are currently two competing models for the frictional strength of the San Andreas Fault in California: the strong-fault model and the weak-fault model. The strong-fault model predicts the maximum horizontal compressive stress axis to be at low angles to the fault, while the relatively weak fault model predicts it to be at high angles. Previous studies have disagreed as to which model is supported by observed stress orientations. We review and compare these studies and present results from several new focal mechanism stress inversions. We find that the observed stress orientations of different studies are generally consistent, implying that the disagreement is one of interpretation. The majority of studies find compressive stress orientations at intermediate angles to the fault, not strictly consistent with either current model. The strong-fault model is acceptable if the San Andreas is assumed to be a nonoptimally orientated fault that fails because optimally oriented, preexisting planes are not present. The relatively weak fault model is not consistent with the stress orientations. We propose two alternative models to better explain the observed intermediate stress orientations: an intermediate-strength San Andreas model and a model in which all major active faults are weak.

California↗

USGS2005

The calcite → aragonite transformation in low-Mg marble: Equilibrium relations, transformations mechanisms, and rates

Experimental transformation of a rather pure natural calcite marble to aragonite marble did not proceed via the expected straightforward polymorphic replacement. Instead, the small amount of Mg in the starting material (0.36 wt %) was excluded from the growing aragonite and diffused preferentially into the remaining calcite grains, producing Mg-rich calcite rods that persisted as relicts. Nucleation of aragonite occurred exclusively on grain boundaries, with aragonite [001] oriented subparallel to calcite [0001]. The aragonite crystals preferentially consumed the calcite crystal on which they nucleated, and the reaction fronts developed preferentially along the {010} and {110} planes of aragonite. Each aragonite neoblast that grew was nearly free of Mg (typically <0.1 wt %). The excess Mg was taken up by the calcite grains in between, stabilizing them and causing a few volume percent rodlike relicts of Mg-enriched calcite (up to 10 wt % MgO) to be left behind by the advancing reaction front. The aragonite growth rates are approximately linear and range from &sim;3 × 10 −11 m s −1 at 600°C to &sim;9 × 10 −9 m s −1 at 850°C, with an apparent activation enthalpy of 166 ± 91 kJ mol −1 . This reaction mechanism and the resultant texture are akin to cellular precipitation reactions in metals. Similar transformation textures have been reported from high-Mg marbles in Japan and China that disproportionated to low-Mg calcite and dolomite.

Journal of Geophysical Research B: Solid Earth↗

USGS2005

Magnitude and location of historical earthquakes in Japan and implications for the 1855 Ansei Edo earthquake

Japan Meteorological Agency (JMA) intensity assignments IJMA are used to derive intensity attenuation models suitable for estimating the location and an intensity magnitude Mjma for historical earthquakes in Japan. The intensity for shallow crustal earthquakes on Honshu is equal to -1.89 + 1.42MJMA - 0.00887?? h - 1.66log??h, where MJMA is the JMA magnitude, ??h = (??2 + h2)1/2, and ?? and h are epicentral distance and focal depth (km), respectively. Four earthquakes located near the Japan Trench were used to develop a subducting plate intensity attenuation model where intensity is equal to -8.33 + 2.19MJMA -0.00550??h - 1.14 log ?? h. The IJMA assignments for the MJMA7.9 great 1923 Kanto earthquake on the Philippine Sea-Eurasian plate interface are consistent with the subducting plate model; Using the subducting plate model and 226 IJMA IV-VI assignments, the location of the intensity center is 25 km north of the epicenter, Mjma is 7.7, and MJMA is 7.3-8.0 at the 1?? confidence level. Intensity assignments and reported aftershock activity for the enigmatic 11 November 1855 Ansei Edo earthquake are consistent with an MJMA 7.2 Philippine Sea-Eurasian interplate source or Philippine Sea intraslab source at about 30 km depth. If the 1855 earthquake was a Philippine Sea-Eurasian interplate event, the intensity center was adjacent to and downdip of the rupture area of the great 1923 Kanto earthquake, suggesting that the 1855 and 1923 events ruptured adjoining sections of the Philippine Sea-Eurasian plate interface.

Journal of Geophysical Research B: Solid Earth↗

USGS2006

Distribution of stress drop, stiffness, and fracture energy over earthquake rupture zones

Using information provided by slip models and the methodology of McGarr and Fletcher (2002), we map static stress drop, stiffness (k = ????/u, where ???? is static stress drop and u is slip), and fracture energy over the slip surface to investigate the earthquake rupture process and energy budget. For the 1994 M6.7 Northridge, 1992 M7.3 Landers, and 1995 M6.9 Kobe earthquakes, the distributions of static stress drop show strong heterogeneity, emphasizing the importance of asperities in the rupture process. Average values of static stress drop are 17, 11, and 4 Mpa for Northridge, Landers, and Kobe, respectively. These values are substantially higher than estimates based on simple crack models, suggesting that the failure process involves the rupture of asperities within the larger fault zone. Stress drop as a function of depth for the Northridge and Landers earthquakes suggests that stress drops are limited by crustal strength. For these two earthquakes, regions of high slip are surrounded by high values of stiffness. Particularly for the Northridge earthquake, the prominent patch of high slip in the central part of the fault is bordered by a ring of high stiffness and is consistent with expectations based on the failure of an asperity loaded at its edge due to exterior slip. Stiffness within an asperity is inversely related to its dimensions. Estimates of fracture energy, based on static stress drop, slip, and rupture speed, were used to investigate the nature of slip weakening at four locations near the hypocenter of the Kobe earthquake for comparison with independent results based on a dynamic model of this earthquake. One subfault updip and to the NE of the hypocenter has a fracture energy of 1.1 MJ/m2 and a slip-weakening distance, Dc, of 0.66 m. Right triangles, whose base and height are Dc and the dynamic stress drop, respectively, approximately overlie the slip-dependent stress given by Ide and Takeo (1997) for the same locations near the hypocenter. The total fracture energy for the Kobe earthquake, 3.7 ?? 1014 J, is about the same as the seismic energy (Ea = 3.2 ?? 1014 J.

Journal of Geophysical Research B: Solid Earth↗

USGS2005

Forecasting the evolution of seismicity in southern California: Animations built on earthquake stress transfer

We develop a forecast model to reproduce the distibution of main shocks, aftershocks and surrounding seismicity observed during 1986-200 in a 300 ?? 310 km area centered on the 1992 M = 7.3 Landers earthquake. To parse the catalog into frames with equal numbers of aftershocks, we animate seismicity in log time increments that lengthen after each main shock; this reveals aftershock zone migration, expansion, and densification. We implement a rate/state algorithm that incorporates the static stress transferred by each M ??? 6 shock and then evolves. Coulomb stress changes amplify the background seismicity, so small stress changes produce large changes in seismicity rate in areas of high background seismicity. Similarly, seismicity rate declines in the stress shadows are evident only in areas with previously high seismicity rates. Thus a key constituent of the model is the background seismicity rate, which we smooth from 1981 to 1986 seismicity. The mean correlation coefficient between observed and predicted M ??? 1.4 shocks (the minimum magnitude of completeness) is 0.52 for 1986-2003 and 0.63 for 1992-2003; a control standard aftershock model yields 0.54 and 0.52 for the same periods. Four M ??? 6.0 shocks struck during the test period; three are located at sites where the expected seismicity rate falls above the 92 percentile, and one is located above the 75 percentile. The model thus reproduces much, but certainly not all, of the observed spatial and temporal seismicity, from which we infer that the decaying effect of stress transferred by successive main shocks influences seismicity for decades. Finally, we offer a M ??? 5 earthquake forecast for 2005-2015, assigning probabilities to 324 10 ?? 10 km cells.

Journal of Geophysical Research B: Solid Earth↗

USGS2005

Vertical motions of the Puerto Rico Trench and Puerto Rico and their cause

The Puerto Rico trench exhibits great water depth, an extremely low gravity anomaly, and a tilted carbonate platform between (reconstructed) elevations of +1300 m and -4000 m. I argue that these features are manifestations of large vertical movements of a segment of the Puerto Rico trench, its forearc, and the island of Puerto Rico that took place 3.3 m.y. ago over a time period as short as 14-40 kyr. I explain these vertical movements by a sudden increase in the slab's descent angle that caused the trench to subside and the island to rise. The increased dip could have been caused by shearing or even by a complete tear of the descending North American slab, although the exact nature of this deformation is unknown. The rapid (14-40 kyr) and uniform tilt along a 250 km long section of the trench is compatible with scales of mantle flow and plate bending. The proposed shear zone or tear is inferred from seismic, morphological, and gravity observations to start at the trench at 64.5??W and trend southwestwardly toward eastern Puerto Rico. The tensile stresses necessary to deform or tear the slab could have been generated by increased curvature of the trench following a counterclockwise rotation of the upper plate and by the subduction of a large seamount.

Journal of Geophysical Research B: Solid Earth↗

USGS2005

Source mechanism of Vulcanian degassing at Popocatépetl Volcano, Mexico, determined from waveform inversions of very long period signals

The source mechanism of very long period (VLP) signals accompanying volcanic degassing bursts at Popocatépetl is analyzed in the 15–70 s band by minimizing the residual error between data and synthetics calculated for a point source embedded in a homogeneous medium. The waveforms of two eruptions (23 April and 23 May 2000) representative of mild Vulcanian activity are well reproduced by our inversion, which takes into account volcano topography. The source centroid is positioned 1500 m below the western perimeter of the summit crater, and the modeled source is composed of a shallow dipping crack (sill with easterly dip of 10°) intersecting a steeply dipping crack (northeast striking dike dipping 83° northwest), whose surface extension bisects the vent. Both cracks undergo a similar sequence of inflation, deflation, and reinflation, reflecting a cycle of pressurization, depressurization, and repressurization within a time interval of 3–5 min. The largest moment release occurs in the sill, showing a maximum volume change of 500–1000 m 3 , pressure drop of 3–5 MPa, and amplitude of recovered pressure equal to 1.2 times the amplitude of the pressure drop. In contrast, the maximum volume change in the dike is less (200–300 m 3 ), with a corresponding pressure drop of 1–2 MPa and pressure recovery equal to the pressure drop. Accompanying these volumetric sources are single-force components with magnitudes of 10 8 N, consistent with melt advection in response to pressure transients. The source time histories of the volumetric components of the source indicate that significant mass movement starts within the sill and triggers a mass movement response in the dike within a few seconds. Such source behavior is consistent with the opening of a pathway for escape of pent-up gases from slow pressurization of the sill driven by magma crystallization. The opening of this pathway and associated rapid evacuation of volcanic gases induces the pressure drop. Pressure recovery in the magma filling the sill is driven by diffusion of gases from the resulting supersaturated melt into bubbles. Assuming a penny-shaped crack at ambient pressure of 40 MPa, the observed pressure and volume variations can be modeled with the following attributes: crack radius (100 m), crack aperture (5 m), bubble number density (10 10 –10 12 m −3 ), initial bubble radius (10 −6 m), final bubble radius (&sim;10 −5 m), and net decrease of gas concentration in the melt (0.01 wt %).

Popocatépetl Volcano↗

USGS2004

Spatial and temporal multiyear sea ice distributions in the Arctic: A neural network analysis of SSM/I data, 1988-2001

Arctic multiyear sea ice concentration maps for January 1988-2001 were generated from SSM/I brightness temperatures (19H, 19V, and 37V) using modified multiple layer perceptron neural networks. Learning data for the neural networks were extracted from ice maps derived from Okean and ERS satellite imagery to capitalize on the stability of active radar multiyear ice signatures. Evaluations of three learning algorithms and several topologies indicated that networks constructed with error back propagation learning and 3-20-1 topology produced the most consistent and physically plausible results. Operational neural networks were developed specifically with January learning data, and then used to estimate daily multiyear ice concentrations from daily-averaged SSM/I brightness temperatures during January. Monthly mean maps were produced for analysis by averaging the respective daily estimates. The 14-year series of January multiyear ice distributions revealed dense and persistent cover in the central Arctic surrounded by expansive regions of highly fluctuating interannual cover. Estimates of total multiyear ice area by the neural network were intermediate to those of other passive microwave algorithms, but annual fluctuations and trends were similar among all algorithms. When compared to Radarsat estimates of multiyear ice concentration in the Beaufort and Chukchi Seas (1997-1999), average discrepancies were small (0.9-2.5%) and spatial coherency was reasonable, indicating the neural network's Okean and ERS learning data facilitated passive microwave inversion that emulated backscatter signatures. During 1988-2001, total January multiyear ice area declined at a significant linear rate of -54.3 x 10 3 km 2 /yr -1 (-1.4%/yr -1 ). The most persistent and extensive decline in multiyear ice concentration (-3.3%/yr -1 ) occurred in the southern Beaufort and Chukchi Seas. In autumn 1996, a large multiyear ice recruitment of over 106 km 2 (mostly in the Siberian Arctic) fully replenished the previous 8-year decline in total area, but it was followed by an accelerated and compensatory decline during the subsequent 4 years. Seventy-five percent of the interannual variation in January multiyear sea ice area was explained by linear regression on two atmospheric parameters: the previous inter's (JFM) Arctic Oscillation index as a proxy to melt duration and the previous year's average sea level pressure gradient across the Fram Strait as a proxy to annual ice export. Consecutive year changes (1994-2001) in January multiyear ice volume were significantly correlated with duration of the intervening melt season (R 2 = 0.73, -80.0 km 3 /d -1 ), emphasizing a large thermodynamic influence on the Arctic's mass sea ice balance during summers with anomalous melt durations.

Journal of Geophysical Research↗

USGS2005

Wave spectral energy variability in the northeast Pacific

The dominant characteristics of wave energy variability in the eastern North Pacific are described from NOAA National Data Buoy Center (NDBC) buoy data collected from 1981 to 2003. Ten buoys at distributed locations were selected for comparison based on record duration and data continuity. Long‐period (LP) [ T > 12] s, intermediate‐period [6 ≤ T ≤ 12] s, and short‐period [ T < 6] s wave spectral energy components are considered separately. Empirical orthogonal function (EOF) analyses of monthly wave energy anomalies reveal that all three wave energy components exhibit similar patterns of spatial variability. The dominant mode represents coherent heightened (or diminished) wave energy along the West Coast from Alaska to southern California, as indicated by composites of the 700 hPa height field. The second EOF mode reveals a distinct El Niño‐Southern Oscillation (ENSO)‐associated spatial distribution of wave energy, which occurs when the North Pacific storm track is extended unusually far south or has receded to the north. Monthly means and principal components (PCs) of wave energy levels indicate that the 1997–1998 El Niño winter had the highest basin‐wide wave energy within this record, substantially higher than the 1982–1983 El Niño. An increasing trend in the dominant PC of LP wave energy suggests that storminess has increased in the northeast Pacific since 1980. This trend is emphasized at central eastern North Pacific locations. Patterns of storminess variability are consistent with increasing activity in the central North Pacific as well as the tendency for more extreme waves in the south during El Niño episodes and in the north during La Niña.

Journal of Geophysical Research C: Oceans↗

USGS2005

Pitted cones and domes on Mars: Observations in Acidalia Planitia and Cydonia Mensae using MOC, THEMIS, and TES data

Domes and cones with summit pits located in Acidalia Planitia and Cydonia Mensae were studied using MOC and THEMIS images and a TES‐derived thermal inertia map. North of 40.5°N latitude, the features have a dome‐like morphology, and south of that latitude, the morphology is more cone‐like. Layering is apparent in the summit craters of fresher looking southern cones, and asymmetric aprons were observed in some instances. Some of the northern domes also display layering in their summit craters, but asymmetric aprons were not observed. The northern domes can also display multiple summit pits or no summit pits at all and can occur in association with higher‐albedo “pancake” features. The northern domes are higher in albedo but have apparent thermal inertias that are lower than the surrounding plains. The apparent thermal inertia values of the southern cones range from values comparable to the surrounding plains to slightly lower. From the TES thermal inertia map, we infer that the thermal inertia values of the pitted cones are between those of basaltic fine dust and sand, while those of the surrounding plains are closer to that of basaltic sand. While a unique interpretation of the origin of the pitted cones is not possible with the available data, we do not find compelling evidence to suggest an origin related to either basaltic volcanism or ground‐ice. Instead, an origin for these features through some combination of mud volcanism and evaporite deposition around geysers and/or springs is most consistent with the observations.

Journal of Geophysical Research E: Planets↗

USGS2005

Three decades of Martian surface changes

The surface of Mars has changed dramatically during the three decades spanned by spacecraft exploration. Comparisons of Mars Global Surveyor images with Viking and Mariner 9 pictures suggest that more than one third of Mars' surface area has brightened or darkened by at least 10%. Such albedo changes could produce significant effects on solar heating and the global circulation of winds across the planet. All of the major changes took place in areas of moderate to high thermal inertia and rock abundance, consistent with burial of rocky surfaces by thin dust layers deposited during dust storms and subsequent exposure of the rocky surfaces by aeolian erosion. Several distinct mechanisms contribute to aeolian erosion on Mars. Prevailing winds dominate erosion at low latitudes, producing diffuse albedo boundaries and elongated wind streaks generally oriented in the direction of southern summer winds. Dust devils darken the mid to high latitudes from 45 to 70 degrees during the summer seasons, forming irregular albedo patterns consisting of dark linear tracks. Dust storms produce regional albedo variations with distinct but irregular margins. Dark sand duties in southern high latitudes appear to be associated with regional darkening that displays diffuse albedo boundaries. No surface changes were observed to repeat regularly on an annual basis, but many of the changes took place in areas that alternate episodically between high- and low-albedo states as thin mantles of dust are deposited and later stripped off. Hence the face of Mars remains recognizable after a century of telescopic observations, in spite of the enormous extent of alteration that has taken place during the era of spacecraft exploration.

Journal of Geophysical Research E: Planets↗

USGS2005

Nitrogen and carbon flow from rock to water: Regulation through soil biogeochemical processes, Mokelumne River watershed, California, and Grand Valley, Colorado

Soil denitrification is an ecologically important nitrogen removal mechanism that releases to the atmosphere the greenhouse gas N2O, an intermediate product from the reduction of NO3- to N 2. In this study we evaluate the relationship between soil carbon and denitrification potential in watersheds with bedrock acting as a nonpoint source of nitrogen, testing the hypothesis that nitrate leaching to stream water is in part regulated by denitrification. Two sites, one in a Mediterranean climate and the other in an arid climate, were investigated to understand the interplay between carbon and denitrification potential. Both sites included carbonaceous bedrock with relatively high nitrogen concentrations (> 1,000 mg N kg-1) and had low background nitrogen concentrations in surface and groundwater. There was a net accumulation of carbon and nitrogen in soil relative to the corresponding bedrock, with the exception of carbonaceous shale from the arid site. There the concentration of carbon in the soil (15,620 mg C kg-1) was less than the shale parent (22,460 mg C kg-1), consistent with the bedrock being a source of soil carbon. Rates of denitrification potential (0.5-83 ??g N kg-1 hr-1) derived from laboratory incubations appeared to be related to the ratio of dissolved organic carbon and nitrate extracted from soils. These data indicate that microbial processes such as denitrification can help maintain background nitrogen concentrations to tens of ??M N in relatively undisturbed ecosystems when nitrogen inputs from weathering bedrock are accompanied by sufficient organic carbon concentrations to promote microbial nitrogen transformations.

Journal of Geophysical Research F: Earth Surface↗