The Ordovician Sierras Pampeanas–Puna basin connection: Basement thinning and basin formation in the Proto-Andean back-arc

The Ordovician Sierras Pampeanas, located in a continental back-arc position at the Proto-Andean margin of southwest Gondwana, experienced substantial mantle heat transfer during the Ordovician Famatina orogeny, converting Neoproterozoic and Early Cambrian metasediments to migmatites and granites. The high-grade metamorphic basement underwent intense extensional shearing during the Early and Middle Ordovician. Contemporaneously, up to 7000 m marine sediments were deposited in extensional back-arc basins covering the pre-Ordovician basement. Extensional Ordovician tectonics were more effective in mid- and lower crustal migmatites than in higher levels of the crust. At a depth of about 13 km the separating boundary between low-strain solid upper and high-strain lower migmatitic crust evolved to an intra-crustal detachment. The detachment zone varies in thickness but does not exceed about 500 m. The formation of anatectic melt at the metamorphic peak, and the resulting drop in shear strength, initiated extensional tectonics which continued along localized ductile shear zones until the migmatitic crust cooled to amphibolite facies P–T conditions. P–T–d–t data in combination with field evidence suggest significant (ca. 52%) crustal thinning below the detachment corresponding to a thinning factor of 2.1. Ductile thinning of the upper crust is estimated to be less than that of the lower crust and might range between 25% and 44%, constituting total crustal thinning factors of 1.7–2.0. While the migmatites experienced retrograde decompression during the Ordovician, rocks along and above the detachment show isobaric cooling. This suggests that the magnitude of upper crustal extension controls the amount of space created for sediments deposited at the surface. Upper crustal extension and thinning is compensated by newly deposited sediments, maintaining constant pressure at detachment level. Thinning of the migmatitic lower crust is compensated by elevation of the crust–mantle boundary. The degree of mechanical coupling between migmatitic lower and solid upper crust across the detachment zone is the main factor controlling upper crustal extension, basin formation, and sediment thickness in the back-arc basin. The initiation of crustal extension in the back-arc, however, crucially depends on the presence of anatectic melt in the middle and lower crust. Consumption of melt and cooling of the lower crust correlate with decreasing deposition rates in the sedimentary basins and decreasing rates of crustal extension.     

Sub-Recent Sexual Populations of Limnocythere inopinata Recorded for the First Time from > 3500 m Altitude on the Tibetan Plateau

