Coupling mechanisms in double sandbar systems. Part 2: impact on alongshore variability of inner‐bar rip channels

Double sandbar systems are common morphological features along sandy, wave‐dominated, micro‐ to meso‐tidal coastlines. In the companion paper, we demonstrated how various alongshore inner‐bar rip‐channel patterns can develop through morphological coupling to an alongshore‐variable outer bar. The simulated coupling patterns are, however, scarcely observed in the field. Instead, inner‐bar rip channels more often possess remarkably smaller and more variable alongshore length scales, suggesting that coupling mechanisms do not play a substantial role in the overall double‐sandbar dynamics. Here we use a numerical model to show that the relative importance of self‐organization and morphological coupling changes in favour of the latter with an increase in waterdepth variability along the outer‐bar crest. Furthermore, we find that the typical alongshore variability in inner‐bar rip‐channel scale is indicative of a mixture of self‐organization and morphological coupling rather than self‐organization alone. Morphological coupling may thus be more important to understanding and predicting the evolution of inner‐bar rip channels than previously envisaged. Copyright © 2010 John Wiley and Sons, Ltd.     

Equivalent linear substructure approximation of soil–foundation–structure interaction: model presentation and validation

An equivalent linear substructure approximation of the soil–foundation–structure interaction is proposed in this paper. Based on the inherent linearity of the approach, the solution of the structural and the soil domain is obtained simultaneously, incorporating the effects of the primary and secondary soil nonlinearities. The proposed approximation is established theoretically and then validated against centrifuge benchmark soil–foundation–structure interaction tests. The equivalent linear substructure approximation is proved to simulate efficiently the effects of the nonlinear soil behavior on the soil–foundation–structure system under a strong earthquake ground motion.     

INTERNATIONAL HYDROLOGICAL DECADE / WORKING GROUP ON “THE INFLUENCE OF MAN ON THE HYDROLOGICAL CYCLE”

Abstract

Floods, as extreme hydrological phenomena, can be described by more than one correlated characteristic, such as peak, volume and duration. These characteristics should be jointly considered since they are generally not independent. For an ungauged site, univariate regional flood frequency analysis (FA) provides a limited assessment of flood events. A recent study proposed a procedure for regional FA in a multivariate framework. This procedure represents a multivariate version of the index-flood model and is based on copulas and multivariate quantiles. The performance of the proposed procedure was evaluated by simulation. However, the model was not tested on a real-world case study data. In the present paper, practical aspects are investigated jointly for flood peak (Q) and volume (V) of a dataset from the Côte-Nord region in the province of Quebec, Canada. The application of the proposed procedure requires the identification of the appropriate marginal distribution, the estimation of the index flood and the selection of an appropriate copula. The results of the case study show that the regional bivariate FA procedure performed well. This performance depends strongly on the performance of the two univariate models and, more specifically, the univariate model of Q. The results show also the impact of the homogeneity of the region on the performance of the univariate and bivariate models.

Editor D. Koutsoyiannis     

Distinguishing the Regional Atmospheric Controls on Precipitation Isotopic Variability in the Central-Southeast Portion of Brazil

Precipitation isotope ratios (O and H) record the history of water phase transitions and fractionation processes during moisture transport and rainfall formation. Here, we evaluated the isotopic composition of precipitation over the central-southeastern region of Brazil at different timescales. Monthly isotopic compositions were associated with classical effects (rainfall amount, seasonality, and continentality), demonstrating the importance of vapor recirculation processes and different regional atmospheric systems (South American Convergence Zone-SACZ and Cold Fronts-CF). While moisture recycling and regional atmospheric processes may also be observed on a daily timescale, classical effects such as the amount effect were not strongly correlated (δ¹⁸O-precipitation rate r ≤ –0.37). Daily variability revealed specific climatic features, such as δ¹⁸O depleted values (~ –6‰ to –8‰) during the wet season were associated with strong convective activity and large moisture availability. Daily isotopic analysis revealed the role of different moisture sources and transport effects. Isotope ratios combined with d-excess explain how atmospheric recirculation processes interact with convective activity during rainfall formation processes. Our findings provide a new understanding of rainfall sampling timescales and highlight the importance of water isotopes to decipher key hydrometeorological processes in a complex spatial and temporal context in central-southeastern Brazil.     

Well placement optimization with the covariance matrix adaptation evolution strategy and meta-models

The amount of hydrocarbon recovered can be considerably increased by finding optimal placement of non-conventional wells. For that purpose, the use of optimization algorithms, where the objective function is evaluated using a reservoir simulator, is needed. Furthermore, for complex reservoir geologies with high heterogeneities, the optimization problem requires algorithms able to cope with the non-regularity of the objective function. In this paper, we propose an optimization methodology for determining optimal well locations and trajectories based on the covariance matrix adaptation evolution strategy (CMA-ES) which is recognized as one of the most powerful derivative-free optimizers for continuous optimization. In addition, to improve the optimization procedure, two new techniques are proposed: (a) adaptive penalization with rejection in order to handle well placement constraints and (b) incorporation of a meta-model, based on locally weighted regression, into CMA-ES, using an approximate stochastic ranking procedure, in order to reduce the number of reservoir simulations required to evaluate the objective function. The approach is applied to the PUNQ-S3 case and compared with a genetic algorithm (GA) incorporating the Genocop III technique for handling constraints. To allow a fair comparison, both algorithms are used without parameter tuning on the problem, and standard settings are used for the GA and default settings for CMA-ES. It is shown that our new approach outperforms the genetic algorithm: It leads in general to both a higher net present value and a significant reduction in the number of reservoir simulations needed to reach a good well configuration. Moreover, coupling CMA-ES with a meta-model leads to further improvement, which was around 20% for the synthetic case in this study.     

