Constraints on Mercury’s Na exosphere: Combined MESSENGER and ground-based data

We have used observations of sodium emission obtained with the McMath-Pierce solar telescope and MESSENGER’s Mercury Atmospheric and Surface Composition Spectrometer (MASCS) to constrain models of Mercury’s sodium exosphere. The distribution of sodium in Mercury’s exosphere during the period January 12-15, 2008, was mapped using the McMath-Pierce solar telescope with the 5″ × 5″ image slicer to observe the D-line emission. On January 14, 2008, the Ultraviolet and Visible Spectrometer (UVVS) channel on MASCS sampled the sodium in Mercury’s anti-sunward tail region. We find that the bound exosphere has an equivalent temperature of 900-1200 K, and that this temperature can be achieved if the sodium is ejected either by photon-stimulated desorption (PSD) with a 1200 K Maxwellian velocity distribution, or by thermal accommodation of a hotter source. We were not able to discriminate between the two assumed velocity distributions of the ejected particles for the PSD, but the velocity distributions require different values of the thermal accommodation coefficient and result in different upper limits on impact vaporization. We were able to place a strong constraint on the impact vaporization rate that results in the release of neutral Na atoms with an upper limit of 2.1 × 10⁶ cm⁻² s⁻¹. The variability of the week-long ground-based observations can be explained by variations in the sources, including both PSD and ion-enhanced PSD, as well as possible temporal enhancements in meteoroid vaporization. Knowledge of both dayside and anti-sunward tail morphologies and radiances are necessary to correctly deduce the exospheric source rates, processes, velocity distribution, and surface interaction.     

Environmental geochemistry of toxic heavy metals in soils around Sarcheshmeh porphyry copper mine smelter plant,Rafsanjan, Kerman,Iran

The Sarcheshmeh copper mine smelter plant is one of the biggest copper producers in Iran. Long-time operation of about 25 years of the smelter plant causes release of potentially toxic heavy metals into the environment. In this paper, geochemical distribution of toxic heavy metals in 28 soil samples was evaluated around the Sarcheshmeh smelter plant. Soils developed over the nonmineralized and uncontaminated areas have an average background concentration of 41.25 mg kg⁻¹ Cu, 26.6 mg kg⁻¹ As, 12.7 mg kg⁻¹ Pb, 0.9 mg kg⁻¹ Sb, 1.9 mg kg⁻¹ Mo, 1.7 mg kg⁻¹ Sn, 0.2 mg kg⁻¹ Cd, 0.15 mg kg⁻¹ Bi, 235 mg kg⁻¹ S and 73.4 mg kg⁻¹ Zn, respectively. As a result of smelting process, the upper soil layers (0–5 cm) were polluted by Cu (>1,397 mg kg⁻¹), Cd (>3.42 mg kg⁻¹), S (>821 mg kg⁻¹), Mo (>10.3 mg kg⁻¹), Sb (>11.7 mg kg⁻¹), As (>120.6 mg kg⁻¹), Pb (>83.8 mg kg⁻¹), Zn (>214.9 mg kg⁻¹), and Sn (>3.7 mg kg⁻¹), respectively. These values are much higher than the normal concentration of the elements in the uncontaminated soil layers. The elemental values decrease with distance travelled away of the smelter plant, especially at minimum wind direction. Furthermore, high contaminated values of Cu (8,430 mg kg⁻¹), As (500 mg kg⁻¹), Pb (331 mg kg⁻¹), Mo (61 mg kg⁻¹), Sb (56.2 mg kg⁻¹), Zn (664 mg kg⁻¹), Cd (17.2 mg kg⁻¹), Bi (13.4 mg kg⁻¹), and S (3,780 mg kg⁻¹) were observed in the upper soil layers close to the smelting waste dumps. Sequential extraction analysis shows that about 270 mg kg⁻¹ Cu, 28 mg kg⁻¹ Pb, 50.33 mg kg⁻¹ Zn, and 47.84 mg kg⁻¹ As were adsorbed by Fe and Mn oxides. The carbonate phases include 151 mg kg⁻¹ Cu, 28 mg kg⁻¹ Pb, 25 mg kg⁻¹ Zn, and 32.99 mg kg⁻¹ As. Organic matter adsorbed 314.6 mg kg⁻¹ Cu and 29.18 mg kg⁻¹ Zn.     

