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.     

Tectonic significance of Triassic mafic rocks in the June Complex,Sanandaj–Sirjan zone,Iran

This study concentrates on the petrological and geochemical investigation of mafic rocks embedded within the voluminous Triassic June Complex of the central Sanandaj–Sirjan zone (Iran), which are crucial to reconstruct the geodynamics of the Neotethyan passive margin. The Triassic mafic rocks are alkaline to sub-alkaline basalts, containing 43.36–49.09 wt% SiO₂, 5.19–20.61 wt% MgO and 0.66–4.59 wt% total alkalis. Based on MgO concentrations, the mafic rocks fall into two groups: cumulates (Mg# = 51.61–58.94) and isotropic basaltic liquids (Mg# = 24.54–42.66). In all samples, the chondrite-normalized REE patterns show enrichment of light REEs with variable (La/Yb)_N ratios ranging from 2.48 to 9.00, which confirm their amalgamated OIB-like and E-MORB-like signatures. Enrichment in large-ion lithophile elements and depletion in high field strength elements (HFSE) relative to the primitive mantle further support this interpretation. No samples point to crustal contamination, all having undergone fractionation of olivine + clinopyroxene + plagioclase. Nevertheless, elemental data suggest that the substantial variations in (La/Sm)_PM and Zr/Nb ratios can be explained by variable degrees of partial melting rather than fractional crystallization from a common parental magma. The high (Nb/Yb)_PM ratio in the alkaline mafic rocks points to the mixing of magmas from enriched and depleted mantle sources. Abundant OIB alkaline basalts and rare E-MORB appear to be linked to the drifting stage on the northern passive margin of the Neotethys Ocean.     

Macroseepage of methane and light alkanes at the La Brea tar pits in Los Angeles

Geologic seepage of methane and light (C₂-C₅) alkanes was measured at the La Brea Tar Pits in Los Angeles. Samples were collected using flux chambers with stainless steel canisters and analyzed using gas chromatography. Average seepage rates from individual seepage sites were 970 ± 330 mg/h of methane, 14.0 ± 5.5 mg/h of ethane, 9.1 ± 3.7 mg/h of propane, 3.7 ± 1.6 mg/h of i-butane, 0.33 ± 0.16 mg/h of n-butane, 260 ± 120 μg/h of i-pentane, and 5.3 ± 1.9 μg/h of n-pentane, while maximum seepage rates exceeded 17 g/h of methane, 270 mg/h of ethane, 190 mg/h of propane, 95 mg/h of i-butane, 10 mg/h of n-butane, 7 mg/h of i-pentane, and 0.1 mg/h of n-pentane. These absolute fluxes have an additional unknown amount of error associated with them due to sampling methodology, and should be taken as the lower limit of emissions. Samples collected revealed generally dry gas, with high methane emissions relative to the light alkanes. Overall emissions from the tar pits were found to come not only from the active geologic seepage, but also from the outgassing of the standing asphalt at the site. Using the gas ratios, which are negligibly affected by errors introduced by sampling methodology, observed in this study, daily emissions of C₂ – C₅ alkanes from the La Brea area were estimated to be 4.7 ± 1.6 Mg, which represents 2–3 % of total emissions in the entire Los Angeles region.     

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.     

Numerical modeling of toxic nonaqueous phase liquid removal from contaminated groundwater systems: mesh effect and discretization error estimation

Numerical modeling has now become an indispensable tool for investigating the fundamental mechanisms of toxic nonaqueous phase liquid (NAPL) removal from contaminated groundwater systems. Because the domain of a contaminated groundwater system may involve irregular shapes in geometry, it is necessary to use general quadrilateral elements, in which two neighbor sides are no longer perpendicular to each other. This can cause numerical errors on the computational simulation results due to mesh discretization effect. After the dimensionless governing equations of NAPL dissolution problems are briefly described, the propagation theory of the mesh discretization error associated with a NAPL dissolution system is first presented for a rectangular domain and then extended to a trapezoidal domain. This leads to the establishment of the finger‐amplitude growing theory that is associated with both the corner effect that takes place just at the entrance of the flow in a trapezoidal domain and the mesh discretization effect that occurs in the whole NAPL dissolution system of the trapezoidal domain. This theory can be used to make the approximate error estimation of the corresponding computational simulation results. The related theoretical analysis and numerical results have demonstrated the following: (1) both the corner effect and the mesh discretization effect can be quantitatively viewed as a kind of small perturbation, which can grow in unstable NAPL dissolution systems, so that they can have some considerable effects on the computational results of such systems; (2) the proposed finger‐amplitude growing theory associated with the corner effect at the entrance of a trapezoidal domain is useful for correctly explaining why the finger at either the top or bottom boundary grows much faster than that within the interior of the trapezoidal domain; (3) the proposed finger‐amplitude growing theory associated with the mesh discretization error in the NAPL dissolution system of a trapezoidal domain can be used for quantitatively assessing the correctness of computational simulations of NAPL dissolution front instability problems in trapezoidal domains, so that we can ensure that the computational simulation results are controlled by the physics of the NAPL dissolution system, rather than by the numerical artifacts. Copyright © 2014 John Wiley & Sons, Ltd.     

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.