ADM1‐Based Robust Interval Observer for Anaerobic Digestion Processes

A robust state estimation scheme is proposed for anaerobic digestion (AD) processes to estimate key variables under the most uncertain scenarios (namely, uncertainties on the process inputs and unknown reaction and specific growth rates). This scheme combines the use of the IWA Anaerobic Digestion Model No. 1 (ADM1), the interval observer theory and a minimum number of measurements to reconstruct the unmeasured process variables within guaranteed lower and upper bounds in which they evolve. The performance of this robust estimation scheme is evaluated via numerical simulations that are carried out under actual operating conditions. It is shown that under some structural and operational conditions, the proposed robust interval observer (RIO) has the property of remaining stable in the face of uncertain process inputs, badly known kinetics and load disturbances. It is also shown that the RIO is indeed a powerful tool for the estimation of biomass (composed of seven different species) from a minimum number of measurements in a system with a total of 32 variables from which 24 correspond to state variables.     

Discrete principal‐monotonicity inference for hydro‐system analysis under irregular nonlinearities,data uncertainties,and multivariate dependencies. Part I: methodology development

Hydrological system analyses are challenged by complexities of irregular nonlinearities, data uncertainties, and multivariate dependencies. Among them, the irregular nonlinearities mainly represent inexistence of regular functions for robustly simulating highly complicated relationships between variables. Few existing studies can enable reliable simulation of hydrological processes under these complexities. This may lead to decreased robustness of the constructed models, unfeasibility of suggestions for human activities, and damages to socio‐economy and eco‐environment. In the first of two companion papers, a discrete principal‐monotonicity inference (DPMI) method is proposed for hydrological systems analysis under these complexities. Normalization of non‐normally distributed samples and invertible restoration of modelling results are enabled through a discrete distribution transformation approach. To mitigate data uncertainties, statistical inference is employed to assess the significance of differences among samples. The irregular nonlinearity between the influencing factors (i.e. predictors) and the hydrological variable of interest (i.e. the predictand) is interpreted as piecewise monotonicity. Monotonicity is further represented as principal monotonicity under multivariate dependencies. Based on stepwise classification and cluster analyses, all paired samples representing the responsive relationship between the predictors and the predictand are discretized as a series of end nodes. A prediction approach is advanced for estimating the predictand value given any combination of predictors. The DPMI method can reveal evolvement rules of hydrological systems under these complexities. Reliance of existing hydro‐system analysis methods on predefined functional forms is removed, avoiding artificial disturbances, e.g. empiricism in selecting model functions under irregular nonlinearities, on the modelling process. Both local and global significances of predictors in driving the evolution of hydrological variables are identified. An analysis of interactions among these complexities is also achieved. The understanding obtained from the DPMI process and associated results can facilitate hydrological prediction, guide water resources management, improve hydro‐system analysis methods, or support hydrological systems analysis in other cases. The effectiveness and advantages of DPMI will be demonstrated through a case study of streamflow simulation in Xingshan Watershed, China, in another paper. Copyright © 2016 John Wiley & Sons, Ltd.     

Treatment of a Trifluraline Effluent by Means of Oxidation‐Coagulation with Fe(VI) and Combined Fenton Processes

