Lava flow hazard at Nyiragongo volcano,D.R.C.

The 2002 eruption of Nyiragongo volcano constitutes the most outstanding case ever of lava flow in a big town. It also represents one of the very rare cases of direct casualties from lava flows, which had high velocities of up to tens of kilometer per hour. As in the 1977 eruption, which is the only other eccentric eruption of the volcano in more than 100 years, lava flows were emitted from several vents along a N–S system of fractures extending for more than 10 km, from which they propagated mostly towards Lake Kivu and Goma, a town of about 500,000 inhabitants. We assessed the lava flow hazard on the entire volcano and in the towns of Goma (D.R.C.) and Gisenyi (Rwanda) through numerical simulations of probable lava flow paths. Lava flow paths are computed based on the steepest descent principle, modified by stochastically perturbing the topography to take into account the capability of lava flows to override topographic obstacles, fill topographic depressions, and spread over the topography. Code calibration and the definition of the expected lava flow length and vent opening probability distributions were done based on the 1977 and 2002 eruptions. The final lava flow hazard map shows that the eastern sector of Goma devastated in 2002 represents the area of highest hazard on the flanks of the volcano. The second highest hazard sector in Goma is the area of propagation of the western lava flow in 2002. The town of Gisenyi is subject to moderate to high hazard due to its proximity to the alignment of fractures active in 1977 and 2002. In a companion paper (Chirico et al., Bull Volcanol, in this issue, 2008) we use numerical simulations to investigate the possibility of reducing lava flow hazard through the construction of protective barriers, and formulate a proposal for the future development of the town of Goma.     

JP3D compression of solar data-cubes: Photospheric imaging and spectropolarimetry

Hyperspectral imaging is an ubiquitous technique in solar physics observations and the recent advances in solar instrumentation enabled us to acquire and record data at an unprecedented rate. The huge amount of data which will be archived in the upcoming solar observatories press us to compress the data in order to reduce the storage space and transfer times. The correlation present over all dimensions, spatial, temporal and spectral, of solar data-sets suggests the use of a 3D base wavelet decomposition, to achieve higher compression rates. In this work, we evaluate the performance of the recent JPEG2000 Part 10 standard, known as JP3D, for the lossless compression of several types of solar data-cubes. We explore the differences in: a) The compressibility of broad-band or narrow-band time-sequence; I or V Stokes profiles in spectropolarimetric data-sets; b) Compressing data in [x,y, λ] packages at different times or data in [x,y,t] packages of different wavelength; c) Compressing a single large data-cube or several smaller data-cubes; d) Compressing data which is under-sampled or super-sampled with respect to the diffraction cut-off.     

Heavy metals in mangrove sediments of the central Arabian Gulf shoreline,Saudi Arabia

To assess heavy metals in mangrove swamps of Sehat and Tarut coastal areas along the Arabian Gulf, 18 sediment samples were collected for Al, V, Cr, Mn, Cu, Zn, Cd, Pb, Hg, Sr, As, Fe, Co, and Ni analysis. The results indicated that the distribution of some metals was largely controlled by anthropogenic inputs, while others were of terrigenous origin and most strongly associated with distribution of aluminum and total organic carbon in sediments. Mangrove sediments were extremely severe enriched with Sr (EF?=?67.59) and very severe enriched with V, Hg, Cd, Cu, As (EF?=?44.28, 37.45, 35.77, 25.97, and 11.53, respectively). Average values of Sr, V, Hg, Cd, Cu, Ni, As, and Cr were mostly higher than the ones recorded from the Mediterranean Sea, the Red Sea, the Gulf of Aqaba, the Caspian Sea, the Arabian and Oman gulfs, coast of Tanzania, sediment quality guidelines, and the background shale and the earth crust. Landfilling due to coastal infrastructure development around mangrove forests, oil spills and petrochemical and desalination effluents from Al-Jubail industrial city to the north were the anthropogenic activities that further enhanced heavy metals in the studied mangrove sediments.     

