Fully coupled elastoplastic hydro-mechanical analysis of unsaturated porous media using a meshfree method

This paper presents the first application of an advanced meshfree method, ie, the edge-based smoothed point interpolation method (ESPIM), in simulation of the coupled hydro-mechanical behaviour of unsaturated porous media. In the proposed technique, the problem domain is spatially discretised using a triangular background mesh, and the polynomial point interpolation method combined with a simple node selection scheme is adopted for creating nodal shape functions. Smoothing domains are formed on top of the background mesh, and a constant smoothed strain, created by applying the smoothing operation over the smoothing domains, is assigned to each smoothing domain. The deformation and flow models are developed based on the equilibrium equation of the mixture, and linear momentum and mass balance equations of the fluid phases, respectively. The effective stress approach is followed to account for the coupling between the flow and deformation models. Further coupling among the phases is captured through a hysteretic soil water retention model that evolves with changes in void ratio. An advanced elastoplastic constitutive model within the context of the bounding surface plasticity theory is employed for predicting the nonlinear behaviour of soil skeleton. Time discretisation is performed by adopting a three-point discretisation method with growing time steps to avoid temporal instabilities. A modified Newton-Raphson framework is designed for dealing with nonlinearities of the discretised system of equations. The performance of the numerical model is examined through a number of numerical examples. The state-of-the-art computational scheme developed is useful for simulation of geotechnical engineering problems involving unsaturated soils.     

Evaluation of the temporal variations of groundwater storage and its interactions with climatic variables using GRACE data and hydrological models: A study from Turkey

