Paleo-stress orientations from calcite twins in the North Pyrenean foreland, determined by the Etchecopar inverse method

During the last decade, it has been shown that accurate analyses of calcite twinning in thin section may allow computation of the paleo-deviatoric stress tensor associated with this deformation. Following from the work of Larroque and Laurent (1988), and Lacombe et al. (1989), this is a new application of the Etchecopar inverse method. Weakly deformed limestones from about 150 to 200 km north of the Pyrenean chain (South France) were studied. The Etchecopar inverse method, which can be used on a microcomputer, is comparable to the method of analysis of striated fault planes proposed by Etchecopar et al. in 1981, as it determines the principal stress orientations and the stress ellipsoid shape ratio φ = (σ₂ − σ₃)/(σ₁ − σ₃). The Etchecopar method using calcite twins has been selected for this study because it is suitable for the determination of polyphase deformation in coarse-grained limestones (Lacombe et al. 1989).

Ten samples of oolitic limestones, cemented with sparry calcite, from the Bajocian (170 Ma) and the Oligocene (30 Ma) were studied. For each of them, only one or two stress tensors were determined: the first is characterized by a principal compressive stress axis σ₁ which is horizontal and approximately N-S. Local departures from this orientation have been observed and are interpreted as local stress perturbations. For the second, which was defined in the Oligocene, but occurred only locally in the Jurassic limestones, the position of σ₁ is nearly vertical; the direction of extension σ₃ is either NE-SW or undefined into the horizontal plane (σ₂ = σ₃); this radial extension is interpreted as the result of two or more superposed extensive deformation phases with vertical σ₁. The N-S compression was never found in the Oligocene limestones, and from its orientation, it is assumed to correspond to the Pyrenean deformation. These results are compared to the more complex deformation pattern obtained from classical microtectonical studies (Arthaud and Choukroune, 1972; Bonijoly and Blès, 1983; Granier and Blès, 1988; Blès et al., 1989). Only the major tectonic events are recorded by calcite twinning. Finally, the suitability of calcite twin analyses for the determination of complex polyphase deformation is discussed.     

Construction and analysis of higher order finite difference schemes for the 1D wave equation

In this article, we propose to: 1. Establish most of the properties conjectured in [2] about the higher order finite difference approximation of the 1D Laplace operator. 2. Generalize to any order the fourth-order accurate scheme in space and time of Shubin and Bell [20] and Cohen [6]. For this new family of 2m–2m schemes, we establish, via elementary mathematics, various stability and dispersion results that are helpful to compare these schemes to the 2–2m schemes of Anné et al. [2]. This revised version was published online in July 2006 with corrections to the Cover Date.     

Expanding the Sémiologie Graphique for contemporary cartography,some ideas from visual semiotics,art history and design

ABSTRACT

In this article, we describe some ideas based on Bertin’s Sémiologie Graphique with a focus on extending it for aesthetic quality issues in cartography. We consider research in semiotics, aesthetics, and graphic design, with some fruitful results about the use of color or composition. We present several tools, oriented toward a pedagogical use and accessibility, helping to examine the color composition of an expressing image or to create meaningful gradients. While Bertin’s concepts still remain relevant and efficient, specific ideas can be developed further, to adapt to current media and gain from scientific advances from related disciplines.     

Toward a generic method for studying water renewal, with application to the epilimnion of Lake Tanganyika

We present a method, based on the concept of age and residence time, to study the water renewal in a semi-enclosed domain. We split the water of this domain into different water types. The initial water is the water initially present in the semi-enclosed domain. The renewing water is defined as the water entering the domain of interest. Several renewing water types may be considered depending on their origin. We present the equations for computing the age and the residence time of a certain water type. These timescales are of use to understand the rate at which the water renewal takes place. Computing these timescales can be achieved at an acceptable extra computer cost.The above-mentioned method is applied to study the renewal of epilimnion (i.e. the surface layer) water in Lake Tanganyika. We have built a finite element reduced-gravity model modified to take into account the water exchange between the epilimnion and the hypolimnion (i.e. the bottom layer), the water supply from precipitation and incoming rivers, and the water loss from evaporation and the only outgoing river. With our water renewal diagnoses, we show that the only significant process in the renewal of epilimnion water in Lake Tanganyika is the water exchange between the epilimnion and the hypolimnion, other phenomena being negligible.     

