Evaluation of spatial and spatiotemporal estimation methods in simulation of precipitation variability patterns

Identification of spatial and spatiotemporal precipitation variations plays an important role in different hydrological applications such as missing data estimation. In this paper, the results of Bayesian maximum entropy (BME) and ordinary kriging (OK) are compared for modeling spatial and spatiotemporal variations of annual precipitation with and without incorporating elevation variations. The study area of this research is Namak Lake watershed located in the central part of Iran with an area of approximately 90,000 km². The BME and OK methods have been used to model the spatial and spatiotemporal variations of precipitation in this watershed, and their performances have been evaluated using cross-validation statistics. The results of the case study have shown the superiority of BME over OK in both spatial and spatiotemporal modes. The results have shown that BME estimates are less biased and more accurate than OK. The improvements in the BME estimates are mostly related to incorporating hard and soft data in the estimation process, which resulted in more detailed and reliable results. Estimation error variance for BME results is less than OK estimations in the study area in both spatial and spatiotemporal modes.     

Hydrogeological aspects of groundwater drainage of the urban areas in Kuwait City

Residential areas in Kuwait City have witnessed a dramatic rise in subsurface water tables over the last three decades. This water rise phenomenon is attributed mainly to over irrigation practices of private gardens along with leakage from domestic and sewage networks. This paper presents a comprehensive study for urban drainage in two selected areas representing the two hydrogeological settings encountered in Kuwait City. In the first area, a vertical drainage scheme was applied successfully over an area of 1 km². The system has been under continuous operation and monitoring for more than 4 years without problems, providing a permanent solution for the water rise problem in this area. The hydrogeological system has approached steady state conditions and the water levels have dropped to about 3·5 m below the ground surface. In the second area a dual drainage scheme, composing of horizontal and vertical elements, is proposed. Horizontal elements are suggested in the areas where the deep groundwater contains hazardous gases that may pose environmental problems. The proposed drainage scheme in the second area has not yet been implemented. Field tests were conducted to assess the aquifer parameters in both areas and a numerical model has been developed to predict the long‐term response of the hydrogeological system in the two areas under consideration. Copyright © 2001 John Wiley & Sons, Ltd.     

Magnetized Kelvin-Helmholtz instability in the presence of a radiation field

The purpose of this study is to analyze the dynamical role of a radiation field on the growth rate of the unstable Kelvin-Helmholtz (KH) perturbations. As a first step toward this purpose, the analyze is done in a general way, irrespective of applying the model to a specific astronomical system. The transition zone between the two layers of the fluid is ignored. Then, we perform a linear analysis and by imposing suitable boundary conditions and considering a radiation field, we obtain appropriate dispersion relation. Unstable modes are studied by solving the dispersion equation numerically, and then growth rates of them are obtained. By analyzing our dispersion relation, we show that for a wide range of the input parameters, the radiation field has a destabilizing effect on KH instability. In eruptions of the galaxies or supermassive stars, the radiation field is dynamically important and because of the enhanced KH growth rates in the presence of the radiation; these eruptions can inject more momentum and energy into their environment and excite more turbulent motions.     

Thin accretion disc with a corona in a central magnetic field

We study the steady-state structure of an accretion disc with a corona surrounding a central, rotating, magnetized star. We assume that the magneto-rotational instability is the dominant mechanism of angular momentum transport inside the disc and is responsible for producing magnetic tubes above the disc. In our model, a fraction of the dissipated energy inside the disc is transported to the corona via these magnetic tubes. This energy exchange from the disc to the corona which depends on the disc physical properties is modified because of the magnetic interaction between the stellar magnetic field and the accretion disc. According to our fully analytical solutions for such a system, the existence of a corona not only increases the surface density but reduces the temperature of the accretion disc. Also, the presence of a corona enhances the ratio of gas pressure to the total pressure. Our solutions show that when the strength of the magnetic field of the central neutron star is large or the star is rotating fast enough, profiles of the physical variables of the disc significantly modify due to the existence of a corona.     

Hot accretion with outflow and thermal conduction

We present self-similar solutions for advection-dominated accretion flows with thermal conduction in the presence of outflows. Possible effects of outflows on the accretion flow are parametrized and a saturated form of thermal conduction, as is appropriate for the weakly-collisional regime of interest, is included in our model. While the cooling effect of outflows is noticeable, thermal conduction provides an extra heating source. In comparison to accretion flows without winds, we show that the disc rotates faster and becomes cooler because of the angular momentum and energy flux which are taking away by the winds. But thermal conduction opposes the effects of winds and not only decreases the rotational velocity, but increases the temperature. However, reduction of the surface density and the enhanced accretion velocity are amplified by both of the winds and the thermal conduction. We find that for stronger outflows, a higher level of saturated thermal conduction is needed to significantly modify the physical profiles of the accretion flow.     

The role of Kelvin–Helmholtz instability in dusty and partially ionized outflows

