Geodynamic evolution of a forearc rift in the southernmost Mariana Arc

The southernmost Mariana forearc stretched to accommodate opening of the Mariana Trough backarc basin in late Neogene time, erupting basalts at 3.7–2.7 Ma that are now exposed in the Southeast Mariana Forearc Rift (SEMFR). Today, SEMFR is a broad zone of extension that formed on hydrated, forearc lithosphere and overlies the shallow subducting slab (slab depth ≤ 30–50 km). It comprises NW–SE trending subparallel deeps, 3–16 km wide, that can be traced ≥ ∼30 km from the trench almost to the backarc spreading center, the Malaguana‐Gadao Ridge (MGR). While forearcs are usually underlain by serpentinized harzburgites too cold to melt, SEMFR crust is mostly composed of Pliocene, low‐K basaltic to basaltic andesite lavas that are compositionally similar to arc lavas and backarc basin (BAB) lavas, and thus defines a forearc region that recently witnessed abundant igneous activity in the form of seafloor spreading. SEMFR igneous rocks have low Na₈, Ti₈, and Fe₈, consistent with extensive melting, at ∼23 ± 6.6 km depth and 1239 ± 40°C, by adiabatic decompression of depleted asthenospheric mantle metasomatized by slab‐derived fluids. Stretching of pre‐existing forearc lithosphere allowed BAB‐like mantle to flow along the SEMFR and melt, forming new oceanic crust. Melts interacted with pre‐existing forearc lithosphere during ascent. The SEMFR is no longer magmatically active and post‐magmatic tectonic activity dominates the rift.     

Parameter and rating curve uncertainty propagation analysis of the SWAT model for two small Mediterranean catchments

Abstract

The SWAT model was tested to simulate the streamflow of two small Mediterranean catchments (the Vène and the Pallas) in southern France. Model calibration and prediction uncertainty were assessed simultaneously by using three different techniques (SUFI-2, GLUE and ParaSol). Initially, a sensitivity analysis was conducted using the LH-OAT method. Subsequent sensitive parameter calibration and SWAT prediction uncertainty were analysed by considering, firstly, deterministic discharge data (assuming no uncertainty in discharge data) and secondly, uncertainty in discharge data through the development of a methodology that accounts explicitly for error in the rating curve (the stage?discharge relationship). To efficiently compare the different uncertainty methods and the effect of the uncertainty of the rating curve on model prediction uncertainty, common criteria were set for the likelihood function, the threshold value and the number of simulations. The results show that model prediction uncertainty is not only case-study specific, but also depends on the selected uncertainty analysis technique. It was also found that the 95% model prediction uncertainty interval is wider and more successful at encompassing the observations when uncertainty in the discharge data is considered explicitly. The latter source of uncertainty adds additional uncertainty to the total model prediction uncertainty.

Editor D. Koutsoyiannis; Associate editor D. Gerten

Citation Sellami, H., La Jeunesse, I., Benabdallah, S., and Vanclooster, M., 2013. Parameter and rating curve uncertainty propagation analysis of the SWAT model for two small Mediterranean watersheds. Hydrological Sciences Journal, 58 (8), 1635?1657.     

Untangling the influence of in-lake productivity and terrestrial organic matter flux on 4,250 years of mercury accumulation in Lake Hambre, Southern Chile

There is ongoing debate about the relative influence of aquatic production, flux, and sedimentation of aquatic and terrestrial organic matter on mercury accumulation in lake sediments. In this study, lake sediments spanning the past 4,250 years, were collected from remote, organic-rich Lake Hambre, Patagonia (53° S) and investigated for changes in the accumulation of pre-anthropogenic mercury and organic matter of aquatic and terrestrial origin. Natural mercury accumulation varied by up to a factor of four, comparable to the recent anthropogenic forcing of the mercury cycle (factor 3–5). Hydrogen and Oxygen indices (HI and OI, Rock–Eval©) and nitrogen/carbon ratios of the organic matter, combined with multi-element sediment data, reveal intense changes in aquatic productivity as well as influx of terrestrial organic matter into the lake. Evaluation of the multi-element dataset using Principal Component Analysis shows clear covariation of mercury with other soil-derived elements such as copper and yttrium. This covariance reflects a common transport mechanism, i.e. leaching of trace-element-bearing organic matter complexes from catchment soils. Correlation between changes in aquatic productivity and mercury concentrations occurs in some sections of the record, but we do not suggest they are linked by a direct causal relationship. Mass balance approaches suggest that mercury scavenging and accumulation in this organic-rich lake is controlled by the supply of mercury from catchment soils rather than the amount of organic material produced within the water column. A common controlling mechanism, i.e. changing climate, however, is thought to independently drive variations in both the flux of terrestrial organic matter mercury complexes and aquatic productivity.     