<正>The central and northwestern Tibetan Plateau region,also called"Qiangtang Plateau"(30°27'25"-35°39'13"N,83o41'14"-95°10'46"E),is the highest and largest arcticalpine area of the Earth covering approximately 600,000km~2 with altitude ranging from 4600 to 5100 m a.s.l.(Song et al.,2012,Fig.1).Ostracod species of     

Predictability of large interannual Arctic sea-ice anomalies

In projections of twenty-first century climate, Arctic sea ice declines and at the same time exhibits strong interannual anomalies. Here, we investigate the potential to predict these strong sea-ice anomalies under a perfect-model assumption, using the Max-Planck-Institute Earth System Model in the same setup as in the Coupled Model Intercomparison Project Phase 5 (CMIP5). We study two cases of strong negative sea-ice anomalies: a 5-year-long anomaly for present-day conditions, and a 10-year-long anomaly for conditions projected for the middle of the twenty-first century. We treat these anomalies in the CMIP5 projections as the truth, and use exactly the same model configuration for predictions of this synthetic truth. We start ensemble predictions at different times during the anomalies, considering lagged-perfect and sea-ice-assimilated initial conditions. We find that the onset and amplitude of the interannual anomalies are not predictable. However, the further deepening of the anomaly can be predicted for typically 1 year lead time if predictions start after the onset but before the maximal amplitude of the anomaly. The magnitude of an extremely low summer sea-ice minimum is hard to predict: the skill of the prediction ensemble is not better than a damped-persistence forecast for lead times of more than a few months, and is not better than a climatology forecast for lead times of two or more years. Predictions of the present-day anomaly are more skillful than predictions of the mid-century anomaly. Predictions using sea-ice-assimilated initial conditions are competitive with those using lagged-perfect initial conditions for lead times of a year or less, but yield degraded skill for longer lead times. The results presented here suggest that there is limited prospect of predicting the large interannual sea-ice anomalies expected to occur throughout the twenty-first century.     

C‐band SAR backscatter characteristics of Arctic sea and land ice during winter

Abstract

Synthetic Aperture Radar (SAR) data has become an important tool for studies of polar regions, due to high spatial resolution even during the polar night and under cloudy skies. We have studied the temporal variation of sea and land ice backscatter of twenty‐four SAR images from the European Remote Sensing satellite (ERS‐1) covering an area in Lady Ann Strait and Jones Sound, Nunavut, from January to March 1992. The presence of fast ice in Jones Sound and glaciers and ice caps on the surrounding islands provides an ideal setting for temporal backscatter studies of ice surfaces. Sample regions for eight different ice types were selected and the temporal backscatter variation was studied. The observed backscatter values for each ice type characterize the radar signatures of the ice surfaces. This time series of twenty‐four SAR images over a 3‐month period provides new insights into the degree of temporal variability of each surface. Ice caps exhibit the highest backscatter value of ‐3.9 dB with high temporal variability. Valley glacier ice backscatter values decrease with decreasing altitude, and are temporally the most stable, with standard deviations of 0.08–0.10 dB over the 90‐day period.

First‐year ice and lead ice show a negative trend in backscatter values in time and a positive correlation of up to 0.59 with air temperature over the 90‐day period. For first‐year ice and lead ice, episodes of large temperature fluctuations (±12°C) are associated with rapid changes in backscatter values (±2 dB). We attribute the backscatter increase to a temperature‐induced increase in brine volume at the base of the snow pack. Multi‐year ice, conglomerate ice and shore ice are relatively stable over the 3‐month period, with a backscatter variation of only a few dBs. An observed lag time of up to three days between backscatter increase/decrease and air temperature can be attributed to the insulation effect of the snow cover over sea ice. The net range of the backscatter values observed on the most temporally stable surface, valley glacier ice, of about 0.30 dB indicates that the ERS‐1 SAR instrument exceeds the 1 dB calibration accuracy specified for the Alaska SAR Facility processor for the three winter months.     

Affine distortion of small GPS networks with large height differences

GPS is a promising tool for real-time monitoring of deformations of slopes or large structures. However, remaining systematic effects in GPS phase observations after double differencing and application of a priori models affect the resulting coordinates. They complicate the proper separation of the actual deformations from pseudo-deformations induced by the systematic effects. This paper shows that for small monitoring networks (baseline lengths <5 km) only affine distortions of the network geometry are generated by the remaining distance dependent systematic effects, e.g. unmodelled tropospheric and ionospheric propagation effects, or satellite orbit errors. Hence, a generic correction model is given by a three-dimensional affine transformation involving a maximum of 12 transformation parameters. For the determination of these parameters, four high quality GPS stations are necessary which are not affected by the actual deformations to be monitored. Based on the analysis of network geometries of synthetic GPS networks with large height differences and considering the physics of the GPS observations it is shown, however, that less than 12 parameters are sufficient for the computation of the corrections. The proposed 8 parameter model was applied to the GPS monitoring network of the Gradenbach landslide. For this small network with large height differences, it was shown that the distortions can be reduced by about 75%.     

Common-Reflection-Surface Stack for Converted Waves

The finite-offset (FO) common-reflection-surface (CRS) stack has been shown to be able to handle not only P-P or S-S but also arbitrarily converted reflections. It can provide different stack sections such as common-offset (CO), common-midpoint (CMP) and common-shot (CS) sections with significantly increased signal-to-noise ratio from the multi-coverage pre-stack seismic data in a data-driven way. It is our purpose in this paper to demonstrate the performance of the FO CRS stack on data involving converted waves in inhomogeneous layered media. In order to do this we apply the FO CRS stack for common-offset to a synthetic seismic data set involving P-P as well as P-S converted primary reflections. We show that the FO CRS stack yields convincing improvement of the image quality in the presence of noisy data and successfully extracts kinematic wavefield attributes useful for further analyses. The extracted emergence angle information is used to achieve a complete separation of the wavefield into its P-P and P-S wave components, given the FO CRS stacked horizontal and vertical component sections.     

A New Equation Solver for Modeling Turbulent Flow in Coupled Matrix‐Conduit Flow Models

Karst aquifers represent dual flow systems consisting of a highly conductive conduit system embedded in a less permeable rock matrix. Hybrid models iteratively coupling both flow systems generally consume much time, especially because of the nonlinearity of turbulent conduit flow. To reduce calculation times compared to those of existing approaches, a new iterative equation solver for the conduit system is developed based on an approximated Newton–Raphson expression and a Gauß–Seidel or successive over‐relaxation scheme with a single iteration step at the innermost level. It is implemented and tested in the research code CAVE but should be easily adaptable to similar models such as the Conduit Flow Process for MODFLOW‐2005. It substantially reduces the computational effort as demonstrated by steady‐state benchmark scenarios as well as by transient karst genesis simulations. Water balance errors are found to be acceptable in most of the test cases. However, the performance and accuracy may deteriorate under unfavorable conditions such as sudden, strong changes of the flow field at some stages of the karst genesis simulations.     

Data Analysis and Modeling to Optimize Thermal Treatment Cost and Performance

The objective of in situ thermal treatment is typically to reduce the contaminant mass or average soil concentration below a specified value. Evaluation of whether the objective has been met is usually made by averaging soil concentrations from a limited number of soil samples. Results from several field sites indicate large performance uncertainty using this approach, even when the number of samples is large. We propose a method to estimate average soil concentration by fitting a log normal probability model to thermal mass recovery data. A statistical approach is presented for making termination decisions from mass recovery data, soil sample data, or both for an entire treatment volume or for subregions that explicitly considers estimation uncertainty which is coupled to a stochastic optimization algorithm to identify monitoring strategies to meet objectives with minimum expected cost. Early termination of heating in regions that reach cleanup targets sooner enables operating costs to be reduced while ensuring a high likelihood of meeting remediation objectives. Results for an example problem demonstrate that significant performance improvement and cost reductions can be achieved using this approach.     

Cyclicity in ridge patterns on the prograding coasts of Estonia

Recent LiDAR surveys have revealed that on postglacially uplifting coasts of Estonia rhythmic coastal landforms (beach ridge sequences and foredune plains) occur to a considerable extent. We studied four of them to reveal age and periodicity in these multiple ridge systems and discussed their genesis in the Subatlantic (semi‐continental) conditions of the Baltic Sea area. Using recent models of Fennoscandian uplift due to glacial isostatic adjustment (GIA), we constructed Holocene apparent sea level curves for the study sites at Õngu, Mänspe, Haldi and Keibu; converted distance–height shore profiles into time series (including corrections on shore profile non‐linearity and variations in GIA‐eustasy balance); and analysed the patterns using spectral analysis. It was suggested that due to non‐tidal conditions, relatively low‐energy hydrodynamic forcing and small aeolian contribution, the ridges mark ancient shorelines. They are relatively modest in height (mostly 0.2–1 m), form regular and extensive (up to 150 ridges) patterns, and date back to ~9000 years before present. We studied <5‐ka‐old sections. The mean ridge spacing varied, depending also on coastal slope, between 21 and 39 m. Both simple counting and spectral analysis involved some specific limitations, yet the estimates for typical spacing were alike, at 32 (±5) years. The regular nature of the low‐ridge patterns originated from relative sea level lowering and gradual sediment accretion/erosion. However, the progradation was rather uplift‐ than accretion‐driven and the stepwise process in ridge formation was probably not autocyclic. It was governed or modulated by quasi‐periodic 25–40 year cyclicity in local wave forcing, relative sea level variations and wind conditions. Being most likely connected to the North Atlantic Oscillation, the quasi‐regular, decadal‐scale, similarly phased variations may magnify each other's effect on the westerly exposed coasts of Estonia. Additionally, some other (e.g. event‐driven) mechanisms may also be present.     

Pleistocene mass-flow deposits of basaltic hyaloclastite on a shallow submarine shelf,South Iceland

A Pleistocene subaqueous, volcanic sequence in South Iceland consists of flows of basaltic hyaloclastite and lava with interbedded sedimentary diamictite units. Emplacement occurred on a distal submarine shelf in drowned valleys along the southern coast of Iceland. The higher sea level was caused by eustatic sea-level change, probably towards the end of a glaciation. This sequence, nearly 700 m thick, rests unconformably on eroded flatlying lavas and sedimentary rocks of likely Tertiary age. A Standard Depositional Unit, describing the flows of hyaloclastite, starts with compact columnar-jointed basalt overlain by cubejointed basalt, and/or pillow lava. This in turn is overlain by thick unstructured hyaloclastite containing aligned basalt lobes, and bedded hyaloclastite at the top. A similar lithofacies succession is valid for proximal to distal locations. The flows were produced by repeated voluminous extrusions of basaltic lava from subaquatic fissures on the Eastern Rift Zone of Iceland. The fissures are assumed to lie in the same general area as the 1783 Laki fissure which produced 12 km³ of basaltic lava. Due to very high extrusion rates, the effective water/melt ratio was low, preventing optimal fragmentation of the melt. The result was a heterogeneous mass of hyaloclastite and fluid melt which moved en masse downslope with the melt at the bottom of the flow and increasingly vesicular hyaloclastite fragments above. The upper and distal parts of the flow moved as low-concentration turbulent suspensions that deposited bedded hyaloclastite.