The constrained normal form algorithm

This paper describes an algorithm which brings a regularizable polynomial perturbation of a three degree of freedom Kepler problem into a normal form which Poisson commutes with the Kepler Hamiltonian. We illustrate the alogrithm with an example: the quadratic Zeeman effect. For other applications of this algorithm see [1],[4], and [5]. The authors have written a program in MAPLE which implements the constrained normal form.     

Stratigraphic analysis and the asymptotic distribution of the coefficient of cross-association

The asymptotic distribution of the coefficient of cross-association between two sequences of nominal data is investigated under the assumption of independent random sampling from a multinomial distribution slightly modified to suit geological applications. This modification is due to merging consecutive strata of the same rock type into one lithological unit, and can be described by a simple Markov chain. A distinction is made between the case in which the theoretical frequencies are known, and the situation where these have to be estimated from the observed data.     

Special Contributions: Pictorial Series of the Manifestations of the Dynamics of the Earth: 5. Kamchatka,Its Early Cartography and Recent Climate Changes

Studia Geophysica et Geodaetica -     

On the Integrability of Limited-Area Numerical Weather Prediction Model ALADIN over Extended Time Periods

This note summarizes results of the first integration of regional numerical weather prediction model ALADIN in a climate mode. The ALADIN model, developed in an international cooperation led by Météo France, is operationally used for weather prediction. The grid step of the model is 12 km; the integration domain covers a major part of Europe. A one-month-long run has been performed with this model on observed boundary conditions (represented by assimilations by the global model ARPEGE). It is demonstrated that no excessive error is generated and accumulated in the model during the integration; hence the model is integrable for extended time periods and may serve a basis for a development towards a regional climate model.     

Geochemistry and mineralogy of Platinum-group elements in the Ransko gabbro–peridotite massif,Bohemian Massif (Czech Republic)