Recent developments of the Middle East catalog

This article summarizes a recent study in the framework of the Global Earth model (GEM) and the Earthquake Model of the Middle East (EMME) project to establish the new catalog of seismicity for the Middle East, using all historical (pre-1900), early and modern instrumental events up to 2006. According to different seismicity, which depends on geophysical, geological, tectonic, and seismicity data, this region is subdivided to nine subregions, consisting of Alborz–Azerbaijan, Afghanistan–Pakistan, Saudi Arabia, Caucasus, Central Iran, Kopeh–Dagh, Makran, Zagros, and Turkey (Eastern Anatolia; after 30° E). After omitting the duplicate events, aftershocks, and foreshocks by using the Gruenthal method, and uniform all magnitude to Mw scale, 28,244 main events remain for the new catalog of Middle East from 1250 B.C. through 2006. The magnitude of completeness (Mc) was determined as 4.9 for five out of nine subregions, where the least values of Mc were found to be 4.2. The threshold of Mc is around 5.5, 5.0, 4.5, and 4.0, for the time after 1950, 1963, 1975, and 2000, respectively. The average of teleseismic depths in all regions is less than 15 km. Totally, majority of depth for Kopeh–Dagh and Central Iran, Zagros, and Alborz–Azerbaijan, approximately, is 15, 13, and 11 km and for Afghanistan–Pakistan, Caucasus, Makran, Turkey (after 30° E), and Saudi Arabia is about 9 km.     

Selection of Less Biased Threshold Angles for SAM Classification Using the Real Value–Area Fractal Technique

The accuracy of classification of the Spectral Angle Mapping (SAM) is warranted by choosing the appropriate threshold angles, which are normally defined by the user. Trial‐and‐error and statistical methods are commonly applied to determine threshold angles. In this paper, we discuss a real value–area (RV–A) technique based on the established concentration–area (C–A) fractal model to determine less biased threshold angles for SAM classification of multispectral images. Short wave infrared (SWIR) bands of the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) images were used over and around the Sar Cheshmeh porphyry Cu deposit and Seridune porphyry Cu prospect. Reference spectra from the known hydrothermal alteration zones in each study area were chosen for producing respective rule images. Segmentation of each rule image resulted in a RV–A curve. Hydrothermal alteration mapping based on threshold values of each RV–A curve showed that the first break in each curve is practical for selection of optimum threshold angles. The hydrothermal alteration maps of the study areas were evaluated by field and laboratory studies including X–ray diffraction analysis, spectral analysis, and thin section study of rock samples. The accuracy of the SAM classification was evaluated by using an error matrix. Overall accuracies of 80.62% and 75.45% were acquired in the Sar Cheshmeh and Seridune areas, respectively. We also used different threshold angles obtained by some statistical techniques to evaluate the efficiency of the proposed RV–A technique. Threshold angles provided by statistical techniques could not enhance the hydrothermal alteration zones around the known deposits, as good as threshold angles obtained by the RV–A technique. Since no arbitrary parameter is defined by the user in the application of the RV‐A technique, its application prevents introduction of human bias to the selection of optimum threshold angle for SAM classification.     

Use of a trench adjacent to a shallow foundation as a mitigation measure for hazards associated with reverse faulting

Acta Geotechnica - Evidence from recent earthquakes reminds us that fault-induced permanent ground displacement has a devastating effect on structures in addition to damage caused by wave...     

Estimation and validation of atmospheric water vapor content using a MODIS NIR band ratio technique based on AIRS water vapor products

