Computation of the elements of eclipsing binary systems in the frequency-domain

The aim of the present paper will be to translate the essential parts of the theory of Fourier analysis of the light changes of eclipsing variables into more practical terms; and describe procedures (illustrated by numerical examples) which should enable their users to obtain the desired results with maximum accuracy and minimum loss of information by processes which can be fully automated.In order to unfold in steps how this can be done, the scope of the present paper-the first of two-will be restricted to an exposition of the analysis of light changes caused by eclipses of spherical stars; while between minima due to this cause the light of the system should remain sensibly constant. An extension of our analysis to incorporate photometric effects arising from mutual distortion of the components of close eclipsing systems between minima as well as within eclipses is being postponed for the second communication.In developing this subject we shall single out for the user's attention only those parts of the whole theory which are of direct relevance to practical work. Their justification can be largely found in sources already published; and new developments essential for our work, not yet made public, will be relegated to several Appendices at the end of the text, in order not to render its text too discursive and deflect the reader's attention from the main theme of its narrative.After a brief outline of the subject given in Section 1, Section 2 will introduce the reader to practical aspects of the Fourier analysis of the light curves; and Section 3 will be devoted to its use to determine the numerical values of the momentsA _2m of the light curves which constitute the cornerstones for all subsequent work. Section 4 will describe an algebraization of the process of determination of the elements for the case of total (annular) eclipses; while Section 5 will do the same for partial eclipses. The concluding Section 6 will be devoted to an error analysis of our problem, and to an outline of the way by which the errors of the individual observations will compound to the uncertainty of the final results. Lastly, Appendices 1–5 concluding the paper will contain additional details of some aspects of our work, or proofs of new processes made use of to obtain our results, whose earlier inclusion would have made the main text too discursive.     

Moments of inertia of the lunar globe,and their bearing on chemical differentiation of its outer layers

It is pointed out that the observed moments of inertia of the Moon, disclosed by its librations, are influenced mainly by the distribution of mass in the outer zone in which the lithostatic pressure is less than 10 kb (i.e., in the outer shell not more than 200 km deep); and a conspicuous departure of such moments from those expected in hydrostatic equilibrium disclosed that these layers could never have been fluid. In the same way, the actual shape of the lunar surface cannot represent a solidified surface of a fluid, petrified at any distance from the Earth.The shape of the Moon, and differences of its moments of inertia must reflect the way in which the initial process of cold accretion fell short of producing a globe with strictly spherically-symmetrical stratification of material; and has nothing to do with tides - present or fossil. Such melting or lava flows as may have occurred at the Moon's surface from time to time must have remained localized, and without much effect on the dynamical properties of the Moon. A global ocean of molten magma some 200 km in depth (postulated sometimes to provide a reservoir in which the differentiation of elements exhibited by surface rocks could have taken place) at any time in the past is incompatible with the dynamical evidence on the motion of the Moon about its center of gravity.Bellcomm, Inc., 955 L'Enfant Plaza North, S.W. Washington, D.C. 20024, U.S.A.     

Special Contribution: Pictorial Series of the Manifestations of the Dynamics of the Earth: 4. Historical Images of Rockbursts' Effects in Mines

Studia Geophysica et Geodaetica -     

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.     

Seismic Activity Around and Under Krakatau Volcano,Sunda Arc: Constraints to the Source Region of Island Arc Volcanics

There is general agreement that calc-alkaline volcanic rocks at convergent plate margins are genetically related to the process of subduction (Ringwood, 1974; Maaloe and Petersen, 1981; Hawkesworth et al., 1997). However, opinions on the mode and site of generation of primary magma for island arc volcanism differ substantially. The site of generation of calc-alkaline magma is thought to be either in the mantle wedge (Plank and Langmuir, 1988; McCulloch and Gamble, 1991) or in the subducting slab (White and Dupré, 1986; Defant and Drummond, 1990; Edwards et al., 1993; Ryan and Langmuir, 1993). We present seismological evidence in favour of the latter concept. A distinctive seismicity pattern around and under the Krakatau volcano was identified during systematic studies of the SE Asian convergent plate margins by means of global seismological data. A column-like cluster of events, probably associated with the dynamics of the volcano, is clearly separated from the events in the Wadati-Benioff zone. The accuracy of hypocentral determinations of the events of the cluster does not differ from the accuracy of the events belonging to the subducting slab. The depths of the cluster events vary from very shallow to about 100 km without any apparent discontinuity. On the other hand, there is a pronounced aseismic gap in the Wadati-Benioff zone directly beneath the volcano at depths between 100-150 km. The Krakatau cluster connects this aseismic gap to the volcano at the surface. The pervasive occurrence of earthquakes in the continental wedge between the subducting slab and the Earth surface bears witness to the brittle character of the continental lithosphere and casts doubt on the existence of large-scale melting of mantle material. The aseismic gap (Hanu and Vank, 1985), interpreted by us as a partially melted domain occurring in subducted slabs in practically all active subduction zones that reach depths greater than 100 km, is here used as evidence for the location of the primary source region of island arc volcanics in the subducting plate.     