A new mammalian fauna from the earliest Eocene (Ilerdian) of the Corbières (Southern France): palaeobiogeographical implications

A new mammal fauna from the earliest Eocene of Le Clot (Corbières, Southern France) is described. Some taxa identified there, such as Corbarimys hottingeri and Paschatherium plaziati, allow a correlation with the previously described Corbières fauna of Fordones. Moreover, the presence at Le Clot of Lessnessina praecipuus, which is defined in Palette (Provence, Southern France) allows correlating both localities. All three of these localities are referred to the MP7 reference level, even if a direct correlation with the type locality of MP7 (Dormaal, Belgium) is not ascertained. A Southern Europe biochronological sequence is proposed for the beginning of the Eocene: Silveirinha, Fordones/Palette/Le Clot, Rians/Fournes. The diagnosis of a new species of a neoplagiaulacid multituberculate (?Ectypodus riansensis nov. sp.) is proposed. From a compilation of data concerning the earliest Eocene faunas from both Northern Europe and Southern Europe, we conclude that two faunas corresponding to two distinct palaeobiogeographical provinces existed in Western Europe. This faunal disparity, probably related to climatic differences, makes correlations particularly difficult between the two provinces.     

Morphological similarities between DBM and a microeconomic model of sprawl

We present a model that simulates the growth of a metropolitan area on a 2D lattice. The model is dynamic and based on microeconomics. Households show preferences for nearby open spaces and neighbourhood density. They compete on the land market. They travel along a road network to access the CBD. A planner ensures the connectedness and maintenance of the road network. The spatial pattern of houses, green spaces and road network self-organises, emerging from agents individualistic decisions. We perform several simulations and vary residential preferences. Our results show morphologies and transition phases that are similar to Dieletric Breakdown Models (DBM). Such similarities were observed earlier by other authors, but we show here that it can be deducted from the functioning of the land market and thus explicitly connected to urban economic theory.     

Instability structures, synsedimentary faults and turbidites, witnesses of a Liassic seismotectonic activity in the Dauphiné Zone (French Alps): A case example in the Lower Pliensbachian at Saint-Michel-en-Beaumont

Instability structures, synsedimentary faults and turbidites have been studied in the Lower Pliensbachian succession of Saint-Michel-en-Beaumont, belonging to the Taillefer block, an ancient half-graben emplaced during the Liassic Tethyan rifting. Geometrical and mechanical analyses demonstrate that the instability structures occurred thanks to movements along spineless synsedimentary normal faults, when the turbiditic and limestone layers were already case-hardened and partly fractured by tension gashes even when the mudstones were still unlithified. Both the tension gashes and the synsedimentary faults are homogeneous in strike with the major regional faults and are in good agreement with the regional direction of extension for this period. The characters of the turbiditic beds, with erosive base, graded bedding, and incomplete Bouma sequence, are in favour of a seismic origin. Instability structures, spineless synsedimentary faults and turbiditic inflows are thus considered as seismites and interpreted as the result of high seismicity periods including some events with M > 5 in the general extensive ambiance of the Liassic Tethyan rifting. The analysis of the geometrical relationships between all these sedimentary features allows to distinguish the successive stage of occurrence of an instability structure, from the sedimentation of alternating marls and limestones, and sudden turbiditic inflows, then early case-hardening of the turbidites, until the important seismotectonic event generating the spineless normal faults, themselves triggering the fall of indurated blocks and locally the forming of breccias. The Ornon Fault, which constitutes the border of the Taillefer block, 15 km eastward, played a major role during the Liassic sedimentation and may represent the major seismic fault related to the seismites occurrence in the Beaumont basin.     

Incremental growth of normal faults: Insights from a laser-equipped analog experiment

