Search USGSSearch

SEARCH · Search USGS

Results for “Developments in Earth Surface Processes”

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 127 records · Page 7Linked to original sources

Towards sustainable groundwater use: Setting long-term goals, backcasting, and managing adaptively

The sustainability of crucial earth resources, such as groundwater, is a critical issue. We consider groundwater sustainability a value‐driven process of intra‐ and intergenerational equity that balances the environment, society, and economy. Synthesizing hydrogeological science and current sustainability concepts, we emphasize three sustainability approaches: setting multigenerational sustainability goals, backcasting, and managing adaptively. As most aquifer problems are long‐term problems, we propose that multigenerational goals (50 to 100 years) for water quantity and quality that acknowledge the connections between groundwater, surface water, and ecosystems be set for many aquifers. The goals should be set by a watershed‐ or aquifer‐based community in an inclusive and participatory manner. Policies for shorter time horizons should be developed by backcasting, and measures implemented through adaptive management to achieve the long‐term goals. Two case histories illustrate the importance and complexity of a multigenerational perspective and adaptive management. These approaches could transform aquifer depletion and contamination to more sustainable groundwater use, providing groundwater for current and future generations while protecting ecological integrity and resilience.

Ground Water

Channel response to sediment wave propagation and movement, Redwood Creek, California, USA

Redwood Creek, north coastal California, USA, has experienced dramatic changes in channel configuration since the 1950s. A series of large floods (in 1955, 1964, 1972 and 1975) combined with the advent of widespread commercial timber harvest and road building resulted in extensive erosion in the basin and contributed high sediment loads to Redwood Creek. Since 1975, no peak flows have exceeded a 5 year recurrence interval. Twenty years of cross-sectional survey data document the downstream movement of a ‘sediment wave’ in the lower 26 km of this gravel-bedded river at a rate of 800 to 1600 m a −1 during this period of moderately low flows. Higher transit rates are associated with reaches of higher unit stream power. The wave was initially deposited at a site with an abrupt decrease in channel gradient and increase in channel width. The amplitude of the wave has attenuated more than 1 m as it moved downstream, and the duration of the wave increased from eight years upstream to more than 20 years downstream. Channel aggradation and subsequent degradation have been accommodated across the entire channel bed. Channel width has not decreased significantly after initial channel widening from large (>25 year recurrence interval) floods. Three sets of longitudinal surveys of the streambed showed the highest increase in pool depths and frequency in a degrading reach, but even the aggrading reach exhibited some pool development through time. The aggraded channel bed switched from functioning as a sediment sink to a significant sediment source as the channel adjusted to high sediment loads. From 1980 to 1990, sediment eroded from temporary channel storage represented about 25 per cent of the total sediment load and 95 per cent of the bedload exported from the basin.

California

Controls on valley spacing in landscapes subject to rapid base-level fall

What controls the architecture of drainage networks is a fundamental question in geomorphology. Recent work has elucidated the mechanisms of drainage network development in steadily uplifting landscapes, but the controls on drainage-network morphology in transient landscapes are relatively unknown. In this paper we exploit natural experiments in drainage network development in incised Plio-Quaternary alluvial fan surfaces in order to understand and quantify drainage network development in highly transient landscapes, i.e. initially unincised low-relief surfaces that experience a pulse of rapid base-level drop followed by relative base-level stasis. Parallel drainage networks formed on incised alluvial-fan surfaces tend to have a drainage spacing that is approximately proportional to the magnitude of the base-level drop. Numerical experiments suggest that this observed relationship between the magnitude of base-level drop and mean drainage spacing is the result of feedbacks among the depth of valley incision, mass wasting and nonlinear increases in the rate of colluvial sediment transport with slope gradient on steep valley side slopes that lead to increasingly wide valleys in cases of larger base-level drop. We identify a threshold magnitude of base-level drop above which side slopes lengthen sufficiently to promote increases in contributing area and fluvial incision rates that lead to branching and encourage drainage networks to transition from systems of first-order valleys to systems of higher-order, branching valleys. The headward growth of these branching tributaries prevents the development of adjacent, ephemeral drainages and promotes a higher mean valley spacing relative to cases in which tributaries do not form. Model results offer additional insights into the response of initially unincised landscapes to rapid base-level drop and provide a preliminary basis for understanding how varying amounts of base-level change influence valley network morphology.