Evaluation of pre-CHAMP gravity models `GRIM5-S1' and `GRIM5-C1' with satellite crossover altimetry

 The new GFZ/GRGS gravity field models GRIM5-S1 and GRIM5-C1, currently used as initial models for the CHAMP mission, have been compared with other recent models (JGM 3, EGM 96) for radial orbit accuracy (by means of latitude lumped coefficients) in computations on altimetry satellite orbits. The bases for accuracy judgements are multi-year averages of crossover sea height differences from Geosat and ERS 1/2 missions. This radially sensitive data is fully independent of the data used to develop these gravity models. There is good agreement between the observed differences in all of the world's oceans and projections of the same errors from the scaled covariance matrix of their harmonic geopotential coefficients. It was found that the tentative scale factor of five for the formal standard deviations of the harmonic coefficients of the new GRIM fields is justified, i.e. the accuracy estimates, provided together with the GRIM geopotential coefficients, are realistic. Received: 20 February 2001 / Accepted: 24 October 2001     

A comparative study of geothermal and meteorological records of climate change in Kamchatka

To reconstruct the recent climate history in Kamchatka, a series of repeated precise temperature logs were performed in a number of boreholes located in a broad east-west strip (between 52 and 54°N) in the central part of Kamchatka west of Petropavlovsk-Kamchatski. Within three years more than 30 temperature logs were performed in 10 holes (one up to six logs per hole) to the depth of up to 400 metres. Measured temperature gradients varied in a broad interval 0 to 60 mK/m and in some holes a sizeable variation in the subsurface temperatures due to advective heat transport by underground water was observed. Measured data were compared with older temperature profiles obtained in the early eighties by Sugrobov and Yanovsky (1993). Even when older data are of poorer precision (accuracy of about 0.1 K), they presented valuable information of the subsurface temperature conditions existing 20–25 years ago. Borehole observations and the inverted ground surface temperature histories (GSTHs) used for the paleoclimate reconstruction were complemented with a detailed survey of meteorological data. Namely, the long-term surface air temperature (SAT) and precipitation records from Petropavlovsk station (in operation since 1890) were used together with similar data from a number of local subsidiary meteo-stations operating in Central Kamchatka since 1950. Regardless of extreme complexity of the local meteorological/climate conditions, diversity of borehole sites and calibration of measuring devices used during the whole campaign, the results of the climate reconstruction supported a general warming of about 1 K characteristic for the 20th century, which followed an inexpressive cooler period typical for the most of the 19th century. In the last three to four decades the warming rate has been locally increasing up to 0.02 K/year. It was also shown that the snow cover played a dominant role in the penetration of the climate “signal” to depth and could considerably smooth down the subsurface response to the changes occurred on the surface.     

Explicit formula for the geoid-quasigeoid separation

The explicit formula for the geoid-to-quasigeoid correction is derived in this paper. On comparing the geoidal height and height anomaly, this correction is found to be a function of the mean value of gravity disturbance along the plumbline within the topography. To evaluate the mean gravity disturbance, the gravity field of the Earth is decomposed into components generated by masses within the geoid, topography and atmosphere. Newton’s integration is then used for the computation of topography-and atmosphere-generated components of the mean gravity, while the combined solution for the downward continuation of gravity anomalies and Stokes’ boundary-value problem is utilized in computing the component of mean gravity disturbance generated by mass irregularities within the geoid. On application of this explicit formulism a theoretical accuracy of a few millimetres can be achieved in evaluation of the geoid-to-quasigeoid correction. However, the real accuracy could be lower due to deficiencies within the numerical methods and to errors within the input data (digital terrain and density models and gravity observations).     

Magnitude scale structure from filtered broadbandP-wave records

A series of Hokkaido events, recorded by the FBV Broadband Seismograph System at the KHC Seismic Station, is used to study the structure of the earthquake magnitude scale on the basis of maximum velocity amplitudesA _vmax of teleseismicPwaves in different period bands. Amplitude-periodband (APB) diagrams are constructed for each event. According to the shape of the APB diagrams the events investigated can be divided into three types: (a) events with largestA _vmax values in the intermediate period range (periods ofA _vmax from 2.2 to 23 sec), (b) events with largestA _vmax values in the short-period range (periods ofA _vmax from 1 to 2 sec), (c) events exhibiting anomalous APB diagrams. Type (a) events seem to represent the process of wave generation that prevails for shallow earthquakes. Type (b) events approach to explosive-like generation of seismic waves. The nature of the exceptionally occurring type (c) events must be clarified in further investigations. The influence of the type of earthquake on the magnitude values estimated on the basis of standard class A and B (short-period and intermediate-band) seismograms is demonstrated. It appears that for estimating correct values of earthquake magnitudes complementary information on the process of seismic wave generation in the focus is necessary. At teleseismic distances this information can be obtained from either APB diagrams or amplitude spectra ofP waves recorded, e.g., by broadband velocity sensing instruments.