Definition and interests of reciprocity and reciprocity gap principles for groundwater flow problems

We introduce the reciprocity and reciprocity gap principles for flow problems in hydrogeology and illustrate their interest in addressing identification problems. The reciprocity principle is derived from mechanics and establishes for flow problems a relationship between different sets of forcing terms, including sources, sinks and boundary conditions, and the resulting head fields. The reciprocity gap principle compares different head fields resulting from the same forcing terms applied to different structures. We give general 2D expressions of the reciprocity and reciprocity gap principles for transient flow problems and give two examples of applications for the identification of transmissivity values and interfaces between different transmissivities. Identification capacities of the reciprocity and reciprocity gap principles yielding direct inversion methods could be used as initial guesses for more advanced inverse problem methodologies.     

Statistical characteristics of flow as indicators of channeling in heterogeneous porous and fractured media

We introduce two new channeling indicators D_ic and D_cc based on the Lagrangian distribution of flow rates. On the basis of the participation ratio, these indicators characterize the extremes of both the flow-tube width distribution and the flow rate variation along flow lines. The participation ratio is an indicator biased toward the larger values of a distribution and is equal to the normalized ratio of the square of the first-order moment to the second-order moment. Compared with other existing indicators, they advantageously provide additional information on the flow channel geometry, are consistently applicable to both porous and fractured media, and are generally less variable for media generated using the same parameters than other indicators. Based on their computation for a broad range of porous and fracture permeability fields, we show that they consistently characterize two different geometric properties of channels. D_ic gives a characteristic scale of low-flow zones in porous media and a characteristic distance between effectively flowing structures in fractured cases. D_cc gives a characteristic scale of the extension of high-flow zones in porous media and a characteristic channel length in fractured media. D_ic is mostly determined by channel density and permeability variability. D_cc is, however, more affected by the nature of the correlation structure like the presence of permeability channels or fractures in porous media and the length distribution in fracture networks.     

Sky luminance distribution in Pamplona (Spain) during the summer period

In this work the outdoor daylight conditions in Pamplona (South Europe) during the summer period have been studied. The selected sky type (from fifteen standards) at a given moment is the one exhibiting the lowest RMSD when comparing the theoretical and experimental luminance distributions in the sky hemisphere. Two year data of luminance distribution registered every 10 min in 145 positions of the sky hemisphere have been used for selecting the sky type. The most frequent sky type in Pamplona is V.5 (cloudless polluted with a broad solar corona), with an occurrence of 29.5%. This result coincides with the one observed in a previous study in Athens.Six types of sky (V.5, IV.4, III.4, III.3, V.4 y II.2) out of the fifteen standards become practically the 80% of all the studied ones. Regarding a possible use in daylight climate studies, the frequency of occurrence of the fifteen types of sky for fourteen solar elevation intervals has been included.     

The roles of surface heat flux and ocean heat transport convergence in determining Atlantic Ocean temperature variability

The temperature variability of the Atlantic Ocean is investigated using an eddy-permitting (1/4°) global ocean model (ORCA-025) forced with historical surface meteorological fields from 1958 to 2001. The simulation of volume-averaged temperature and the vertical structure of the zonally averaged temperature trends are compared with those from observations. In regions with a high number of observations, in particular above a depth of 500 m and between 22° N and 65° N, the model simulation and the dataset are in good agreement. The relative contribution of variability in ocean heat transport (OHT) convergence and net surface heat flux to changes in ocean heat content is investigated with a focus on three regions: the subpolar and subtropical gyres and the tropics. The surface heat flux plays a relatively minor role in year-to-year changes in the subpolar and subtropical regions, but in the tropical North Atlantic, its role is of similar significance to the ocean heat transport convergence. The strongest signal during the study period is a cooling of the subpolar gyre between 1970 and 1990, which subsequently reversed as the mid-latitude OHT convergence transitioned from an anomalously weak to an anomalously strong state. We also explore whether model OHT anomalies can be linked to surface flux anomalies through a Hovmöller analysis of the Atlantic sector. At low latitudes, increased ocean heat gain coincides with anomalously strong northward transport, whereas at mid-high latitudes, reduced ocean heat loss is associated with anomalously weak heat transport.     