Earth Surface Processes and Landforms

A Next Generation Lunar Orbiter mission

The Moon is the scientific foundation for our knowledge of the early evolution and impact history of the terrestrial planets. Over the last decades the lunar science community has made significant progress in addressing key lunar science and exploration goals, while defining many new high-priority scientific questions regarding the formation and evolution of the Moon. On a broad scale, the last Planetary Decadal Survey defined three scientific objectives to guide studies of the inner planets and Moon. These objectives address the origins and diversity of terrestrial planets, the evolution of life on terrestrial planets, and climate processes on Earth-like planets [1]. Many of these objectives carry forward into exploration goals of a renewed lunar human exploration program, and urgently addressing these objectives will enable the rapid development of exploration plans. A large (Flagship or New Frontiers class) Next Generation Lunar Orbiter (NGLO) mission would address the last Planetary Decadal Survey objective of understanding the origin and diversity of terrestrial planets by studying the geochemistry and geology of the Moon at an unparalleled resolution compared to other lunar mission datasets. Also, NGLO would address the objective of studying the evolution of life on terrestrial planets by furthering knowledge about the composition and distribution of volatile elements on the lunar surface and better characterizing the past and present-day impact rates in the inner Solar System in order to better understand the original delivery of water to Earth. Key exploration goals, including identifying the nature and distribution of lunar volatiles (i.e., water, ice), mapping and characterizing potentially valuable lunar resources, and establishing a human presence on the Moon, also would be addressed by NGLO.

Bulletin of the AAS

Field-trip guide to Mount Hood, Oregon, highlighting eruptive history and hazards

This guidebook describes stops of interest for a geological field trip around Mount Hood volcano. It was developed for the 2017 International Association of Volcanology and Chemistry of the Earth’s Interior (IAVCEI) Scientific Assembly in Portland, Oregon. The intent of this guidebook and accompanying contributions is to provide an overview of Mount Hood, including its chief geologic processes, magmatic system, eruptive history, local tectonics, and hazards, by visiting a variety of readily accessible localities. We also describe coeval, largely monogenetic, volcanoes in the region. Accompanying the field-trip guidebook are separately authored contributions that discuss in detail the Mount Hood magmatic system and its products and behavior (Kent and Koleszar, this volume); Mount Hood earthquakes and their relation to regional tectonics and the volcanic system (Thelen and Moran, this volume); and young surface faults cutting the broader Mount Hood area whose extent has come to light after acquisition of regional light detection and ranging coverage (Madin and others, this volume). The trip makes an approximately 175-mile (280-kilometer) clockwise loop around Mount Hood, starting and ending in Portland. The route heads east on Interstate 84 through the Columbia River Gorge National Scenic Area. The guidebook points out only a few conspicuous features of note in the gorge, but many other guides to the gorge are available. The route continues south on the Mount Hood National Scenic Byway on Oregon Route 35 following Hood River, and returns to Portland on U.S. Highway 26 following Sandy River. The route traverses rocks as old as the early Miocene Eagle Creek Formation and overlying Columbia River Basalt Group of middle Miocene age, but chiefly lava flows and clastic products of arc volcanism of late Miocene to Holocene age.

Oregon

A geomorphic-process-based cellular automata model of colluvial wedge morphology and stratigraphy

