Transtensional deformation in the Lake Tahoe region, California and Nevada, USA

Dextral transtensional deformation is occurring along the Sierra Nevada–Great Basin boundary zone (SNGBBZ) at the eastern edge of the Sierra Nevada microplate. In the Lake Tahoe region of the SNGBBZ, transtension is partitioned spatially and temporally into domains of north–south striking normal faults and transitional domains with conjugate strike-slip faults. The normal fault domains, which have had large Holocene earthquakes but account only for background seismicity in the historic period, primarily accommodate east–west extension, while the transitional domains, which have had moderate Holocene and historic earthquakes and are currently seismically active, primarily record north–south shortening. Through partitioned slip, the upper crust in this region undergoes overall constrictional strain.Major fault zones within the Lake Tahoe basin include two normal fault zones: the northwest-trending Tahoe–Sierra frontal fault zone (TSFFZ) and the north-trending West Tahoe–Dollar Point fault zone. Most faults in these zones show eastside down displacements. Both of these fault zones show evidence of Holocene earthquakes but are relatively quiet seismically through the historic record. The northeast-trending North Tahoe–Incline Village fault zone is a major normal to sinistral-oblique fault zone. This fault zone shows evidence for large Holocene earthquakes and based on the historic record is seismically active at the microearthquake level. The zone forms the boundary between the Lake Tahoe normal fault domain to the south and the Truckee transition zone to the north.Several lines of evidence, including both geology and historic seismicity, indicate that the seismically active Truckee and Gardnerville transition zones, north and southeast of Lake Tahoe basin, respectively, are undergoing north–south shortening. In addition, the central Carson Range, a major north-trending range block between two large normal fault zones, shows internal fault patterns that suggest the range is undergoing north–south shortening in addition to east–west extension.A model capable of explaining the spatial and temporal partitioning of slip suggests that seismic behavior in the region alternates between two modes, one mode characterized by an east–west minimum principal stress and a north–south maximum principal stress as at present. In this mode, seismicity and small-scale faulting reflecting north–south shortening concentrate in mechanically weak transition zones with primarily strike-slip faulting in relatively small-magnitude events, and domains with major normal faults are relatively quiet. A second mode occurs after sufficient north–south shortening reduces the north–south S_hmax in magnitude until it is less than S_v, at which point S_v becomes the maximum principal stress. This second mode is then characterized by large earthquakes on major normal faults in the large normal fault domains, which dominate the overall moment release in the region, producing significant east–west extension.     

Successful application of multiscale methods in a real reservoir simulator environment

For the past 10 years or so, a number of so-called multiscale methods have been developed as an alternative approach to upscaling and to accelerate reservoir simulation. The key idea of all these methods is to construct a set of prolongation operators that map between unknowns associated with cells in a fine grid holding the petrophysical properties of the geological reservoir model and unknowns on a coarser grid used for dynamic simulation. The prolongation operators are computed numerically by solving localized flow problems, much in the same way as for flow-based upscaling methods, and can be used to construct a reduced coarse-scale system of flow equations that describe the macro-scale displacement driven by global forces. Unlike effective parameters, the multiscale basis functions have subscale resolution, which ensures that fine-scale heterogeneity is correctly accounted for in a systematic manner. Among all multiscale formulations discussed in the literature, the multiscale restriction-smoothed basis (MsRSB) method has proved to be particularly promising. This method has been implemented in a commercially available simulator and has three main advantages. First, the input grid and its coarse partition can have general polyhedral geometry and unstructured topology. Secondly, MsRSB is accurate and robust when used as an approximate solver and converges relatively fast when used as an iterative fine-scale solver. Finally, the method is formulated on top of a cell-centered, conservative, finite-volume method and is applicable to any flow model for which one can isolate a pressure equation. We discuss numerical challenges posed by contemporary geomodels and report a number of validation cases showing that the MsRSB method is an efficient, robust, and versatile method for simulating complex models of real reservoirs.     

Constraining sequence stratigraphy in north Australian basins: SHRIMP U–Pb zircon geochronology between Mt Isa and McArthur River