Capturing the residence time boundary layer—application to the Scheldt Estuary

At high Peclet number, the residence time exhibits a boundary layer adjacent to incoming open boundaries. In a Eulerian model, not resolving this boundary layer can generate spurious oscillations that can propagate into the area of interest. However, resolving this boundary layer would require an unacceptably high spatial resolution. Therefore, alternative methods are needed in which no grid refinement is required to capture the key aspects of the physics of the residence time boundary layer. An extended finite element method representation and a boundary layer parameterisation are presented and tested herein. It is also explained how to preserve local consistency in reversed time simulations so as to avoid the generation of spurious residence time extrema. Finally, the boundary layer parameterisation is applied to the computation of the residence time in the Scheldt Estuary (Belgium/The Netherlands). This timescale is simulated by means of a depth-integrated, finite element, unstructured mesh model, with a high space–time resolution. It is seen that the residence time temporal variations are mainly affected by the semi-diurnal tides. However, the spring–neap variability also impacts the residence time, particularly in the sandbank and shallow areas. Seasonal variability is also observed, which is induced by the fluctuations over the year of the upstream flows. In general, the residence time is an increasing function of the distance to the mouth of the estuary. However, smaller-scale fluctuations are also present: they are caused by local bathymetric features and their impact on the hydrodynamics.     

Development of a hybrid simulation controller for full‐scale experimental investigation of seismic retrofits for soft‐story woodframe buildings

Hybrid simulations of a full‐scale soft‐story woodframe building specimen with various retrofits were carried out as part of the Network for Earthquake Engineering Simulation Research project – NEES‐Soft: seismic risk reduction for soft‐story woodframe buildings. The test structure in the hybrid simulation was a three‐story woodframe building that was divided into a numerical substructure of the first story with various retrofits and a full‐scale physical substructure of the upper two stories. Four long‐stroke actuators, two at the second floor and two at the roof diaphragm, were attached to the physical substructure to impose the simulated seismic responses including both translation and in‐plane rotation. Challenges associated with this first implementation of a full‐scale hybrid simulation on a woodframe building were identified. This paper presents the development and validation of a scalable and robust hybrid simulation controller for efficient test site deployment. The development consisted of three incremental validation phases ranging from small‐scale, mid‐scale, to full‐scale tests conducted at three laboratories. Experimental setup, procedure, and results of each phase of the controller development are discussed, demonstrating the effectiveness and efficiency of the incremental controller development approach for large‐scale hybrid simulation programs with complex test setup. Copyright © 2016 John Wiley & Sons, Ltd.     

Impact of draining hilly lands on runoff and on‐site erosion: a case study from humid Ethiopia

The use of drainage ditches on farmland has an impact on erosion processes both on‐site and off‐site, though their environmental impacts are not unequivocal. Here we study the runoff response and related rill erosion after installing drainage ditches and assess the effects of stone bunds in north Ethiopia. Three different land management systems were studied in 10 cropland catchments around Wanzaye during the rainy season of 2013: (1) the exclusive use of drainage ditches (locally called feses), (2) the exclusive use of stone bunds, and (3) a mixture of both systems. Stone bunds are an effective soil and water conservation technique, making the land more resistant against on‐site erosion, and allowing feses to be installed at a larger angle with the contour. The mean rill volumes for the 10 studied cropland catchments during the rainy season of 2013 was 3.73 ± 4.20 m³ ha^?1 corresponding to a soil loss of 5.72 ± 6.30 ton ha^?1. The establishment of feses causes larger rill volumes (R = 0.59, N = 10), although feses are perceived as the best way to avoid soil erosion when no stone bunds are present. The use of feses increases event‐based runoff coefficients (RCs) on cropland from c. 5% to values up to 39%. Also, a combination of low stone bund density and high feses density results in a higher RC, whereas catchments with a high stone bund density and low feses density have a lower RC. Peak runoff discharges decrease when stone bund density increases, whereas feses density is positively related to the peak runoff discharge. A multiple linear relation in which both feses and stone bund densities are used as explanatory variable, performs best in explaining runoff hydrograph peakedness (R² = 83%). Copyright © 2015 John Wiley & Sons, Ltd.     