We conducted a laser-equipped analog experiment aimed at quasi-continuously monitoring the growth of a dense population of normal faults in homogeneous conditions. To further understand the way geological faults progressively gain in slip and length as they accumulate more strain, we measured with great precision the incremental slip and length changes that the analog faults sustain as they grow. These measurements show that the analog faults share common features with the natural ones. In particular, during their growth, the faults develop and maintain cumulative slip profiles that are generally triangular and asymmetric. The growth takes place through two distinct phases: an initial, short period of rapid lateral lengthening, followed by a longer phase of slip accumulation with little or no lengthening. The incremental slip is found to be highly variable in both space (along the faults) and time, resulting in variable slip rates. In particular, ‘short- and long-term’ slip rates are markedly different. We also find that slip measurements at local points on fault traces do not contain clear information on the slip increment repeat mode. Finally, while the fault growth process is highly heterogeneous when considered at the scale of a few slip events, it appears homogeneous and self-similar at longer time scales which integrate many slip increments. This is likely to be the result of a feedback between stress heterogeneities and slip development. The long-term scale homogeneity also implies that the long-term faulting process is primarily insensitive to the short-term heterogeneities that are rapidly smoothed or redistributed. We propose a new conceptual scenario of fault growth that integrates the above observations and we suggest that faults grow in a bimodal way as a result of a self-driven and self-sustaining process.     

Patos Lagoon outflow within the Río de la Plata plume using an airborne salinity mapper: Observing an embedded plume

Major river systems discharging into continental shelf waters frequently form buoyant coastal currents that propagate along the continental shelf in the direction of coastal trapped wave propagation (with the coast on the right/left, in the northern/southern hemisphere). The combined flow of the Uruguay and Paraná Rivers, which discharges freshwater into the Río de la Plata estuary (Lat. ∼36°S), often gives rise to a buoyant coastal current (the ‘Plata plume’) that extends northward along the continental shelf off Uruguay and Southern Brazil. Depending upon the prevailing rainfall, wind and tidal conditions, the Patos/Mirim Lagoon complex (Lat. ∼32°S) may also produce a freshwater outflow plume that expands across the inner continental shelf. Under these circumstances the Patos outflow plume can be embedded in temperature, salinity and current fields that are strongly influenced by the larger Plata plume. The purpose of this paper is to present observations of such an embedded plume structure and to determine the dynamical characteristics of the ambient and embedded plumes.     

Improving modelling and inversion in refraction seismics with a first-order Eikonal solver

A first-order Eikonal solver is applied to modelling and inversion in refraction seismics. The method calculates the traveltime of the fastest wave at any point of a regular grid, including head waves as used in refraction. The efficiency, robustness and flexibility of the method give a very powerful modelling tool to find both traveltimes and raypaths. Comparisons with finite-difference data show the validity of the results. Any arbitrarily complex model can be studied, including the exact topography of the surface, thus avoiding static corrections. Later arrivals are also obtained by applying high-slowness masks over the high-velocity zones. Such an efficient modelling tool may be used interactively to invert for the model, but a better method is to apply the refractor-imaging principle of Hagedoorn to obtain the refractors from the picked traveltime curves. The application of this principle has already been tried successfully by previous authors, but they used a less well-adapted Eikonal solver. Some of their traveltimes were not correct in the presence of strong velocity variations, and the refractor-imaging principle was restricted to receiver lines along a plane surface. With the first-order Eikonal solver chosen, any topography of the receiving surface can be considered and there is no restriction on the velocity contrast. Based on synthetic examples, the Hagedoorn principle appears to be robust even in the case of first arrivals associated with waves diving under the refractor. The velocities below the refractor can also be easily estimated, parallel to the imaging process. In this way, the model can be built up successively layer by layer, the refractor-imaging and velocity-mapping processes being performed for each identified refractor at a time. The inverted model could then be used in tomographic inversions because the calculated traveltimes are very close to the observed traveltimes and the raypaths are available.     

海带(Laminaria japonica)幼孢子体生长和光合作用的N需求

根据室内和围隔实验中海带幼孢子体在不同硝态氮浓度下的生长情况和光合作用速率(Pmax),分析了幼孢子体的N需求,得到其最大生长率(μm)为0.12d-1,维持生存的最低组织N含量(NQ)为16.8μg/mgDW,以最大生长率生长所必需组织N的临界值(NC)为20.4μg/mgDW,每天以最大速率生长的N需求(Nreq)为2.45μg/mgDWd-1。同时,不同处理组的初始NO3-N浓度越高,海带幼孢子体吸收速率和组织N的累加速率越高,且呈明显的线性相关关系(R2分别为0.8393和0.7793,P<0.05)。现场围隔实验中,叶绿素a含量(R2=0.7907,P<0.05)和组织N含量(R2=0.9147,P<0.01)与Pmax也呈明显的线性相关关系。同时,根据海带幼孢子体N的需求和营养吸收状况,分析认为,海带幼孢子体存在受到N限制的风险,但凭其营养吸收能力有适应N限制的能力。还根据海带的这种生理特征,探讨了大型海藻的养殖对富营养化海水的生态调控。