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.     

Determination of gravity anomaly at sea level from inversion of satellite gravity gradiometric data

Gravity gradients can be used to determine the local gravity field of the Earth. This paper investigates downward continuation of all elements of the disturbing gravitational tensor at satellite level using the second-order partial derivatives of the extended Stokes formula in the local-north oriented frame to determine the gravity anomaly at sea level. It considers the inversion of each gradient separately as well as their joint inversion. Numerical studies show that the gradients T_zz, T_xx, T_yy and T_xz have similar capability of being continued downward to sea level in the presence of white noise, while the gradient T_yz is considerably worse than the others. The bias-corrected joint inversion process shows the possibility of recovering the gravity anomaly with 1 mGal accuracy. Variance component estimation is also tested to update the observation weights in the joint inversion.     

Declustering experimental variograms by global estimation with fourth order moments

The variogram is a key parameter for geostatistical estimation and simulation. Preferential sampling may bias the spatial structure and often leads to noisy and unreliable variograms. A novel technique is proposed to weight variogram pairs in order to compensate for preferential or clustered sampling . Weighting the variogram pairs by global kriging of the quadratic differences between the tail and head values gives each pair the appropriate weight, removes noise and minimizes artifacts in the experimental variogram. Moreover, variogram uncertainty could be computed by this technique. The required covariance between the pairs going into variogram calculation, is a fourth order covariance that must be calculated by second order moments. This introduces some circularity in the calculation whereby an initial variogram must be assumed before calculating how the pairs should be weighted for the experimental variogram. The methodology is assessed by synthetic and realistic examples. For synthetic example, a comparison between the traditional and declustered variograms shows that the declustered variograms are better estimates of the true underlying variograms. The realistic example also shows that the declustered sample variogram is closer to the true variogram.     

Measuring the effect of an ongoing urbanization process on biodiversity conservation suitability index: integrating scenario-based urban growth modelling with Conservation Assessment and Prioritization System (CAPS)

The present study adopts an integrative modelling methodology, which combines the strengths of the SLEUTH model and the Conservation Assessment and Prioritization System (CAPS) method. By developing a scenario-based geographic information system simulation environment for Hashtpar City, Iran, the manageability of the landscape under each urban growth scenario is analysed. In addition, the CAPS approach was used for biodiversity conservation suitability mapping. The SLEUTH model was implemented to generate predictive urban layers of the years 2020, 2030, 2040 and 2050 for each scenario (dynamic factors for conservation suitability mapping). Accordingly, conservation suitability surface of the area is updated for each time point and under each urban development storyline. Two-way analysis of variance and Duncan’s new multiple range tests were employed to compare the functionality of the three scenarios. Based on results, the managed urban growth scenario depicted better results for manageability of the landscape and less negative impact on conservation suitability values.     

Near-field probabilistic seismic hazard analysis with characteristic earthquake effects

In conventional seismic hazard analysis, uniform distribution over area and magnitude range is assumed for the evaluation of source seismicity which is not able to capture peculiar characteristic of near-fault ground motion well. For near-field hazard analysis, two important factors need to be considered: (1) rupture directivity effects and (2) occurrence of scenario characteristic ruptures in the nearby sources. This study proposed a simple framework to consider these two effects by modifying the predictions from the conventional ground motion model based on pulse occurrence probability and adjustment of the magnitude frequency distribution to account for the rupture characteristic of the fault. The results of proposed approach are compared with those of deterministic and probabilistic seismic hazard analyses. The results indicate that characteristic earthquake and directivity consideration both have significant effects on seismic hazard analysis estimates. The implemented approach leads to results close to deterministic seismic hazard analysis in the short period ranges (T < 1.0 s) and follows probabilistic seismic hazard analysis results in the long period ranges (T > 1.0 s). Finally, seismic hazard maps based on the proposed method could be developed and compared with other methods.     

Strong ground motion record selection for the reliable prediction of the mean seismic collapse capacity of a structure group

How to select a limited number of strong ground motion records (SGMRs) is an important challenge for the seismic collapse capacity assessment of structures. The collapse capacity is considered as the ground motion intensity measure corresponding to the drift‐related dynamic instability in the structural system. The goal of this paper is to select, from a general set of SGMRs, a small number of subsets such that each can be used for the reliable prediction of the mean collapse capacity of a particular group of structures, i.e. of single degree‐of‐freedom systems with a typical behaviour range. In order to achieve this goal, multivariate statistical analysis is first applied, to determine what degree of similarity exists between each selected small subset and the general set of SGMRs. Principal Component analysis is applied to identify the best way to group structures, resulting in a minimum number of SGMRs in a proposed subset. The structures were classified into six groups, and for each group a subset of eight SGMRs has been proposed. The methodology has been validated by analysing a first‐mode‐dominated three‐storey‐reinforced concrete structure by means of the proposed subsets, as well as the general set of SGMRs. The results of this analysis show that the mean seismic collapse capacity can be predicted by the proposed subsets with less dispersion than by the recently developed improved approach, which is based on scaling the response spectra of the records to match the conditional mean spectrum. Copyright © 2010 John Wiley & Sons, Ltd.