Dynamical History of the Earth–Moon System

Differential equations describing the tidal evolution of the earth's rotation and of the lunar orbital motion are presented in a simple close form. The equations differ in form for orbits fixed to the terrestrial equator and for orbits with the nodes precessing along the ecliptic due to solar perturbations. Analytical considerations show that if the contemporary lunar orbit were equatorial the evolution would develop from an unstable geosynchronous orbit of the period about 4.42 h (in the past) to a stable geosynchronous orbit of the period about 44.8 days (in the future). It is also demonstrated that at the contemporary epoch the orbital plane of the fictitious equatorial moon would be unstable in the Liapunov's sense, being asymptotically stable at early stages of the evolution. Evolution of the currently near-ecliptical lunar orbit and of the terrestrial rotation is traced backward in time by numerical integration of the evolutional equations. It is confirmed that about 1.8 billion years ago a critical phase of the evolution took place when the equatorial inclination of the moon reached small values and the moon was in a near vicinity of the earth. Before the critical epoch t _cr two types of the evolution are possible, which at present cannot be unambiguously distinguished with the help of the purely dynamical considerations. In the scenario that seems to be the most realistic from the physical point of view, the evolution also has started from a geosynchronous equatorial lunar orbit of the period 4.19 h. At t < t _cr the lunar orbit has been fixed to the precessing terrestrial equator by strong perturbations from the earth's flattening and by tidal effects; at the critical epoch the solar perturbations begin to dominate and transfer the moon to its contemporary near-ecliptical orbit which evolves now to the stable geosynchronous state. Probably this scenario is in favour of the Darwin's hypothesis about originating the moon by its separation from the earth. Too much short time scale of the evolution in this model might be enlarged if the dissipative Q factor had somewhat larger values in the past than in the present epoch. Values of the length of day and the length of month, estimated from paleontological data, are confronted with the results of the developed model.     

Recharge and salinization in the Madrid Basin

The purpose of this work was to reinvestigate the existing hydrogeological conceptual model of the basin of Madrid, Spain. A cumulative chemical isotopic diagram which enabled the distinction between different groups of water as well as calculation of the mode of their blending was used for this investigation. It was found that the groups of discharge were lighter in their isotopic composition than that of recharge. The previous explanation of this fact, backed by carbon-14 dating, was the long residence time due to flow lines going down to depths of more than 1000 m. This flow model assumes homogenous conditions to these depths. This assumption can not be supported by evidence from deep wells. Thus a modified model is suggested which maintains homogenous conditions only to about 300 m and a deep confined aquifer below containing paleowater. The higher degree of depletion of this water has been explained by a colder climate on top of an altitude effect. Another interesting observation was the correlation between the isotopic composition of the rains, the month of the rain event and the composition of the recharge group groundwater. It could be seen that the winter rains resemble the groundwater composition, which shows that practically all the spring and summer rains were evapotranspirated.     

Numerical simulations of high-speed solar wind streams within 1 AU and their signatures at 1 AU

A parametric study of the evolution within, and signatures at, 1 AU of high-speed streams is performed with the use of a MHD, 21/2-D, time-dependent model. This study is an extension of an earlier one by Smith and Dryer (1990) who examined the ecliptic plane consequences of relatively short-duration, energetic solar disturbances. The present study examines both the erupting and corotating parts of long-duration, high-speed streams characteristic of coronal hole flows. By examining the variation of the simulated plasma velocity, density, temperature, and magnetic field at 1 AU, as well as the location of the solar coronal hole sources relative to the observer at 1 AU, we are able to provide some insight into the identification of the solar sources of interplanetary disturbances. We present and discuss two definitions for angle locating the solar source of interplanetary disturbances at 1 AU.We apply our results to the suggestion by Hewish (1988) that low-latitude coronal holes are suitably positioned to be the sources of major geomagnetic storms when the holes are in the eastern half of the solar hemisphere at the time of the commencement of the storm. Our results indicate that, for these cases, the streams emanating from within the hole must be very fast, greater than 1000 km s^–1, or very wide, greater than 60°, at the inner boundary of 18 solar radii in our simulation.     

The stable isotope signature of kilometre-scale fracturedominated metamorphic fluid pathways, Mary Kathleen, Australia

