Parameter estimation in semi‐distributed hydrological catchment modelling using a multi‐criteria objective function

Output generated by hydrologic simulation models is traditionally calibrated and validated using split‐samples of observed time series of total water flow, measured at the drainage outlet of the river basin. Although this approach might yield an optimal set of model parameters, capable of reproducing the total flow, it has been observed that the flow components making up the total flow are often poorly reproduced. Previous research suggests that notwithstanding the underlying physical processes are often poorly mimicked through calibration of a set of parameters hydrologic models most of the time acceptably estimates the total flow. The objective of this study was to calibrate and validate a computer‐based hydrologic model with respect to the total and slow flow. The quick flow component used in this study was taken as the difference between the total and slow flow. Model calibrations were pursued on the basis of comparing the simulated output with the observed total and slow flow using qualitative (graphical) assessments and quantitative (statistical) indicators. The study was conducted using the Soil and Water Assessment Tool (SWAT) model and a 10‐year historical record (1986–1995) of the daily flow components of the Grote Nete River basin (Belgium). The data of the period 1986–1989 were used for model calibration and data of the period 1990–1995 for model validation. The predicted daily average total flow matched the observed values with a Nash–Sutcliff coefficient of 0·67 during calibration and 0·66 during validation. The Nash–Sutcliff coefficient for slow flow was 0·72 during calibration and 0·61 during validation. Analysis of high and low flows indicated that the model is unbiased. A sensitivity analysis revealed that for the modelling of the daily total flow, accurate estimation of all 10 calibration parameters in the SWAT model is justified, while for the slow flow processes only 4 out of the set of 10 parameters were identified as most sensitive. Copyright © 2007 John Wiley & Sons, Ltd.     

Modal analysis of non-classically damped linear systems

A critical, textbook-like review of the generalized modal superposition method of evaluating the dynamic response of nonclassically damped linear systems is presented, which it is hoped will increase the attractiveness of the method to structural engineers and its application in structural engineering practice and research. Special attention is given to identifying the physical significance of the various elements of the solution and to simplifying its implementation. It is shown that the displacements of a non-classically damped n-degree-of-freedom system may be expressed as a linear combination of the displacements and velocities of n similarly excited single-degree-of-freedom systems, and that once the natural frequencies of vibration of the system have been determined, its response to an arbitrary excitation may be computed with only minimal computational effort beyond that required for the analysis of a classically damped system of the same size. The concepts involved are illustrated by a series of examples, and comprehensive numerical data for a three-degree-of-freedom system are presented which elucidate the effects of several important parameters. The exact solutions for the system are also compared over a wide range of conditions with those computed approximately considering the system to be classically damped, and the interrelationship of two sets of solutions is discussed.     

Seismic risk assessment for the mega-city of Istanbul: Ductility, strength and maximum interstory drift demands

This study investigates the seismic demands due to the catastrophic 1999 M_w=7.4 Kocaeli and the M_w=7.1 Duzce earthquakes. The inelastic response spectra for the Kocaeli and Duzce earthquakes are investigated for systems with known strength and ductility. An analytical fiber element model is developed for a typical reinforced concrete building in Turkey. The interstory drifts are calculated from nonlinear dynamic analysis using 26 recorded strong-motion data from the 1999 Kocaeli and the Duzce earthquakes. In the dynamic analysis, the structural members are modeled by employing distributed plasticity fiber elements and both geometrical as well as material nonlinearities are taken into account. This study shows that the ductility and interstory drift demands due to the Kocaeli and the Duzce earthquakes were very severe (well above the code prescribed values) even for moderately inelastic structures. It is apparent from the results that the forward directivity effect is the most influential factor on the interstory drift demand. Both the distance to the fault rupture and the site conditions affect seismic demands, but the site conditions and the local topography are more influential than the distance from the fault rupture. This study shows that substantial damage should be expected in a future earthquake at all districts of Istanbul, but especially at Avcilar, Cekmece, Fatih, Bakirkoy and the Zeytinburnu districts. It is also shown that base isolation may substantially improve the performance of a structure in the inelastic domain and base isolated structures may be designed for lower minimum lateral strengths and higher strength reduction factors.     

