Palynological evidence for early Holocene aridity in the southern Sierra Nevada, California

Sediments of Balsam Meadow have produced a 11,000-yr pollen record from the southern Sierra Nevada of California. The Balsam Meadow diagram is divided into three zones. (1) The Artemisia zone (11,000–7000 yr B.P.) is characterized by percentages of sagebrush (Artemisia) and other nonarboreal pollen higher than can be found in the modern local vegetation. Vegetation during this interval was probably similar to the modern vegetation on the east slope of the Sierra Nevada and the climate was drier than that of today. (2) Pinus pollen exceeded 80% from 7000 to 3000 yr B.P. in the Pinus zone. The climate was moister than during the Artemisia zone. (3) Fir (Abies, Cupressaceae, and oak (Quercus) percentages increased after 3000 yr B.P. in the Abies zone as the modern vegetation at the site developed and the present cool-moist climatic regime was established. Decreased fire frequency after 1200 yr B.P. is reflected in decreased abundance of macroscopic charcoal and increased concentration of Abies magnifica and Pinus murrayana needles.     

Turbulence mechanisms at an agricultural site

An extensive set of turbulence data from the 3- and 12-m heights taken over an agricultural site (Marsta, Sweden) are analyzed and compared with data from ideal sites.In unstable air, Monin-Obukhov similarity is found to be valid for the non-dimensional gradients of wind, _m , temperature, _h , and humidity, _e , for (only a few data), for _T /|T _*|,/ _E /|E _*| and for the non-dimensionalized inertial subrange spectra of temperature and humidity. Where comparison is possible, the unstable data also agree with those found in the Kansas study, with one remarkable exception, the inertial subrange constant of the temperature spectrum, ₁, being only 0.39, compared to the value 0.80 found at the Kansas site.On the stable side, most similarity predictions break down, with most of the data differing systematically from the corresponding Kansas results, the only exception being . The inertial subrange constants for temperature, ₁, and for humidity, ₁ are found to have the same values, 0.39 and 0.30, respectively, as they do on the unstable side. Remarkable similarity is found for the shape of the stable u- and - and e-spectra. In addition, this shape is found to be identical with that found in Kansas. The peak wavelength of the stable u-, and -spectra is found to be about four times larger than it is for the corresponding Kansas spectra. This is interpreted to be a result of the increased macro-roughness at the Marsta site as compared with that at the Kansas site. A possible explanation for the low ₁-value is discussed, suggesting that ₁ is not a universal constant, but instead dependent on the turbulent structure.     

Fabrication of synthetic calcite–muscovite rocks with variable texture — An analogue to cataclasite fabrics?

A series of large diameter calcite–muscovite aggregates has been prepared from calcite and muscovite powders, in order to gain a better understanding of how texture develops in impure carbonate rocks. The development of the microstructure and the crystallographic preferred orientation (CPO, texture) during the preparation process is described. The synthetic rocks have been fabricated from powders of calcite and muscovite by uniaxial cold-pressing at loads up to 400 MPa and subsequent hot isostatic pressing (HIPping) at pressures of 150 to 170 MPa and a temperature of 670 °C. The resulting textures and microstructures are homogeneous throughout the samples. The calcite CPO is generated by rigid body rotation and twinning during cold-pressing and is not significantly altered by recrystallization during HIPping. Grain growth during HIPping is observed in pure calcite samples, but is inhibited through high porosity and the presence of muscovite in the mixed aggregates. The preferred orientation of the calcite c-axes is found to increase with increasing uniaxial cold pressure, and to be independent of the muscovite content. The magnetic bulk susceptibility of the starting material has been changed by the formation of ferromagnetic impurities during fabrication. Comparison of the samples to natural calcite fabrics from fault zones show the potential of the experiments and fabric analyses presented to analyze and to better understand the deformation mechanisms of fault zones.     

Mortar coupling and upscaling of pore-scale models

Pore-scale models are becoming increasingly useful as predictive tools for modeling flow and transport in porous media. These models can accurately represent the 3D pore-structure of real media. Currently first-principles modeling methods are being employed for obtaining qualitative and quantitative behavior. Generally, artificial, simple boundary conditions are imposed on a model that is used as a stand-alone tool for extracting macroscopic parameters. However, realistic boundary conditions, reflecting flow and transport in surrounding media, may be necessary for behavior that occurs over larger length scales or including pore-scale models in a multiscale setting. Here, pore-scale network models are coupled to adjacent media (additional pore-scale or continuum-scale models) using mortars. Mortars are 2D finite-element spaces employed to couple independent subdomains by enforcing continuity of pressure and flux at shared boundary interfaces. While mortars have been used in the past to couple subdomains of different models, physics, and meshes, they are extended here for the first time to pore-scale models. The approach is demonstrated by modeling single-phase flow in coupled pore-scale models, but the methodology can be utilized to model dynamic processes and perform multiscale modeling in 3D continuum simulators for flow and transport.     

Topographic analyses of K*lauea Volcano, Hawai'i, from interferometric airborne radar

 We analyze digital topographic data collected in September 1993 over a ∼500-km² portion of K*lauea Volcano, Hawai'i, by the C-band (5.6-cm wavelength) topographic synthetic aperture radar (TOPSAR) airborne interferometric radar. Field surveys covering an ∼1-km² area of the summit caldera and the distal end of an ∼8-m-thick 'a'* flow indicate that the 10-m spatial resolution TOPSAR data have a vertical accuracy of 1–2 m over a variety of volcanic surfaces. After conversion to a common datum, TOPSAR data agree favorably with a digital elevation model (DEM) produced by the U.S. Geological Survey (USGS), with the important exception of the region of the ongoing eruption (which postdates the USGS DEM). This DEM comparison gives us confidence that subtracting the USGS data from TOPSAR data will produce a reasonable estimate of the erupted volume as of September 1993. This subtraction produces dense rock equivalent (DRE) volumes of 392, 439, and 90×10⁶ m³ for the Pu'u '*'*, K*pa'ianah*, and episode 50–53 stages of the eruption, respectively. These are 124, 89, and 94% of the volumes calculated by staff of the Hawaiian Volcano Observatory (HVO) but do not include lava of K*pa'ianah* and episodes 50–53 that flowed into the ocean and are thus invisible to TOPSAR. Accounting for this lava increases the TOPSAR volumes to 124, 159, and 129% of the HVO volumes. Including the ±2-m uncertainty derived from the field surveys produces TOPSAR-derived volumes for the eruption as a whole that range between 81 and 125% of the USGS-derived values. The vesicularity- and ocean-corrected TOPSAR volumes yield volumetric eruption rates of 4.5, 4.5, and 2.7 m³/s for the three stages of the eruption, which compare with HVO-derived values of 3.6, 2.8, and 2.1 m³/s, respectively. Our analysis shows that care must be taken when vertically registering the TOPSAR and USGS DEMs to a common datum because C-band TOPSAR penetrates only partially into thick forest and therefore produces a DEM within the tree canopy, whereas the USGS DEM is adjusted for vegetation. Received: 28 April 1998 / Accepted: 1 February 1999     

Simulation of Water‐Table Aquifers Using Specified Saturated Thickness

Simulating groundwater flow in a water‐table (unconfined) aquifer can be difficult because the saturated thickness available for flow depends on model‐calculated hydraulic heads. It is often possible to realize substantial time savings and still obtain accurate head and flow solutions by specifying an approximate saturated thickness a priori, thus linearizing this aspect of the model. This specified‐thickness approximation often relies on the use of the “confined” option in numerical models, which has led to confusion and criticism of the method. This article reviews the theoretical basis for the specified‐thickness approximation, derives an error analysis for relatively ideal problems, and illustrates the utility of the approximation with a complex test problem. In the transient version of our complex test problem, the specified‐thickness approximation produced maximum errors in computed drawdown of about 4% of initial aquifer saturated thickness even when maximum drawdowns were nearly 20% of initial saturated thickness. In the final steady‐state version, the approximation produced maximum errors in computed drawdown of about 20% of initial aquifer saturated thickness (mean errors of about 5%) when maximum drawdowns were about 35% of initial saturated thickness. In early phases of model development, such as during initial model calibration efforts, the specified‐thickness approximation can be a very effective tool to facilitate convergence. The reduced execution time and increased stability obtained through the approximation can be especially useful when many model runs are required, such as during inverse model calibration, sensitivity and uncertainty analyses, multimodel analysis, and development of optimal resource management scenarios.     

Soil carbon storage following road removal and timber harvesting in redwood forests

Soil carbon storage plays a key role in the global carbon cycle and is important for sustaining forest productivity. Removal of unpaved forest roads has the potential for increasing carbon storage in soils on forested terrain as treated sites revegetate and soil properties improve on the previously compacted road surfaces. We compared soil organic carbon (SOC) content at several depths on treated roads to SOC in adjacent second‐growth forests and old‐growth redwood forests in California, determined whether SOC in the upper 50 cm of soil varies with the type of road treatment, and assessed the relative importance of site‐scale and landscape‐scale variables in predicting SOC accumulation in treated road prisms and second‐growth redwood forests. Soils were sampled at 5, 20, and 50 cm depths on roads treated by two methods (decommissioning and full recontouring), and in adjacent second‐growth and old‐growth forests in north coastal California. Road treatments spanned a period of 32 years, and covered a range of geomorphic and vegetative conditions. SOC decreased with depth at all sites. Treated roads on convex sites exhibited higher SOC than on concave sites, and north aspect sites had higher SOC than south aspect sites. SOC at 5, 20, and 50 cm depths did not differ significantly between decommissioned roads (treated 18–32 years previous) and fully recontoured roads (treated 2–12 years previous). Nevertheless, stepwise multiple regression models project higher SOC developing on fully recontoured roads in the next few decades. The best predictors for SOC on treated roads and in second‐growth forest incorporated aspect, vegetation type, soil depth, lithology, distance from the ocean, years since road treatment (for the road model) and years since harvest (for the forest model). The road model explained 48% of the variation in SOC in the upper 50 cm of mineral soils and the forest model, 54%. Copyright © 2015 John Wiley & Sons, Ltd.     

North polar region of Mars: Advances in stratigraphy, structure, and erosional modification

We have remapped the geology of the north polar plateau on Mars, Planum Boreum, and the surrounding plains of Vastitas Borealis using altimetry and image data along with thematic maps resulting from observations made by the Mars Global Surveyor, Mars Odyssey, Mars Express, and Mars Reconnaissance Orbiter spacecraft. New and revised geographic and geologic terminologies assist with effectively discussing the various features of this region. We identify 7 geologic units making up Planum Boreum and at least 3 for the circumpolar plains, which collectively span the entire Amazonian Period. The Planum Boreum units resolve at least 6 distinct depositional and 5 erosional episodes. The first major stage of activity includes the Early Amazonian (∼3 to 1 Ga) deposition (and subsequent erosion) of the thick (locally exceeding 1000 m) and evenly-layered Rupes Tenuis unit (Abrt), which ultimately formed approximately half of the base of Planum Boreum. As previously suggested, this unit may be sourced by materials derived from the nearby Scandia region, and we interpret that it may correlate with the deposits that regionally underlie pedestal craters in the surrounding lowland plains. The second major episode of activity during the Middle to Late Amazonian (1 Ga) began with a section of dark, sand-rich and light-toned ice-rich irregularly-bedded sequences (Planum Boreum cavi unit, Abb_c) along with deposition of evenly-bedded light-toned ice- and moderate-toned dust-rich layers (Planum Boreum 1 unit, Abb₁). These units have transgressive and gradational stratigraphic relationships. Materials in Olympia Planum underlying the dunes of Olympia Undae are interpreted to consist mostly of the Planum Boreum cavi unit (Abb_c). Planum Boreum materials were then deeply eroded to form spiral troughs, Chasma Boreale, and marginal scarps that define the major aspects of the polar plateau's current regional topography. Locally- to regionally-extensive (though vertically minor) episodes of deposition of evenly-bedded, light- and dark-toned layered materials and subsequent erosion of these materials persisted throughout the Late Amazonian. Sand saltation, including dune migration, is likely to account for much of the erosion of Planum Boreum, particularly at its margin, alluding to the lengthy sedimentological history of the circum-polar dune fields. Such erosion has been controlled largely by topographic effects on wind patterns and the variable resistance to erosion of materials (fresh and altered) and physiographic features. Some present-day dune fields may be hundreds of kilometers removed from possible sources along the margins of Planum Boreum, and dark materials, comprised of sand sheets, extend even farther downwind. These deposits also attest to the lengthy period of erosion following emplacement of the Planum Boreum 1 unit. We find no evidence for extensive glacial flow, topographic relaxation, or basal melting of Planum Boreum materials. However, minor development of normal faults and wrinkle ridges may suggest differential compaction of materials across buried scarps. Timing relations are poorly-defined mostly because resurfacing and other uncertainties prohibit precise determinations of surface impact crater densities. The majority of the stratigraphic record may predate the recent (<20 Ma) part of the orbitally-driven climate record that can be reliably calculated. Given the strong stratigraphic but loose temporal constraints of the north polar geologic record, a comparison of north and south polar stratigraphy permits a speculative scenario in which major Amazonian depositional and erosional episodes driven by global climate activity is plausible.     

New Ways of Using Well Construction Reports for Hydrostratigraphic Analyses

Hydrostratigraphic models require subsurface data. Those data can be gleaned from Well Construction Reports (WCRs) if the reports are properly reviewed and processed. In the United States, WCRs or similar documents are generally required to be submitted to a state agency. The quality of those reports varies considerably and depends on the water well contractor's level of geologic knowledge. Nevertheless, these types of reports are usually the largest available subsurface data set. We present two examples from Wisconsin, each of which had approximately 2000 WCRs, to demonstrate a variety of new applications of WCRs including obtaining information about the depositional environment; selecting sites to perform surface geophysical measurements; and creating three‐dimensional hydrostratigraphic models for groundwater modeling.