The Ransko gabbro-peridotite massif in Eastern Bohemia is a strongly differentiated intrusive complex of Lower Cambrian age. The complex hosts low grade Ni-Cu ores mainly developed close to the contact of olivine-rich rocks with gabbros, in troctolites and, to a much lesser extent, in both pyroxene and olivine gabbros and plagioclase-rich peridotites. The ore zone is characterized by strong serpentinization and uralitization. The total Ni + Cu locally reaches up to 4 wt%. Anomalous concentrations of platinum-group elements (PGE's) (maximum 532 ppb Pd, 182 ppb Pt, 53 ppb Rh, 15 ppb Ru, 41 ppb Ir) were detected in samples of Cu-Ni and Ni-Cu ores (maximum 2.63 wt% Ni and 2.31 wt% Cu) from the Jezírka orebody. The main ore paragenesis includes pyrrhotite, pentlandite, chalcopyrite, cubanite, pyrite, magnetite, mackinawite, valleriite, ilmenite and sphalerite. During this work, michenerite, froodite, sperrylite, gold, native bismuth, altaite, tsumoite, hessite, an unnamed Bi-Ni telluride, cobaltite-gersdorffite and galena were newly identified. The host rocks originated through partial melting of a slightly depleted mantle source with noble metals scavenged from this primitive magma prior to the development of these rocks.     

Topex/Poseidon Altimetry and Dynamics of the Ocean-Atmosphere System

The T/P altimeter data 1993 – 1997 (cycles 11 – 194) has been analyzed with emphases on seasonal variations in sea surface topography (SST). The amplitude of the annual variations amounted to (5.9±0.3) mm when inverted barometer (IB) corrections were applied and (2.0±0.4) mm without any IB corrections. The amplitude of the semi-annual variations in SST was small with IB corrections applied: (0.6±0.3) mm. However, when no IB corrections were applied, it was (1.8±0.4) mm, i.e. the semiannual variations are at the same level as the annual variations with no IB corrections. The phase angle offset of the annual term has shifted by about 180° when IB correction was applied. The dynamics of the ocean-atmosphere system is discussed and it is concluded that it could, at least partly, be responsible for the above observed effects.     

Assessment of gully erosion rates through interviews and measurements: a case study from northern Ethiopia

Gullying has been widespread in the Ethiopian Highlands during the 20th century. It threatens the soil resource, lowers crop yields in intergully areas through enhanced drainage and desiccation, and aggravates flooding and reservoir siltation. Knowing the age and rates of gully development during the last few decades will help explain the reasons for current land degradation. In the absence of historical written or photographic documentation, the AGERTIM method (Assessment of Gully Erosion Rates Through Interviews and Measurements) has been developed. It comprises measurements of contemporary gully volumes, monitoring of gully evolution over several years and semi‐structured interview techniques. Gully erosion rates in the Dogu'a Tembien District, Tigray, Ethiopia, were estimated in three representative case‐study areas. In Dingilet, gullying started around 1965 after gradual environmental changes (removal of vegetation from cropland in the catchment and eucalyptus plantation in the valley bottom); rill‐like incisions grew into a gully, which increased rapidly in the drier period between 1977 and 1990. The estimated evolution of the total gully volume in the other areas show patterns similar to those of the Dingilet gully. Average gully erosion rate over the last 50 years is 6·2 t ha^?1 a^?1. Since 1995, no new gullies have developed in the study area. Area‐specific short‐term gully erosion rates are now on average 1·1 t ha^?1 a^?1. The successful application of the AGERTIM method requires an understanding of the geomorphology of the study area and an integration of the researchers with the rural society. It reveals that rapid gully development in the study area is some 50 years old and is mainly caused by human‐induced environmental degradation. Under the present‐day conditions of ‘normal’ rain and catchment‐wide soil and water conservation, gully erosion rates are decreasing. Copyright © 2006 John Wiley & Sons, Ltd.     

Time-lapse seismic analysis of pressure depletion in the Southern Gas Basin

In the Southern Gas Basin (SGB) of the North Sea there are many mature gas fields where time‐lapse monitoring could be very beneficial in extending production life. However, the conditions are not immediately attractive for time‐lapse seismic assessment. This is primarily because the main production effect to be assessed is a pore pressure reduction and frame stiffening because of gas production in tight sandstone reservoirs that also have no real seismic direct hydrocarbon indicators. Modelling, based on laboratory measurements, has shown that such an effect would be small and difficult to detect in seismic data. This paper makes two main contributions. Firstly, this is, to our knowledge, the first time‐lapse study in the SGB and involves a real‐data assessment of the viability for detecting production in such an environment. Secondly, the feasibility of using markedly different legacies of data in such a study is addressed, including an assessment of the factors influencing the crossmatching. From the latter, it is found that significant, spatially varying time shifts need to be, and are successfully, resolved through 3‐D warping. After the warping, the primary factors limiting the crossmatching appear to be residual local phase variations, possibly induced by the differing migration strategies, structure, reverberations and different coherencies of the volumes, caused by differences in acquisition‐structure azimuth and acquisition fold. Despite these differences, a time‐lapse amplitude signature is observed that is attributable to production. The character of the 4‐D amplitude anomalies may also indicate variations in stress sensitivity, e.g. because of zones of fracturing. Additionally, warping‐derived time attributes have been highlighted as a potential additional avenue for detection of pressure depletion in such reservoirs. Although the effects are subtle, they may indicate changes in stress/pressure in and around the reservoir because of production. However, to fully resolve the subtle time‐lapse effects in such a reservoir, the data differences need to be better addressed, which may be possible by full re‐processing and pre‐stack analysis, but more likely dedicated 4‐D acquisition would be required.     

Remotely sensed estimation of forest canopy density: A comparison of the performance of four methods

In recent years, a number of alternative methods have been proposed to predict forest canopy density from remotely sensed data. To date, however, it remains difficult to decide which method to use, since their relative performance has never been evaluated. In this study the performance of: (1) an artificial neural network, (2) a multiple linear regression, (3) the forest canopy density mapper and (4) a maximum likelihood classification method was compared for prediction of forest canopy density using a Landsat ETM+ image. Comparison of confusion matrices revealed that the regression model performed significantly worse than the three other methods. These results were based on a z-test for comparison of weighted kappa statistics, which is an appropriate statistic for analysis of ranked categories. About 89% of the variance of the observed canopy density was explained by the artificial neural networks, which outperformed the other three methods in this respect. Moreover, the artificial neural networks gave an unbiased prediction, while other methods systematically under or over predicted forest canopy density. The choice of biased method could have a high impact on canopy density inventories.