Efficient computation of volume fractions for multi-material cell complexes in a grid by slicing

The paper presents a novel slicing based method for computation of volume fractions in multi-material solids given as a B-rep whose faces are triangulated and shared by either one or two materials. Such objects occur naturally in geoscience applications and the said computation is necessary for property estimation problems and iterative forward modeling. Each facet in the model is cut by the planes delineating the given grid structure or grid cells. The method, instead of classifying the points or cells with respect to the solid, exploits the convexity of triangles and the simple axis-oriented disposition of the cutting surfaces to construct a novel intermediate space enumeration representation called slice-representation, from which both the cell containment test and the volume-fraction computation are done easily. Cartesian and cylindrical grids with uniform and non-uniform spacings have been dealt with in this paper. After slicing, each triangle contributes polygonal facets, with potential elliptical edges, to the grid cells through which it passes. The volume fractions of different materials in a grid cell that is in interaction with the material interfaces are obtained by accumulating the volume contributions computed from each facet in the grid cell. The method is fast, accurate, robust and memory efficient. Examples illustrating the method and performance are included in the paper.

Applicability of Statistical Learning Algorithms for Spatial Variability of Rock Depth

Two algorithms are outlined, each of which has interesting features for modeling of spatial variability of rock depth. In this paper, reduced level of rock at Bangalore, India, is arrived from the 652 boreholes data in the area covering 220 sqa <.km. Support vector machine (SVM) and relevance vector machine (RVM) have been utilized to predict the reduced level of rock in the subsurface of Bangalore and to study the spatial variability of the rock depth. The support vector machine (SVM) that is firmly based on the theory of statistical learning theory uses regression technique by introducing epsilon-insensitive loss function has been adopted. RVM is a probabilistic model similar to the widespread SVM, but where the training takes place in a Bayesian framework. Prediction results show the ability of learning machine to build accurate models for spatial variability of rock depth with strong predictive capabilities. The paper also highlights the capability ofRVM over the SVM model.

Dust aerosols over India and adjacent continents retrieved using METEOSAT infrared radiance Part II: quantification of wind dependence and estimation of radiative forcing

Long-range transport of continental dust makes these particles a significant constituent even at locations far from their sources. It is important to study the temporal variations in dust loading over desert regions and the role of meteorology, in order to assess its radiative impact. In this paper, infrared radiance (10.5-12.5 mu m), acquired by the METEOSAT-5 satellite (similar to 5-km resolution) during 1999 and 2003 was used to quantify wind dependence of dust aerosols and to estimate the radiative forcing. Our analysis shows that the frequency of occurrence of dust events was higher during 2003 compared to 1999. Since the dust production function depends mainly on the surface wind speed over regions which are dry and without vegetation, the role of surface wind on IDDI was examined in detail. It was found that an increase of IDDI with wind speed was nearly linear and the rate of increase in IDDI with surface wind was higher during 2003 compared to 1999. It was also observed that over the Indian desert, when wind speed was the highest during monsoon months (June to August), the dust production rate was lower because of higher soil moisture (due to monsoon rainfall). Over the Arabian deserts, when the wind speed is the highest during June to August, the dust production rate is also highest, as soil moisture is lowest during this season. Even though nothing can be said precisely on the reason why 2003 had a greater number of dust events, examination of monthly mean soil moisture at source regions indicates that the occurrence of high winds simultaneous with high soil moisture could be the reason for the decreased dust production efficiency in 1999. It appears that the deserts of Northwest India are more efficient dust sources compared to the deserts of Saudi Arabia and Northeast Africa (excluding Sahara). The radiative impact of dust over various source regions is estimated, and the regionally and annually averaged top of the atmosphere dust radiative forcing (short wave, clear-sky and over land) over the entire study region (0-35 degrees N; 30 degrees-100 degrees E) was in the range of -0.9 to +4.5 W m(-2). The corresponding values at the surface were in the range of -10 to -25 W m(-2). Our studies demonstrate that neglecting the diurnal variation of dust can cause errors in the estimation of long wave dust forcing by as much as 50 to 100%, and nighttime retrieval of dust can significantly reduce the uncertainties. A method to retrieve dust aerosols during nighttime is proposed. The regionally and annually averaged long wave dust radiative forcing was +3.4 +/- 1.6 W m(-2).

Changes in seasonal snow cover in Hindu Kush-Himalayan region

The changes in seasonal snow covered area in the Hindu Kush-Himalayan (HKH) region have been examined using Moderate – resolution Imaging Spectroradiometer (MODIS) 8-day standard snow products. The average snow covered area of the HKH region based on satellite data from 2000 to 2010 is 0.76 million km2 which is 18.23% of the total geographical area of the region. The linear trend in annual snow cover from 2000 to 2010 is −1.25±1.13%. This is in consistent with earlier reported decline of the decade from 1990 to 2001. A similar trend for western, central and eastern HKH region is 8.55±1.70%, +1.66% ± 2.26% and 0.82±2.50%, respectively. The snow covered area in spring for HKH region indicates a declining trend (−1.04±0.97%). The amount of annual snowfall is correlated with annual seasonal snow cover for the western Himalaya, indicating that changes in snow cover are primarily due to interannual variations in circulation patterns. Snow cover trends over a decade were also found to vary across seasonally and the region. Snow cover trends for western HKH are positive for all seasons. In central HKH the trend is positive (+15.53±5.69%) in autumn and negative (−03.68±3.01) in winter. In eastern HKH the trend is positive in summer (+3.35±1.62%) and autumn (+7.74±5.84%). The eastern and western region of HKH has an increasing trend of 10% to 12%, while the central region has a declining trend of 12% to 14% in the decade between 2000 and 2010. Snow cover depletion curve plotted for the hydrological year 2000–2001 reveal peaks in the month of February with subsidiary peaks observed in November and December in all three regions of the HKH.