Mn-oxides and sequestration of heavy metals in a suburban catchment basin of the Chesapeake Bay watershed

The Chesapeake Bay is greatly impacted by numerous pollutants including heavy metals and understanding the controls on the distribution of heavy metals in the watershed is critical to mitigation and remediation efforts in controlling this type of pollution. Clasts from a stormwater catchment basin draining a subdivision near George Mason University, Fairfax VA (38°50.090°N 78°19.204°W) were investigated using X-ray diffraction (XRD), Scanning electron microcopy (SEM) and energy dispersive spectroscopy (EDS) to determine the nature of Mn-oxide coatings and relationship to bound heavy metals. Mn-oxides are poorly crystalline and occur as subhedral to anhedral platy particles and more rarely as euhedral plates. Micronodules are a commonly observed texture. Chemical compositions of coatings are variable with average major constituent concentrations being Mn (33.38 wt%), Fe (11.88 wt%), Si (7.33 wt%), Al (5.03 wt%), and Ba (0.90 wt%). Heavy metals are found in the coatings with Zn being most prevalent, occurring in approximately 58% of analyses with an average concentration of (0.66 wt%). Minor amounts of Co, Ni, Pb, and Cl are observed. Heavy metals and Cl are interpreted as being derived from road pollution. Mn-oxides can serve as a sequestration mechanism for pollution but may also release heavy metals. Field and laboratory observations indicate Mn-oxides occurring on the surface of the clasts can be mechanically mobilized. This is a mechanism for transporting heavy metals into the Chesapeake Bay watershed. Deicing agents may serve as a mechanism to release heavy metals through cation exchange and increased ionic strength. This is the first detailed mineralogical investigation of Mn-oxides and the roles they may play in pollution in the Chesapeake Bay.     

Position and velocity perturbations for the determination of geopotential from space geodetic measurements

Although space geodetic observing systems have been advanced recently to such a revolutionary level that low Earth Orbiting (LEO) satellites can now be tracked almost continuously and at the unprecedented high accuracy, none of the three basic methods for mapping the Earth’s gravity field, namely, Kaula linear perturbation, the numerical integration method and the orbit energy-based method, could meet the demand of these challenging data. Some theoretical effort has been made in order to establish comparable mathematical modellings for these measurements, notably by Mayer-Gürr et al. (J Geod 78:462–480, 2005). Although the numerical integration method has been routinely used to produce models of the Earth’s gravity field, for example, from recent satellite gravity missions CHAMP and GRACE, the modelling error of the method increases with the increase of the length of an arc. In order to best exploit the almost continuity and unprecedented high accuracy provided by modern space observing technology for the determination of the Earth’s gravity field, we propose using measured orbits as approximate values and derive the corresponding coordinate and velocity perturbations. The perturbations derived are quasi-linear, linear and of second-order approximation. Unlike conventional perturbation techniques which are only valid in the vicinity of reference mean values, our coordinate and velocity perturbations are mathematically valid uniformly through a whole orbital arc of any length. In particular, the derived coordinate and velocity perturbations are free of singularity due to the critical inclination and resonance inherent in the solution of artificial satellite motion by using various types of orbital elements. We then transform the coordinate and velocity perturbations into those of the six Keplerian orbital elements. For completeness, we also briefly outline how to use the derived coordinate and velocity perturbations to establish observation equations of space geodetic measurements for the determination of geopotential.     

Geothermal energy potential of Tulsishyam thermal springs of Gujarat,India

