The SuperCOSMOS Sky Survey – III. Astrometry

In this, the third in a series of three papers concerning the SuperCOSMOS Sky Survey, we describe the astrometric properties of the data base. We describe the algorithms employed in the derivation of the astrometric parameters of the data, and demonstrate their accuracies by comparison with external data sets using the first release of data, the South Galactic Cap survey. We show that the celestial coordinates, which are tied to the International Celestial Reference Frame via the Tycho–2 reference catalogue, are accurate to better than ±0.2 arcsec at J , R ∼19,18 , rising to ±0.3 arcsec at J , R ∼22,21 , with positional-dependent systematic effects from bright to faint magnitudes at the ∼0.1-arcsec level. The proper motion measurements are shown to be accurate to typically ±10 mas yr⁻¹ at J , R ∼19,18 , rising to ±50 mas yr⁻¹ at J , R ∼22,21 , and are tied to zero using the extragalactic reference frame. We show that the zero-point errors in the proper motions are ≤1 mas yr⁻¹ for R >17 , and are no larger than ∼10 mas yr⁻¹ for R <17 mas yr⁻¹ .     

Parallax distortion by the weak microlensing effect

Parallax measurements allow distances to celestial objects to be determined. Together with measurements of their position on the celestial sphere, they give a full three-dimensional picture of the location of the objects relative to the observer. The distortion of the parallax value of a distant source affected by weak microlensing is considered. This means that the weak microlensing leads to distortion of the distance scale. The gravitational deflection causes a retrograde apparent motion of the image, which is revealed as a negative parallax. It is shown that the distortions may significantly change the parallax values when they amount to several microseconds of arc. In particular, at this level many measured values of parallaxes should be negative.     

APPLICATION OF MAP AND FILE INFORMATIONVISUALIZATION SYSTEM TO COMPREHENSIVEDIVISION OF NATURAL DISASTERS- Takingthe ChangjiangValleg as an Example

In recent years，scholars at home and abroad have method，and that vanous disasters would be analyzed assynthetically studied natural dlsaste‘theoretically and an Integral．SHI Peilun（1991）putforwad a scientificmethodologlcally，as well as its cases analys。s（CND，term—regional disaster system，which Indicates that1987； MA et al，1990； MARBLE，1990； NE et al，the situation of a disaster（calamity loss）results from1999； PATAK et al，1982； SHI，1991；VAN et al，h…     

Report of the International Astronomical Union Division I Working Group on Precession and the Ecliptic

The IAU Working Group on Precession and the Equinox looked at several solutions for replacing the precession part of the IAU 2000A precession–nutation model, which is not consistent with dynamical theory. These comparisons show that the (Capitaine et al., Astron. Astrophys., 412, 2003a) precession theory, P03, is both consistent with dynamical theory and the solution most compatible with the IAU 2000A nutation model. Thus, the working group recommends the adoption of the P03 precession theory for use with the IAU 2000A nutation. The two greatest sources of uncertainty in the precession theory are the rate of change of the Earth’s dynamical flattening, ΔJ₂, and the precession rates (i.e. the constants of integration used in deriving the precession). The combined uncertainties limit the accuracy in the precession theory to approximately 2 mas cent⁻². Given that there are difficulties with the traditional angles used to parameterize the precession, z_A, ζ_A, and θ_A, the working group has decided that the choice of parameters should be left to the user. We provide a consistent set of parameters that may be used with either the traditional rotation matrix, or those rotation matrices described in (Capitaine et al., Astron. Astrophys., 412, 2003a) and (Fukushima Astron. J., 126, 2003). We recommend that the ecliptic pole be explicitly defined by the mean orbital angular momentum vector of the Earth–Moon barycenter in the Barycentric Celestial Reference System (BCRS), and explicitly state that this definition is being used to avoid confusion with previous definitions of the ecliptic. Finally, we recommend that the terms precession of the equator and precession of the ecliptic replace the terms lunisolar precession and planetary precession, respectively.     

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.     