Relative effect of slope and equilibrium line altitude on the retreat of Himalayan glaciers

It has been observed that a majority of glaciers in the Himalayas have been retreating. In this paper, we show that there are two major factors which control the advance/retreat of the Himalayan glaciers. They are the slope of the glacier and changes in the equilibrium line altitude. While it is well known, that these factors are important, we propose a new way of combining them and use it to predict retreat. The functional form of this model has been derived from numerical simulations using an ice-flow code. The model has been successfully applied to the movement of eight Himalayan glaciers during the past 25 years. It explains why the Gangotri glacier is retreating while Zemu of nearly the same length is stationary, even if they are subject to similar environmental changes. The model has also been applied to a larger set of glaciers in the Parbati basin, for which retreat based on satellite data is available, though over a shorter time period.

Influence of anthropogenic contamination on groundwater chemistry in Mulbagal town, Kolar District, India

Geochemical studies are performed to examine the impact of leachate infiltration from on-site sewage disposal systems on the groundwater chemistry in Mulbagal town, Kolar District, Karnataka State, India. The leachate infiltration imposed nitrate concentrations ranging from 4 mg/L to 388 mg/L in the groundwater samples; it was observed that 79% of the samples exhibited nitrate concentrations in excess of drinking water permissible limit (45 mg/L). The average (of 43 measurements) E. coli levels in the groundwater samples corresponded to 189 MPN/100 mL and 55% of the samples tested exhibit pathogen contamination. Results also showed that the groundwater in the study area is characterized by acidic pH, large calcium + magnesium ion and Na/Cl ratios of < unity causing majority of the ground water samples to classify as Ca-Mg-Cl type and Na-Cl type. Saturation index (SI) computation using Visual MINTEQ program showed that the groundwater samples are under-saturated with respect to calcite. The theoretical SI values (of calcite) however suggested that the groundwater samples ought to be over-saturated with calcite. Under-saturation of the groundwater samples with calcite is attributed to increased dissolution of the mineral in the acidic environment of the groundwater.

Coupling geomorphology parameters with lithology for micro-level groundwater resource assessment: a case study from semi-arid hard rock terrain in Tamil Nadu, India

The standard procedure of groundwater resource estimation in India till date is based on the specific yield parameters of each rock type (lithology) derived through pumping test analysis. Using the change in groundwater level, specific yield, and area of influence, groundwater storage change could be estimated. However, terrain conditions in the form of geomorphological variations have an important bearing on the net groundwater recharge. In this study, an attempt was made to use both lithology and geomorphology as input variables to estimate the recharge from different sources in each lithology unit influenced by the geomorphic conditions (lith-geom), season wise separately. The study provided a methodological approach for an evaluation of groundwater in a semi-arid hard rock terrain in Tirunelveli, Tamil Nadu, India. While characterizing the gneissic rock, it was found that the geomorphologic variations in the gneissic rock due to weathering and deposition behaved differently with respect to aquifer recharge. The three different geomorphic units identified in gneissic rock (pediplain shallow weathered (PPS), pediplain moderate weathered (PPM), and buried pediplain moderate (BPM)) showed a significant variation in recharge conditions among themselves. It was found from the study that Peninsular gneiss gives a net recharge value of 0.13 m/year/unit area when considered as a single unit w.r.t. lithology, whereas the same area considered with lith-geom classes gives recharge values between 0.1 and 0.41 m/year presenting a different assessment. It is also found from this study that the stage of development (SOD) for each lith-geom unit in Peninsular gneiss varies from 168 to 230 %, whereas the SOD is 223 % for the lithology as a single unit.

Seismic intensity map of South India for estimated future earthquakes

In this study, an attempt has been made to prepare the seismic intensity map for south India considering the probable earthquakes in the region. Anbazhagan et al. (Nat Hazards 60:1325-1345, 2012) have identified eight probable future earthquake zones in south India based on rupture-based seismic hazard analysis. Anbazhagan et al. (Eng Geol 171:81-95, 2014) has estimated the maximum future earthquake magnitude at these eight zones using regional rupture character. In this study, the whole south India is divided into several grids of size 1(o) x 1(o) and the intensity at each grid point is calculated using the regional intensity model for the maximum earthquake magnitude at each of the eight zones. The intensity due to earthquakes at these zones is mapped and thus eight seismic intensity maps are prepared. The final seismic intensity map of south India is obtained by considering the maximum intensity at each grid point due to the estimated earthquakes. By looking at the seismic intensity map, one can expect slight to heavy damage due to the probable earthquake magnitudes. Heavy damage may happen close to the probable earthquake zones.