Atmospheric water vapor validation needs simultaneous, well-defined, and independent information which are not easily available causing limitations in the development of remote sensing water vapor retrieval algorithms. This study is concerned with the retrieval of total atmospheric water vapor content and its validation. A band ratio method has been used to estimate the water vapor content based on Moderate Resolution Imaging Spectroradiometer (MODIS) Near InfraRed (NIR) data. The method uses MODIS bands 17, 18, and 19 as NIR bands and band 2 to remove the land cover reflectance. Furthermore, the Atmospheric Infrared Sounder (AIRS) has been used for both algorithm development and analysis of the results. The method has been modified to take into account the dry condition of the central parts of Iran. Using some various datasets, the method is implemented and evaluated quantitatively. The validation of the water vapor estimates has been undertaken by an analysis of AIRS data. The validation results shows error as low as 9 % for the estimated water vapor using the MODIS NIR band ratio method.     

Estimation from incomplete data files based on magnitude-frequency relation for earthquakes

The problem of incorporating the available seismological information provided by the major events of the historical catalog with those for the short period of instrumental data is investigated. Assuming that the frequency-magnitude law of Gutenberg and Richter is valid for both periods, an estimation procedure for the main parameter of this law and the rate of earthquake occurrence for historical period is presented. Application of the proposed method is demonstrated, using both real and simulated data. An extension to allow for variable quality of complete data with different magnitude values is also included.     

3D Elastic Solutions for Laterally Loaded Discs: Generalised Brazilian and Point Load Tests

This paper investigates the application of a double Fourier series technique to the construction of an elastic stress field in a cylindrical bar subject to lateral boundary loads. The lateral loads, including the constant load boundary conditions, are represented by two Fourier series: one on the perimeter of the circular section (r ₀, θ) and the other on the longitudinal curved surface parallel to the bar axis (z). The technique invokes acceptable potential functions of the Papkovich–Neuber displacement field, satisfying the governing partial differential equations, to assign appropriate odd and even trigonometric Fourier terms in cylindrical coordinates (r, θ, z). The generic solution decomposes the problem of interest to a state of stress caused by two independent boundary conditions along the z axis and θ-polar angle, both superimposed on a solution for which these potentials are the product of the trigonometric terms of the independent variables (θ, z). Constants appearing in the resultant second-order partial differential equations are determined from the generally mixed (tractions and/or displacements) boundary conditions. While the solutions are satisfied exactly at the ends of an infinite bar, they are satisfied weakly on average, in the light of Saint Venant’s approximation at the two ends of a finite bar. The application of the proposed analysis is verified against available elastic solutions for axisymmetric and non-axisymmetric engineering problems such as the indirect Brazilian Tensile Strength and Point Load Strength tests.     

Integrating seismic attributes in the accurate modeling of geological structures and determining the storage of the gas reservoir in Gorgan Plain (North of Iran)

Three dimensional seismic operation of Gorgan Plain was studied around a well, which is situated in North of Iran following the hitting of a thin overpressure gas layer (thickness of 9.6 m), with the purpose of the accurate modeling of geological structures and determining the approximate gas storages. The geological structures of the reservoir were modeled using the seismic attributes (coherence, instantaneous amplitude and spectral decomposition (FFT)). The obtained results clearly demonstrated the shape and volume of the existing structural traps in the studied area. In order to estimate the thickness of gas layer in the 3D seismic volume and determining the gas storage, the thickness changes based on the seismic amplitudes were used because its thickness was less than the critical resolution thickness for this layer. However, due to its low thickness, the lack of indicator peak in seismic sections and strong faults of area, it was difficult to pursue this layer in the seismic volume and map its exact amplitude. Considering this issue, a new method with integrating of seismic attributes was recommended. First, the instantaneous amplitude attribute of the thin reservoir layer reflector in computed synthetic seismogram were fabricated and then the frequency regarding the highest amount (dominant frequency) was chosen by Fourier Transform. Finally, spectral decomposition (FFT) with the resulting frequency was gained over the cross-section of the layer's instantaneous amplitude attribute in the 3D seismic volume choosing a proper time window. In such a situation, an increase of its thickness was seen as its amplitude increase and the minimum gas storage of this reservoir was calculated using the area of the restricted part of high thickness (over 9.6 m).