Identification of physical processes inherent in artificial neural network rainfall runoff models

The emergence of artificial neural network (ANN) technology has provided many promising results in the field of hydrology and water resources simulation. However, one of the major criticisms of ANN hydrologic models is that they do not consider/explain the underlying physical processes in a watershed, resulting in them being labelled as black‐box models. This paper discusses a research study conducted in order to examine whether or not the physical processes in a watershed are inherent in a trained ANN rainfall‐runoff model. The investigation is based on analysing definite statistical measures of strength of relationship between the disintegrated hidden neuron responses of an ANN model and its input variables, as well as various deterministic components of a conceptual rainfall‐runoff model. The approach is illustrated by presenting a case study for the Kentucky River watershed. The results suggest that the distributed structure of the ANN is able to capture certain physical behaviour of the rainfall‐runoff process. The results demonstrate that the hidden neurons in the ANN rainfall‐runoff model approximate various components of the hydrologic system, such as infiltration, base flow, and delayed and quick surface flow, etc., and represent the rising limb and different portions of the falling limb of a flow hydrograph. Copyright © 2004 John Wiley & Sons, Ltd.     

Vertical stress distribution due to an arbitrarily shaped foundation using triangulation and an analytical expression for a triangular loaded region

Using the basic Boussinesq's equation, the expression for the vertical stress distribution (σ_z) underneath any point on the ground surface due to a general triangular loaded region in a preferred orientation with a linearly varied loading has been successfully derived. When the triangle is not in a preferred orientation, a simple axis transformation is required and the expression will be equally applicable. Based on this expression, σ_z due to an arbitrarily shaped loaded foundation can simply be determined by first triangulating the loaded area and summing up the contributions from each generated triangular region. The procedures for triangulating and calculating the stress distribution can be simply automated through computer programs.     

Response of circular footings and anchor plates in non-homogeneous elastic soils

The axisymmetric elastic response of circular footings and anchor plates in a linearly non-homogeneous elastic soil is analysed. It is assumed that footings/anchors are flexible and subjected to axisymmetric vertical loads. The response of the footing/anchor is modelled by using the classical Poisson–Kirchhoff thin plate theory. A variational technique is used to analyse the interaction problem. A representation for the contact stress is established by using a fundamental solution corresponding to a unit vertical pressure acting over an annular region in the interior of the non-homogeneous soil. The fundamental solution can be derived by using rigorous analytical procedures. The influence of the footing flexibility and the degree of soil non-homogeneity on the displacements, bending moments and contact stresses of a surface footing is examined over a wide range of governing parameters. In the case of anchor plates the influence of depth of embedment, degree of soil non-homogeneity and anchor flexibility on the anchor displacement is investigated.     

Small angle neutron scattering by opals

Micro- and non-crystalline opals, chalcedony and flint show diffuse small angle neutron scattering (SANS). Precious opals give rise to two additional intensity maxima at very small scattering angles which are due to Bragg reflections from the closest packed non-crystalline silica spheres. A small angle texture diagram reveals that the closest packing is faulty. Synthetic non-crystalline opals yield much less intense small angle scattering due to lower contrast between silica spheres and interstitial cement or particles; in this case intensity maxima were not observed. The outer part of the scattering curves of opal-CT and microcrystalline quartz deviates from Porod's law. The specific surface of natural non-crystalline opals ranges from 0.006 to 0.018 nm^–1. In microcrystalline opals, the specific surface is about 10 times larger than in non-crystalline opals.     

Comparison of results from several AOGCMs for global and regional sea-level change 1900–2100

 Sea-level rise is an important aspect of climate change because of its impact on society and ecosystems. Here we present an intercomparison of results from ten coupled atmosphere-ocean general circulation models (AOGCMs) for sea-level changes simulated for the twentieth century and projected to occur during the twenty first century in experiments following scenario IS92a for greenhouse gases and sulphate aerosols. The model results suggest that the rate of sea-level rise due to thermal expansion of sea water has increased during the twentieth century, but the small set of tide gauges with long records might not be adequate to detect this acceleration. The rate of sea-level rise due to thermal expansion continues to increase throughout the twenty first century, and the projected total is consequently larger than in the twentieth century; for 1990–2090 it amounts to 0.20–0.37 m. This wide range results from systematic uncertainty in modelling of climate change and of heat uptake by the ocean. The AOGCMs agree that sea-level rise is expected to be geographically non-uniform, with some regions experiencing as much as twice the global average, and others practically zero, but they do not agree about the geographical pattern. The lack of agreement indicates that we cannot currently have confidence in projections of local sea-level changes, and reveals a need for detailed analysis and intercomparison in order to understand and reduce the disagreements. Received: 1 September 2000 / Accepted: 20 April 2001     

Observational Selection of Comet Discoveries

The distribution of perihelion distances of nearly periodic comets was studied using different assumptions about comet discovery probabilities. It was shown that the concentration of perihelia at small heliocentric distances cannot be explained solely by the influence of visibility conditions.     

Grain boundary diffusion in enstatite

We have investigated grain boundary diffusion rates in enstatite by heating single crystals of quartz packed in powdered San Carlos olivine (Mg_0.90Fe_0.10)₂SiO₄ at controlled oxygen fugacities in the range 10^?5.7 to 10^?8.7?atm and temperatures from 1350° to 1450?°C for times from 5 to 100?h at 1?atm total pressure. Following the experiments, the thickness of the coherent polycrystalline reaction rim of pyroxene that had formed between the quartz and olivine was measured using backscatter scanning imaging in the electron microprobe. Quantitative microprobe analysis indicated that the composition of this reaction phase is (Mg_0.92Fe_0.08)₂Si₂O₆. The rate of growth of the pyroxene increases with increasing temperature, is independent of the oxygen fugacity, and is consistent with a parabolic rate law, indicating that the growth rate is controlled by ionic diffusion through the pyroxene rim. Microstructural observations and platinum marker experiments suggest that the reaction phase is formed at the olivine-pyroxene interface, and is therefore controlled by the diffusion of silicon and oxygen. The parabolic rate constants determined from the experiments were analyzed in terms of the oxide activity gradient across the rim to yield mean effective diffusivities for the rate-limiting ionic species, assuming bulk transport through the pyroxene layer. These effective diffusivities are faster than the lattice diffusivities for the slowest species (silicon) calculated from creep experiments, but slower than measured lattice diffusivities for oxygen in enstatite. Thus, silicon grain boundary diffusion is most likely to be the rate-limiting process in the growth of the pyroxene rims. Also, as oxygen transport through the pyroxene rims must be faster than silicon transport, diffusion of oxygen along the grain boundaries must be faster than through the lattice. The grain boundary diffusivity for silicon in orthopyroxenite is then given by D¯gbSiδ=(3.3±3.0)×10^?9f0.0O₂e^?400±65/RT?m³s^?1, where the activation energy for diffusion is in kJ/mol, and δ is the grain boundary width in m. Calculated growth rates for enstatite under these conditions are significantly slower than predicted by an extrapolation from similar experiments performed at 1000?°C under high pressure (hydrous) conditions by Yund and Tullis (1992), perhaps due to water-enhancement of diffusion in their experiments.