Surface uplift and climate change - the geomorphic evolution of the western escarpment of the Andes of northern Chile between the Miocene and present

The Western Escarpment of the Andes at 18.30°S (Arica area, northern Chile) is a classical example for a transient state in landscape evolution. This part of the Andes is characterized by the presence of >10,000 km2 plains that formed between the Miocene and the present, and >1500 m deeply incised valleys. Although processes in these valleys scale the rates of landscape evolution, determinations of ages of incision, and more importantly, interpretations of possible controls on valley formation have been controversial. This paper uses morphometric data and observations, stratigraphic information, and estimates of sediment yields for the time interval between ca. 7.5 Ma and present to illustrate that the formation of these valleys was driven by two probably unrelated components. The first component is a phase of base-level lowering with magnitudes of∼300–500 m in the Coastal Cordillera. This period of base-level change in the Arica area, that started at ca. 7.5 Ma according to stratigraphic data, caused the trunk streams to dissect headward into the plains. The headward erosion interpretation is based on the presence of well-defined knickzones in stream profiles and the decrease in valley widths from the coast toward these knickzones. The second component is a change in paleoclimate. This interpretation is based on (1) the increase in the size of the largest alluvial boulders (from dm to m scale) with distal sources during the last 7.5 m.y., and (2) the calculated increase in minimum fluvial incision rates of ∼0.2 mm/yr between ca. 7.5 Ma and 3 Ma to ∼0.3 mm/yr subsequently. These trends suggest an increase in effective water discharge for systems sourced in the Western Cordillera (distal source). During the same time, however, valleys with headwaters in the coastal region (local source) lack any evidence of fluvial incision. This implies that the Coastal Cordillera became hyperarid sometime after 7.5 Ma. Furthermore, between 7.5 Ma and present, the sediment yields have been consistently higher in the catchments with distal sources (∼15 m/m.y.) than in the headwaters of rivers with local sources (<7 m/m.y.). The positive correlation between sediment yields and the altitude of the headwaters (distal versus local sources) seems to reflect the effect of orographic precipitation on surface erosion. It appears that base-level change in the coastal region, in combination with an increase in the orographic effect of precipitation, has controlled the topographic evolution of the northern Chilean Andes.

Evaluating stress patterns for Plio-Quaternary brittle deformation along the Calama-Olacapato-El Toro fault zone, Central Andes

Understanding the geometry and kinematics of the major structures of an orogen is important to elucidate its style of deformation, as well as its tectonic evolution. We describe the temporal and spatial changes in the state of stress of the trans-orogen area of the Calama-Olacapato-El Toro (COT) Fault Zone in the Central Andes, at about 24°S within the northern portion of the Puna Plateau between the Argentina-Chile border. The importance of the COT derives principally from the Quaternary-Holocene activity recognized on some segments, which may shed new light on its possible control on Quaternary volcanism and on the seismic hazard evaluation of the area. Field geological surveys along with kinematic analysis and numerical inversion of ∼140 new fault-slip measurements have revealed that this portion of the COT zone, previously considered a continuous, long-lived lineament, in reality has been subjected to three different kinematic regimes: 1) a Miocene transpressional phase with the maximum principal stress (σ1) chiefly trending NNE-SSW; 2) an extensional phase that started by 9 Ma, with a horizontal NW-SE-striking minimum principal stress (σ3) – permutations between σ2 and σ3 axes have been recognized at two sites – and 3) a left-lateral strike-slip phase with a horizontal ∼E-W &sigma1 and ∼N-S σ3 dating to the Late Pliocene-Quaternary. Spatially, in the Quaternary, the left-lateral component decreases toward the westernmost tip of the COT, where it transitions to extension; this produced to a N-S horst and graben structure. Hence, even if transcurrence is still active in the eastern portion of the COT, as focal mechanisms of crustal earthquakes indicate, our study demonstrates that extension is becoming the predominant structural style of deformation, at least in the western region. These major temporal and spatial changes in the tectonic regimes are attributed in part to changes in the magnitude of the boundary forces due to subduction processes. The overall orogen-perpendicular extension might be the result of vertical stress larger than both the horizontal stresses induced by gravitational effect of a thickened crust.