The development of colluvial wedges at the base of fault scarps following normal-faulting earthquakes serves as a sedimentary record of paleoearthquakes and is thus crucial in assessing seismic hazard. Although there is a large body of observations of colluvial wedge development, connecting this knowledge to the physics of sediment transport can open new frontiers in our understanding. To explore theoretical colluvial wedge evolution, we develop a cellular automata model driven by the production and disturbance (e.g., bioturbative reworking) of mobile regolith and fault-scarp collapse. We consider both 90 and 60 ∘ dipping faults and allow the colluvial wedges to develop over 2000 model years. By tracking sediment transport time, velocity, and provenance, we classify cells into analogs for the debris and wash sedimentary facies commonly described in paleoseismic studies. High values of mobile regolith production and disturbance rates produce relatively larger and more wash-facies-dominated wedges, whereas lower values produced relatively smaller, debris-facies-dominated wedges. Higher lateral collapse rates lead to more debris facies relative to wash facies. Many of the modeled colluvial wedges fully developed within 2000 model years after the earthquake, with many being much faster when process rates are high. Finally, for scenarios with the same amount of vertical displacement, differently sized colluvial wedges developed depending on the rates of geomorphic processes and fault dip. A change in these variables, say by environmental change such as precipitation rates, could theoretically result in different colluvial wedge facies assemblages for the same characteristic earthquake rupture scenario. Finally, the stochastic nature of collapse events, when coupled with high disturbance, illustrates that multiple phases of colluvial deposition are theoretically possible for a single earthquake event.

Earth Surface Dynamics

Volcanology: Lessons learned from Synthetic Aperture Radar imagery

Twenty years of continuous Earth observation by satellite SAR have resulted in numerous new insights into active volcanism, including a better understanding of subsurface magma storage and transport, deposition of volcanic materials on the surface, and the structure and development of volcanic edifices. This massive archive of data has resulted in fundamental leaps in our understanding of how volcanoes work – for example, identifying magma accumulation at supposedly quiescent volcanoes, even in remote areas or in the absence of ground-based data. In addition, global compilations of volcanic activity facilitate comparison of deformation behavior between different volcanic arcs and statistical evaluation of the strong link between deformation and eruption. SAR data are also increasingly used in timely hazard evaluation thanks to decreases in data latency and growth in processing and analysis techniques. The existing archive of SAR imagery is on the cusp of being enhanced by a new generation of satellite SAR missions, in addition to ground-based and airborne SAR systems, which will provide enhanced temporal and spatial resolution, broader geographic coverage, and improved availability of data to the scientific community. Now is therefore an opportune time to review the contributions of SAR imagery to volcano science, monitoring, and hazard mitigation, and to explore the future potential for SAR in volcanology. Provided that the ever-growing volume of SAR data can be managed effectively, we expect the future application of SAR data to expand from being a research tool for analyzing volcanic activity after the fact, to being a monitoring and research tool capable of imaging a wide variety of processes on different temporal and spatial scales as those processes are occurring. These data can then be used to develop new models of how volcanoes work and to improve quantitative forecasts of volcanic activity as a means of mitigating risk from future eruptions.

Journal of Volcanology and Geothermal Research

Wind-driven particle mobility on Mars: Insights from Mars Exploration Rover observations at "El Dorado" and surroundings at Gusev Crater

The ripple field known as 'El Dorado' was a unique stop on Spirit's traverse where dust-raising, active mafic sand ripples and larger inactive coarse-grained ripples interact, illuminating several long-standing issues of Martian dust mobility, sand mobility, and the origin of transverse aeolian ridges. Strong regional wind events endured by Spirit caused perceptible migration of ripple crests in deposits SSE of El Dorado, erasure of tracks in sandy areas, and changes to dust mantling the site. Localized thermal vortices swept across El Dorado, leaving paths of reduced dust but without perceptibly damaging nearly cohesionless sandy ripple crests. From orbit, winds responsible for frequently raising clay-sized dust into the atmosphere do not seem to significantly affect dunes composed of (more easily entrained) sand-sized particles, a long-standing paradox. This disparity between dust mobilization and sand mobilization on Mars is due largely to two factors: (1) dust occurs on the surface as fragile, low-density, sand-sized aggregates that are easily entrained and disrupted, compared with clay-sized air fall particles; and (2) induration of regolith is pervasive. Light-toned bed forms investigated at Gusev are coarse-grained ripples, an interpretation we propose for many of the smallest linear, light-toned bed forms of uncertain origin seen in high-resolution orbital images across Mars. On Earth, wind can organize bimodal or poorly sorted loose sediment into coarse-grained ripples. Coarse-grained ripples could be relatively common on Mars because development of durable, well-sorted sediments analogous to terrestrial aeolian quartz sand deposits is restricted by the lack of free quartz and limited hydraulic sediment processing. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research E: Planets