Fifty‐five new SHRIMP U–Pb zircon ages from samples of northern Australian ‘basement’ and its overlying Proterozoic successions are used to refine and, in places, significantly change previous lithostratigraphic correlations. In conjunction with sequence‐stratigraphic studies, the 1800–1580 Ma rock record between Mt Isa and the Roper River is now classified into three superbasin phases—the Leichhardt, Calvert and Isa. These three major depositional episodes are separated by ～20 million years gaps. The Isa Superbasin can be further subdivided into seven supersequences each 10–15 million years in duration. Gaps in the geological record between these supersequences are variable; they approach several million years in basin‐margin positions, but are much smaller in the depocentres. Arguments based on field setting, petrography, zircon morphology, and U–Pb systematics are used to interpret these U–Pb zircon ages and in most cases to demonstrate that the ages obtained are depositional. In some instances, zircon crystals are reworked and give maximum depositional ages. These give useful provenance information as they fingerprint the source(s) of basin fill. Six new ‘Barramundi’ basement ages (around 1850 Ma) were obtained from crystalline units in the Murphy Inlier (Nicholson Granite and Cliffdale Volcanics), the Urapunga Tectonic Ridge (‘Mt Reid Volcanics’ and ‘Urapunga Granite’), and the central McArthur Basin (Scrutton Volcanics). New ages were also obtained from units assigned to the Calvert Superbasin (ca 1740–1690 Ma). SHRIMP results show that the Wollogorang Formation is not one continuous unit, but two different sequences, one deposited around 1730 Ma and a younger unit deposited around 1722 Ma. Further documentation is given of a regional 1725 Ma felsic event adjacent to the Murphy Inlier (Peters Creek Volcanics and Packsaddle Microgranite) and in the Carrara Range. A younger ca 1710 Ma felsic event is indicated in the southwestern McArthur Basin (Tanumbirini Rhyolite and overlying Nyanantu Formation). Four of the seven supersequences in the Isa Superbasin (ca 1670–1580 Ma) are reasonably well‐constrained by the new SHRIMP results: the Gun Supersequence (ca 1670–1655 Ma) by Paradise Creek Formation, Moondarra Siltstone, Breakaway Shale and Urquhart Shale ages grouped between 1668 and 1652 Ma; the Loretta Supersequence (ca 1655–1645 Ma) by results from the Lady Loretta Formation, Walford Dolomite, the upper part of the Mallapunyah Formation and the Tatoola Sandstone between ca 1653 and 1647 Ma; the River Supersequence (ca 1645–1630 Ma) by ages from the Teena Dolomite, Mt Les and Riversleigh Siltstones, and Barney Creek, Lynott, St Vidgeon and Nagi Formations clustering around 1640 Ma; and the Term Supersequence (ca 1630–1615 Ma) by ages from the Stretton Sandstone, lower Doomadgee Formation and lower part of the Lawn Hill Formation, mostly around 1630–1620 Ma. The next two younger supersequences are less well‐constrained geochronologically, but comprise the Lawn Supersequence (ca 1615–1600 Ma) with ages from the lower Balbirini Dolomite, and lower Doomadgee, Amos and middle Lawn Hill Formations, clustered around 1615–1610 Ma; and the Wide Supersequence (ca 1600–1585 Ma) with only two ages around 1590 Ma, one from the upper Balbirini Dolomite and the other from the upper Lawn Hill Formation. The Doom Supersequence (<1585 Ma) at the top of the Isa Superbasin is essentially unconstrained. The integration of high‐precision SHRIMP dating from continuously analysed stratigraphic sections, within a sequence stratigraphic context, provides an enhanced chronostratigraphic framework leading to more reliable interpretations of basin architecture and evolution.     

Re-interpretation of a mixed petrographic province — The Sintra intrusive complex (Portugal) and related rocks