USING AUTO- AND CROSS-CORRELATIONS FROM SEISMIC NOISE TO MONITOR VELOCITY CHANGES AT VILLARRICA VOLCANO, CHILE

We used the Green's functions from auto-correlations and cross-correlations of seismic ambient noise to monitor temporal velocity changes in the subsurface at Villarrica volcano in the Southern Andes of Chile. Campaigns were conducted from March to October 2010 and February to April 2011 with 8 broadband and 6 short-period stations, respectively. We prepared the data by removing the instrument response, normalizing with a root-mean-square method, whitening the spectra, and filtering from 1 to 10 Hz. This frequency band was chosen based on the relatively high background noise level in that range. Hour-long auto- and cross-correlations were computed and the Green's functions stacked by day and total time. To track the temporal velocity changes we stretched a 24 hour moving window of correlation functions from 90% to 110% of the original and cross correlated them with the total stack. All of the stations' auto-correlations detected what is interpreted as an increase in velocity in 2010, with an average increase of 0.13%. Cross-correlations from station V01, near the summit, to the other stations show comparable changes that are also interpreted as increases in velocity. We attribute this change to the closing of cracks in the subsurface due either to seasonal snow loading or regional tectonics. In addition to the common increase in velocity across the stations, there are excursions in velocity on the same order lasting several days. Amplitude decreases as the station's distance from the vent increases suggesting these excursions may be attributed to changes within the volcanic edifice. In at least two occurrences the amplitudes at stations V06 and V07, the stations farthest from the vent, are smaller. Similar short temporal excursions were seen in the auto-correlations from 2011, however, there was little to no increase in the overall velocity.

Great Falls Tectonic Zone: A Perspective from Idaho

An overview of geologic constraints on the age and extent of the Great Falls tectonic zone, a northeast-trending set of faults identified in west-central Montana.

The Geology of Central Chile's Long-Lived Volcanic Arc

Discussion of subduction zones that are associated with volcanic arcs and major chemical exchanges between the Earth's surface and underlying mantle.

Dating brittle tectonic movements with cleft monazite: Fluid-rock interaction and formation of REE minerals

[1] Two millimeter-sized hydrothermal monazites from an open fissure (cleft) that developed late during a dextral transpressional deformation event in the Aar Massif, Switzerland, have been investigated using electron microprobe and ion probe. The monazites are characterized by high Th/U ratios typical of other hydrothermal monazites. Deformation events in the area have been subdivided into three phases: (D1) main thrusting including formation of a new schistosity, (D2) dextral transpression, and (D3) local crenulation including development of a new schistosity. The two younger deformational structures are related to a subvertically oriented intermediate stress axis, which is characteristic for strike slip deformation. The inferred stress environment is consistent with observed kinematics and the opening of such clefts. Therefore, the investigated monazite-bearing cleft formed at the end of D2 and/or D3, and during dextral movements along NNW dipping planes. Interaction of cleft-filling hydrothermal fluid with wall rock results in rare earth element (REE) mineral formation and alteration of the wall rock. The main newly formed REE minerals are Y-Si, Y-Nb-Ti minerals, and monazite. Despite these mineralogical changes, the bulk chemistry of the system remains constant and thus these mineralogical changes require redistribution of elements via a fluid over short distances (centimeter). Low-grade alteration enables local redistribution of REE, related to the stability of the accessory phases. This allows high precision isotope dating of cleft monazite. 232Th/208Pb ages are not affected by excess Pb and yield growth domain ages between 8.03 ± 0.22 and 6.25 ± 0.60 Ma. Monazite crystallization in brittle structures is coeval or younger than 8 Ma zircon fission track data and hence occurred below 280°C.

Deciphering the driving forces of erosion rates on millennial to million-year timescales in glacially impacted landscapes: An example from the Western Alps