On the state of lithospheric stress in the absence of applied tectonic forces

Numerous published analyses of the nontectonic state of stress are based on Hooke's law and the boundary condition of zero horizontal deformation. This approach has been used to determine the gravitational stress state as well as the effects of processes such as erosion and temperature changes on the state of lithospheric stress. The major disadvantage of these analyses involves the assumption of lateral constraint which seems unrealistic in view of the observational fact that the crust can deform horizontally in response to applied loads. If the same problems are addressed by assuming that the remote stress state is constant, instead of the condition of zero horizontal deformation, then the resulting stress states are entirely different and in good accord with observations. In the absence of applied tectonic forces the only likely gravitational stress states are those for which all three principal stresses are nearly equal. To the contrary, the gravitational stress states developed on the basis of the lateral constraint assumption can be ruled out. The processes of erosion and sedimentation have slight tendencies to increase and decrease, respectively, the state of deviatoric stress. In particular, for initial stress states in the range of slightly extensional to compressional, erosion has the effect of enhancing the ratio of average horizontal to vertical stress, which may explain, at least in part, the common observation of high near-surface horizontal stresses. Temperature changes have only minor effects on the stress state, as averaged over the thickness of the lithosphere.

Journal of Geophysical Research Solid Earth

Coincidence and spatial variability of geology, soils, and vegetation, Mill Run watershed, Virginia

The Mill Run watershed is a structurally‐controlled synclinal basin on the eastern limb of the Massanutten Mountain complex of northwestern Virginia. Bedrock contacts are obscured by coarse sandstone debris from exposures near basin divides. Colluvium blankets more than half the basin, masking geomorphic surfaces, affecting vegetation patterns, and contributing to the convexity of the alluvial, terrace, pediment and erosion surfaces. Examination of the bedrock geology, geomorphology, soils and vegetation shows distinct distributional correspondences. Vegetation is strongly interdependent with geomorphology, bedrock geology, and soils. On convex colluvial slopes, mixed hardwood forests are most common. In concave coves and deep gorges, mixed hardwoods are replaced by conifers. In thin colluvium, in poorly developed soils, and on blockfields, chestnut oak is singularly prevalent. Conifers dominate shaley bedrock areas. Soils and surficial sediments have a major effect on near‐surface hydrology. During wet seasons, cemented horizons in the subsurface cause temporary saturation in the superjacent horizons; lateral movement of soil‐water effectively eliminates a vertical component of ground‐water recharge. Vegetation is strongly dependent on water availability and thus reflects the distribution of subsurface barriers and sediment‐soil fades changes.

Virginia

Interferometric synthetic aperture radar: Building tomorrow's tools today

A synthetic aperture radar (SAR) system transmits electromagnetic (EM) waves at a wavelength that can range from a few millimeters to tens of centimeters. The radar wave propagates through the atmosphere and interacts with the Earth’s surface. Part of the energy is reflected back to the SAR system and recorded. Using a sophisticated image processing technique, called SAR processing (Curlander and McDonough, 1991), both the intensity and phase of the reflected (or backscattered) signal of each ground resolution element (a few meters to tens of meters) can be calculated in the form of a complex-valued SAR image representing the reflectivity of the ground surface. The amplitude or intensity of the SAR image is determined primarily by terrain slope, surface roughness, and dielectric constants, whereas the phase of the SAR image is determined primarily by the distance between the satellite antenna and the ground targets, slowing of the signal by the atmosphere, and the interaction of EM waves with ground surface. Interferometric SAR (InSAR) imaging, a recently developed remote sensing technique, utilizes the interaction of EM waves, referred to as interference, to measure precise distances. Very simply, InSAR involves the use of two or more SAR images of the same area to extract landscape topography and its deformation patterns.

Alaska Satellite Facility News and Notes

Estimating sand concentrations using ADCP‐based acoustic inversion in a large fluvial system characterized by bi‐modal suspended‐sediment distributions