Sub-alkaline to calc-alkaline granite dominates the high-level plutonic massif of the Serra de Sintra, while two overlapping ring complexes recognised within it consist of four distinct alkaline to sub-alkaline quartz-bearing syenitic units, together with intrusions of a remarkable kaersutite- and (titan?)biotite-bearing theralite heteromorph (Mafraite), contaminated to varying degrees by sialic admixture. Large volumes of sub-volcanic breccia testify to considerable late explosive activity in and near the ring centres. Minor intrusives are approximate compositional counterparts of the plutonic rocks, but more closely resemble lavas forming the extensive “basaltic mantle” of the Lisbon region, where compositions range from olivine basalt, basanite and theralite (luscladite) to subordinate quartz-bearing trachyandesite and trachyte. Undersaturated salic rocks are absent throughout the province. Magmatic differentiation probably played a minimal part in petrogenesis of the Sintra intrusives. The granite is interpreted as a direct product of basement syntexis, the other rocks as products of selective palingenetic reactions between mafraitic magma and crystalline basement (or granite). Evolution of the volcanic and hypabyssal rocks may have been controlled by similar reactions, but evidence favouring differentiation of this suite from sialcontaminated basic parents is almost as strong. There are some points of resemblance between this province and that of the Auvergne region in southern France.     

Paleoclimatic changes in southern Africa during the intensification of Northern Hemisphere glaciation: evidence from ODP Leg 175 Site 1085

Pliocene age sediments from Ocean Drilling Program Leg 175, Site 1085-A and B in the Cape Basin were analyzed to investigate the impact of the intensification of Northern Hemisphere glaciation (INHG) on the South Atlantic Benguela Current system from 4 to 2 Ma. Proxies for productivity (concentrations and mass accumulation rates of total organic carbon, carbon to nitrogen ratios, percent calcium carbonate, and percent biogenic silica) as well as weight percent sand (a proxy for preservation or winnowing) peak at 3.2, 3.0, 2.4, and 2.25 Ma. Normative calculations of allied trace and major elemental determinations indicate synchronous increases in productivity peaks, as well as high concentrations and accumulations of terrigenous sediments. Coeval increases in hemipelagic sedimentation and productivity indicators could be the result of enhanced eolian sedimentation resulting from strengthened winds, leading to elevated rates of upwelling and enhanced productivity. However, rapid burial, as indicated by high sedimentation rates, could also enhance preservation. The very high concentrations (>30%) and accumulations (up to 60 g/cm²/kyr) limit the likelihood that eolian sedimentation was the only transport mechanism, invoking an additional fluvial source. Rapid burial by either eolian or fluvial transport links these intervals of enhanced preservation and productivity with continental climate changes resulting from (1) increased winds and/or dust availability due to higher aridity in the Namibia/northern South Africa region; (2) lowered sea-level related to increased ice volume; (3) increased sediment load due to wetter conditions in the continental interior; or (4) some combination. Peaks at 3.2, 2.4 and 2.25 Ma are coincident with maximum precession, suggesting a link between hemipelagic sedimentation and enhanced monsoonal circulation over southern Africa. The Site 1085 sedimentary record during the INHG seems to be controlled by low-latitude processes linked to precession rather than hig-latitude processes.     

The cosmic ray energy spectrum between 10 and 10 eV

The energy spectrum of cosmic rays with primary energies between 10¹⁴ eV and 10¹⁶ eV has been studied with the CASA-MIA air shower array. The measured differential energy spectrum is a power law (dj/dE ∝ E^−y) with spectral indices γ of 2.66±0.02 below approximately 10¹⁵ eV and 3.00±0.05 above. A new method is used for measuring primary energy derived from ground-based data in a compositionally insensitive way. In contrast with some previous reports, the “knee” of the energy spectrum does not appear sharp, but rather a smooth transition over energies from 10¹⁵ eV to 3.0 × 10¹⁵ eV.     

Lithological and Structural Controls on Regional 3-D Fluid Flow Patterns during Greenschist Facies Metamorphism of the Dalradian of the SW Scottish Highlands