Strategies to mitigate aliasing of loading signals while estimating GPS frame parameters

Although GNSS techniques are theoretically sensitive to the Earth center of mass, it is often preferable to remove intrinsic origin and scale information from the estimated station positions since they are known to be affected by systematic errors. This is usually done by estimating the parameters of a linearized similarity transformation which relates the quasi-instantaneous frames to a long-term frame such as the International Terrestrial Reference Frame (ITRF). It is well known that non-linear station motions can partially alias into these parameters. We discuss in this paper some procedures that may allow reducing these aliasing effects in the case of the GPS techniques. The options include the use of well-distributed sub-networks for the frame transformation estimation, the use of site loading corrections, a modification of the stochastic model by downweighting heights, or the joint estimation of the low degrees of the deformation field. We confirm that the standard approach consisting of estimating the transformation over the whole network is particularly harmful for the loading signals if the network is not well distributed. Downweighting the height component, using a uniform sub-network, or estimating the deformation field perform similarly in drastically reducing the amplitude of the aliasing effect. The application of these methods to reprocessed GPS terrestrial frames permits an assessment of the level of agreement between GPS and our loading model, which is found to be about 1.5 mm WRMS in height and 0.8 mm WRMS in the horizontal at the annual frequency. Aliased loading signals are not the main source of discrepancies between loading displacement models and GPS position time series.     

Ocean heat transport into the Arctic in the twentieth and twenty-first century in EC-Earth

The ocean heat transport into the Arctic and the heat budget of the Barents Sea are analyzed in an ensemble of historical and future climate simulations performed with the global coupled climate model EC-Earth. The zonally integrated northward heat flux in the ocean at 70°N is strongly enhanced and compensates for a reduction of its atmospheric counterpart in the twenty first century. Although an increase in the northward heat transport occurs through all of Fram Strait, Canadian Archipelago, Bering Strait and Barents Sea Opening, it is the latter which dominates the increase in ocean heat transport into the Arctic. Increased temperature of the northward transported Atlantic water masses are the main reason for the enhancement of the ocean heat transport. The natural variability in the heat transport into the Barents Sea is caused to the same extent by variations in temperature and volume transport. Large ocean heat transports lead to reduced ice and higher atmospheric temperature in the Barents Sea area and are related to the positive phase of the North Atlantic Oscillation. The net ocean heat transport into the Barents Sea grows until about year 2050. Thereafter, both heat and volume fluxes out of the Barents Sea through the section between Franz Josef Land and Novaya Zemlya are strongly enhanced and compensate for all further increase in the inflow through the Barents Sea Opening. Most of the heat transported by the ocean into the Barents Sea is passed to the atmosphere and contributes to warming of the atmosphere and Arctic temperature amplification. Latent and sensible heat fluxes are enhanced. Net surface long-wave and solar radiation are enhanced upward and downward, respectively and are almost compensating each other. We find that the changes in the surface heat fluxes are mainly caused by the vanishing sea ice in the twenty first century. The increasing ocean heat transport leads to enhanced bottom ice melt and to an extension of the area with bottom ice melt further northward. However, no indication for a substantial impact of the increased heat transport on ice melt in the Central Arctic is found. Most of the heat that is not passed to the atmosphere in the Barents Sea is stored in the Arctic intermediate layer of Atlantic water, which is increasingly pronounced in the twenty first century.     

Late Pliocene to Early Quaternary extensional detachment in the La Paz–El Cabo area (Baja California Sur, Mexico): implications on the opening of the Gulf of California and the mechanics of oblique rifting