Quantifying sediment flux within rivers is a challenge for many disciplines due, mainly, to difficulties inherent to traditional sediment sampling methods. These methods are operationally complex, high cost, and high risk. Additionally, the resulting data provide a low spatial and temporal resolution estimate of the total sediment flux, which has impeded advances in the understanding of the hydro‐geomorphic characteristics of rivers. Acoustic technologies have been recognized as a leading tool for increasing the resolution of sediment data by relating their echo intensity level measurements to suspended sediment. Further effort is required to robustly test and develop these techniques across a wide range of conditions found in natural river systems. This article aims to evaluate the application of acoustic inversion techniques using commercially available, down‐looking acoustic Doppler current profilers (ADCPs) in quantifying suspended sediment in a large sand bed river with varying bi‐modal particle size distributions, wash load and suspended‐sand ratios, and water stages. To achieve this objective, suspended sediment was physically sampled along the Paraná River, Argentina, under various hydro‐sedimentological regimes. Two ADCPs emitting different sound frequencies were used to simultaneously profile echo intensity level within the water column. Using the sonar equation, calibrations were determined between suspended‐sand concentrations and acoustic backscatter to solve the inverse problem. The study also analyzed the roles played by each term of the sonar equation, such as ADCP frequency, power supply, instrument constants, and particle size distributions typically found in sand bed rivers, on sediment attenuation and backscatter. Calibrations were successfully developed between corrected backscatter and suspended‐sand concentrations for all sites and ADCP frequencies, resulting in mean suspended‐sand concentration estimates within about 40% of the mean sampled concentrations. Noise values, calculated using the sonar equation and sediment sample characteristics, were fairly constant across evaluations, suggesting that they could be applied to other sand bed rivers.

Earth Surface Processes and Landforms

Radarclinometry of the Earth and Venus from space-shuttle and Venera-15 imagery

The project to develop a line-integral approach to 2-dimensional radarclinometry and to bring it to the status of producing topographic maps from real radar images has been concluded. The final developments of the theory itself have involved a trial-and-error resolution of the curvature decision process at each integration step over range as follows: (1) Locally Indeterminate Azimuth-Azimuth Curvature is invoked if the range-directed path of integration is within 1 ??? in angle of the tangent to a local characteristic curve of the partial differential equation of radarclinometry (equivalent to a lapse in the necessity for an auxiliary curvature assumption); (2) Local Cylindricity is invoked if the local image isophote has a radius-of-curvature greater than 50 pixels; (3) Least-Squared Local Sphericity is invoked if the characteristic curve trends at greater than 70 ??? to the range direction (the auxiliary curvature assumption is becoming a sufficiently strong influence as to warrant the overconstraint), and (4) the default hypothesis, which is invoked most often, is the localization through the Euler/Lagrange equation from the calculus of variations of the global principle of minimization of the surface area of the terrain. The development of the set of line integrals into a 2-dimensional topographic surface is not practically achieved by branching the line integral at the range threshold, because the radarclinometry equations are too frequently coupled but weakly to the slope component in the direction of radar-azimuth, and under circumstances for which the powerfully influential auxiliary curvature assumption is too unrealistic. In other words, a line integration in radar-azimuth is far more frequently directed orthogonally to the local characteristic curve than is one carried out over range. Such orthogonality results in stepping the strike under the exclusive control of the curvature assumption. Instead, a quasi-surface-integration step is taken by modeling the dependence on initial strike of the gravitational potential energy of the vertical slab of terrain under the range-profile. The adopted starting strike for the range integral is the one which minimizes the gravitational potential energy. This radarclinometric method, in combination with my recently published method for determining an effective radar back-scattering function from one-dimensional slope statistics and image pixel-signal statistics, was applied to three images. First, to separate theoretical difficulties from experimental impediments, an artificial radar image was generated from a topographic map of the Lake Champlain West quadrangle in the Adirondack Mountains. Except for the regional trend in elevation, to which radarclinometry is insensitive by design, the agreement between the original and derived topography appears good. The morphologies agree and the range of relief is the same to within 4%. As an example of data of the highest quality available from space-borne radar at the present time, a SIR-B image of very rugged terrain in the coastal mountains of Oregon was similarly processed. The result, after filtering to redistribute photoclinometric errors about the two-dimensional spatial spectrum, agrees with ground truth almost as well. As an example of the worst possible data, in terms of signal-to-noise ratio and radar incidence angle (no detraction from the praise due the first high resolution space-borne radar-imaging of Venus intended), a Venera-15 image segment in Sedna Planitia just north-east of Sapho was processed, using Venera altimetry and Pioneer roughness data for slope statistics, in spite of the resolution mis-match. Considerably more trial-and-error filtering was required. The result appears plausible, but an error check is, of course, impossible. ?? 1990 Kluwer Academic Publishers.