The amination water (AW) effluent stream from the industrial production of the trifluraline herbicide was submitted to an oxidation‐coagulation treatment with potassium ferrate, combined with advanced oxidation processes. The experimental results obtained by analysis of variance (ANOVA) for the oxidation‐coagulation‐Fenton process, evaluating the variables pH (A), Fe(VI) concentration (B), and H₂O₂ concentration (C), demonstrated that the regression equation resulting from the Response Surface Methodology (RSM) experimental design, for the quadratic model, was η_Abs (%) = 36.9– 21.58A + 8.37A² + 1.36B + 0.92B² + 1.08C + 1.52C² + 1.27AB – 1.34AC + 1.33BC. The maximum absorptiometric color reduction occurred at pH 3, with corresponding maximum amounts of iron and hydrogen peroxide. The absorptiometric color and COD reduction were 96% and 57%, respectively. For the oxidation‐coagulation‐photo‐Fenton process, the analyzed variables were pH (A), Fe(VI) concentration (B), H₂O₂ concentration (C), and temperature (D). The regression equation resulting from the quadratic model was η_Abs (%) = 38.3 – 20.2A + 8.12A² – 0.27B + 3.73B² + 0.3C + 3.6C² + 1.67D + 3.1D² + 1.72AB + 0.51AC – 1.82AD + 0.74BC – 1.11BD + 0.03CD. The ANOVA response showed that the highest absorptiometric color reduction occurred at pH 3, with respective maximum amounts of iron and hydrogen peroxide at 60°C. The maximum efficiencies achieved by the proposed treatment process for the trifluraline effluent stream were 95% and 85%, for absorptiometric color and COD reduction, respectively.     

Robust Regulation of Alkalinity in Highly Uncertain Continuous Anaerobic Digestion Processes

In this paper, a multiple‐input multiple‐output (MIMO) model‐based robust control scheme is proposed for the indirect control of both total alkalinity and the ratio (intermediate alkalinity)/(total alkalinity) by regulating volatile fatty acid concentrations and strong ions concentration, while guaranteeing the so‐called operational stability in anaerobic digestion (AD) processes. The proposed MIMO regulator is an adaptive controller derived from an AD model which incorporates the physicochemical equilibrium of the system as well as the use of a robust interval observer to estimate key process bounds that are used in the computation of the control efforts. Numerical simulations were carried out for a number of operating conditions under the most uncertain scenarios. Results showed that the proposed multivariable control law is able to recover the system stability around a pre‐determined set point in the face of parameter uncertainty and load disturbances.     

Modelling microseismicity of a producing reservoir from coupled fluid‐flow and geomechanical simulation

In this paper, we investigate production induced microseismicity based on modelling material failure from coupled fluid‐flow and geomechanical simulation. The model is a graben style reservoir characterized by two normal faults subdividing a sandstone reservoir into three compartments. The results are analysed in terms of spatial and temporal variations in distribution of material failure. We observe that material failure and hence potentially microseismicity is sensitive to not only fault movement but also fluid movement across faults. For sealing faults, failure is confined to the volume in and around the well compartment, with shear failure localized along the boundaries of the compartment and shear‐enhanced compaction failure widespread throughout the reservoir compartment. For non‐sealing faults, failure is observed within and surrounding all three reservoir compartments as well as a significant distribution located near the surface of the overburden. All shear‐enhanced compaction failures are localized within the reservoir compartments. Fault movement leads to an increase in shear‐enhanced compaction events within the reservoir as well as shear events located within the side‐burden adjacent to the fault. We also evaluate the associated moment tensor mechanisms to estimate the pseudo scalar seismic moment of failure based on the assumption that failure is not aseismic. The shear‐enhanced compaction events display a relatively normal and tight pseudo scalar seismic moment distribution centred about 10⁶ Pa, whereas the shear events have pseudo scalar seismic moments that vary over three orders of magnitude. Overall, the results from the study indicate that it may be possible to identify compartment boundaries based on the results of microseismic monitoring.     

Basement structural architecture and hydrocarbon conduit potential of polygonal faults in the Cooper‐Eromanga Basin,Australia