Abstract Large calcite veins and pods in the Proterozoic Corella Formation of the Mount Isa Inlier provide evidence for kilometre-scale fluid transport during amphibolite facies metamorphism. These 10- to 100-m-scale podiform veins and their surrounding alteration zones have similar oxygen and carbon isotopic ratios throughout the 200 × 10-km Mary Kathleen Fold Belt, despite the isotopic heterogeneity of the surrounding wallrocks. The fluids that formed the pods and veins were not in isotopic equilibrium with the immediately adjacent rocks. The pods have δ¹³C_calcite values of –2 to –7% and δ¹⁸O_calcite values of 10.5 to 12.5%. Away from the pods, metadolerite wallrocks have δ¹⁸O_whole-rock values of 3.5 to 7%. and unaltered banded calc-silicate and marble wallrocks have δ¹³C_calcite of –1.6 to –0.6%, and δ¹⁸O_calcite of 18 to 21%. In the alteration zones adjacent to the pods, the δ¹⁸O values of both metadolerite and calc-silicate rocks approach those of the pods. Large calcite pods hosted entirely in calc-silicates show little difference in isotopic composition from pods hosted entirely in metadolerite. Thus, 100- to 500-m-scale isotopic exchange with the surrounding metadolerites and calc-silicates does not explain the observation that the δ¹⁸O values of the pods are intermediate between these two rock types. Pods hosted in felsic metavolcanics and metasiltstones are also isotopically indistinguishable from those hosted in the dominant metadolerites and calc-silicates. These data suggest the veins are the product of infiltration of isotopically homogeneous fluids that were not derived from within the Corella Formation at the presently exposed crustal level, although some of the spread in the data may be due to a relatively small contribution from devolatilization reactions in the calc-silicates, or thermal fluctuations attending deformation and metamorphism. The overall L-shaped trend of the data on plots of δ¹³C vs. δ¹⁸O is most consistent with mixing of large volumes of externally derived fluids with small volumes of locally derived fluid produced by devolatilization of calc-silicate rocks. Localization of the vein systems in dilatant sites around metadolerite/calc-silicate boundaries indicates a strong structural control on fluid flow, and the stable isotope data suggest fluid migration must have occurred at scales greater than at least 1 km. The ultimate source for the external fluid is uncertain, but is probably fluid released from crystallizing melts derived from the lower crust or upper mantle. Intrusion of magmas below the exposed crustal level would also explain the high geothermal gradient calculated for the regional metamorphism.     

A Critical review of theoretical models of negatively polarized light scattered by atmosphereless solar system bodies

About a dozen physical mechanisms and models aspire to explain the negative polarization of light scattered by atmosphereless celestial bodies. This is too large a number for the reliable interpretation of observational data. Through a comparative analysis of the models, our main goal is to answer the question: Does any one model have an advantage over the others? Our analysis is based on new laboratory polarimetric and photometric data as well as on theoretical results. We show that the widely used models due to Hopfield and Wolff cannot realistically explain the phase-angle dependence of the degree of polarization observed at small phase angles. The so-called interference or coherent backscattering mechanism is the most promising model. Models based on that mechanism use well-defined physical parameters to explain both negative polarization and the opposition effect. They are supported by laboratory experiments, particularly those showing enhancement of negative polarization with decreasing particle size down to the wavelength of light. According to the interference mechanism, pronounced negative branches of polarization, like those of C-class asteroids, may indicate a high degree of optical inhomogeneity of light-scattering surfaces at small scales. The mechanism also seems appropriate for treating the negative polarization and opposition effects of cometary dust comae, planetary rings, and the zodiacal light.     

Three-dimensional calculations of high- and low-mass planets embedded in protoplanetary discs

We analyse the non-linear, three-dimensional response of a gaseous, viscous protoplanetary disc to the presence of a planet of mass ranging from 1 Earth mass (1 M_⊕) to 1 Jupiter mass (1 M_J) by using the zeus hydrodynamics code. We determine the gas flow pattern, and the accretion and migration rates of the planet. The planet is assumed to be in a fixed circular orbit about the central star. It is also assumed to be able to accrete gas without expansion on the scale of its Roche radius. Only planets with masses   M _p≳ 0.1 M_J  produce significant perturbations in the surface density of the disc. The flow within the Roche lobe of the planet is fully three-dimensional. Gas streams generally enter the Roche lobe close to the disc mid-plane, but produce much weaker shocks than the streams in two-dimensional models. The streams supply material to a circumplanetary disc that rotates in the same sense as the orbit of the planet. Much of the mass supply to the circumplanetary disc comes from non-coplanar flow. The accretion rate peaks with a planet mass of approximately 0.1 M_J and is highly efficient, occurring at the local viscous rate. The migration time-scales for planets of mass less than 0.1 M_J, based on torques from disc material outside the Roche lobes of the planets, are in excellent agreement with the linear theory of type I (non-gap) migration for three-dimensional discs. The transition from type I to type II (gap) migration is smooth, with changes in migration times of about a factor of 2. Starting with a core which can undergo runaway growth, a planet can gain up to a few M_J with little migration. Planets with final masses of the order of 10 M_J would undergo large migration, which makes formation and survival difficult.     

Star-Driven Wind and Jet Models

We outline some main results from recent analytical modelling of axisymmetric jets from the coronae of young stars and compare them to disk-wind and X-wind models. We emphazise the roles of the magnetic rotator and the disk in the formation and the evolution of the jet. We conjecture that with time both the efficiency of the magnetic rotator and the role of the disk as a primary source for the wind decline.