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.     

Outflow channel sources, reactivation, and chaos formation, Xanthe Terra, Mars

The undulating, warped, and densely fractured surfaces of highland regions east of Valles Marineris (located north of the eastern Aureum Chaos, east of the Hydraotes Chaos, and south of the Hydaspis Chaos) resulted from extensional surface warping related to ground subsidence, caused when pressurized water confined in subterranean caverns was released to the surface. Water emanations formed crater lakes and resulted in channeling episodes involved in the excavation of Ares, Tiu, and Simud Valles of the eastern part of the circum-Chryse outflow channel system. Progressive surface subsidence and associated reduction of the subsurface cavernous volume, and/or episodes of magmatic-driven activity, led to increases of the hydrostatic pressure, resulting in reactivation of both catastrophic and non-catastrophic outflow activity. Ancient cratered highland and basin materials that underwent large-scale subsidence grade into densely fractured terrains. Collapse of rock materials in these regions resulted in the formation of chaotic terrains, which occur in and near the headwaters of the eastern circum-Chryse outflow channels. The deepest chaotic terrain in the Hydaspis Chaos region resulted from the collapse of pre-existing outflow channel floors. The release of volatiles and related collapse may have included water emanations not necessarily linked to catastrophic outflow. Basal warming related to dike intrusions, thermokarst activity involving wet sediments and/or dissected ice-enriched country rock, permafrost exposed to the atmosphere by extensional tectonism and channel incision, and/or the injection of water into porous floor material, may have enhanced outflow channel floor instability and subsequent collapse. In addition to the possible genetic linkage to outflow channel development dating back to at least the Late Noachian, clear disruption of impact craters with pristine ejecta blankets and rims, as well as preservation of fine tectonic fabrics, suggest that plateau subsidence and chaos formation may have continued well into the Amazonian Period. The geologic and paleohydrologic histories presented here have important implications, as new mechanisms for outflow channel formation and other fluvial activity are described, and new reactivation mechanisms are proposed for the origin of chaotic terrain as contributors to flooding. Detailed geomorphic analysis indicates that subterranean caverns may have been exposed during chaos formation, and thus chaotic terrains mark prime locations for future geologic, hydrologic, and possible astrobiologic exploration.     

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.