A thorough and complete understanding of the structural geology and evolution of the Cooper‐Eromanga Basin has been hampered by low‐resolution seismic data that becomes particularly difficult to interpret below the thick Permian coal measures. As a result, researchers are tentative to interpret the basement fault architecture within the basin, which is largely undefined. To provide a better understanding of the basement fault geometry, all available two‐dimensional seismic lines together with 12 three‐dimensional seismic surveys were structurally interpreted with assistance from seismic attribute analysis. The Upper Cretaceous Cadna‐owie Formation and top Permian reflectors were analysed using a common seismic attribute technique (incoherency) that was used to infer the presence of faults that may have otherwise been overlooked. Detailed basement fault maps for each seismic survey were constructed and used in conjunction with two‐dimensional seismic data interpretation to produce a regional basement fault map. Large north‐northeast–south‐southwest‐striking sinistral strike–slip faults were identified within the Patchawarra Trough appearing to splay from the main northeast–southwest‐striking ridge. These sinistral north‐northeast–south‐southwest‐striking faults, together with field‐scale southeast–northwest‐striking dextral strike–slip faults, are optimally oriented to have potentially developed as a conjugated fault set under a south‐southeast–north‐northwest‐oriented strike–slip stress regime. Geomechanical modelling for a regionally extensive system of Cretaceous polygonal faults was performed to calculate the Leakage Factor and Dilation Tendency of individual faults. Faults that extend into Lower Cretaceous oil‐rich reservoirs with strikes of between 060°N and 140°N and a high to near‐vertical dip angle were identified to most likely be acting as conduits for the tertiary migration of hydrocarbons from known Lower Cretaceous hydrocarbon reservoirs into shallow Cretaceous sediments. This research provides valuable information on the regional basement fault architecture and a more detailed exploration target for the Cooper‐Eromanga Basin, which were previously not available in literature.     

Quantifying different riverbank erosion processes during an extreme flood event

Riverbank erosion is a major contributor to catchment sediment budgets. At large spatial scales data is often restricted to planform channel change, with little information on process distributions and their sediment contribution. This study demonstrates how multi‐temporal LiDAR and high resolution aerial imagery can be used to determine processes and volumes of riverbank erosion at a catchment scale. Remotely sensed data captured before and after an extreme flood event, enabled a digital elevation model of difference (DoD) to be constructed for the channel and floodplain. This meant that: the spatial area that could be assessed was extensive; three‐dimensional forms of bank failures could be mapped at a resolution that enabled process inference; and the volume and rates of different bank erosion processes over time could be assessed. A classification of riverbank mass failures, integrating form and process, identified a total of 437 mass failure polygons throughout the study area. These were interpreted as wet flow mass failures based on the presence of a well defined scarp wall and the absence of failed blocks on the failure floor. The failures appeared to be the result of: bank exfiltration, antecedent moisture conditions preceding the event, and the historic development of the channel. Using one‐dimensional hydraulic modelling to delineate geomorphic features within the main boundary of the macrochannel, an estimated 1 466 322 m² of erosion was interpreted as fluvial entrainment, occurring across catchment areas from 30 to 1668 km². Only 8% of the whole riverbank planform area was occupied by mass failures, whilst fluvial entrainment covered 33%. A third of the volume of material eroded came from mass failures, even though they occupied 19% of the eroded bank area. The availability of repeat LiDAR surveys, combined with high‐resolution aerial photography, was very effective in erosion process determination and quantification at a large spatial scale. Copyright © 2013 John Wiley & Sons, Ltd.     

Large ice‐contact slope movements: glacial buttressing,deformation and erosion

Glaciers and slope movements may act simultaneously to erode and modify glaciated slopes. Undercutting by glaciers can destabilize slopes but the extent to which slope failure may progress prior to subsequent glacier withdrawal has not hitherto been considered. The traditional view has been that the buttressing effect of ice prevents slope movement. The problem with this view is that ice is one‐third the density of rock and flows under low applied stress. Consequently, failed slopes may move into the glacier if they exert a stress in excess of the resistance provided by the glacier. Slope movement rate depends on ice rheology and other factors influencing driving and resisting stresses. Simple viscous equations are used to investigate these variables. The equations predict that small (<125 000 m³) ice‐contact rockslides can deform ice at several mm/year, increasing to several m/year for very large (>10⁸ m³) rockslides. To test these estimates, field evidence is presented of slope movements in glaciated valleys of New Zealand; narrowing or squeezing of glaciers adjacent to unstable rock slopes is demonstrated and considered to be the result of slope movement. For one site, geomorphic mapping and slope movement monitoring data show that movement rates are of similar order of magnitude to those predicted by the viscous equations; closer agreement could be achieved with the application of modelling techniques that can more realistically model the complex slope geometries and stability factors encountered, or by obtaining additional empirical data to calibrate the models. This research implies that, while the concept of glacial debuttressing – the reduction of slope support from withdrawal of glaciers – is valid, complete debuttressing is not a prerequisite for the movement of ice‐contact rock slopes. These slope movements may contribute to the erosional processes of glaciers and the evolution of glaciated slopes in a previously unrecognized way. Copyright © 2012 John Wiley & Sons, Ltd.     

