The exponential decline in saturated hydraulic conductivity with depth: a novel method for exploring its effect on water flow paths and transit time distribution

The strong vertical gradient in soil and subsoil saturated hydraulic conductivity is characteristic feature of the hydrology of catchments. Despite the potential importance of these strong gradients, they have proven difficult to model using robust physically based schemes. This has hampered the testing of hypotheses about the implications of such vertical gradients for subsurface flow paths, residence times and transit time distribution. Here we present a general semi‐analytical solution for the simulation of 2D steady‐state saturated‐unsaturated flow in hillslopes with saturated hydraulic conductivity that declines exponentially with depth. The grid‐free solution satisfies mass balance exactly over the entire saturated and unsaturated zones. The new method provides continuous solutions for head, flow and velocity in both saturated and unsaturated zones without any interpolation process as is common in discrete numerical schemes. This solution efficiently generates flow pathlines and transit time distributions in hillslopes with the assumption of depth‐varying saturated hydraulic conductivity. The model outputs reveal the pronounced effect that changing the strength of the exponential decline in saturated hydraulic conductivity has on the flow pathlines, residence time and transit time distribution. This new steady‐state model may be useful to others for posing hypotheses about how different depth functions for hydraulic conductivity influence catchment hydrological response. Copyright © 2016 John Wiley & Sons, Ltd.     

Countermeasures for breakwater foundation subjected to foreshocks and main shock of earthquake loading

Coastal protective structures, such as composite breakwaters, are generally vulnerable to earthquake. It was observed that breakwaters damage mainly due to failure of their foundations. However, the seismically induced failure process of breakwater foundation has not been well understood. This study describes failure mechanism of breakwater foundation as well as a newly developed reinforcing model for breakwater foundation that can render resiliency to breakwater against earthquake-related disasters. Steel sheet piles and gabions were used as reinforcing materials for foundation. The experimental program consisted of a series of shaking table tests for conventional and reinforced foundation of breakwater. Numerical analyses were conducted using finite difference method, and it was observed that the numerical models were capable to elucidate the seismic behavior of soil–reinforcement–breakwater system. This paper presents an overview of the results of experimental and numerical studies of the seismic response of breakwater foundation. Overall, the results of these studies show the effectiveness of the reinforced foundation in mitigating the earthquake-induced damage to the breakwater. Moreover, numerical simulation was used for parametric study to determine the effect of different embedment depths of sheet piles on the performance of breakwater foundation subjected to seismic loading.     

Fracture Permeability Alteration due to Chemical and Mechanical Processes: A Coupled High-Resolution Model

Reactive fluid-flow experiments in fractures subjected to normal stress suggest the potential for either increased or decreased permeability resulting from fracture-surface dissolution. We present a computational model that couples mechanical deformation and chemical alteration of fractures subjected to constant normal stress and reactive fluid flow. The model explicitly represents micro-scale roughness of the fracture surfaces and calculates elastic deformation of the rough surfaces using a semi-analytical approach that ensures the surfaces remain in static equilibrium. A depth-averaged reactive transport model calculates chemical alteration of the surfaces, which leads to alteration of the contacting fracture surfaces. The mechanical deformation and chemical alteration calculations are explicitly coupled, which is justified by the disparate timescales required for equilibration of mechanical stresses and reactive transport processes. An idealized analytical representation of dissolution from a single contacting asperity shows that under reaction-limited conditions, contacting asperities can dissolve faster than the open regions of the fracture. Computational simulations in fractures with hundreds of contacting asperities show that the transition from transport-limited conditions (low flow rates) to reaction-rate-limited conditions (high flow rates) causes a shift from monotonically increasing permeability to a more complicated process in which permeability initially decreases and then increases as contacting asperities begin to dissolve. These results are qualitatively consistent with a number of experimental observations reported in the literature and suggest the potential importance of the relative magnitude of mass transport and reaction kinetics on the evolution of fracture permeability in fractures subjected to combined normal stress and reactive fluid flow.     

The formation of an eccentric gap in a gas disc by a planet in an eccentric orbit

We investigate the effect of a planet on an eccentric orbit on a two-dimensional low-mass gaseous disc. At a planet eccentricity above the planet's Hill radius divided by its semimajor axis, we find that the disc morphology differs from that exhibited by a disc containing a planet in a circular orbit. An eccentric gap is created with eccentricity that can exceed the planet's eccentricity and precesses with respect to the planet's orbit. We find that a more massive planet is required to open a gap when the planet is on an eccentric orbit. We attribute this behaviour to spiral density waves excited at corotation resonances by the eccentric planet. These act to increase the disc's eccentricity and exert a torque opposite in sign to that exerted by the Lindblad resonances. The reduced torque makes it more difficult for waves driven by the planet to overcome viscous inflow in the disc.     

On the direction of rotation of the spirals in galaxies