In many regions, tectonic uplift is the main driver of erosion over million-year (Myr) timescales, but climate changes can markedly affect the link between tectonics and erosion, causing transient variations in erosion rates. Here we study the driving forces of millennial to Myr-scale erosion rates in the French Western Alps, as estimated from in situ produced cosmogenic 10Be and a newly developed approach integrating detrital and bedrock apatite fission-track thermochronology. Millennial erosion rates from 10Be analyses vary between ~0.27 and ~1.33 m/kyr, similar to rates measured in adjacent areas of the Alps. Significant positive correlations of millennial erosion rates with geomorphic measures, in particular with the LGM ice thickness, reveal a strong transient morphological and erosional perturbation caused by repeated Quaternary glaciations. The perturbation appears independent of Myr-scale uplift and erosion gradients, with the effect that millennial erosion rates exceed Myr-scale erosion rates only in the internal Alps where the latter are low (<0.4 km/Myr). These areas, moreover, exhibit channels that clearly plot above a general linear positive relation between Myr-scale erosion rates and normalized steepness index. Glacial erosion acts irrespective of rock uplift and thus not only leads to an overall increase in erosion rates but also regulates landscape morphology and erosion rates in regions with considerable spatial gradients in Myr-scale tectonic uplift. Our study demonstrates that climate change, e.g., through occurrence of major glaciations, can markedly perturb landscape morphology and related millennial erosion rate patterns, even in regions where Myr-scale erosion rates are dominantly controlled by tectonics.

Upper crustal deformation in continent-continent collision: A case study from the Bernard nappe complex (Valais, Switzerland)

The Penninic nappes in the Swiss Alps formed during continental collision between the Adriatic and European plates in Cenozoic times. Although intensely studied, the finite geometry of the basement-bearing Penninic nappes in western Switzerland has remained a matter of debate for decades (e.g., “Siviez-Mischabel dilemma”) and the paleogeographic origin of various nappes has been disputed. Here, we present new structural data for the central part of the Penninic Bernard nappe complex, which contains pre-Permian basement and Permo-Mesozoic metasedimentary units. Our lithological and structural observations indicate that the discrepancy between the different structural models proposed for the Bernard nappe complex can be explained by a lateral discontinuity. In the west, the presence of a Permian graben caused complex isoclinal folding, whereas in the east, the absence of such a graben resulted mainly in imbricate thrusting. The overall geometry of the Bernard nappe complex is the result of three main deformation phases: (1) detachment of Mesozoic cover sediments along Triassic evaporites (Evolène phase) during the early stages of collision, (2) Eocene top-to-the-N(NW) nappe stacking (Anniviers phase), and (3) subsequent backfolding and backshearing (Mischabel phase). The southward localized backshearing is key to understand the structural position and paleogeographic origin of units, such as the Frilihorn and Cimes Blanches “nappes” and the Antrona ophiolites. Based on these observations, we present a new tectonic model for the entire Penninic region of western Switzerland and discuss this model in terms of continental collision zone processes.

Softening the Lower Crust: Modes of Syn-Transport Transposition Around and Adjacent to a Deep Crustal Granulite Nappe, Parry Sound Domain, Grenville Province, Ontario, Canada