Late Holocene slip rate and average recurrence interval of great earthquakes on the Shangzhi segment of the Yilan‐Yitong Fault Zone,northeastern China: Constraints from paleo‐earthquakes and historical written records

The Yilan‐Yitong Fault Zone (YYFZ) is considered to be the key branch of the Tancheng‐Lujiang Fault Zone (TLFZ) in northeastern China. Although the Mesozoic and early Cenozoic deformation of the YYFZ has been studied intensively over the past century, few estimates of slip rate and recurrence interval of large earthquakes in the late Quaternary, which are the two most important parameters for understanding the potential seismic hazard of this crucial structure, were obtained. Based on integrated interpretations of high resolution satellite images and detailed geologic and geomorphic mapping, linear landforms were identified, including fault scarps and troughs, along the Shangzhi segment of the YYFZ, which exceeds 25 km in length. Synthesized results of trench excavations and differential GPS measurements of terrace surfaces indicate that two events (E1, E2) occurred along the Shangzhi segment during the late Holocene, which resulted in 3.2 ±0.1 m of total vertical co‐seismic displacement with clear features of thrust motion. ¹⁴C dating of samples suggests that event E1 occurred between 440 ±30 years BP and 180 ±30 years BP and that event E2 occurred between 4 090 ±30 years BP and 3 880 ±30 years BP, which indicates that the minimum vertical slip rate of the Shangzhi segment of the YYFZ has been approximately 0.8 ±0.03 mm/year during the late Holocene. Constraints from paleo events and the slip rate suggest that the average recurrence interval of major earthquakes on the YYFZ is 3 800 ±200 years. Historical documents in Korea show that event E1 possibly corresponds to the earthquake that occurred in AD 1810 (the Qing Dynasty in Chinese history) in the Ningguta area, which had surface‐wave magnitude (Ms) of 6.8–7.5. Studies of kinematics show that the right‐lateral strike‐slip with a reverse component has been dominant along the YYFZ during the late Holocene.     

A process‐based conversion model for caesium‐137 derived erosion rates on agricultural land: an integrated spatial approach

There is increasing recognition that ¹³⁷Cs data remain one of the few sources of spatially distributed information concerning soil erosion. However, many of the conversion models that have been used to convert ¹³⁷Cs data into soil redistribution rates failed to account for some of the key factors affecting the redistribution of ¹³⁷Cs in agricultural landscapes. The conversion model presented in this paper aims to overcome some of the limitations associated with existing models and therefore to provide more realistic estimates of soil erosion rates on agricultural land. The conversion model aims at coupling soil redistribution processes directly with ¹³⁷Cs redistribution. Emphasis is placed on the spatial representation of soil redistribution processes and the adequate simulation of tillage processes. The benefits of the presented model arise from the two‐dimensional spatial integration of mass balance models with soil erosion models. No a priori assumptions about the intensity of any soil redistribution process are necessary and the level of agreement between observed and simulated ¹³⁷Cs inventories enables us to evaluate the performance of the model. The spatial implementation and the use of fuzzy parameter sets also allow us to assess the uncertainties associated with soil erosion estimates. It was shown that an adequate simulation of tillage processes is necessary and that simplified tillage models may lead to erroneous estimates of soil redistribution. The model was successfully applied to a study site in the Belgian Loam Belt and the results indicated that tillage is the dominant process. Furthermore, the uncertainties associated with the estimation of water erosion rates were much higher than those associated with tillage, especially for depositional areas. Copyright © 2003 John Wiley & Sons, Ltd.     