Monitoring of the fluctuations of groundwater storage is particularly important in arid and semi-arid regions where water scarcity brings about various challenges. Remote sensing data and techniques play a preponderant role in developing solutions to environmental problems. The launch of Gravity Recovery and Climate Experiment (GRACE) satellites has eased the remote monitoring and evaluation of groundwater resources with an unprecedented precision over large scales. Within the scope of the current study, the latest release (RL06) of GRACE mass concentrations (Mascons) from Jet Propulsion Laboratory (JPL) dataset as well as Global Land Data Assimilation System (GLDAS) models of Noah and Catchment Land Surface Model (CLSM) were used to provide Groundwater Storage Anomalies (GWSA) over Turkey. The temporal interactions of the estimated GWSA with the climatic variables of precipitation and temperature (derived from the reanalysis datasets of CHELSA [Climatologies at High resolution for the Earth's Land Surface Areas] and FLDAS [the Famine Early Warning Systems Network Land Data Assimilation System], respectively) were investigated statistically. The results suggest that there is a descending trend (from 2003 to 2016) for Terrestrial Water Storage Anomalies (TWSA) and GWSA over Turkey with a total loss of 11 and 6 cm of water, respectively. The statistical analysis results also indicate that the monthly variations of GWSA over Turkey are highly correlated with precipitation and temperature at 2-month lag. The analysis of the climatology (long-term) values of monthly GWSA, precipitation and temperature also revealed high agreement between the variables.     

Geostatistical features of streambed vertical hydraulic conductivities in Frenchman Creek Watershed in Western Nebraska

This study evaluated the spatial variability of streambed vertical hydraulic conductivity (K_v) in different stream morphologies in the Frenchman Creek Watershed, Western Nebraska, using different variogram models. Streambed K_v values were determined in situ using permeameter tests at 10 sites in Frenchman, Stinking Water and Spring Creeks during the dry season at baseflow conditions. Measurements were taken both in straight and meandering stream channels during a 5 day period at similar flow conditions. Each test site comprised of at least three transects and each transect comprised of at least three K_v measurements. Linear, Gaussian, exponential and spherical variogram models were used with Kriging gridding method for the 10 sites. As a goodness-of-fit statistic for the variogram models, cross-validation results showed differences in the median absolute deviation and the standard deviation of the cross-validation residuals. Results show that using the geometric means of the 10 sites for gridding performs better than using either all the K_v values from the 93 permeameter tests or 10 K_v values from the middle transects and centre permeameters. Incorporating both the spatial variability and the uncertainty involved in the measurement at a reach segment can yield more accurate grid results that can be useful in calibrating K_v at watershed or sub-watershed scales in distributed hydrological models.     

A comparison of chemical and chemodynamical models

Representative results from a comparison of the chemical evolution of spherical collapse models without and with a intercloud medium are presented. The hot metal-rich gas distributes quickly the metals produced in supernovae throughout the galaxy, thus leading to a more homogeneous chemical evolution and flatter metallicity gradients in the gas and the stars. The stellar population is somewhat less concentrated towards the centre. The strong outflow results in a substantial loss of metals from the galaxy to its surroundings, and a lower effective yield in the galaxy. This revised version was published online in August 2006 with corrections to the Cover Date.     

Transient‐state groundwater flow in various geometries with wells at arbitrary positions: analytical methods

We have developed a method for analytically solving the porous medium flow equation in many different geometries for horizontal (two‐dimensional), homogeneous and isotropic aquifers containing impermeable boundaries and any number of pumping or injection wells located at arbitrary positions within the system. Solutions and results are presented for rectangular and circular aquifers but the method presented here is easily extendible to many geometries. Results are also presented for systems where constant head boundary conditions can be emulated internal to the aquifer boundary. Recommendations for extensions of the present work are briefly discussed. Copyright © 2003 John Wiley & Sons, Ltd.     

Chaos and elliptical galaxies

Recent results on chaos in triaxial galaxy models are reviewed. Central mass concentrations like those observed in early-type galaxies - either stellar cusps, or massive black holes — render most of the box orbits in a triaxial potential stochastic. Typical Liapunov times are 3–5 crossing times, and ensembles of stochastic orbits undergo mixing on timescales that are roughly an order of magnitude longer. The replacement of the regular orbits by stochastic orbits reduces the freedom to construct self-consistent equilibria, and strong triaxiality can be ruled out for galaxies with sufficiently high central mass concentrations.     

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.     

Tidal effects on temperature front in the Yellow Sea

Temperature front (TF) is one of the important features in the Yellow Sea, which forms in spring, thrives in summer, and fades in autumn as thermocline declines. TF intensity ⋎S _T ⋎ is defined to describe the distribution of TF. Based on the MASNUM wave-tide-circulation coupled model, temperature distribution in the Yellow Sea was simulated with and without tidal effects. Along 36°N, distribution of TF from the simulated results are compared with the observations, and a quantitative analysis is introduced to evaluate the tidal effects on the forming and maintaining processes of the TF. Tidal mixing and the circulation structure adapting to it are the main causes of the TF. Supported by the National Basic Research Program of China (No. G1999043809) and the National Science Foundation of China (No. 49736190).     

Generalized strain probing of constitutive models

Advanced material constitutive models are used to describe complex soil behaviour. These models are often used in the solution of boundary value problems under general loading conditions. Users and developers of constitutive models need to methodically investigate the represented soil response under a wide range of loading conditions. This paper presents a systematic procedure for probing constitutive models. A general incremental strain probe, 6D hyperspherical strain probe (HSP), is introduced to examine rate‐independent model response under all possible strain loading conditions. Two special cases of HSP, the true triaxial strain probe (TTSP) and the plane‐strain strain probe (PSSP), are used to generate 3‐D objects that represent model stress response to probing. The TTSP, PSSP and general HSP procedures are demonstrated using elasto‐plastic models. The objects resulting from the probing procedure readily highlight important model characteristics including anisotropy, yielding, hardening, softening and failure. The PSSP procedure is applied to a Neural Network (NN) based constitutive model. It shows that this probing is especially useful in understanding NN constitutive models, which do not contain explicit functions for yield surface, hardening, or anisotropy. Copyright © 2004 John Wiley & Sons, Ltd.     

Numerical modelling of permafrost in bedrock in northern Fennoscandia during the Holocene

The occurrence of permafrost in bedrock in northern Fennoscandia and its dependence on past and presently ongoing climatic variations was investigated with one- (1D) and two-dimensional (2D) numerical models by solving the transient heat conduction equation with latent heat effects included. The study area is characterized by discontinuous permafrost occurrences such as palsa mires and local mountain permafrost. The ground temperature changes during the Holocene were constructed using climatic proxy data. This variation was used as a forcing function at the ground surface in the calculations. Several versions of the present ground temperature were applied, resulting in different subsurface freezing–thawing conditions in the past depending on the assumed porosity and geothermal conditions.Our results suggest that in high altitude areas with a cold climate (present mean annual ground temperature between 0°C and −3°C), there may have been considerable variations in permafrost thickness (ranging from 0 to 150 m), as well as periods of no permafrost at all. The higher is the porosity of bedrock filled with ice, the stronger is the retarding effect of permafrost against climatic variations.Two-dimensional models including topographic effects with altitude-dependent ground temperatures and slope orientation and inclination dependent solar radiation were applied to a case of mountain permafrost in Ylläs, western Finnish Lapland, where bedrock permafrost is known to occur in boreholes to a depth of about 60 m. Modelling suggests complicated changes in permafrost thickness with time as well as contrasting situations on southern and northern slopes of the mountain.Extrapolating the climatic warming of the last 200 years to the end of the next century when the anticipated increase in the annual average air temperature is expected to be about 2 K indicates that the permafrost occurrences in bedrock in northern Fennoscandia would be thawing rapidly in low-porosity formations. However, already a porosity of 5% filled with ice would retard the thawing considerably.