T -distributed Random Fields: A Parametric Model for Heavy-tailedWell-log Data1

Histograms of observations from spatial phenomena are often found to be more heavy-tailed than Gaussian distributions, which makes the Gaussian random field model unsuited. A T-distributed random field model with heavy-tailed marginal probability density functions is defined. The model is a generalization of the familiar Student-T distribution, and it may be given a Bayesian interpretation. The increased variability appears cross-realizations, contrary to in-realizations, since all realizations are Gaussian-like with varying variance between realizations. The T-distributed random field model is analytically tractable and the conditional model is developed, which provides algorithms for conditional simulation and prediction, so-called T-kriging. The model compares favourably with most previously defined random field models. The Gaussian random field model appears as a special, limiting case of the T-distributed random field model. The model is particularly useful whenever multiple, sparsely sampled realizations of the random field are available, and is clearly favourable to the Gaussian model in this case. The properties of the T-distributed random field model is demonstrated on well log observations from the Gullfaks field in the North Sea. The predictions correspond to traditional kriging predictions, while the associated prediction variances are more representative, as they are layer specific and include uncertainty caused by using variance estimates.     

Cyclic variations in the differential rotation of the solar corona

The rotation of the solar corona is analyzed using the original database on the brightness of the FeXIV 530.3 nm coronal green line covering six recent activity cycles. The rate of the differential rotation of the corona depends on the cycle phase. In decay phases, there are only small differences in the rotation, which are similar to that of a rigid body. The differences are more significant (though less pronounced than in the photosphere) during rise phases, just before maxima, and sometimes at maxima. The total rate of the coronal rotation is represented as a superposition of two, i.e., fast and slow modes. The synodic period of the fast mode is approximately 27 days at the equator and varies slightly with time. This mode displays weak differences in rotation and is most pronounced in the middle of decay phases. The slow mode is manifested only at high latitudes during the rise phases of activity, and displays a mean period of 31 days. The relative contribution of each mode to the total rotational rate is determined as a function of time and heliographic latitude. These results indicate that the structure of the velocity field in the convective zone must also vary with time. This conclusion can be verified by helioseismology measurements in the near future.     

The Removal of Dissolved Metals by Hydroxysulphate Precipitates during Oxidation and Neutralization of Acid Mine Waters, Iberian Pyrite Belt

This study examines the removal of dissolved metals during the oxidation and neutralization of five acid mine drainage (AMD) waters from La Zarza, Lomero, Esperanza, Corta Atalaya and Poderosa mines (Iberian Pyrite Belt, Huelva, Spain). These waters were selected to cover the spectrum of pH (2.2–3.5) and chemical composition (e.g., 319–2,103 mg/L Fe; 2.85–33.3 g/L SO₄⁼) of the IPB mine waters. The experiments were conducted in the laboratory to simulate the geochemical evolution previously recognized in the field. This evolution includes two stages: (1) oxidation of dissolved Fe(II) followed by hydrolysis and precipitation of Fe(III), and (2) progressive pH increase during mixing with fresh waters. Fe(III) precipitates at pH < 3.5 (stages 1 and 2) in the form of schwertmannite, whereas Al precipitates during stage 2 at pH 5.0 in the form of several hydroxysulphates of variable composition (hydrobasaluminite, basaluminite, aluminite). During these stages, trace elements are totally or partially sorbed and/or coprecipitated at different rates depending basically on pH, as well as on the activity of the SO₄⁼ anion (which determines the speciation of metals). The general trend for the metals which are chiefly present as aqueous free cations (Pb²⁺, Zn²⁺, Cu²⁺, Cd²⁺, Mn²⁺, Co²⁺, Ni²⁺) is a progressive sorption at increasing pH. On the other hand, As and V (mainly present as anionic species) are completely scavenged during the oxidation stage at pH < 3.5. In waters with high activities (> 10⁻¹) of the SO ₄⁼ ion, some elements like Al, Zn, Cd, Pb and U can also form anionic bisulphate complexes and be significantly sorbed at pH < 5. The removal rates at pH 7.0 range from around 100% for As, V, Cu and U, and 60–80% for Pb, to less than 20% for Zn, Co, Ni and Mn. These processes of metal removal represent a significant mechanism of natural attenuation in the IPB.     

Avalanche Risk During Backcountry Skiing – An Analysis of Risk Factors

Skier-triggered avalanches are the main cause of avalanche accidents in backcountry skiing. The risk of accidents during backcountry skiing was analysed statistically and related to factors such as elevation level, aspect, stability rating and the time of the year. The analysis is based on a database about terrain usage and avalanche accidents from a large heli-skiing operator in Canada, which makes it possible to study the conditional probability of accidents given the recorded pattern of terrain usage. This study shows that the historical risk of accidentally triggering an avalanche greater than size 1 depends highly on the stability rating, with the highest risk occurring during “poor” stability. The risk is greater at high elevations, and it is lower during the late season than earlier on. Skier risk does not depend as much on aspect as may be indicated from avalanche data alone. However, it is relatively high in the N–NE–E sector. These factors are not independent of each other and therefore analyses of combined factors were also performed. Questionnaires and interviews were used to gain knowledge about the terrain selection of professional mountain guides. These results indicate that when selecting terrain, guides first look at the overall shape and size of the terrain, but avalanche history of terrain and inclination are also important factors. Finally, remarks in avalanche reports were analysed, and common human factors identified.     