Range imaging: a new method for high‐resolution topographic measurements in small‐ and medium‐scale field sites

Topographic measurements are essential for the study of earth surface processes. Three‐dimensional data have been conventionally obtained through terrestrial laser scanning or photogrammetric methods. However, particularly in steep and rough terrain, high‐resolution field measurements remain challenging and often require new creative approaches. In this paper, range imaging is evaluated as an alternative method for obtaining surface data in such complex environments. Range imaging is an emerging time‐of‐flight technology, using phase shift measurements on a multi‐pixel sensor to generate a distance image of a surface. Its suitability for field measurements has yet not been tested. We found ambient light and surface reflectivity to be the main factors affecting error in distance measurements. Low‐reflectivity surfaces and strong illumination contrasts under direct exposure to sunlight lead to noisy distance measurements. However, regardless of lighting conditions, the accuracy of range imaging was markedly improved by averaging multiple images of the same scene. For medium ambient lighting (shade) and a light‐coloured surface the measurement uncertainty was approximately 9 mm. To further test the suitability of range imaging for field applications we measured a reach of a steep mountain stream with a horizontal resolution of approximately 1 cm (in the focal plane of the camera), allowing for the interpolation of a digital elevation model on a 2 cm grid. Comparison with an elevation model obtained from terrestrial laser scanning for the same site revealed that both models show similar degrees of topographic detail. Despite limitations in measurement range and accuracy, particularly at bright ambient lighting, range imaging offers three‐dimensional data in real time and video mode without the need of post‐processing. Therefore, range imaging is a useful complement or alternative to existing methods for high‐resolution measurements in small‐ to medium‐scale field sites. Copyright © 2012 John Wiley & Sons, Ltd.     

Spatial distribution and content of soil organic matter in an agricultural field in eastern Canada,as estimated from geostatistical tools

Soil erosion induces soil redistribution within the landscape and thus contributes to the spatial variability of soil quality. This study complements a previous experimentation initiated by the authors focusing on soil redistribution as a result of soil erosion, as indicated by caesium‐137 (¹³⁷Cs) measurements, in a small agricultural field in Canada. The spatial variability of soil organic matter (SOM) was characterized using geostatistics, which consider the randomized and structured nature of spatial variables and the spatial distribution of the samples. The spatial correlation of SOM (in percentages) patterns in the topsoil was established taking into account the spatial structure present in the data. A significant autocorrelation and reliable variograms were found with a R² ≥ 0·9, thus demonstrating a strong spatial dependence. Ordinary Kriging (OK) interpolation provided the best cross validation (r² = 0·35). OK and inverse distance weighting power two (IDW2) interpolation approaches produced similar estimates of the total SOM content of the topsoil (0–20 cm) of the experimental field, i.e. 211 and 213 tonnes, respectively. However, the two approaches produced differences in the spatial distribution patterns and the relative magnitude of some SOM content classes. The spatialization of SOM and soil redistribution variability – as evidenced by ¹³⁷Cs measurements – is a first step towards the assessment of the impact of soil erosion on SOM losses to recommend conservation measures. Copyright © 2009 John Wiley & Sons, Ltd.     

Continuous field estimation of dissolved organic carbon concentration and biochemical oxygen demand using dual‐wavelength fluorescence,turbidity and temperature

