Degradation of Geopotential Recovery from Short Repeat-Cycle Orbits: Application to GRACE Monthly Fields

Throughout 2004 the GRACE (Gravity Recovery And Climate Experiment) orbit contracted slowly to yield a sparse repeat track of 61 revolutions every 4 days on 19 September 2004. As a result, we show from linear perturbation theory that geopotential information previously available to fully resolve a gravity field every month of 120× 120 (degree by order) in spherical harmonics was compressed then into about one-fourth of the necessary observation space. We estimate from this theory that the ideal gravity field resolution in September 2004 was only about 30 × 30. More generally, we show that any repeat-cycle mission for geopotential recovery with full resolution L × L requires the number of orbit-revolutions-to-repeat to be greater than 2L.     

First results of DORIS data analysis at Geodetic Observatory Pecný

In a cooperation between the Astronomical Institute, University of Bern (AIUB), the Geodetic Observatory Pecný (GOPE), and the Institut Géographique National (IGN), DORIS data analysis capabilities were implemented into a development version of the Bernese GPS software. The DORIS Doppler observables are reformulated such that they are similar to global navigation satellite system (GNSS) carrier-phase observations, allowing the use of the same observation models and algorithms as for GNSS carrier-phase data analysis with only minor software modifications. As such, the same algorithms may be used to process DORIS carrier-phase observations. First results from the analysis of 3 weeks of DORIS data (September 2004, five DORIS-equipped satellites) at GOPE are promising and are presented here. They include the comparison of station coordinates with coordinate estimates derived by the Laboratoire d’Etudes en Géophysique et Océanographie Spatiale/Collecte Localisation Satellites analysis centre (LCA) and the Institut Géographique National/Jet Propulsion Laboratory (IGN/JPL), and the comparison of Earth orientation parameters (EOPs) with the International Earth Rotation and Reference Frames Service (IERS) C04 model. The modified Bernese results are of a slightly lower, but comparable, quality than corresponding solutions routinely computed within the IDS (International DORIS Service). The weekly coordinate repeatability RMS is of the order of 2–3 cm for each 3D station coordinate. Comparison with corresponding estimates of station coordinates from current IDS analysis centers demonstrates similar precision. Daily pole component estimates show a mean difference from IERS-C04 of 0.6  mas in X _p and  ? 0.5  mas in Y _p and a RMS of 0.8  mas in X _p and 0.9  mas in Y _p (mean removed). An automatic analysis procedure is under development at GOPE, and routine DORIS data processing will be implemented in the near future.     

Low-field variation of magnetic susceptibility measured by the KLY-4S Kappabridge and KLF-4A magnetic susceptibility meter: Accuracy and interpretational programme

The KLY-4S Kappabridge and KLF-4A Magnetic Susceptibility Meter enable automated measurement of susceptibility variation with field in the ranges of 2–450 A/m and 5–300 A/m (in effective values), respectively. Unfortunately, the measurement accuracy decreases with decreasing field and it is not easy to decide whether the susceptibility variation at the lowest fields is natural phenomenon or results from measuring errors. To overcome this problem, the accuracies of both the above instruments were investigated experimentally using artificial specimens (mixture of pure magnetite and plaster of Paris) with variable susceptibilities ranging from 1 × 10⁻⁵ to 5 × 10⁻². The complete curve of the field variation of susceptibility of each specimen was measured 10 times and the relative error was calculated for each field. In the KLY-4S Kappabridge, in specimens with susceptibilities higher than 100 × 10⁻⁶, the relative errors are lower than 3% in all fields and lower than 1% in the fields stronger than 10 A/m. In the KLF-4A Magnetic Susceptibility Meter, in relatively strongly magnetic specimens with susceptibilities 5 × 10⁻⁴ to 5 × 10⁻², the relative error is less than 1.5% in the entire field range. While the former instrument is convenient for investigating almost all rock types, the latter instrument is convenient for measuring moderately and strongly magnetic specimens. To facilitate work with field variation of susceptibility curves, showing variable accuracies with field, the programme FieldVar was written. One of its options is plotting the measured data with corresponding field-variable error bars. In this way, a tool is offered for interpreting such susceptibility changes that are sound and reasonable from the point of view of measuring accuracy.     

Tin–polymetallic sulfide deposits in the eastern part of the Dachang tin field (South China) and the role of black shales in their origin

