A parallel computing approach to viewshed analysis of large terrain data using graphics processing units

Viewshed analysis, often supported by geographic information system, is widely used in many application domains. However, as terrain data continue to become increasingly large and available at high resolutions, data-intensive viewshed analysis poses significant computational challenges. General-purpose computation on graphics processing units (GPUs) provides a promising means to address such challenges. This article describes a parallel computing approach to data-intensive viewshed analysis of large terrain data using GPUs. Our approach exploits the high-bandwidth memory of GPUs and the parallelism of massive spatial data to enable memory-intensive and computation-intensive tasks while central processing units are used to achieve efficient input/output (I/O) management. Furthermore, a two-level spatial domain decomposition strategy has been developed to mitigate a performance bottleneck caused by data transfer in the memory hierarchy of GPU-based architecture. Computational experiments were designed to evaluate computational performance of the approach. The experiments demonstrate significant performance improvement over a well-known sequential computing method, and an enhanced ability of analyzing sizable datasets that the sequential computing method cannot handle.     

The structure of rapidly rotating relaxed globular clusters

A model of rapidly rotating globular clusters is constructed assuming uniform angular velocity and a truncated Maxwellian distribution in the stellar velocities. Since the first-order theory developed by Woolley and Dickens becomes inaccurate for rapid rotation, a small-mass envelope is fitted to their models by the method of Monaghan and Roxburgh. A comparison is made of the critical values derived by the two methods.     

3D Subsurface Characterization of Banded Iron Formation Mineralization using Large-Scale Gravity Data: A Case Study in Parts of Bharatpur,Dausa and Karauli Districts of Rajasthan,India

Our study interprets large-scale gravity data to delineate concealed banded iron formation (BIF) iron mineralization in India's Rajasthan province. The study area belongs to the Bharatpur, Dausa, and Karauli districts of Rajasthan. We measured 1462 gravity readings to understand the rock types, depth and geometry of the different rock formations in the proposed study area. We also collected representative lithologies from more than 100 locations in the study area and calculated their density values. The measured gravity datasets are investigated via qualitative (e.g., Bouguer anomaly, first derivative and second derivative) and quantitative (radially averaged power spectrum, 3D Euler deconvolution, and 3D inversion) approach. The qualitative methods suggest a general NE–SW orientation of the BIFs, controlled by the general trend of the study area's structural setting. The lithological contact between the Bhilwara and Vindhyan Supergroups is demarcated by a NE–SW trending steep gravity gradient zone. In this area, representative lithologies yield high densities (about 3.746 gm/cc), and the samples identified as BIF represent exploration targets for iron ore. We have also developed our own in-house 3D gravity inversion code in this study. A model space inversion algorithm is converted into a data space using the identity relationship. It makes inversion algorithm very user-friendly on conventional desktop computers. The outcomes from the 3D inversion suggest that the concealed iron ore thickens to the west. This interpretation is also in good correlation with Euler 3D deconvolution of the gravity data.

     

Application of Global Particle Swarm Optimization for Inversion of Residual Gravity Anomalies Over Geological Bodies with Idealized Geometries

A global particle swarm optimization (GPSO) technique is developed and applied to the inversion of residual gravity anomalies caused by buried bodies with simple geometry (spheres, horizontal, and vertical cylinders). Inversion parameters, such as density contrast of geometries, radius of body, depth of body, location of anomaly, and shape factor, were optimized. The GPSO algorithm was tested on noise-free synthetic data, synthetic data with 10% Gaussian noise, and five field examples from different parts of the world. The present study shows that the GPSO method is able to determine all the model parameters accurately even when shape factor is allowed to change in the optimization problem. However, the shape was fixed a priori in order to obtain the most consistent appraisal of various model parameters. For synthetic data without noise or with 10% Gaussian noise, estimates of different parameters were very close to the actual model parameters. For the field examples, the inversion results showed excellent agreement with results from previous studies that used other inverse techniques. The computation time for the GPSO procedure is very short (less than 1 s) for a swarm size of less than 50. The advantage of the GPSO method is that it is extremely fast and does not require assumptions about the shape of the source of the residual gravity anomaly.     

