Flow dynamics at the continental scale: Streamflow correlation and hydrological similarity

Streamflow variability in space and time critically affects anthropic water uses and ecosystem services. Unfortunately, spatiotemporal patterns of flow regimes are often unknown, as discharge measurements are usually recorded at a limited number of hydrometric stations unevenly distributed along river networks. Advances in understanding the physical processes that control the spatial patterns of river flows are therefore necessary to predict water availability at ungauged locations or to extrapolate pointwise streamflow observations. This work explores the use of the spatial correlation of river flows as a metric to quantify the similarity between hydrological responses of two catchments. Following a stochastic framework, 340,000 cross‐correlations between pairs of daily streamflows time series are predicted at a seasonal timescale across the contiguous United States using 413 catchments of the MOPEX dataset. Model predictions of streamflow correlation obtained in absence of run‐off information are successfully used to identify catchment outlets sharing similar discharge dynamics and flow regimes across a broad range of geomorphoclimatic conditions, without relying on calibration. The selection of reference streamgauges based on predicted streamflow correlation generally outperforms the selection based on spatial proximity, especially as the density of available gauged sections decreases. Interestingly, correlated outlets share a broad spectrum of hydrological signatures (mean discharge, flow variability, and recession properties), suggesting that catchments forced by analogous frequency and intensity of effective rainfall events might exhibit common geomorphoecological traits leading to similar hydrological responses. The proposed framework provides a physical basis to assist the regionalization of flow dynamics and to interpret the spatial variability of flow regimes along stream networks.     

Proposal of design rules for ductile X-CBFS in the framework of EUROCODE 8

Cross concentrically braced frames (X-CBFs) are commonly used as primary seismic resisting system, owing to their large lateral stiffness, simplicity of design, and relatively low constructional cost. Current EN 1998-1 provides design rules theoretically aiming at developing ductile global plastic mechanism, namely enforcing plastic deformations in the diagonal members, while the remaining structural members and connections should elastically behave. However, as widely demonstrated by many existing studies, the design and the corresponding seismic performance of EC8-compliant X-CBFs are generally affected by several criticisms, eg, difficulties in sizing of diagonal members, massive and non-economical structures, and poor seismic behavior. In light of these considerations, the research activity presented in this paper is addressed to revise the design rules and requirements given EN 1998-1 for X-CBFs to simplify the design process and to improve the ductility and the dissipative capacity of this structural system. Hence, design rules are proposed for the next version of EN 1998-1 and numerically validated by means of nonlinear dynamic analyses.     

A reduced-order model for Monte Carlo simulations of stochastic groundwater flow

We explore the ability of the greedy algorithm to serve as an effective tool for the construction of reduced-order models for the solution of fully saturated groundwater flow in the presence of randomly distributed transmissivities. The use of a reduced model is particularly appealing in the context of numerical Monte Carlo (MC) simulations that are typically performed, e.g., within environmental risk assessment protocols. In this context, model order reduction techniques enable one to construct a surrogate model to reduce the computational burden associated with the solution of the partial differential equation governing the evolution of the system. These techniques approximate the model solution with a linear combination of spatially distributed basis functions calculated from a small set of full model simulations. The number and the spatial behavior of these basis functions determine the computational efficiency of the reduced model and the accuracy of the approximated solution. The greedy algorithm provides a deterministic procedure to select the basis functions and build the reduced-order model. Starting from a single basis function, the algorithm enriches the set of basis functions until the largest error between the full and the reduced model solutions is lower than a predefined tolerance. The comparison between the standard MC and the reduced-order approach is performed through a two-dimensional steady-state groundwater flow scenario in the presence of a uniform (in the mean) hydraulic head gradient. The natural logarithm of the aquifer transmissivity is modeled as a second-order stationary Gaussian random field. The accuracy of the reduced basis model is assessed as a function of the correlation scale and variance of the log-transmissivity. We explore the performance of the reduced model in terms of the number of iterations of the greedy algorithm and selected metrics quantifying the discrepancy between the sample distributions of hydraulic heads computed with the full and the reduced model. Our results show that the reduced model is accurate and is highly efficient in the presence of a small variance and/or a large correlation length of the log-transmissivity field. The flow scenarios associated with large variances and small correlation lengths require an increased number of basis functions to accurately describe the collection of the MC solutions, thus reducing significantly the computational advantages associated with the reduced model.     

The scattering of SH waves by a finite crack with a superposition-based diffraction technique

The problem of diffraction of cylindrical and plane horizontally polarized shear waves (SH waves) by a finite crack embedded in a plane bidimensional elastic full-space is revisited. Particularly, we construct an approximate solution by the addition of independent diffracted terms. In our method the derivation of the fundamental case of a semi-infinite crack obtained as a degenerate case of a generalized wedge is first considered. This result is then used as a building block to compute the diffraction of the main incident waves. The interaction between the opposite edges of the crack is later considered in terms of a series, one term at a time until a desired tolerance is reached. Moreover, we propose a procedure to determine the number of required interactions as a function of frequency. The solution derived with the superposition technique is shown to be effective at low and high frequencies and as shown by comparisons with a direct boundary element method software, highly accurate solutions are obtained after retaining just a few terms of the infinite series.     

On the estimation of the maximum depth of investigation of transient electromagnetic soundings: the case of the Vizcaino transect,Mexico

