Studies in the application of recurrence relations to special perturbation methods

The numerical integration of the differential equations describing dynamical systems has been shown in previous papers of this series to be most effectively accomplished by an explicit Taylor series method.In this paper we show that one explicit Taylor series method, developed earlier in this series and which appears to possess a high degree of versatility, yields considerable gains in efficiency over classical single-step and multi-step methods. (In this context efficiency is a measure of the time taken to carry out a calculation of a specific accuracy).For a given accuracy criterion governing the local truncation error (LTE) it is found that the Taylor series method is generallytwice as fast as the classical multi-step method and up totwenty times faster than the classical single-step method.     

Optimizing urban stormwater control strategies and assessing aquifer recharge through drywells in an urban watershed

Hydrogeology Journal - A coupled simulation-optimization model (SOM) is developed in this work that links the US Environmental Protection Agency’s Storm Water Management Model (SWMM) with a...     

Analysis of water control in an underground mine under strong karst media influence (Vazante mine,Brazil)

This work presents analysis of groundwater flow conditions and groundwater control measures for Vazante underground mine located in the state of Minas Gerais, Brazil. According to field observations, groundwater flow processes in this mine are highly influenced by the presence of karst features located in the near-surface terrain next to Santa Catarina River. The karstic features, such as caves, sinkholes, dolines and conduits, have direct contact with the aquifer and tend to increase water flow into the mine. These effects are more acute in areas under the influence of groundwater-level drawdown by pumping. Numerical analyses of this condition were carried out using the computer program FEFLOW. This program represents karstic features as one-dimensional discrete flow conduits inside a three-dimensional finite element structure representing the geologic medium following a combined discrete-continuum approach for representing the karst system. These features create preferential flow paths between the river and mine; their incorporation into the model is able to more realistically represent the hydrogeological environment of the mine surroundings. In order to mitigate the water-inflow problems, impermeabilization of the river through construction of a reinforced concrete channel was incorporated in the developed hydrogeological model. Different scenarios for channelization lengths for the most critical zones along the river were studied. Obtained results were able to compare effectiveness of different river channelization scenarios. It was also possible to determine whether the use of these impermeabilization measures would be able to reduce, in large part, the elevated costs of pumping inside the mine.     

A first classification scheme of flow-bed interaction for clay-laden density currents and soft substrates

Many aquatic environments exhibit soft, muddy substrates, but this important property has largely been ignored in process-based models of Earth-surface flow. Novel laboratory experiments were carried out to shed light on the feedback processes that occur when particulate density currents (turbidity currents) move over a soft mud substrate. These experiments revealed multiple types of flow-bed interaction and large variations in bed deformation and bed erosion, which are interpreted to be related to the interplay between the shear forces of the current and the stabilising forces in the bed. Changes in this force balance were simulated by varying the clay concentrations in the flow and in the bed. Five different interaction types are described, and dimensional and non-dimensional phase diagrams for flow-bed interaction are presented.     

Comparative efficiency analysis of different nonlinear modelling strategies to simulate the biaxial response of RC columns

The performance of different nonlinear modelling strategies to simulate the response of RC columns subjected to axial load combined with cyclic biaxial horizontal loading is compared. The models studied are classified into two categories according to the nonlinearity distribution assumed in the elements: lumped-plasticity and distributed inelasticity. For this study, results of tests on 24 columns subjected to cyclic uniaxial and biaxial lateral displacements were numerically reproduced. The analyses show that the global envelope response is satisfactorily represented with the three modelling strategies, but significant differences were found in the strength degradation for higher drift demands and energy dissipation.     

Lithology and fluid prediction from prestack seismic data using a Bayesian model with Markov process prior

We invert prestack seismic amplitude data to find rock properties of a vertical profile of the earth. In particular we focus on lithology, porosity and fluid. Our model includes vertical dependencies of the rock properties. This allows us to compute quantities valid for the full profile such as the probability that the vertical profile contains hydrocarbons and volume distributions of hydrocarbons. In a standard point wise approach, these quantities can not be assessed. We formulate the problem in a Bayesian framework, and model the vertical dependency using spatial statistics. The relation between rock properties and elastic parameters is established through a stochastic rock model, and a convolutional model links the reflectivity to the seismic. A Markov chain Monte Carlo (MCMC) algorithm is used to generate multiple realizations that honours both the seismic data and the prior beliefs and respects the additional constraints imposed by the vertical dependencies. Convergence plots are used to provide quality check of the algorithm and to compare it with a similar method. The implementation has been tested on three different data sets offshore Norway, among these one profile has well control. For all test cases the MCMC algorithm provides reliable estimates with uncertainty quantification within three hours. The inversion result is consistent with the observed well data. In the case example we show that the seismic amplitudes make a significant impact on the inversion result even if the data have a moderate well tie, and that this is due to the vertical dependency imposed on the lithology fluid classes in our model. The vertical correlation in elastic parameters mainly influences the upside potential of the volume distribution. The approach is best suited to evaluate a few selected vertical profiles since the MCMC algorithm is computer demanding.     

Sea water intrusion by sea-level rise: scenarios for the 21st century

This study presents a method to assess the contributions of 21st-century sea-level rise and groundwater extraction to sea water intrusion in coastal aquifers. Sea water intrusion is represented by the landward advance of the 10,000 mg/L iso-salinity line, a concentration of dissolved salts that renders groundwater unsuitable for human use. A mathematical formulation of the resolution of sea water intrusion among its causes was quantified via numerical simulation under scenarios of change in groundwater extraction and sea-level rise in the 21st century. The developed method is illustrated with simulations of sea water intrusion in the Seaside Area sub-basin near the City of Monterey, California (USA), where predictions of mean sea-level rise through the early 21st century range from 0.10 to 0.90 m due to increasing global mean surface temperature. The modeling simulation was carried out with a state-of-the-art numerical model that accounts for the effects of salinity on groundwater density and can approximate hydrostratigraphic geometry closely. Simulations of sea water intrusion corresponding to various combinations of groundwater extraction and sea-level rise established that groundwater extraction is the predominant driver of sea water intrusion in the study aquifer. The method presented in this work is applicable to coastal aquifers under a variety of other scenarios of change not considered in this work. For example, one could resolve what changes in groundwater extraction and/or sea level would cause specified levels of groundwater salinization at strategic locations and times.     

Seasonal patterns of benthic arthropods in a subtropical Andean basin

We studied the seasonal variation of arthropod assemblages and the response to high discharge events in a subtropical Andean basin. Using abundance of common taxa and taxon richness of each sampled site, we examined the temporal variation of arthropods and related these changes with 12 environmental variables. Seasonal patterns of benthic arthropods were confirmed by uni- and multivariate techniques. Benthic abundance and taxon richness peaked on autumn-winter months, and declined abruptly on the month of peak discharge (summer season). Distinct taxonomic groups were dominant in each season influenced by discharge and seven chemical variables.Assemblage resistance to disturbance by spates was low. Rapid recovery of benthic arthropods after floods reflected high resilience of the system. The dominance of groups of short life cycles such as Ephemeroptera played an important role in this process.