Root properties of vegetation communities and their impact on the erosion resistance of river dikes

Predicted climate change and the associated sea level rise poses an increased threat of flooding due to wave overtopping events at sea and river dikes. To safeguard the land from flooding it is important to keep the soil erosion resistance at the dikes high. As plant roots can be very effective in reducing soil erosion rates by concentrated flow, the main goal of this study is to explore the variability in root system characteristics of five dike vegetation communities along the Scheldt River (Belgium) and to assess their effectiveness in controlling soil erosion rates during concentrated flow. This study is the first one to investigate systematically the erosion‐reducing potential of the root properties of representative dike vegetation communities in a temperate humid climate. Results show that the presence of Urtica dioica resulted in large differences in root length density (RLD) among dike vegetation communities. Observed RLD values in the topsoil ranged from 129 to 235 km m^‐3 for dike vegetation communities without U. dioica, while smaller values ranging from 22 to 58 km m^?3 were found for vegetation communities with U. dioica. The erosion‐reducing effect of the dike vegetation communities was estimated based on a global Hill curve model, linking the RLD to the soil detachment ratio (SDR; i.e. the ratio of the soil detachment rate for root‐permeated topsoils to the soil detachment rate for root‐free topsoils). Concentrated flow erosion rates are likely to be reduced to 13–16% of the erosion rates for root‐free topsoils if U. dioica is absent compared to 22–30% for vegetation communities with U. dioica. Hence, to maintain a high resistance of the soil against concentrated flow erosion it is important to avoid the overgrowth of grassland by U. dioica through an effective vegetation management. Copyright © 2016 John Wiley & Sons, Ltd.     

The Burren: a glacial,karstic and biokarstic expression of a limestone plateau in western Ireland

The term glaciokarst describes a landscape where both glacial and karstic processes have contributed to geomorphological evolution and has been applied to a range of environments from the high arctic to the alpine Mediterranean. Nevertheless, glaciokarstic environments globally often exhibit significant variation in geomorphological processes and landforms due to these geographical differences. The Burren, County Clare, Republic of Ireland, is often quoted as a quintessential glaciokarstic landscape. However, the Burren and other similar environments would appear to lean towards one end of the glaciokarst spectrum, where solutional and biological processes have been dominant throughout the Holocene, in contrast to ice‐dominated glaciokarsts where karstic and biokarstic processes are temperature‐limited and cryospheric processes remain the principal geomorphological agents. Holocene landscape evolution and the development of a range of meso‐, micro‐ and nano‐scale karren features on limestone surfaces appears to be largely biokarstic in origin. Karstification of many glacially scoured limestone pavements would have begun under acidic soil cover, with biological soil processes contributing to smooth, rounded cryptokarstic surface forms. Holocene soil erosion is attributed to anthropogenic activity, climatic fluctuations and the evolution of the karstic groundwater system leading to vertical soil loss through widening grikes. Exposed limestone pavements subject to subaerial conditions often exhibit extensive lichen colonization which has been shown to influence the overall rate of karstification and contribute to the development of micro‐ and meso‐scale bioweathering features. Where cryptokarstic features have been exhumed from beneath soil cover, their evolution under subaerial conditions leads to intermediate, polygenetic karren features. In light of our current understanding of the Burren landscape, it is proposed that the term glaciobiokarstic may be a better expression to encompass the biological processes that have played a fundamental role in the evolution of the Burren and similar landscapes, without neglecting the contribution of glacial and inorganic karstic processes. Copyright © 2016 John Wiley & Sons, Ltd.