Effect of sampling and linkage on fault length and length–displacement relationship

The power-law exponent (n) in the equation: D=cL ⁿ , with D = maximum displacement and L = fault length, would be affected by deviations of fault trace length. (1) Assuming n=1, numerical simulations on the effect of sampling and linkage on fault length and length–displacement relationship are done in this paper. The results show that: (a) uniform relative deviations, which means all faults within a dataset have the same relative deviation, do not affect the value of n; (b) deviations of the fault length due to unresolved fault tip decrease the values of n and the deviations of n increase with the increasing length deviations; (c) fault linkage and observed dimensions either increase or decrease the value of n depending on the distribution of deviations within a dataset; (d) mixed deviations of the fault lengths are either negative or positive and cause the values of n to either decrease or increase; (e) a dataset combined from two or more datasets with different values of c and orders of magnitude also cause the values of n to deviate. (2) Data including 19 datasets and spanning more than eight orders of fault length magnitudes (10⁻²–10⁵ m) collected from the published literature indicate that the values of n range from 0.55 to 1.5, the average value being 1.0813, and the peak value of n _d (double regression) is 1.0–1.1. Based on above results from the simulations and published data, we propose that the relationship between the maximum displacement and fault length in a single tectonic environment with uniform mechanical properties is linear, and the value of n deviated from 1 is mainly caused by the sampling and linkage effects.     

韧性剪切变形对岩石地球化学行为的制约——以北阿尔金巴什考供韧性剪切带为例

韧性剪切变形对岩石地球化学行为的制约一直是地质学家们探讨的课题。本文以构成北阿尔金红柳沟——拉配泉俯冲碰撞杂岩带与北阿尔金地块边界的巴什考供斜向逆冲型韧性剪切带为例，通过对韧性剪切带内花岗岩变形前后不同变形强度构造岩的地球化学组成进行对比，确定等比线斜率，探讨韧性变形对岩石体积和成分变异的影响。计算结果表明，在糜棱岩化过程中，糜棱岩化花岗岩体积亏损21％，花岗质糜棱岩体积亏损31％。质量平衡计算结果和等比线图表明，韧；陛剪切作用导致SiO2，流失量最大，A12O3、K2O及Ba、Rb、Sr等都有不同程度的丢失，显示出较强的活动性，MnO、P2O5、Sc位于等比线上或附近，表现出相对的稳定性。岩石中活动组分的变异是流体渗滤作用引起的，不活动组分的变异是体．积亏损造成的。     

Geometry and accuracy of reflecting points in bistatic satellite altimetry

The analysis of the time and space distribution of specular (reflecting) points in bistatic altimetry between GPS and CHAMP satellites or SAC-C (taken as examples) is extended from Wagner and Klokočník (2003 J. Geod 77: 128–138). We demonstrate a significantly higher number and density of reflecting points in bistatic altimetry in comparison with traditional monostatic altimetry. After an outline of our older accuracy assessment for the vertical position of the reflecting point, we add a new independent derivation and compare both approaches. We account for orbit errors of both the transmitters (GPS) and receiver (CHAMP) satellites, and the measurement (delay) error. We found that the accuracy of the vertical position of the reflecting point decreases only slowly with increasing off-nadir angle and that the orbit errors must be accounted for if decimeter and better accuracy is required. In this paper, we do not study errors such as state of the ocean, technical parameters of the receiving system, and atmospheric corrections.     

Daily Meteorological Observations in Cádiz – San Fernando. Analysis of the Documentary Sources and the Instrumental Data Content (1786–1996)

Meteorological observations in the city of Cádiz are acknowledged as having been made from the middle of the 18th century onwards although they were only recorded and preserved in documentary form from 1789 onwards. Data readings were taken at the new Naval Observatory in San Fernando, ten kilometres from Cádiz, from 1797 onwards. Continuous series for temperature and atmospheric pressure at a daily resolution have been compiled and constructed from 1817 until 1996. The data series is composed of thrice daily observations made at Cádiz (1821–1880) by local observers and hourly data from Naval Observatory at San Fernando (1870–1996).     

Simultaneous use of geochemical and geophysical methods to characterise abandoned landfills

. Groundwater sampling and geophysical methods determined a serious contamination problem associated with refilled exploitation sites at the Cal Dimoni area, Llobregat delta, Barcelona, Spain. To characterise this process, hydrogeochemical analyses were performed and showed the following modifications to groundwater chemical composition: increasing pH values, changing redox conditions, significant increases in total organic carbon (TOC) and certain trace elements, and high groundwater conductivity values. Major ion content accumulations were found under the refilled area. In contrast, elements involved in the oxidation–reduction processes, such as iron, manganese and nitrates, clearly diminished. Electromagnetic prospecting methods were also performed and delineated the contamination plume extent. These methods also showed separate sources of contamination, one clearly related to the groundwater–refilled zone leachate interaction, another as a consequence of the manure–accumulation surface site. Geochemical and geophysical methods have shown similar results for locating groundwater contamination sources, and for determining leachate generation mechanisms and flow paths.