The Dafulou and Huile vein and stratabound cassiterite-sulfide deposits and sheeted ore veins at the Kangma cassiterite-sulfide deposit are located in the eastern part of the Dachang tin field. These deposits are hosted in a sedimentary sequence containing significant concentrations of organic matter in the form of Lower Devonian calcareous black shales and hornfels. These rocks together with the younger intrusion of Longxianggai granite (91DŽ Ma) actively participated in the formation of Sn-polymetallic deposits. The following three major stages have been distinguished in stratiform and vein-type orebodies at Dafulou, Huile and Kangma: stage I (cassiterite, pyrrhotite, arsenopyrite, tourmaline, carbonate), stage II - main sulfide stage (quartz, cassiterite, arsenopyrite, pyrrhotite, sphalerite, stannite, pyrite, carbonates) and stage III (native Bi, galena, electrum, sulfosalts). Stage IV (post-ore), recognized at Huile is represented by barren carbonates and zeolites. Whole rock geochemistry has revealed that at Dafulou, Bi and Cu correlate strongly with S, whereas V and Pb correlate well with C_org (organic carbon). The similar distribution patterns of selected elements in average slightly mineralized low-Ca black shales indicate a fluid composition similar for all deposits studied. Studies of graphitization of the organic matter in black shales adjacent to orebodies indicate that d₍₀₀₂₎ and FWHM (full width in half maximum)/peak height values gradually decrease in the following sequence: Dafulou deposit M Kangma deposit M Huile deposit. The pyrolysate of wall rocks at the Dafulou deposit is relatively enriched in asphaltenes and maltenes (55.6-72.0% of the pyrolysate) comparable with pyrolysate obtained from more distal black shales (19.2-28.5%). Typical GC-MS spectra of pyrolysate from distal black shales are dominated by alkanes in the n-C15 to n-C25 range, aromatic molecules being represented mostly by alkyl-naphthalenes. In contrast, only traces of aliphatic hydrocarbons in the n-C14 to n-C18 range and elemental sulfur were identified in pyrolysates from pyrrhotitized wall rocks. The earliest fluid inclusions of the studied system occur in the quartz-tourmaline-cassiterite assemblage of stage I at Dafulou. These inclusions are H₂O-CO₂-CH₄-rich, with 10 to 20 vol% of aqueous phase. P-T conditions of the trapping of inclusions are estimated to be up to 400 °C and 1.3 to 2.0 kbar (between 5.0 and 7.5 km under lithostatic pressure). In contrast, the presence of a low density gaseous CO₂-CH₄ phase indicates relatively low pressures during the formation of the breccia-type quartz-calcite-cassiterite-sulfide mineralization (stage II), when P-T conditions probably reached approx. 380 to 400 °C and 0.6 kbar (up to 6 km under hydrostatic pressure). Fluid inclusion data and oxygen isotope thermometry indicate that cassiterite-sulfide ores of the main sulfide stage (stage II) formed from aqueous-carbonic fluid (CO₂/CH₄ =ᄺ) at temperatures of up to 390 °C at Dafulou and in a temperature range of 250 to 360 °C at Huile and 260 to 370 °C at Kangma. The '³⁴S values of sulfides from Dafulou range mostly between -1 and -6‰, whereas sulfides from the Kangma and Huile deposits are characterized by more negative '³⁴S values (between -8 and -11‰, and between -9 and -12‰, respectively). These data suggest that bacteriogenic sulfides of black shales were a dominant source of reduced sulfur for epigenetic (vein and replacement) mineralization. Oxygen isotopic compositions of five quartz-cassiterite pairs from Dafulou and Huile show a relatively narrow range of calculated oxygen isotope temperatures (250-320 °C, using the equation of Alderton 1989) and high '¹⁸O_fluid values between +8 and +10‰ (SMOW), which are in agreement with fluid derivation from and/or high temperature equilibration with the Longxianggai granite. The carbon and oxygen isotope composition of carbonates reflects variable carbon sources. Stage I calcite is characterized by narrow ranges of '¹³C (-7.0 to -9.5‰ PDB) and '¹⁸O (+15.0 to +17.5‰ SMOW). This calcite shows ubiquitous deformation, evidenced by intense development of twins. Fluid compositions calculated at 330 °C for the Dafulou and Huile deposits and at 270-300 °C for the Kangma deposit ('¹⁸O_fluid between +10.0 and +11.5‰ SMOW, '¹³C_fluid between -5.5 and -7.5‰ PDB), agree with fluid derivation from and/or equilibration with the peraluminous, high-'¹⁸O Longxianggai granite and suggest a significant influence of contact metasedimentary sequences (carbon derived from decomposition and/or alteration of organic matter of calcareous black shales). The '¹³ C values of organic matter from the Lower to Upper Devonian host rocks at the Dafulou deposit (-24.0 and -28.0‰) fit with a marine origin from algae. However, organic matter adjacent to the host rock-ore contact displays a slight enrichment in ¹³C. The organic carbon from the Huile and Kangma deposits is even more ¹³C enriched (-24.6 to -23.5‰). The most heavy '¹³ C values (-16.5‰) were detected in hornfels sampled at the contact of the Upper Devonian sediments with the Longxianggai granite. The '¹³C data broadly correlate with the degree of structural ordering (degree of graphitization) of organic matter, which indicates that both variables are related to thermal overprint.     

Ice cover variability in the Caspian and Aral seas from active and passive microwave satellite data