Greenschist facies metabasite sills enclosed within phyllitesand psammites of the Dalradian Supergroup of the SW ScottishHighlands have been used as quantitative sensors of fluid flow,from which ‘fossil’ fluid flow-paths have been mapped.Infiltration of H₂O–CO₂ fluid caused carbonation of themetabasite assemblage. Consequently, metabasite sills have developeda characteristic asymmetric zonation of reacted (carbonate-bearing)margins and unreacted (carbonate-absent) interiors. Reactionfront advection theory has been used to constrain the time-integratedfluid flux and flow direction in one dimension, for each sill,from the widths of the reacted margins. Three-dimensional fluxesand flow directions have been constrained from one-dimensionalfluxes by geometrical computation, and regional fluid flow patternshave been deduced. On a regional scale, fluid flow was channelled through phyllitesand focused towards the axial zone of major antiformal foldstructures. Fluid flow was channelled through the 6 km thickArdrishaig Phyllites, within which fluid fluxes were found toincrease exponentially with proximity to the axial surface ofthe Ardrishaig Anticline. Fluid flow through the psammites waslimited to narrow zones of intense deformation associated withmetabasic or phyllitic intercalations. Although time-integratedfluid fluxes were locally high (4 10² m³/m²), average fluxes(100 m¹/m²) are compatible with availability of fluid from underlyingdevolatilizing rocks. Complete isotopic homogenization of phyllites,psammites and metabasites within the axial region of the ArdrishaigAnticline is predicted by these high measured fluxes. KEY WORDS: 3-D flow; Dalradian; SW Scotland; greenschist facies; metabasites ^*Corresponding author. Telephone: 0131-650-5885 Fax: 0131-668-3184 E-mail: askelton{at}glg.ed.ac.uk     

Climate Change Impacts for the Conterminous USA: An Integrated Assessment

Here we simulate dryland agriculture in the United States in order to assess potential future agricultural production under a set of general circulation model (GCM)-based climate change scenarios. The total national production of three major grain crops—corn, soybeans, and winter wheat—and two forage crops—alfalfa and clover hay—is calculated for the actual present day core production area (CPA) of each of these crops. In general, higher global mean temperature (GMT) reduces production and higher atmospheric carbon dioxide concentration ([CO₂]) increases production. Depending on the climatic change scenarios employed overall national production of the crops studied changes by up to plus or minus 25% from present-day levels. Impacts are more significant regionally, with crop production varying by greater than ±50% from baseline levels. Analysis of currently possible production areas (CPPAs) for each crop indicates that the regions most likely to be affected by climate change are those on the margins of the areas in which they are currently grown. Crop yield variability was found to be primarily influenced by local weather and geographic features rather than by large-scale changes in climate patterns and atmospheric composition. Future US agronomic potential will be significantly affected by the changes in climate projected here. The nature of the crop response will depend primarily on to what extent precipitation patterns change and also on the degree of warming experienced.     

An analysis of the components of the luminescence signals of selected polymineral and quartz samples from loess in western China and southern Tajikistan,and their suitability for optical dating

Previous luminescence dating studies on loess from China and Tajikistan have focused on the establishment of the regional chronology of the loess sequences. In order to improve the precision and accuracy of optical ages derived from the loess of the last glacial period in these regions, we have examined the components of luminescence signals in three loess samples from western China and southern Tajikistan. Our results show that the polymineral IRSL and post-IR OSL, and quartz OSL signals from loess of the two regions are represented by three components, which display different bleaching and growth characteristics. While the composition of the polymineral IRSL signals is similar between samples with the same age from the two regions, in the case of quartz there is significant discrepancy in the proportion of the fast and medium components of the OSL signals. Greater difference is observed in the composition of the polymineral post-IR OSL signals for the loess from the two regions. The three components of polymineral IRSL signals yield almost identical equivalent dose values as that derived from the total IRSL signal. An apparent agreement in equivalent dose is observed between the fast component of the polymineral post-IR OSL and the quartz OSL for the loess of western China but not in the loess of the same age from southern Tajikistan. The fast component of the quartz OSL yields an equivalent dose 25% higher than that based on the total signal for the sample from the base of the Late Pleistocene loess in southern Tajikistan. This demonstrates the importance of signal selection for an accurate luminescence dating of Central Asian loess.