We investigate the linear theory of Kelvin–Helmholtz instability at the interface between a partially ionized dusty outflow and the ambient material analytically. We model the interaction as a multifluid system in a planar geometry. The unstable modes are independent from the charge polarity of the dust particles. Although our results show a stabilizing effect for charged dust particles, the growth time-scale of the growing modes gradually becomes independent of the mass or charge of the dust particles when the magnetic-field strength increases. We show that growth time-scale decreases with increasing the magnetic field. Also, as the mass of the dust particles increases, the growth time-scale of the unstable mode increases.     

How deep can forest vegetation cover extend their hydrological reinforcing contribution?

An experimental campaign was set up to quantify the contribution of evapotranspiration fluxes on hillslope hydrology and stability for different forest vegetation cover types. Three adjacent hillslopes, respectively, covered by hardwood, softwood, and grass were instrumented with nine access tubes each to monitor soil water dynamics at the three depths of 30, 60, and 100 cm, using a PR2/6 profile probe (Delta‐T Devices Ltd) for about 6 months including wet periods. Soil was drier under softwood and wetter under grass at all the three depths during most of the monitoring period. Matric suction derived via the soil moisture measurements was more responsive to changes in the atmospheric conditions and also recovered faster at the 30 cm depth. Results showed no significant differences between mean matric suction under hardwood (101.6 kPa) with that under either softwood or grass cover. However, a significant difference was found between mean matric suction under softwood (137.5 kPa) and grass (84.3 kPa). Results revealed that, during the wettest period, the hydrological effects from all three vegetation covers were substantial at the 30 cm depth, whereas the contribution from grass cover at 60 cm (2.0 kPa) and 100 cm (1.1 kPa) depths and from hardwood trees at 100 cm depth (1.2 kPa) was negligible. It is surmised that potential instability would have occurred at these larger depths along hillslopes where shallow hillslope failures are most likely to occur in the region. The hydrological effects from softwood trees, 8.1 and 3.9 kPa, were significant as the corresponding factor of safety values showed stable conditions at both depths of 60 and 100 cm, respectively. Therefore, the considerable hydrological reinforcing effects from softwood trees to the 100 cm depth suggest that a hillslope stability analysis would show that hillslopes with softwood trees will be stable even during the wet season.     

Soil water dynamics under different forest vegetation cover: Implications for hillslope stability

Though it is well known that vegetation affects the water balance of soils through canopy interception and evapotranspiration, its hydrological contribution to soil hydrology and stability is yet to be fully quantified. To improve understanding of this hydrological process, soil water dynamics have been monitored at three adjacent hillslopes with different vegetation covers (deciduous tree cover, coniferous tree cover, and grass cover), for nine months from December 2014 to September 2015. The monitored soil moisture values were translated into soil matric suction (SMS) values to facilitate the analysis of hillslope stability. Our observations showed significant seasonal variations in SMS for each vegetation cover condition. However, a significant difference between different vegetation covers was only evident during the winter season where the mean SMS under coniferous tree cover (83.6 kPa) was significantly greater than that under grass cover (41 kPa). The hydrological reinforcing contribution due to matric suction was highest for the hillslope with coniferous tree cover, while the hillslope with deciduous tree cover was second and the hillslope with grass cover was third. The greatest contributions for all cover types were during the summer season. During the winter season, the wettest period of the monitoring study, the additional hydrological reinforcing contributions provided by the deciduous tree cover (1.5 to 6.5 kPa) or the grass cover (0.9 to 5.4 kPa) were insufficient to avoid potential slope failure conditions. However, the additional hydrological reinforcing contribution from the coniferous tree cover (5.8 to 10.4 kPa) was sufficient to provide potentially stable hillslope conditions during the winter season. Our study clearly suggests that during the winter season the hydrological effects from both deciduous tree and grass covers are insufficient to promote slope stability, while the hydrological reinforcing effects from the coniferous tree cover are sufficient even during the winter season. Copyright © 2018 John Wiley & Sons, Ltd.     

Random vibration of nonlinear structures with stiffness and strength deterioration by modified tail equivalent linearization method

In this research the tail equivalent linearization method(TELM) has been extended to study structures with degrading materials. The responses of such structures to excitations are non-stationary, even if the excitations are stationary. Non-stationary behavior of the system cannot be considered by conventional TELM. Applying the conventional TELM, the only distinction in the design point excitation for two stationary excitations with different durations is in the addition of a zero value part at the beginning of the design point of the longer excitation. This means that the failure probability is the same for the non-stationary systems under excitations with different durations. Therefore, this solution cannot be correct. In this study, in using TELM for systems with degrading materials, hysteretic energy is replaced by average hysteretic energy, calculated by averaging the obtained hysteretic energy of the structure subjected to a few random sample load realization. In this way, the degradation parameters under design point coincide with those under sample load realizations. Since the average of the hysteretic energy is converges very fast, the modified TELM only requires about tens to hundreds solutions of the response in addition to the ordinary calculations of conventional TELM.