Modelling the effect of changing precipitation inputs on deep soil water utilization

Forests in the Southeastern United States are predicted to experience future changes in seasonal patterns of precipitation inputs as well as more variable precipitation events. These climate change‐induced alterations could increase drought and lower soil water availability. Drought could alter rooting patterns and increase the importance of deep roots that access subsurface water resources. To address plant response to drought in both deep rooting and soil water utilization as well as soil drainage, we utilize a throughfall reduction experiment in a loblolly pine plantation of the Southeastern United States to calibrate and validate a hydrological model. The model was accurately calibrated against field measured soil moisture data under ambient rainfall and validated using 30% throughfall reduction data. Using this model, we then tested these scenarios: (a) evenly reduced precipitation; (b) less precipitation in summer, more in winter; (c) same total amount of precipitation with less frequent but heavier storms; and (d) shallower rooting depth under the above 3 scenarios. When less precipitation was received, drainage decreased proportionally much faster than evapotranspiration implying plants will acquire water first to the detriment of drainage. When precipitation was reduced by more than 30%, plants relied on stored soil water to satisfy evapotranspiration suggesting 30% may be a threshold that if sustained over the long term would deplete plant available soil water. Under the third scenario, evapotranspiration and drainage decreased, whereas surface run‐off increased. Changes in root biomass measured before and 4 years after the throughfall reduction experiment were not detected among treatments. Model simulations, however, indicated gains in evapotranspiration with deeper roots under evenly reduced precipitation and seasonal precipitation redistribution scenarios but not when precipitation frequency was adjusted. Deep soil and deep rooting can provide an important buffer capacity when precipitation alone cannot satisfy the evapotranspirational demand of forests. How this buffering capacity will persist in the face of changing precipitation inputs, however, will depend less on seasonal redistribution than on the magnitude of reductions and changes in rainfall frequency.     

Three-dimensional finite element simulation of fracture propagation in rock specimens with pre-existing fissure(s) under compression and their strength analysis

A FORTRAN program, consistent with the commercially available finite element (FE) code ABAQUS, is developed based on a three-dimensional (3D) linear elastic brittle damage constitutive model with two damage criteria. To consider the heterogeneity of rock, the developed FORTRAN program is used to set the stiffness and strength properties of each element of the FE model following a Weibull distribution function. The reliability of the program is assessed against available experimental results for granite cylindrical specimens with a throughgoing, flat and inclined fissure. The calibration procedure of the material parameters is explained in detail, and it is shown that the compressive to tensile strength ratio can have a substantial influence on the failure response of the specimens. Numerical simulations are conducted for models with different levels of heterogeneity. The results show a smaller load bearing capacity for models with less homogeneity, representing gradual coalescence of fully damaged elements forming throughout the models during loading. The maximum load bearing capacity is studied for various combinations of inclination angles of two centrally aligned, throughgoing and flat fissures of equal length embedded in cylindrical models under uniaxial and multiaxial loading conditions. The key role of the compressive to tensile strength ratio is highlighted by repeating certain simulations with a lower compressive to tensile strength ratio. It is proven that the peak loads of the rock models with sufficiently small compressive to tensile strength ratios containing two throughgoing fissures of equal length are similar, provided that the minimum inclination angles of the models are the same. The results are presented and discussed with respect to the existing experimental findings in the literature, suggesting that the numerical model applied in this study can provide useful insight into the failure behaviour of rock-like materials.     

Effects of the geometry of two‐dimensional fractures on their hydraulic aperture and on the validity of the local cubic law

Flow through rough fractures is investigated numerically in order to assess the validity of the local cubic law for different fracture geometries. Two‐dimensional channels with sinusoidal walls having different geometrical properties defined by the aperture, the amplitude, and the wavelength of the walls' corrugations, the corrugations asymmetry, and the phase shift between the two walls are considered to represent different fracture geometries. First, it is analytically shown that the hydraulic aperture clearly deviates from the mean aperture when the walls' roughness, the phase shift, and/or the asymmetry between the fracture walls are relatively high. The continuity and the Navier–Stokes equations are then solved by means of the finite element method and the numerical solutions compared to the theoretical predictions of the local cubic law. Reynolds numbers ranging from 0.066 to 66.66 are investigated so as to focus more particularly on the effect of flow inertial effects on the validity of the local cubic law. For low Reynolds number, typically less than 15, the local cubic law properly describes the fracture flow, especially when the fracture walls have small corrugation amplitudes. For Reynolds numbers higher than 15, the local cubic law is valid under the conditions that the fracture presents a low aspect ratio, small corrugation amplitudes, and a moderate phase lag between its walls.     