Dissolved organic matter (DOM) quality and quantity is not measured routinely in‐situ limiting our ability to quantify DOM process dynamics. This is problematic given legislative obligations to determine event based variability; however, recent advances in field deployable optical sensing technology provide the opportunity to address this problem. In this paper, we outline a new approach for in‐situ quantification of DOM quantity (Dissolved Organic Carbon: DOC) and a component of quality (Biochemical Oxygen Demand: BOD) using a multi‐wavelength, through‐flow fluorescence sensor. The sensor measured tryptophan‐like (Peak T) and humic‐like (Peak C) fluorescence, alongside water temperature and turbidity. Laboratory derived coefficients were developed to compensate for thermal quenching and turbidity interference (i.e., light attenuation and scattering). Field tests were undertaken on an urban river with ageing wastewater and stormwater infrastructure (Bourn Brook; Birmingham, UK). Sensor output was validated against laboratory determinations of DOC and BOD collected by discrete grab sampling during baseflow and stormflow conditions. Data driven regression models were then compared to laboratory correction methods. A combination of temperature and turbidity compensated Peak T and Peak C was found to be a good predictor of DOC concentration (R² = 0.92). Conversely, using temperature and turbidity correction coefficients provided low predictive power for BOD (R² = 0.46 and R² = 0.51, for Peak C and T, respectively). For this study system, turbidity appeared to be a reasonable proxy for BOD, R² = 0.86. However, a linear mixed effect model with temperature compensated Peak T and turbidity provided a robust BOD prediction (R² = 0.95). These findings indicate that with careful initial calibration, multi‐wavelength fluorescence, coupled with turbidity, and temperature provides a feasible proxy for continuous, in‐situ measurement of DOC concentration and BOD. This approach represents a cost effective monitoring solution, particularly when compared to UV – absorbance sensors and DOC analysers, and could be readily adopted for research and industrial applications.     

Projecting climate change impacts on stream flow regimes with tracer‐aided runoff models ‐ preliminary assessment of heterogeneity at the mesoscale

The northern mid‐high latitudes form a region that is sensitive to climate change, and many areas already have seen – or are projected to see – marked changes in hydroclimatic drivers on catchment hydrological function. In this paper, we use tracer‐aided conceptual runoff models to investigate such impacts in a mesoscale (749 km²) catchment in northern Scotland. The catchment encompasses both sub‐arctic montane sub‐catchments with high precipitation and significant snow influence and drier, warmer lowland sub‐catchments. We used downscaled HadCM3 General Circulation Model outputs through the UKCP09 stochastic weather generator to project the future climate. This was based on synthetic precipitation and temperature time series generated from three climate change scenarios under low, medium and high greenhouse gas emissions. Within an uncertainty framework, we examined the impact of climate change at the monthly, seasonal and annual scales and projected impacts on flow regimes in upland and lowland sub‐catchments using hydrological models with appropriate process conceptualization for each landscape unit. The results reveal landscape‐specific sensitivity to climate change. In the uplands, higher temperatures result in diminishing snow influence which increases winter flows, with a concomitant decline in spring flows as melt reduces. In the lowlands, increases in air temperatures and re‐distribution of precipitation towards autumn and winter lead to strongly reduced summer flows despite increasing annual precipitation. The integration at the catchment outlet moderates these seasonal extremes expected in the headwaters. This highlights the intimate connection between hydrological dynamics and catchment characteristics which reflect landscape evolution. It also indicates that spatial variability of changes in climatic forcing combined with differential landscape sensitivity in large heterogeneous catchments can lead to higher resilience of the integrated runoff response. Copyright © 2012 John Wiley & Sons, Ltd.     

An ArcGIS® Approach to Include Tectonic Structures in Point Data Regionalization

Point data derived from drilling logs must often be regionalized. However, aquifers may show discontinuous surface structures, such as the offset of an aquitard caused by tectonic faults. One main challenge has been to incorporate these structures into the regionalization process of point data. We combined ordinary kriging and inverse distance weighted (IDW) interpolation to account for neotectonic structures in the regionalization process. The study area chosen to test this approach is the largest porous aquifer in Austria. It consists of three basins formed by neotectonic events and delimited by steep faults with a vertical offset of the aquitard up to 70 m within very short distances. First, ordinary kriging was used to incorporate the characteristic spatial variability of the aquitard location by means of a variogram. The tectonic faults could be included into the regionalization process by using breaklines with buffer zones. All data points inside the buffer were deleted. Last, IDW was performed, resulting in an aquitard map representing the discontinuous surface structures. This approach enables one to account for such surfaces using the standard software package ArcGIS^®; therefore, it could be adopted in many practical applications.     