Fluid infiltration during contact metamorphism of interbedded marble and calc-silicate hornfels, Twin Lakes area, central Sierra Nevada, California

Abstract In the Twin Lakes area, central Sierra Nevada, California, most contact metamorphosed marbles contain calcite + dolomite + forsterite ± diopside ± phlogopite ± tremolite, and most calc-silicate hornfelses contain calcite + diopside + wollastonite + quartz ± anorthite ± K-feldspar ± grossular ± titanite. Mineral-fluid equilibria involving calcite + dolomite + tremolite + diopside + forsterite in two marble samples and wollastonite + anorthite + quartz + grossular in three hornfels samples record P± 3 kbar and T± 630° C. Various isobaric univariant assemblages record CO₂-H₂O fluid compositions of χCO₂= 0.61–0.74 in the marbles and χCO₂= 0.11 in the hornfelses. Assuming a siliceous dolomitic limestone protolith consisting of dolomite + quartz ° Calcite ± K-feldspar ± muscovite ± rutile, all plausible prograde reaction pathways were deduced for marble and hornfels on isobaric T-XCO₂ diagrams in the model system K₂O-CaO-MgO-Al₂O₃-SiO₂-H₂O-CO₂. Progress of the prograde reactions was estimated from measured modes and mass-balance calculations. Time-integrated fluxes of reactive fluid which infiltrated samples were computed for a temperature gradient of 150 °C/km along the fluid flow path, calculated fluid compositions, and estimated reaction progress using the mass-continuity equation. Marbles and hornfelses record values in the range 0.1–3.6 × 10⁴ cm³/cm² and 4.8–12.9 × 10⁴ cm³/cm², respectively. For an estimated duration of metamorphism of 10⁵ years, average in situ metamorphic rock permeabilities, calculated from Darcy's Law, are 0.1–8 × 10^?6 D in the marbles and 10–27 × 10^?6 D in the hornfelses. Reactive metamorphic fluids flowed up-temperature, and were preferentially channellized in hornfelses relative to the marbles. These results appear to give a general characterization of hydrothermal activity during contact metamorphism of small pendants and screens (dimensions ± 1 km or less) associated with emplacement of the Sierra Nevada batholith.     

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.     

Electric Discharges in the Protoplanetary Nebula as a Source of Impulse Magnetic Fields to Promote Dust Aggregation

Lightning discharge generated in the protoplanetary nebula is viewed as a temporally isolated surge in the flow of electrically charged particles, similar to that of terrestrial lightning. If the current is intense enough, a powerful circular impulse magnetic field is generated around the instantaneous virtual electric conductor. Such magnetic field is capable of magnetizing dust grains containing ferromagnetic components present in its vicinity to their saturation levels. As a result, dust grains attract one another, forming the aggregates. This magnetically driven attraction suggests an important process possibly operational at an early stage of the planetary accretion. Based on both a classical model for electric conductor, and the theory of Lienard–Wiechert electromagnetic potentials, our calculations show that the magnetic impulse due to a discharge channel of a few cm in diameter transferring a charge of about 10⁴ electrons reaches as high as 10 T. At these magnetic fields, the ferromagnetic dust grains, and possibly the already-formed larger aggregates as well, are easily magnetized to the saturation levels, producing compact clusters exhibiting permanent magnetic moments.     

Optimal Compensation and Implementation for Adaptive Optics Systems

This paper develops a compensation algorithm based on Linear–Quadratic–Gaussian (LQG) control system design whose parameters are determined (in part) by a model of the atmosphere. The model for the atmosphere is based on the open-loop statistics of the atmosphere as observed by the wavefront sensor, and is identified from these using an auto-regressive, moving average (ARMA) model. The (LQG) control design is compared with an existing compensation algorithm for a simulation developed at ESO that represents the operation of MACAO adaptive optics system on the 8.2 m telescopes at Paranal, Chile. This revised version was published online in July 2006 with corrections to the Cover Date.