Thermal evolution of the moon

The thermal history and current state of the lunar interior are investigated using constraints imposed by recent geological and physical data. Theoretical temperature models are computed taking into account different initial conditions, heat sources, differentiation and simulated convection. To account for the early formation of the lunar highlands, the time duration of magmatism and presentday temperatures estimated from lunar electrical conductivity profiles, it is necessary to restrict initial temperatures and abundances of radioactivie elements. Successful models require that the outer half of the Moon initially heated to melting temperatures, probably due to rapid accretion. Differentiation of radioactive heat sources toward the lunar surface occurred during the first 1.6 billion years. Temperatures in the outer 500 km are currently low, while the deep interior (radius less than 700 to 1000 km) is warmer than 1000°C, and is of primordial material. In some models there is a partially melted core. The calculated surface heat flux is between 25 and 30 erg/cm² s.Presently at the Research Triangle Institute, Research Triangle, North Carolina 27709, U.S.A.     

Meteoric diagenesis during sea-level lowstands: Evidence from modern hydrochemical studies on northern Guam

Hydrochemical data for meteoric waters of the uplifted carbonate platform of northern Guam show that, contrary to recent models of lowstand diagenesis, phreatic dissolution is active beneath a 60–180 m thick vadose zone. Overland flow during high intensity rainfall events is focussed into surface detention ponds, which drain very rapidly via the epikarst and vertical fissures to the freshwater lens. We estimate that these waters contribute 13% of dissolved calcium in samples from pumping wells but may also deliver aggressive and/or organic-rich waters to the lens. Fast-flow vadose waters mix with lens-top waters to form fresher cap that discharges rapidly coastwards via cavernous porosity. Slower vadose percolation, sampled as cave drips, equilibrates with bedrock within the upper 30 m of the vadose zone, accounting for some 46% of dissolved calcium in the lens waters. The remaining 41% calcium is generated by net dissolution within the lens.     

Indirect approach to invariant point determination for SLR and VLBI systems: an assessment

We assess the accuracy of some indirect approaches to invariant point (IVP), or system reference point, determination of satellite laser ranging (SLR) and very long baseline interferometry (VLBI) systems using both observed and simulated survey data sets. Indirect IVP determination involves the observation of targets located on these systems during specific rotational sequences and by application of geometrical models that describe the target motion during these sequences. Of concern is that most SLR and VLBI systems have limited rotational freedom thereby placing constraint on the reliability of parameter estimation, including the IVP position. We assess two current approaches to IVP analysis using survey data observed at the Yarragadee (Australia) SLR and the Medicina (Italy) VLBI sites and also simulated data of a large rotationally constrained (azimuth-elevation) VLBI system. To improve reliability we introduce and assess some new geometric conditions, including inter-axis, inter-circle and inter-target conditions, to existing IVP analysis strategies. The error component of a local tie specifically associated with the indirect determination of SLR and VLBI IVP is less than 0.5 mm. For systems with significant rotational limits we find that the inter-axis and inter-circle conditions are critical to the computation of unbiased IVP coordinates at the sub-millimetre level. When the inter-axis and inter-circle geometric conditions are not imposed, we retrieve biased vertical coordinates of the IVP (in our simulated VLBI system) in the range of 1.2–3.4 mm. Using the new geometric conditions we also find that the axis-offset estimates can be recovered at the sub- millimetre accuracy (0.5 mm).     

Variability of the Atlantic thermohaline circulation described by three-dimensional empirical orthogonal functions

We describe the use of bivariate three-dimensional empirical orthogonal functions (EOFs) in characterising low frequency variability of the Atlantic thermohaline circulation (THC) in the Hadley Centre global climate model, HadCM3. We find that the leading two modes are well correlated with an index of the meridional overturning circulation (MOC) on decadal timescales, with the leading mode alone accounting for 54% of the decadal variance. Episodes of coherent oscillations in the sub-space of the leading EOFs are identified; these episodes are of great interest for the predictability of the THC, and could indicate the existence of different regimes of natural variability. The mechanism identified for the multi-decadal variability is an internal ocean mode, dominated by changes in convection in the Nordic Seas, which lead the changes in the MOC by a few years. Variations in salinity transports from the Arctic and from the North Atlantic are the main feedbacks which control the oscillation. This mode has a weak feedback onto the atmosphere and hence a surface climatic influence. Interestingly, some of these climate impacts lead the changes in the overturning. There are also similarities to observed multi-decadal climate variability.     

Instantaneous photochemical rates in the global stratosphere

Starting with the average actual distribution of ozone (Dütsch [15]) and temperature in the stratosphere, we have calculated the solar intensity as a function of wavelength and the instantaneous rates (molecules cm^–3 sec^–1) for each Chapman reaction and for each of several reactions of the oxides of nitrogen. The calculation is similar to that ofBrewer andWilson [5]. These reaction rates were calculated independently in each volume element in spherical polar coordinates defined by R=1 km from zero to 50, =5° latitude, and ø=15° longitude (thus including day and night conditions). Calculations were made for two times: summer-winter (January 15) and spring-fall (March 22). As input data we take observed solar intensities (Ackerman [1]) and observed, critically evaluated. constants for elementary chemical and photochemical reactions; no adjustable parameters are employed. (These are not photochemical equilibrium calculations.) According to the Chapman model, the instantaneous, integrated, world-wide rate of formation of ozone from sunlight is about five times faster than the rate of ozone destruction, and locally (lower tropical stratosphere) the rate of ozone formation exceeds the rate of destruction by a factors as great as 1000. The global rates of increase of ozone are more than 50 times faster thanBrewer andWilson's [5] estimate of the average annual transfer rate of ozone to the troposphere. The rate constants of the Chapman reactions are believed to be well-enough known that it is highly improbable that these discrepancies are, due to erroneous rate constants. It is concluded that something else besides neutral oxygen species is very important in stratospheric ozone photochemistry. The inclusion of a uniform concentration of the oxides of nitrogen (NO_x as, NO and NO₂) averaging 6.6×10^–9 mole fraction gives a balance between global ozone formation and destruction rates. The inclusion of a uniform mole fraction of NO_x at 28×10^–9 also gives a global balance. These calculations support the hypethesis (Crutzen [10],Johnston [24]) that the oxides of nitrogen are the most important factor in the global, natural ozone balance. Several authors have recently evaluated the natural source strength of NO_x in the stratosphere; the projected fleets of supersonic transports would constitute an artificial source of NO_x about equal to the natural value, thus promising more or less to double an active natural stratospheric ingredient.