Implementation of a coupled plastic damage distinct lattice spring model for dynamic crack propagation in geomaterials

This paper studies dynamic crack propagation by employing the distinct lattice spring model (DLSM) and 3‐dimensional (3D) printing technique. A damage‐plasticity model was developed and implemented in a 2D DLSM. Applicability of the damage‐plasticity DLSM was verified against analytical elastic solutions and experimental results for crack propagation. As a physical analogy, dynamic fracturing tests were conducted on 3D printed specimens using the split Hopkinson pressure bar. The dynamic stress intensity factors were recorded, and crack paths were captured by a high‐speed camera. A parametric study was conducted to find the influences of the parameters on cracking behaviors, including initial and peak fracture toughness, crack speed, and crack patterns. Finally, selection of parameters for the damage‐plasticity model was determined through the comparison of numerical predictions and the experimentally observed cracking features.     

Overview of In - Situ Biodegradation and Enhancement

Microbial degradation technologies have been developed to restore ground water quality in aquifers polluted by organic contaminants effectively in recent years. However, in course of the degradation, the formation of biofilms in ground water remediation technology can be detrimental to the effectiveness of a ground water remediation project. Several alternatives are available to a remedial design engineer, such as Permeable Reactive Barriers (PRBs) and in -situ bioremediation, Hydrogen Releasing Compounds (HRCs) barrier, Oxygen Releasing Compounds (ORCs) barrier etc. which are efficient and cost- effective technologies. Excessive biomass formation renders a barrier ineffective in degrading the contaminants, Efforts are made to develop kinetics models which accurately determine bio - fouling and bio - filn formation and to control excessive biomass formation.     

Discovery of 13 New Variable Stars in the Field of the Open Cluster NGC 2168 （M35）

A wide-field time-series CCD photometric survey of variable stars in the field of the open cluster NGC 2168 was carried out using the BATC Schmidt telescope. In total 13 new variable stars are discovered with three W UMa systems, one EA type and two EB type eclipsing binaries (one of them could be a W UMa system), and seven pulsating stars including three candidates of δScuti stars.     

Multi-wavelength Diagnostics of Bombardment by Non-thermal Particles in Solar Flares

Chromospheric lines, including Ha, Lyα, Lyβand CaⅡK, CaⅡ8542, are systemically and quantitatively investigated with respect to the non-thermal excitation and ionization due to particle beam bombardment for a series of solar semi-empirical atmospheric models. As a result we propose to use the contrast in the integrated intensity of hydrogen lines to estimate the total energy flux of the bombarding beam during the solar flare impulsive phase. Partial frequency redistribution is considered in the Lyαline calculation and a smaller intensity enhancement in the Ha line-centers is found than in the previous results of Fang et al.     

Influence of the Choice of Core-Envelope Transition Point on the Binary Merger of Two Main-sequence Components

We have studied the influence of different choices of core-envelope transition point on the final merger of contact binaries with two main-sequence components. A binary of 1.00 + 0.90M⊙ with an initial orbital period of 0.35d is examined. The mass fraction of the primary mixed with the matter of the secondary, qmix, determined by the chosen core-envelope transition point, ranges from 0.04 to 1.00 in our analysis. If as qmix< 0.8, none of the     

The Difference between the α-disks of Seyfert 1 Galaxiesand Quasars

In a previous paper, it was suggested that contamination of the nuclear luminosity by the host galaxy plays an important role in determining the parameters of the standard a disk of AGNs. Using the nuclear absolute B band magnitude instead of the total absolute B band magnitude, we have recalculated the central black hole masses, accretion rates and disk inclinations for 20 Seyfert 1 galaxies and 17 Palomar-Green (PG) quasars. It is found that a small value of a is needed for the Seyfert 1 galaxies than for the PG quasars. This difference in a possibly leads to the different properties of Seyfert 1 galaxies and quasars. Furthermore, we find most of the objects in this sample are not accreting at super-Eddington rates if we adopt the nuclear optical luminosity in our calculation.