Lunar regolith breccia Dhofar 287B: A record of lunar volcanism

Abstract— Dhofar 287 (Dho 287), a recently found lunar meteorite, consists in large part (95%) of low‐Ti mare basalt (Dho 287A) and a minor, attached portion (?5%) of regolith breccia (Dho 287B). The present study is directed mainly at the breccia portion of this meteorite. This breccia consists of a variety of lithic clasts and mineral fragments set in a fine‐grained matrix and minor impact melt. The majority of clasts and minerals appear to have been mainly derived from the low‐Ti basalt suite, similar to that of Dho 287A. Very low‐Ti (VLT) basalts are a minor lithology of the breccia. These are significantly lower in Mg# and slightly higher in Ti compared to Luna 24 and Apollo 17 VLT basalts. Picritic glasses constitute another minor component of the breccia and are compositionally similar to Apollo 15 green glasses. Dho 287B also contains abundant fragments of Mg‐rich pyroxene and anorthite‐rich plagioclase grains that are absent in the lithic clasts. Such fragments appear to have been derived from a coarse‐grained, Mg#‐rich, Na‐poor lithology. A KREEP component is apparent in chemistry, but no highlands lithologies were identified. The Dho 287 basaltic lithologies cannot be explained by near‐surface fractionation of a single parental magma. Instead, magma compositions are represented by a picritic glass; a low‐Ti, Na‐poor glass; and a low‐Ti, Na‐enriched source (similar to the Dho 287A parental melt). Compositional differences among parent melts could reflect inhomogeneity of the lunar mantle. Alternatively, the low‐Ti, Na‐poor, and Dho 287A parent melts could be of hybrid compositions, resulting from assimilation of KREEP by picritic magma. Thus, the Dho 287B breccia contains lithologies from multiple magmatic eruptions, which differed in composition, formational conditions, and cooling histories. Based on this study, the Dho 287 is inferred to have been ejected from a region located distal to highlands terrains, possibly from the western limb of the lunar nearside, dominated by mare basalts and KREEP‐rich lithologies.     

Prioritization of subwatersheds based on quantitative morphometric analysis in lower Bhavani basin,Tamil Nadu,India using DEM and GIS techniques

Watershed prioritization is one of the most important processes in natural resource management system especially in areas of sustainable watershed development and planning. Morphometric characteristics are the viable entity to understand the hydrological behavior of the subwatershed. For prioritization of subwatershed, morphometric analysis was utilized by using the linear, areal, and relief aspects of the drainage basin. In this context, remote sensing and GIS has been proved to be an efficient tool to identify the morphological features. The Survey of India (SOI) topographical maps, satellite data IRS-LISS III, and Cartosat DEM data were utilized to understand the drainage pattern and also for prioritization of subwatershed areas. The prioritization of subwatershed has been attempted using novel and quantitative approaches based on compound parameter ranking for soil erosion. Lower compound factors were chosen as the most feasible for soil erosion. Based on the observation, eight subwatersheds with a higher degree of the slope were severely prone to soil erosion and remaining 21 subwatersheds occur in low-lying areas that can be developed as sustainable watersheds. The identified subwatershed requires immediate soil remediation and water conservation measures for efficient watershed planning and management. The proposed study might be helpful for resource planners, government agencies, private sectors, and other stake holders to take up soil conservation measures and fixation of water-harvesting structures for better decision making.     

Locating oceanic Tuna resources in the Eastern Arabian sea using remote sensing