Earth, Moon and Planets

Understanding natural systems; a perspective for land-use planning in Appalachian Kentucky

An eight-county area at the headwaters of the Kentucky River has been designated the Kentucky River Area Development District (KRADD) by the Appalachian Regional Commission. The objective of the project described in this report has been to provide materials to KRADD planners in a format and containing terminology usable by local people untrained in earth science. Experimental maps (not included in this report) have been prepared largely from preexisting data. Time and cost limitations required a regional analysis as well as somewhat more detailed examples of selected localities. Most of the maps produced to meet these needs show the abundance and distribution of naturally occurring materials and the areas affected by various geomorphic processes. Three types of maps, showing current land use, slope, and flood-prone areas, present both basic and derived data directly applicable to specific land-use decisions. Basic map information on quality and quantity of surface and ground water, bedrock and surficial geology, and mineral fuels can be interpreted for a wide variety of current and potential uses. Texts accompanying the maps explain bedrock control of geomorphic processes, distribution and significance of surficial deposits, and hydrologic characteristics of the intricately dissected eastern Kentucky terrain. Within this conceptual framework, geomorphic processes and the landscape may be evaluated in humanly significant terms of low to high potential risk, thereby indicating both opportunities and limitations for land use.

Kentucky

Gap analysis: Concepts, methods, and recent results

Rapid progress is being made in the conceptual, technical, and organizational requirements for generating synoptic multi-scale views of the earth's surface and its biological content. Using the spatially comprehensive data that are now available, researchers, land managers, and land-use planners can, for the first time, quantitatively place landscape units - from general categories such as 'Forests' or 'Cold-Deciduous Shrubland Formation' to more categories such as 'Picea glauca-Abies balsamea-Populus spp. Forest Alliance' - in their large-area contexts. The National Gap Analysis Program (GAP) has developed the technical and organizational capabilities necessary for the regular production and analysis of such information. This paper provides a brief overview of concepts and methods as well as some recent results from the GAP projects. Clearly, new frameworks for biogeographic information and organizational cooperation are needed if we are to have any hope of documenting the full range of species occurrences and ecological processes in ways meaningful to their management. The GAP experience provides one model for achieving these new frameworks.

Landscape Ecology

An overview of the Landsat Data Continuity Mission

The Landsat Data Continuity Mission (LDCM) is the follow-on mission to Landsat 7 and will be the eighth mission in the Landsat series. The mission is in development via an interagency partnership between the National Aeronautics and Space Administration (NASA) and the Department of Interior (DOI) / United States Geological Survey (USGS). The LDCM satellite will carry two earth-observing sensors, the Operational Land Imager (OLI) to collect image data for nine spectral bands in the reflective portion of the spectrum and the Thermal Infrared Sensor (TIRS) to collect coincident image data for two thermal spectral bands. The LDCM ground segment will control the satellite and will receive, process, archive, and distribute the science data collected by the OLI and TIRS instruments. The USGS Earth Resources Observation & Science Center (EROS) will distribute LDCM data products at no cost to requestors. The mission objective is to continues the Landsat program's collection, archive, and distribution of multispectral imagery affording global, synoptic, and repetitive coverage of the Earth's land surfaces at a scale where natural and human-induced changes can be detected, differentiated, characterized, and monitored over time. The LDCM launch readiness date is currently December, 2012.

Conference Paper

SPEAR: The next generation GFDL modeling system for seasonal to multidecadal prediction and projection