The paper discusses time and space variations of ice extent in the Caspian and Aral seas during the last decade (1992-2002). It uses synergy of data from active (radar altimeter) and passive (radiometer) microwave nadir-looking instruments onboard the TOPEX/Poseidon satellite. The proposed approach is substantiated and validated using both in situ and satellite imagery data for the Caspian Sea. The results indicate significant spatial and temporal variability of ice conditions, with a significant decrease of both the duration of ice season and ice extent during the last four winters (1998-2002). The TOPEX/Poseidon-derived time series of sea ice extent are very valuable in view of the fragmentary and mostly unpublished data on ice conditions on the Caspian and Aral seas since the mid-1980s.     

Interaction of Weak Shock Waves with Current Sheets in an Active Region to Produce Nanoflares and Chains of Type I Radio Bursts

Interaction of weak shock waves with a current sheet is investigated by a two-dimensional numerical magnetohydrodynamic model. In accordance with solar coronal conditions, a ratio of thermal to magnetic pressures of 0.1 and a shock Alfvén Mach number slightly above 1 are considered. It is found that even weak shock waves trigger magnetic field reconnection in current sheets. Based on this result, it is suggested that drifting chains of type I radio bursts are radio manifestations of the interactions of weakly super-Alfvénic shock waves with pre-existing current sheets distributed in an active region. This model of type I noise storms is then discussed in connection with the concept of nanoflares (localized reconnections) and the heating of the solar corona.     

Modelling climate change impacts on maize growth and development in the Czech Republic

Summary The crop growth model CERES-Maize is used to estimate the direct (through enhanced fertilisation effect of ambient CO₂) and indirect (through changed climate conditions) effects of increased concentration of atmospheric CO₂ on maize yields. The analysis is based on multi-year crop model simulations run with daily weather series obtained alternatively by a direct modification of observed weather series and by a stochastic weather generator. The crop model is run in two settings: stressed yields are simulated in water and nutrient limited conditions, potential yields in water and nutrient unlimited conditions. The climate change scenario was constructed using the output from the ECHAM3/T42 model (temperature), regression relationships between temperature and solar radiation, and an expert judgement (precipitation). Results: (i) After omitting the two most extreme misfits, the standard error between the observed and modelled yields is 11%. (ii) The direct effect of doubled CO₂: The stressed yields would increase by 36–41% in the present climate and by 61–66% in the 2 × CO₂ climate. The potential yields would increase only by 9–10% as the improved water use efficiency does not apply. (iii) The indirect effect of doubled CO₂: The stressed yields would decrease by 27–29% (14–16%) at present (doubled) ambient CO₂ concentration. The increased temperature shortens the phenological phases and does not allow for the optimal development of the crop. The simultaneous decrease of precipitation and increase of temperature and solar radiation deepen the water stress, thereby reducing the yields. The reduction of the potential yields is significantly smaller as the effect of the increased water stress does not apply. (iv) If both direct and indirect effects of doubled CO₂ are considered, the stressed yields should increase by 17–18%, and the potential yields by 5–14%. (v) The decrease of the stressed yields due to the indirect effect may be reduced by applying earlier planting dates. Received March 9, 2001 Revised September 25, 2001     

The 3-D Truncation Filtering Methodology Defined for Planar and Spherical Models: Interpreting Gravity Data Generated by Point Masses

This paper focuses on one particular way of linear filtering the gravity data to facilitate gravity inversion or interpretation. With the use of integral transforms the gravity anomalies are transformed into new quantities that allow an easier interpretation with the help of pattern recognition. As the integral transforms are in fact filters, and as the regions of integration are caps with a variable radius, which can be systematically changed as a free parameter, we refer to such methodology as the truncation filtering. Such filters may be understood as weighted spherical windows moving over the surface, on which the gravity anomaly is defined, the kernel of the transform being the weight function. The objective of this paper is to define and deploy the truncation filtering for a planar model, i.e. for a homogenous horizontally infinite layer with embedded anomalous masses, and for a spherical model, i.e., for a homogenous massive sphere with embedded anomalous masses. Instead of the original gravity anomaly, the quantities resulting from the truncation filtering are interpreted/inverted. As we shall see, this approach has certain benefits. The fundamental concept of the truncation filtering methodology is demonstrated here in terms of the model consisting of one point mass anomaly.The relationship between the depth of the point mass and the instant of the onset of the dimple pattern observed in sequences produced by truncation filtering the synthetic gravity data generated by point masses is, for both the planar and spherical models, compiled by computer simulations, as well as derived analytically. It is shown, that the dimple pattern is a consequence of truncating the domain of the filter and is free of the choice of the kernel of the filter. It is shown, that in terms of the mean earth and depths of point masses no greater than some 100 km the spherical model may be replaced by a planar model from the perspective of the truncation filtering methodology. It is also shown, that from the viewpoint of the truncation filtering methodology the rigorous gravity anomaly may be approximated by the vertical component of the gravity disturbance. The relationship between the instant of the dimple onset and the depth of the point mass thus becomes linear and independent of the magnitude (mass) of the point mass.     

Induction vector estimates in the Polish-Czechoslovak part of the Carpathians

u¶rt;um ¶rt;a ma ua m u¶rt;uu n mam amu n aum auau ¶rt;a a n -a am anam.