We test an approach proposed in the literature for estimating the maximum depth of investigation (MDI) of in-loop transient electromagnetic soundings (TEM) with data from the Vizcaino transect, a profile of 38 TEM soundings crossing the Baja California peninsula. The validity of this approach is confirmed by comparing the MDI with the deepest interface of 16 stratified models. In these soundings, all located over the Vizcaino basin, the measured data are not affected by induced polarization. The models indicate the presence of a dipping conductor interpreted as a zone of seawater intrusion with a large lateral extension of over 70 km. The remaining 22 soundings, located over igneous and metamorphic rocks, show reversals in the voltage polarity, indicating the presence of induced polarization effects. The layered Cole-Cole models for these soundings suggest a significant decrease in the MDI. This is confirmed by analyzing the depth behavior of the subsurface current densities. We further analyze the noise level of a data set comprising close to 2000 late-time voltages of about 400 TEM sites acquired in northwestern Mexico. No difference was found between the stationary noise levels of winter and summer, presumably because near thunderstorms are practically absent in this part of Mexico.     

6th April 2009 L��Aquila earthquake, Italy: reinforced concrete building performance

On 6th April 2009 an earthquake of magnitude M _w =  6.3 occurred in the Abruzzo region; the epicentre was very close to the city of L’Aquila (about 6 km away). The event produced casualties and damage to buildings, lifelines and other infrastructures. An analysis of the main damage that reinforced concrete (RC) structures showed after the event is presented in this study. In order to isolate the main causes of structural and non-structural damage, the seismological characteristics of the event are examined, followed by an analysis of the existing RC building stock in the area. The latter issue came under scrutiny after the release of official data about structural types and times of construction, combined with a detailed review of the most important seismic codes in force in the last 100 years in Italy. Comparison of the current design provisions of the Italian and European codes with previous standards allows the main weaknesses of the existing building stock to be determined. Damage to structural and non-structural elements is finally analyzed thanks to photographic material collected in the first week after the event; the main causes of damage are then inferred.     

Intra‐event scale bar–bank interactions and their role in channel widening

Channel bars and banks strongly affect the morphology of both braided and meandering rivers. Accordingly, bar formation and bank erosion processes have been greatly explored. There is, however, a lack of investigations addressing the interactions between bed and bank morphodynamics, especially over short timescales. One major implication of this gap is that the processes leading to the repeated accretion of mid‐channel bars and associated widenings remain unsolved. In a restored section of the Drau River, a gravel‐bed river in Austria, mid‐channel bars have developed in a widening channel. During mean flow conditions, the bars divert the flow towards the banks. One channel section exhibited both an actively retreating bank and an expanding mid‐channel bar, and was selected to investigate the morphodynamic processes involved in bar accretion and channel widening at the intra‐event timescale. We repeatedly surveyed riverbed and riverbank topography, monitored riverbank hydrology and mounted a time‐lapse camera for continuous observation of riverbank erosion processes during four flow events. The mid‐channel bar was shown to accrete when it was submerged during flood events, which at the subsequent flow diversion during lower discharges narrowed the branch along the bank and increased the water surface elevation upstream from the riffle, which constituted the inlet into the branch. These changes of bed topography accelerated the flow along the bank and triggered bank failures up to 20 days after the flood events. Four analysed flow events exhibited a total bar expansion from initially 126 m² to 295 m², while bank retreat was 6 m at the apex of the branch. The results revealed the forcing role of bar accretion in channel widening and highlighted the importance of intra‐event scale bed morphodynamics for bank erosion, which were summarized in a conceptual model of the observed bar–bank interactions. Copyright © 2015 John Wiley & Sons, Ltd.     

Indentation of a continent with a built-in thickness change: experiment and nature

Orogens oblique to the direction of plate convergence are currently attributed to obliquity between the margins of one or both of the sutured continents to their direction convergence. We use a single analogue experiment and natural examples to illustrate a potential additional factor: variations in strength of the indented continent at a high angle to the convergence direction. The wavelengths of structures in laterally shortened lithosphere depend on the strength of the most competent layers. Lateral variations in crustal thickness must therefore lead to structures oblique to any applied lateral compression.

An analogue experiment was performed to explore this phenomenon. A two-layer ‘indented continent’ was modelled by a brittle upper crust of sand above a lower crust of high-viscosity polymer floating on a single layer of low-viscosity syrup representing the mantle. The well-known strike-slip structures allowing lateral escape to distant weak boundaries were hindered by lateral boundaries in front of the indenter. This allowed us to focus on the effects of a thickness change built into the ‘indented continent’ along a zone parallel to the direction in which a vertical rigid wall advancing at a steady rate represented the indenter. Vertical escape led to an ‘orogenic belt’ oblique to the advancing wall; this obliquity influences subsequent lateral escape. Model scaling and interpretations are based on Extended Thin Sheet Approximation (ETSA) and standard theories of faulting.

Four sectors of the Alpine–Himalayan orogen (Iran, Tunisia, the Eastern Alps and the Himalaya) are oblique to the continental convergence direction, and we point to thickness changes at high angles to the suture that may account for this geometry. As crustal thicknesses north of oblique sectors of the Himalayas are not yet known, we speculate on them.

We infer from the main difference between our experiment and all our examples chosen from nature that vertical orogenic escape was oblique to our model suture but can be parallel to natural sutures.