Tulsishyam thermal springs are located in the Saurashtra region of Gujarat, India with discharge temperatures varying from 39 to 42 °C. The pH of these thermal springs varies from 7.1 to 7.4, indicating neutral character. Though these thermal springs propagate through the near surface layer of Deccan basalt, detailed geochemical analysis of the thermal waters using Piper diagram suggests that the water is interacting with the granitic basement rock. Silica and cation geothermometry estimates have reservoir temperature in the range of 138 to 207 °C categorizing it into a low to moderate enthalpy geothermal system. Furthermore, the area has high heat flow values of 53–90 mW/m² because of shallow Moho depth. The prevailing conditions suggest that the geothermal energy can potentially be exploited through an enhanced geothermal system (EGS). The study also indicates different mineral phases that may precipitate out of water during exploitation of geothermal energy and it should be taken into account while designing an EGS for the area.     

Assessment of H<Subscript>2</Subscript>S emission hazards into tunnels: the Nosoud tunnel case study from Iran

The presence of a gassy ground condition is an important problem in tunneling. In this study, the effects of groundwater H₂S and CH₄ emissions are investigated and characterized together with the factors that created these conditions in Nosoud tunnel in Iran. Through the geological investigations, the presence of these gasses was not detected prior to the construction of the tunnel. Groundwater sampling indicated that about 1 L of H₂S is released per 100 L of the water inflow into the Nosoud tunnel under normal conditions. However, the volume of the released gas was varying with the changes in the groundwater discharge rate. Thus, estimation of groundwater inflow into the tunnel is necessary for predicting the volume of gas emission. Based on the experience of the Nosoud tunnel excavations, there are several geological and hydrogeological factors that must be considered as the indicators of gas emissions during tunneling. Considering the importance of ground water gas emission into the tunnels located in gassy conditions, the present work was conducted to predict the H₂S seepage before the excavation using geological and hydrogeological indicators.     

Hydrogeological assessment of the Tabarteh fault zone by physicochemical analysis and multivariate statistical methods

Hydrogeologically, faults may impede, conduit, exert no influence, or may play a combination of these roles on groundwater flow. The object of this paper is to study the hydrogeological role of the Tabarteh fault, which is located on the border of Zagros and Central Iran tectonic zones in an alluvial aquifer. The recorded data of water table levels, chemical parameters, and discharge rate of wells, in addition to geological maps and geophysical results, were collected and evaluated. The outcrop of travertine in limited areas and the emergence of a few small springs within the alluvium show a barrier role of the fault in the groundwater flow. The spatial analysis of chemical components, head time series, and groundwater flow direction assessment demonstrated that the fault acts as both a barrier and a non-barrier in different sections. The multivariate statistical methods of cluster and discriminant analyses also confirm the dual role of the fault.     

Effects of historical and projected land use/cover change on runoff and sediment yield in the Netravati river basin,Western Ghats,India

In this study, the effects of changes in historical and projected land use land cover (LULC) on monthly streamflow and sediment yield for the Netravati river basin in the Western Ghats of India are explored using land use maps from six time periods (1972, 1979, 1991, 2000, 2012, and 2030) and the soil and water assessment tool (SWAT). The LULC for 2030 is projected using the land change modeller with the assumption of normal growth. The sensitivity analysis, model calibration, and validation indicated that the SWAT model could reasonably simulate streamflow and sediment yield in the river basin. The results showed that the spatial extent of the LULC classes of urban (1.80–9.96%), agriculture (31.38–55.75%), and water bodies (1.48–2.66%) increased, whereas that of forest (53.04–27.03%), grassland (11.17–4.41%), and bare land (1.09–0.16%) decreased from 1972 to 2030. The streamflow increased steadily (7.88%) with changes in LULC, whereas the average annual sediment yield decreased (0.028%) between 1972 and 1991 and increased later (0.029%) until 2012. However, it may increase by 0.43% from 2012 to 2030. The results indicate that LULC changes in urbanization and agricultural intensification have contributed to the increase in runoff, amounting to 428.65 and 58.67 mm, respectively, and sediment yield, amounting to 348 and 43 ton/km², respectively, in the catchment area from 1972 to 2030. The proposed methodology can be applied to other river basins for which temporal digital LULC maps are available for better water resource management plans.