The Parry Sound domain is a granulite nappe-stack transported cratonward during reactivation of the ductile lower and middle crust in the late convergence of the Mesoproterozoic Grenville orogeny. Field observations suggest the following with respect to the ductile sheath: (1) Formation of a carapace of transposed amphibolite facies gneiss derived from and enveloping the western extremity of the Parry Sound domain and separating it from high-strain gneiss of adjacent allochthons. This ductile sheath formed dynamically around the moving granulite nappe through the development of systems of progressively linked shear zones. (2) Transposition initiated by hydration (amphibolization) of granulite facies gneiss by introduction of fluid along cracks accompanying pegmatite emplacement. Shear zones nucleated along pegmatite margins and subsequently linked and rotated. The source of the pegmatites was most likely subjacent migmatitic and pegmatite-rich units or units over which Parry Sound domain was transported. Comparison of gneisses of the ductile sheath with high-strain layered gneiss of adjacent allochthons show the mode of transposition of penetratively layered gneiss depended on whether or not the gneiss protoliths were amphibolite or granulite facies tectonites before initiation of transposition, resulting in, e.g., folding before shearing, no folding before shearing, respectively. Meter-scale truncation along high-strain gradients at the margins of both types of transposition-related shear zones observed within and marginal to Parry Sound domain mimic features at kilometer scales, implying that apparent truncation by transposition originating in a manner similar to the ductile sheath may be a common feature of deep crustal ductile reworking. Citation: Culshaw, N., C. Gerbi, and J. Marsh (2010), Softening the lower crust: Modes of syn-transport transposition around and adjacent to a deep crustal granulite nappe, Parry Sound domain, Grenville Province, Ontario, Canada, Tectonics, 29, TC5013, doi:10.1029/2009TC002537.

Along-Strike Differences in the Southern Alps of New Zealand: Consequences of Inherited Variation in Rheology

Long- and short-term strain variations along the Australian-Pacific plate boundary through the South Island of New Zealand, including a 300% increase in orogen width, coexistence of oblique thrusting on orthogonal structures, and variability in the locus of orogenic gold deposits, coincide with rheologically relevant geological variation. Our model investigates the consequences of thin, strong lower crust in the north and thick, weak lower crust in the south. Solution of the full 3-D mechanical equations reproduces the larger wavelength strain patterns of the orogen. A 3-D perturbation-based analytical solution leads to the identification of the sensitivity of displacement type to minor stress changes. Transition from boundary-normal thrusting to boundary-parallel thrusting occurs at the transition from strong to weak lower crust and is related to an increase in either tau(yz) (shear stress in the yz plane) or the ratio of the coordinate normal stresses, (sigma(yy)/sigma(xx)), where x and y are in the horizontal and z is vertical. Both mechanisms are compatible with the geologically dependent rheological variation employed in our model. Citation: Upton, P., P. O. Koons, D. Craw, C. M. Henderson, and R. Enlow (2009), Along-strike differences in the Southern Alps of New Zealand: Consequences of inherited variation in rheology, Tectonics, 28, TC2007, doi:10.1029/2008TC002353.

Three-Dimensional Mechanics of Yakutat Convergence in the Southern Alaskan Plate Corner

Three-dimensional numerical models are used to investigate the mechanical evolution of the southern Alaskan plate corner where the Yakutat and the Pacific plates converge on the North American plate. The evolving model plate boundary consists of Convergent, Lateral, and Subduction subboundaries with flow separation of incoming material into upward or downward trajectories forming dual, nonlinear advective thermal/mechanical anomalies that fix the position of major subaerial mountain belts. The model convergent subboundary evolves into two teleconnected orogens: Inlet and Outlet orogens form at locations that correspond with the St. Elias and the Central Alaska Range, respectively, linked to the East by the Lateral boundary. Basins form parallel to the orogens in response to the downward component of velocity associated with subduction. Strain along the Lateral subboundary varies as a function of orogen rheology and magnitude and distribution of erosion. Strain-dependent shear resistance of the plate boundary associated with the shallow subduction zone controls the position of the Inlet orogen. The linkages among these plate boundaries display maximum shear strain rates in the horizontal and vertical planes where the Lateral subboundary joins the Inlet and Outlet orogens. The location of the strain maxima shifts with time as the separation of the Inlet and Outlet orogens increases. The spatiotemporal predictions of the model are consistent with observed exhumation histories deduced from thermochronology, as well as stratigraphic studies of synorogenic deposits. In addition, the complex structural evolution of the St Elias region is broadly consistent with the predicted strain field evolution. Citation: Koons, P. O., B. P. Hooks, T. Pavlis, P. Upton, and A. D. Barker (2010), Three-dimensional mechanics of Yakutat convergence in the southern Alaskan plate corner, Tectonics, 29, TC4008, doi: 10.1029/2009TC002463.