The detection of trailing arms in twenty spiral galaxies, completed by the 60s, allowed one to suggest the overwhelming prevalence of this sense of rotation. Pasha and Dambis (1982) have recently determined the direction of rotation in a further 54 galaxies and found two leading spirals, NGC 4490 and NGC 5395. The current sample of galaxies with an established sense of rotation contains 79 trailing and two leading systems. The two spirals found to be leading belong to M51 type galaxies.In this paper we present a more detailed investigation of the system NGC 4490/85 which confirms that the arms in NGC 4490 are leading. Statistics are invoked to estimate a probable percentageX _p and an upper limitX _up of galaxies with leading arms. It gives, for 79 trailing and two leading spirals, the valuesX _p 3.5% andX _up 8.2%.     

Satellite image based quantification of invasion and patch dynamics of mesquite (Prosopis juliflora) in Great Rann of Kachchh,Kachchh Biosphere Reserve,Gujarat, India

The invasion of alien species is a significant threat to global biodiversity and the top driver of climate change. The present study was conducted in the Great Rann of Kachchh, part of Kachchh Biosphere Reserve, Gujarat, India, which has been severely affected by invasion of Prosopis juliflora. The invasive weed infestation has been identified using multi-temporal remote sensing datasets of 1977, 1990, 1999, 2005 and 2011. Spatial analyses of the transition matrix, extent of invasive colonies, patchiness, coalescence and rate of spread were carried out. During the study period of three and half decades, almost 295 km² of the natural land cover was converted into Prosopis cover. This study has shown an increment of 42.9% of area under Prosopis cover in the Great Rann of Kachchh, part of the Kachchh Biosphere Reserve during 1977 to 2011. Spatial analysis indicates high occupancy of Prosopis cover with most of the invasion (95.9%) occurring in the grasslands and only 4.1% in other land cover types. The process of Prosopis invasion shows high patch initiation, followed by coalescence, indicating aggressive colonization of species. The number of patches within an area of < 1 km² increased from 1977 to 2011, indicating the formation of new Prosopis habitats by replacing the grasslands. The largest patch of Prosopis cover increased from 144 km² in 1977 to 430 km² in 2011. The estimated mean patch size was 7.8 km² in 1977. The mean patch size was largest during 2011, i.e., 9 km². The annual spread rate for Prosopis has been estimated as 2.1% during 2005–2011. The present work has investigated the long term changes in Prosopis cover in the Great Rann of Kachchh, part of Kachchh Biosphere Reserve. The spatial database generated will be useful in preparing strategies for the management of Prosopis juliflora.     

Assessment of PM<Subscript>2.5</Subscript> chemical compositions in Delhi: primary vs secondary emissions and contribution to light extinction coefficient and visibility degradation

Haze-fog conditions over northern India are associated with visibility degradation and severe attenuation of solar radiation by airborne particles with various chemical compositions. PM_2.5 samples have been collected in Delhi, India from December 2011 to November 2012 and analyzed for carbonaceous and inorganic species. PM₁₀ measurements were made simultaneously such that PM_10–2.5 could be estimated by difference. This study analyzes the temporal variation of PM_2.5 and carbonaceous particles (CP), focusing on identification of the primary and secondary aerosol emissions, estimations of light extinction coefficient (b_ext) and the contributions by the major PM_2.5 chemical components. The annual mean concentrations of PM_2.5, organic carbon (OC), elemental carbon (EC) and PM_10–2.5 were found to be 153.6 ± 59.8, 33.5 ± 15.9, 6.9 ± 3.9 and 91.1 ± 99.9 μg m^?3, respectively. Total CP, secondary organic aerosols and major anions (e.g., SO₄ ^2? and NO₃ ^?) maximize during the post-monsoon and winter due to fossil fuel combustion and biomass burning. PM_10–2.5 is more abundant during the pre-monsoon and post-monsoon. The OC/EC varies from 2.45 to 9.26 (mean of 5.18 ± 1.47), indicating the influence of multiple combustion sources. The b_ext exhibits highest values (910 ± 280 and 1221 ± 371 Mm^?1) in post-monsoon and winter and lowest in monsoon (363 ± 110 and 457 ± 133 Mm^?1) as estimated via the original and revised IMPROVE algorithms, respectively. Organic matter (OM =1.6 × OC) accounts for ~39 % and ~48 % of the b_ext, followed by (NH₄)₂SO₄ (~21 % and ~24 %) and EC (~13 % and ~10 %), according to the original and revised algorithms, respectively. The b_ext estimates via the two IMPROVE versions are highly correlated (R² = 0.95, root mean square error = 38 % and mean bias error = 28 %) and are strongly related to visibility impairment (r = ?0.72), mostly associated with anthropogenic rather than natural PM contributions. Therefore, reduction of CP and precursor gas emissions represents an urgent opportunity for air quality improvement across Delhi.