We demonstrate that Pliocene to Early Quaternary sedimentary formations in Baja California Sur (Mexico) were deposited syn-tectonically over a major detachment associated with the exhumation of Mesozoic crust. The detachment dips to the ENE and is associated with E–W stretching. This large extensional structure strikes almost parallel to the general trend of the Gulf of California and extension is oblique to the East-Pacific seafloor-spreading direction. Crustal-scale stretching in this area was still active after the beginning of seafloor spreading c.  3.6 Ma ago. The detachment is capped by Late Pleistocene–Holocene alluvial sediments the deposition of which seems to be partly syn-tectonic and controlled by minor stretching subparallel to the present-day North American–Pacific kinematic vector. We discuss the implications of our observations on strain partitioning during opening of the California Gulf as well as on the structure of the Gulf of California margin.     

Why chromatic imaging matters

During the last two decades, the first generation of beam combiners at the Very Large Telescope Interferometer has proved the importance of optical interferometry for high-angular resolution astrophysical studies in the near- and mid-infrared. With the advent of 4-beam combiners at the VLTI, the u ? v coverage per pointing increases significantly, providing an opportunity to use reconstructed images as powerful scientific tools. Therefore, interferometric imaging is already a key feature of the new generation of VLTI instruments, as well as for other interferometric facilities like CHARA and JWST. It is thus imperative to account for the current image reconstruction capabilities and their expected evolutions in the coming years. Here, we present a general overview of the current situation of optical interferometric image reconstruction with a focus on new wavelength-dependent information, highlighting its main advantages and limitations. As an Appendix we include several cookbooks describing the usage and installation of several state-of-the art image reconstruction packages. To illustrate the current capabilities of the software available to the community, we recovered chromatic images, from simulated MATISSE data, using the MCMC software SQUEEZE. With these images, we aim at showing the importance of selecting good regularization functions and their impact on the reconstruction.     

Decoupling of arsenic and iron release from ferrihydrite suspension under reducing conditions: a biogeochemical model

High levels of arsenic in groundwater and drinking water are a major health problem. Although the processes controlling the release of As are still not well known, the reductive dissolution of As-rich Fe oxyhydroxides has so far been a favorite hypothesis. Decoupling between arsenic and iron redox transformations has been experimentally demonstrated, but not quantitatively interpreted. Here, we report on incubation batch experiments run with As(V) sorbed on, or co-precipitated with, 2-line ferrihydrite. The biotic and abiotic processes of As release were investigated by using wet chemistry, X-ray diffraction, X-ray absorption and genomic techniques. The incubation experiments were carried out with a phosphate-rich growth medium and a community of Fe(III)-reducing bacteria under strict anoxic conditions for two months. During the first month, the release of Fe(II) in the aqueous phase amounted to only 3% to 10% of the total initial solid Fe concentration, whilst the total aqueous As remained almost constant after an initial exchange with phosphate ions. During the second month, the aqueous Fe(II) concentration remained constant, or even decreased, whereas the total quantity of As released to the solution accounted for 14% to 45% of the total initial solid As concentration. At the end of the incubation, the aqueous-phase arsenic was present predominately as As(III) whilst X-ray absorption spectroscopy indicated that more than 70% of the solid-phase arsenic was present as As(V). X-ray diffraction revealed vivianite Fe(II)₃(PO₄)₂.8H₂O in some of the experiments. A biogeochemical model was then developed to simulate these aqueous- and solid-phase results. The two main conclusions drawn from the model are that (1) As(V) is not reduced during the first incubation month with high Eh values, but rather re-adsorbed onto the ferrihydrite surface, and this state remains until arsenic reduction is energetically more favorable than iron reduction, and (2) the release of As during the second month is due to its reduction to the more weakly adsorbed As(III) which cannot compete against carbonate ions for sorption onto ferrihydrite. The model was also successfully applied to recent experimental results on the release of arsenic from Bengal delta sediments.