Diagenesis of a light, tight-oil chert reservoir at the Ediacaran/Cambrian boundary, Sultanate of Oman

The Al Shomou Silicilyte Member (Athel Formation) in the South Oman Salt Basin shares many of the characteristics of a light, tight-oil (LTO) reservoir: it is a prolifi c source rock mature for light oil, it produces light oil from a very tight matrix and reservoir, and hydraulic fracking technology is required to produce the oil. What is intriguing about the Al Shomou Silicilyte, and different from other LTO reservoirs, is its position related to the Precambrian/Cambrian Boundary (PCB) and the fact that it is a ‘laminated chert‘ rather than a shale. In an integrated diagenetic study we applied microstructural analyses (SEM, BSE) combined with state-of-the-art stable isotope and trace element analysis of the silicilyte matrix and fractures. Fluid inclusion microthermometry was applied to record the salinity and minimum trapping temperatures. The microstructural investigations reveal a fi ne lamination of the silicilyte matrix with a mean lamina thickness of ca. 20 μm consisting of predominantly organic matter-rich and fi nely crystalline quartz-rich layers, respectively. Authigenic, micron-sized idiomorphic quartz crystals are the main matrix components of the silicilyte. Other diagenetic phases are pyrite, apatite, dolomite, magnesite and barite cements. Porosity values based on neutron density logs and core plug data indicate porosity in the silicilyte ranges from less than 2% to almost to 40%. The majority of the pore space in the silicilyte is related to (primary) inter-crystalline pores, with locally important oversized secondary pores. Pore casts of the silica matrix show that pores are extremely irregular in three dimensions, and are generally interconnected by a complex web or meshwork of fi ne elongate pore throats. Mercury injection capillary data are in line with the microstructural observations suggesting two populations of pore throats, with an effective average modal diameter of 0.4 μm. The acquired geochemical data support the interpretation that the primary source of the silica is the ambient seawater rather than hydrothermal or biogenic. A maximum temperature of ca. 45°C for the formation of microcrystalline quartz in the silicilyte is good evidence that the lithifi cation and crystallization of quartz occurred in the fi rst 5 Ma after deposition. Several phases of brittle fracturing and mineralization occurred in response to salt tectonics during burial. The sequences of fracture-fi lling mineral phases (dolomite - layered chalcedony – quartz – apatite - magnesite I+II - barite – halite) indicates a complex fl uid evolution after silicilyte lithifi cation. Primary, all-liquid fl uid inclusions in the fracturefi lling quartz are good evidence of growth beginning at low temperatures, i.e. ≤ 50ºC. Continuous precipitation during increasing temperature and burial is documented by primary two-phase fl uid inclusions in quartz cements that show brines at 50°C and fi rst hydrocarbons at ca. 70°C. The absolute timing of each mineral phase can be constrained based on U-Pb geochronometry, and basin modelling. Secondary fl uid inclusions in quartz, magnesite and barite indicate reactivation of the fracture system after peak burial temperature during the major cooling event, i.e. uplift, between 450 and 310 Ma. A number of fi rst-order trends in porosity and reservoir-quality distribution are observed which are strongly related to the diagenetic and fl uid history of the reservoir: the early in-situ generation of hydrocarbons and overpressure development arrests diagenesis and preserves matrix porosity. Chemical compaction by pressure dissolution in the fl ank areas could be a valid hypothesis to explain the porosity variations in the silicilitye slabs resulting in lower porosity and poorer connectivity on the fl anks of the reservoir. Most of the hydrocarbon storage and production comes from intervals characterized by Amthor et al. 114488 preserved micropores, not hydrocarbon storage in a fracture system. The absence of oil expulsion results in present-day high oil saturations. The main diagenetic modifi cations of the silicilyte occurred and were completed relatively early in its history, i.e. before 300 Ma. An instrumental factor for preserving matrix porosity is the diffi culty for a given slab to evacuate all the fl uids (water and hydrocarbons), or in other words, the very good sealing capacity of the salt embedding the slab.