Levels and spatial distribution of polycyclic aromatic hydrocarbons (PAHs) in sediments from Lenga Estuary, central Chile

The Lenga Estuary is a small brackish wetland located southwest of San Vicente Bay, Region VIII, Chile. Surface sediment from nine sites in the estuary were analysed for PAHs and compared to Sediment Quality Guidelines (SQG). Sediment samples were freeze dried and soxhlet extracted for 16 h using DCM. Identification and quantification was carried out by HPLC. Organic carbon was also determined. Results showed total PAH concentrations ranged from 290 to 6118 (2025 ± 1975) ng g⁻¹ d.w. (2025 ± 1975). Results for organic carbon percentages ranged from 1% to 7%. Statistical analysis showed a significant positive correlation (Pearson test) between organic carbon percentage PAHs. Comparison of contaminant levels and international Sediment Quality Guidelines (SQG) (ERL and ER) suggested that sediment of the Lenga estuary did not show any ecotoxicologial risk for benthic organisms where high levels of PAHs were detected. Monitoring of this and other contaminants is recommended in Chile.     

Levels and spatial distribution of polycyclic aromatic hydrocarbons (PAHs) in superficial sediment from 15 Italian marine protected areas (MPA)

Surface sediment from 15 Italian marine protected areas (MPA) were analysed for polycyclic aromatic hydrocarbons (PAHs). The organic carbon percentage was also determined. Total PAH concentrations (ng g⁻¹ d.w.) ranged from 0.71 (Penisola del Sinis) to 1550 (Miramare). Individual PAH analysis showed that three and four rings PAHs were the most frequently detected isomers and accounted for 60-70% of the PAH total concentrations. PAH ratio analysis showed a prevalence of pyrolytic PAH origin at most of the MPAs with exception of Porto Cesareo and Ustica where a petrogenic origin was detected. Results for organic carbon percentages ranged from 0.3% (Capo Rizzuto) to 2% (Punta Campanella). These results are comparable to other Mediterranean marine environments. However, our results shows that some MPAs, such as Miramare, Porto Cesareo, Isola Capo Rizzuto and Punta Campanella, are subject to strong pressure from urban and industrial activities where high PAH levels were detected.     

Kinematic significance of mingling-rolling structures in lava flows: a case study from Porri Volcano (Salina, Southern Tyrrhenian Sea)

 A basaltic andesite lava flow from Porri Volcano (Salina, Southern Tyrrhenian Sea) is composed of two different magmas. Magma A (51 vol.% of crystals) has a dacitic glass composition, and magma B (18 vol.% of crystals), a basaltic glass composition. Magma B is hosted in A and consists of sub-spherical enclaves and boudin-like, banding and rolling structures (RS). Four types of RS have been recognized: σ–type;δ–type; complex σ-δ–types and transitional structures between sub-spherical enclaves and rolling structures. An analysis of the RS has been performed in order to reconstruct the flow kinematics and the mechanism of flow emplacement. Rolling structures have been selected in three sites located at different distances from the vent. In all sites most RS show the same sense of shear. Kinematic analysis of RS allows the degree of flow non-coaxiality to be determined. The non-coaxiality is expressed by the kinematic vorticity number Wk, a measure of the ratio Sr between pure shear strain rate and simple shear strain rate. The values of Wk calculated from the measured shapes of microscopic RS increase with increasing distance from the vent, from approximately 0.5 to 0.9. Results of the structural analysis reveal that the RS formed during the early–intermediate stage of flow emplacement. They represent originally sub-spherical enclaves deformed at low shear strain. At higher strain, RS deformed to give boudin-like and stretched banding structures. Results of the kinematic analysis suggest that high viscosity lava flows are heterogeneous non-ideal shear flows in which the degree of non-coaxiality increases with the distance from the vent. In the vent area, deformation is intermediate between simple shear and pure shear. Farther from the vent, deformation approaches ideal simple shear. Lateral extension processes occur only in the near-vent zone, where they develop in response to the lateral push of magma extruded from the vent. Lateral shortening processes develop in the distal zone and record the gravity-driven movement of the lava. The lava flow advanced by two main mechanisms, lateral translation and rolling motion. Lateral translation equals rolling near the vent, while rolling motion prevailed in the distal zones. Received: 6 November 1997 / Accepted: 29 November 1997     