In Situ Determination of Au and Cu in Natural Pyrite by Near‐Infrared Femtosecond Laser Ablation‐Inductively Coupled Plasma‐Quadrupole Mass Spectrometry: No Evidence for Matrix Effects

Gold and copper concentrations were determined in natural pyrite by near‐infrared femtosecond LA‐ICP‐QMS, using both sulfide reference materials (pyrrhotite Po‐726 and in‐house natural chalcopyrite Cpy‐RM) and NIST SRM 610 as external calibrators. Firstly, using NIST SRM 610 as the external calibrator, we calculated the Au concentration in Po‐726 and the Cu concentration in Cpy‐RM. The calculated concentration averages for Au and Cu were similar to the values published for Po‐726 and Cpy‐RM, respectively. Secondly, we calculated Au and Cu concentrations taking NIST SRM 610 as an unknown sample and using Po‐726 and Cpy‐RM as external calibrators. Again, the average values obtained closely reflected the preferred concentrations for NIST SRM 610. Finally, we calculated Au and Cu concentrations in natural pyrite using sulfide and silicate reference materials as external calibrators. In both cases, calculated concentrations were very similar, independent of the external calibrator used. The aforementioned data, plus the fact that we obtained very small differences in relative sensitivity values (percentage differences are between 5% and 17% for ⁵⁷Fe, ⁶³Cu and ¹⁹⁷Au) on analyses of silicate and sulfide RMs, indicate that there were no matrix effects related to the differences in material composition. Thus, it is possible to determine Au and Cu in natural sulfides using NIST silicate glasses as an external calibrator.     

Regional frequency analysis of autumnal floods in the province of Quebec, Canada

A large number of models have been proposed over the last years for regional flood frequency analysis in northern regions. However, these models dealt generally with snowmelt-caused spring floods. This paper deals with the adaptation, application, and comparison of two regional frequency analysis methods, canonical correlation analysis (CCA) and universal canonical kriging (UCK), on autumnal floods of 29 stations from the C?te-Nord region (QC, Canada). Three possible periods during which autumnal floods can take place are tested. The absolute and specific flood peak and volume quantiles are also studied. A jack-knife resampling procedure is applied to compare the performance of each model according to the selected period and the type of quantile. The period of September 1st to December 15th is found to be optimal to represent autumnal floods and specific quantiles were shown to lead to better results than absolute quantiles. Variables that explain best the autumnal floods are the basin area, the fraction of the area covered with lakes, and the average of mean July, August, and September maximal temperatures. The CCA model performs slightly better than UCK.     

HP–UHP metamorphism as an indicator of slab dip variations in the Alpine arc

HP/UHP and LT metamorphic units that commonly occur in the inner parts of mountain belts result from the subduction of continental and oceanic material, most often exhumed prior to continental collision. The prograde pressure–temperature history of HP–UHP rocks strongly depends on the convergence rate and on the subduction zone geometry. The maximum pressure recorded provides a proxy for the depth of shearing off and stacking of HP metamorphic nappes. A 2-D thermal model of continental subduction at lithospheric scale is used to compute the length and pressure peak of detached HP metamorphic units as a function of the slab dip angle and the convergence rate. Model results are applied to the metamorphic nappe pile of the inner Alps. A mean convergence rate of 1 cm/year during the subduction of the Briançonnais terrane is indicated by the paleogeographic reconstructions between 46 and 38 Ma. On this basis, the available petrological data and lengths of metamorphic units are used to compute the variations of the slab dip angle. The slab dip angle is shown to increase, from the northeast to the southwest, along the Alpine arc with estimated values of 20° for Suretta, 30–45° for Monte Rosa and Gran Paradiso, and 60° for Dora Maira. From Eocene to Oligocene times, the increase in slab dip angle is controlled by changes of buoyancy, due to the spatial configuration of the Valaisan trough and the incoming of crustal material within the subduction zone.     

Anisotropic Closure in Squeezing Rocks: The Example of Saint-Martin-la-Porte Access Gallery

In this paper, we analyse the convergence measurements recorded in a gallery excavated in severely squeezing ground. The procedure consists in a preliminary geometrical treatment of the raw data to evaluate the principal axes of deformation by assuming an elliptic deformation of the walls of the gallery. Then the convergence law proposed by Sulem et al. (Int J Rock Mech Min Sci Geomech Abstr 24(3):145–154, 1987a), which is extended to account for anisotropic closure, is fitted on the displacements along the two axes of the obtained ellipse. This procedure is more robust and relevant than fitting the convergence recorded on the most deforming segment or fitting the average value of the convergence along the various segments. An attempt is made to correlate the amount and the direction of anisotropic deformation with the lithology and some geological features described by the dominant discontinuity families.