Tuna fishery resources are currently under exploited. The resource potential of tunas in the Indian Exclusive Economic Zone (EEZ) beyond 50 m depths is around 2.09 lakh tonnes as estimated by Fishery Survey of India. The distribution and availability of the tuna are governed by environmental factors like temperature, thermocline depth, availability of prey, visibility etc. Remote sensing provides synoptic information on productivity in terms of chlorophyll and Sea Surface Temperature (SST). In the present paper, satellite remote sensing data from Indian Remote Sensing Satellite IRS- P4 Ocean Colour Monitor (OCM) sensor for chlorophyll-a and diffuse attenuation coefficient (K) and National Oceanic and Atmospheric Administration (NOAA) - Advanced Very High Resolution Radiometer (AVHRR) sensor data for sea surface temperature were analysed and correlated within situ catch data of oceanic tunas, Skipjack(Katsuwonus pelamis) and Yellowfin tuna(Thunnus albacares), off Maharashtra coast. Higher catches were found to be associated with moderate to good primary productivity and in the vicinity of thermal fronts. Relationship between Hooking rate and SST has shown that SST of 28–30°C range is optimum for skipjack and 28–31°C for yellowfin tuna. Besides satellite derived chlorophyll and SST for identification of potential tuna fishing zones, role of diffuse attenuation coefficient (K) for visibility factor is also discussed.     

Evaluation of vertical impact factor coefficients for continuous and integral railway bridges under high-speed moving loads

In the railway bridge analysis and design method,dynamic train loads are regarded as static loads enhanced by an impact factor(IF).The IF coefficients for various railway bridges have been reported as a function of span length or frequency of the bridges in Eurocode(2003).However,these IF coefficient values neglect the effects of very high speeds(>200 km/h)and soil-structure interaction(SSI).In this work,a comprehensive study to assess the impact factor coefficients of mid-span vertical displacements for continuous and integral railway bridges subjected to high-speed moving loads is reported.Three different configurations,each for the three-dimensional(3D)continuous and integral bridge,are considered.Also,single-track(1-T)and two-track(2-T)“real train”loading cases for both these bridge types are considered.Subsequently,finite element analysis of the full-scale 3D bridge models,to identify their IF values,considering the effects of SSI for three different soil conditions,is conducted.The IF values obtained from the study for both bridge types are comparable and are greater than the values recommended by Eurocode(2003).The results reveal that with a loss of soil stiffness,the IF value reduces;thus,it confirms the importance of SSI analysis.     

Analysis of magnetic anomalies due to some two-dimensional bodies using the Mellin transform

An interpretation technique using the Mellin transform is suggested for the analysis of magnetic anomalies due to some two-dimensional structures namely (i) a vertical sheet of both finite and infinite depth extent, (ii) a thick dyke and (iii) a horizontal circular cylinder. The Mellin transformed magnetic anomalies resemble gamma functions which are amenable to an easy interpretation. This procedure is illustrated with a small number of synthetic examples in each case. The practicality of the method is exemplified with the well-known vertical magnetic anomalies of Kursk (USSR) in the case of an infinite sheet model and Karimnagar magnetic anomaly (India) in the case of a horizontal circular cylinder. The results are compared with the techniques already available and found to be reliable.     

Growth and decay of the Arabian Sea mini warm pool during May 2000: Observations and simulations

Oceanographic surveys were carried out in the southeastern Arabian Sea in two phases during 14–17 May and 23–27 May 2000 to study the evolution of the Arabian Sea mini warm pool. The hydrographic data collected during these experiments, satellite imagery of sea surface temperature for the corresponding period and ship drift data were utilized for this purpose. Both satellite and in situ observations indicated the presence of a mini warm pool up to 14 May and its dissipation thereafter. The Arabian Sea mini warm pool region was characterized by low-salinity Bay of Bengal water (<22 sigma-t). This watermass was advected into the southeastern Arabian Sea by the East India Coastal Current mainly during November–January and recirculated in this region in an anticyclonic eddy. Immediately after the dissipation of the mini warm pool, this watermass could not be traced in this region. Simulation of mini warm pool characteristics using the Princeton Ocean Model (POM) agreed well with in situ and satellite observations. Moreover, the model showed strengthening of southerly currents in the warm pool region during its dissipation period, which caused the disappearance of the Bay of Bengal watermass. Specific experiments with the POM revealed the significance of heat flux components, low wind speed and salinity stratification in the formation and sustenance of the mini warm pool.