Impact of coupled perturbations of atmospheric trace gases on Earth's climate and ozone

We have studied the effects on the ozone concentration and surface temperature, of perturbations in the atmospheric content of nitrous oxide, methane, carbon dioxide and chlorofluorocarbons (CFC). The sensitivity study has been carried out with a radiative-convective-photochemical model. The doubling of carbon dioxide concentration has the effect of warming the troposphere and cooling the stratosphere. As a result of this cooling, the change of ozone columnar density produced by 10 ppb of chlorine amount to 9.3% as compared to –10.9% obtained without temperature feedback. Perturbation in nitrous oxide correspond to an increase in NO _x of the stratosphere with consequent ozone reduction while doubling the methane concentration correspond to a slight increase in columnar density. The effect of the increased methane concentration in the stratosphere contributes to reduce the effect of CFC due to the enhanced formation of HCl. The perturbation of these two minor constituents appreciably increase the greenhouse effect to 2.30 from 1.67°, obtained when carbon dioxide alone is considered.     

Google Earth as a data source for investigating river forms and processes: Discriminating river types using form-based process indicators

Google Earth provides potential for exploiting an enormous reservoir of freely-available remotely sensed data to support river science and management. In this paper, we consider how the platform can support investigation of river physical forms and processes by developing an empirically-based reach-scale classification of semi-natural European single thread to transitional rivers. Using strict reach and image selection criteria, we identified 194 reaches of 68 rivers for analysis. Measurements of channel dimensions and counts of in-channel and floodplain features, standardised for reach length and channel width where necessary, were used to derive a series of geomorphologically-relevant process indicators. A suite of multivariate analyses were then applied to this data set, resulting in the discrimination of five river types: laterally stable, laterally active sinuous-meandering; transitional (near-braided); bedrock; and cascade/step dominated. The results of the classification were tested by examining the characteristics and distribution of the river classes in relation to known independent controls of river form including reach-scale energy and valley confinement conditions. Our results show that if methods of data extraction are carefully developed, physically meaningful river reach discrimination can be achieved using Google Earth. Although there are limits to the types of information that can be extracted such that field investigations cannot always be avoided, there is enormous potential to mine Google Earth across different space and time scales, supporting the assembly of large, reliable data sets relevant to river forms and processes in a very cost-effective way. © 2019 John Wiley & Sons, Ltd.     

River,Coastal and Estuarine Morphodynamics Selected papers from the 10th anniversary of the RCEM Symposium

This Special Issue collects 17 selected contributions from participants to the 10th edition of the RCEM (River, Coastal and Estuarine Morphodynamics) Symposium, organized in Padova-Trento (Italy) in September 2017. The series of biennial RCEM symposia has the key goal of enhancing interaction and promoting integration among the scientific communities focused on the morphological dynamics of river, coastal and estuarine environments, through various combinations of theoretical, observational, experimental and modelling approaches. The 17 contributions to this Special Issue contain four state-of-science reviews and overall offer a broad view of the cross-cutting perspective adopted when addressing morphodynamics. Such a perspective accounts for the mutual interplay between morphology, fluid dynamics and other environmental factors, and has presently become a widespread paradigm to address landscape evolution. © 2020 John Wiley & Sons, Ltd.