Comparison of scenario‐neutral approaches for estimation of climate change impacts on flood characteristics

Scenario‐neutral assessments of climate change impact on floods analyse the sensitivity of a catchment to a range of changes in selected meteorological variables such as temperature and precipitation. The key challenges of the approach are the choice of the meteorological variables and statistics thereof and how to generate time series representing altered climatologies of the selected variables. Different methods have been proposed to achieve this, and it remains unclear if and to which extent they result in comparable flood change projections. Here, we compare projections of annual maximum floods (AMFs) derived from three different scenario‐neutral methods for a prealpine study catchment. The methods chosen use different types of meteorological data, namely, observations, regional climate model output, and weather generator data. The different time series account for projected changes in the seasonality of temperature and precipitation, in the occurrence statistics of precipitation, and of daily precipitation extremes. Resulting change in mean AMF peak magnitudes and volumes differs in sign between the methods (range of ?6% to +7% for flood peak magnitudes and ?11% to +14% for flood volumes). Moreover, variability of projected peak magnitudes and flood volumes depends on method with one approach leading to a generally larger spread. The differences between the methods vary depending on whether peak magnitudes or flood volumes are considered and different relationships between peak magnitude and volume change result. These findings can be linked to differing flood regime changes among the three approaches. The study highlights that considering selected aspects of climate change only when performing scenario‐neutral studies may lead to differing representations of flood generating processes by the approaches and thus different quantifications of flood change. As each method comes with its own strengths and weaknesses, it is recommended to combine several scenario‐neutral approaches to obtain more robust results.     

Optimal Design of Multitype Groundwater Monitoring Networks Using Easily Accessible Tools

Monitoring networks are expensive to establish and to maintain. In this paper, we extend an existing data‐worth estimation method from the suite of PEST utilities with a global optimization method for optimal sensor placement (called optimal design) in groundwater monitoring networks. Design optimization can include multiple simultaneous sensor locations and multiple sensor types. Both location and sensor type are treated simultaneously as decision variables. Our method combines linear uncertainty quantification and a modified genetic algorithm for discrete multilocation, multitype search. The efficiency of the global optimization is enhanced by an archive of past samples and parallel computing. We demonstrate our methodology for a groundwater monitoring network at the Steinlach experimental site, south‐western Germany, which has been established to monitor river—groundwater exchange processes. The target of optimization is the best possible exploration for minimum variance in predicting the mean travel time of the hyporheic exchange. Our results demonstrate that the information gain of monitoring network designs can be explored efficiently and with easily accessible tools prior to taking new field measurements or installing additional measurement points. The proposed methods proved to be efficient and can be applied for model‐based optimal design of any type of monitoring network in approximately linear systems. Our key contributions are (1) the use of easy‐to‐implement tools for an otherwise complex task and (2) yet to consider data‐worth interdependencies in simultaneous optimization of multiple sensor locations and sensor types.     

Formation,uplift and dissection of planation surfaces at passive continental margins – a new approach

The usefulness of large‐scale, low‐relief, high‐level landscapes as markers of uplift events has become a subject of disagreement among geomorphologists. We argue that the formation of low‐relief surfaces over areas of large extent and cutting across bedrock of different age and resistance must have been guided by distinct base levels. In the absence of other options the most likely base level is sea level. We have analysed West Greenland landscapes in a recent study by combining the cooling history from apatite fission‐track analysis (AFTA) data with the denudation history from landscape analysis and the stratigraphic record. An important difference between our approach and that of classical geomorphology is that we now have the ability to document when thick sections of rocks have been deposited and then removed. The present‐day high‐level plateau in West Greenland is the remnant of a planation surface that was formed by denudation that lasted c. 20 million years during which up to 1 km of cover was removed after maximum burial at the Eocene–Oligocene transition. Here we present additional AFTA data to show that the planation surface is the end‐product of Cenozoic denudation even in basement areas and argue that Phanerozoic sediments – most likely of Cretaceous–Palaeogene age – must have been present prior to denudation. The planation surface was offset by reactivated faults and uplifted to present‐day altitudes of up to 2 km. The uplift occurred in two late Neogene phases that caused incision of valleys below the planation surface and their subsequent uplift. We therefore find that the elevated and deeply dissected plateau is evidence of episodic post‐rift uplift that took place millions of years after cessation of sea‐floor spreading west of Greenland. We suggest that other margins with similar morphology may also be characterized by episodic post‐rift uplift unrelated to the processes of rifting and continental separation, rather than being permanently uplifted since the time of rifting, as is commonly assumed. Copyright © 2009 John Wiley & Sons, Ltd.     