We document the development and simulation characteristics of the next generation modeling system for seasonal to decadal prediction and projection at the Geophysical Fluid Dynamics Laboratory (GFDL). SPEAR ( S eamless System for P rediction and EA rth System R esearch) is built from component models recently developed at GFDL—the AM4 atmosphere model, MOM6 ocean code, LM4 land model, and SIS2 sea ice model. The SPEAR models are specifically designed with attributes needed for a prediction model for seasonal to decadal time scales, including the ability to run large ensembles of simulations with available computational resources. For computational speed SPEAR uses a coarse ocean resolution of approximately 1.0° (with tropical refinement). SPEAR can use differing atmospheric horizontal resolutions ranging from 1° to 0.25°. The higher atmospheric resolution facilitates improved simulation of regional climate and extremes. SPEAR is built from the same components as the GFDL CM4 and ESM4 models but with design choices geared toward seasonal to multidecadal physical climate prediction and projection. We document simulation characteristics for the time mean climate, aspects of internal variability, and the response to both idealized and realistic radiative forcing change. We describe in greater detail one focus of the model development process that was motivated by the importance of the Southern Ocean to the global climate system. We present sensitivity tests that document the influence of the Antarctic surface heat budget on Southern Ocean ventilation and deep global ocean circulation. These findings were also useful in the development processes for the GFDL CM4 and ESM4 models.

Journal of Advances in Modeling Earth Systems

Petrographic and geochemical evidence for the formation of primary, bacterially induced lacustrine dolomite: La Roda 'white earth' (Pliocene, Central Spain)

Upper Pliocene dolomites ('white earth') from La Roda, Spain, offer a good opportunity to evaluate the process of dolomite formation in lakes. The relatively young nature of the deposits could allow a link between dolomites precipitated in modern lake systems and those present in older lacustrine formations. The La Roda Mg-carbonates (dolomite unit) occur as a 3??5- to 4-m- thick package of poorly indurated, white, massive dolomite beds with interbedded thin deposits of porous carbonate displaying root and desiccation traces as well as local lenticular gypsum moulds. The massive dolomite beds consist mainly of loosely packed 1- to 2-??m-sized aggregates of dolomite crystals exhibiting poorly developed faces, which usually results in a subrounded morphology of the crystals. Minute rhombs of dolomite are sparse within the aggregates. Both knobbly textures and clumps of spherical bodies covering the crystal surfaces indicate that bacteria were involved in the formation of the dolomites. In addition, aggregates of euhedral dolomite crystals are usually present in some more clayey (sepiolite) interbeds. The thin porous carbonate (mostly dolomite) beds exhibit both euhedral and subrounded, bacterially induced dolomite crystals. The carbonate is mainly Ca-dolomite (51-54 mol% CaCO3), showing a low degree of ordering (degree of ordering ranges from 0??27 to 0??48). Calcite is present as a subordinate mineral in some samples. Sr, Mn and Fe contents show very low correlation coefficients with Mg/Ca ratios, whereas SiO2 and K contents are highly correlated. ??18O- and ??13C-values in dolomites range from -3??07??? to 5??40??? PDB (mean = 0??06, ?? = 1??75) and from -6??34??? to -0??39??? PDB (mean = -3??55, ?? = 1??33) respectively. Samples containing significant amounts of both dolomite and calcite do not in general show significant enrichment or depletion in 18O and 13C between the two minerals. The correlation coefficient between ??18O and ??13C for dolomite is extremely low and negative (r = -0??05), whereas it is higher and positive (r = 0??47) for calcite. The lacustrine dolomite deposit from La Roda is interpreted mainly as a result of primary precipitation of dolomite in a shallow, hydrologically closed perennial lake. The lake was supplied by highly saturated HCO3-/CO32- groundwater that leached dolomitic Mesozoic formations. Precipitation of dolomite from alkaline lake waters took place under a semi-arid to arid climate. However, according to our isotopic data, strong evaporative conditions were not required for the formation of the La Roda dolomite. A significant contribution by bacteria to the formation of the dolomites is assumed in view of both petrographic and geochemical evidence.

Sedimentology