MHYDAS‐Erosion: a distributed single‐storm water erosion model for agricultural catchments

In this paper, we present MHYDAS‐Erosion, a dynamic and distributed single‐storm water erosion model developed as a module of the existing hydrological MHYDAS model. As with many catchment erosion models, MHYDAS‐Erosion is able to simulate sediment transport, erosion and deposition by rill and interrill processes. Its originality stems from its capacity to integrate the impact of land management practices (LMP) as key elements controlling the sedimentological connectivity in agricultural catchments. To this end, the water‐sediment pathways are first determined by a specific process‐oriented procedure defined and controlled by the user, which makes the integration of LMP easier. The LMP dynamic behaviours are then integrated into the model as a time‐dependent function of hydrological variables and LMP characteristics. The first version of the model was implemented for vegetative filters and tested using water and sediment discharge measurements at three nested scales of a densely instrumented catchment (Roujan, OMERE Observatory, southern France). The results of discharge and soil loss for simulated rainfall events have been found to acceptably compare with available data. The average R² values for water and sediment discharge are 0·82 and 0·83, respectively. The sensitivity of the model to changes in the proportion of LMP was assessed for a single rain event by considering three scenarios of the Roujan catchment management with vegetative filters: 0% (Scenario 1), 18% (Scenario 2, real case) and 100% (Scenario 3). Compared to Scenario 2 (real case), soil losses decreased for Scenario 3 by 65% on the agricultural plot scale, 62% on the sub‐catchment scale and 45% at the outlet of the catchment and increased for Scenario 1 by 0% on the plot scale, 26% on the sub‐catchment scale and 18% at the outlet of the catchment. Copyright © 2010 John Wiley & Sons, Ltd.     

SOFIE project – 3D shaking table test on a seven‐storey full‐scale cross‐laminated timber building

Multi‐storey buildings made of cross‐laminated timber panels (X‐lam) are becoming a stronger and economically valid alternative in Europe compared with traditional masonry or concrete buildings. During the design process of these multi‐storey buildings, also their earthquake behaviour has to be addressed, especially in seismic‐prone areas such as Italy. However, limited knowledge on the seismic performance is available for this innovative massive timber product. On the basis of extensive testing series comprising monotonic and reversed cyclic tests on X‐lam panels, a pseudodynamic test on a one‐storey X‐lam specimen and 1D shaking table tests on a full‐scale three‐storey specimen, a full‐scale seven‐storey building was designed according to the European seismic standard Eurocode 8 and subjected to earthquake loading on a 3D shaking table. The building was designed with a preliminary action reduction factor of three that had been derived from the experimental results on the three‐storey building. The outcomes of this comprehensive research project called ‘SOFIE – Sistema Costruttivo Fiemme’ proved the suitability of multi‐storey X‐lam structures for earthquake‐prone regions. The buildings demonstrated self‐centring capabilities and high stiffness combined with sufficient ductility to avoid brittle failures. The tests provided useful information for the seismic design with force‐based methods as defined in Eurocode 8, that is, a preliminary experimentally based action reduction factor of three was confirmed. Valid, ductile joint assemblies were developed, and their importance for the energy dissipation in buildings with rigid X‐lam panels became evident. The seven‐storey building showed relatively high accelerations in the upper storeys, which could lead to secondary damage and which have to be addressed in future research. Copyright © 2013 John Wiley & Sons, Ltd.