ALBERT

Description: We have constructed a model that predicts the evolution of CO2 on Mars from the end of the heavy bombardment period to the present. The model draws on published estimates of the main processes believed to affect the fate of CO2 during this period: chemical weathering, regolith uptake, polar cap formation, and atmospheric escape. Except for escape, the rate at which these processes act is controlled by surface temperatures which we calculate using a modified version of the Gierasch and Toon energy balance model. The modifications account for the change in solar luminosity with time, the greenhouse effect, and a polar and equatorial energy budget. Using published estimates for the main parameters, we find no evolutionary scenario in which CO2 is capable of producing a warm (global mean temperatures greater than 250 K) and wet (surface pressures greater than 30 mbar) early climate, and then evolves to present conditions with approximately 7 mbar in the atmosphere, less than 300 mbar in the regolith, and less than 5 mbar in the caps. Such scenarios would only exist if the early sun were brighter than standard solar models suggest, if greenhouse gases other than CO2 were present in the early atmosphere, or if the polar albedo were significantly lower than 0.75. However, these scenarios generally require the storage of large amounts of CO2 (greater than 1 bar) in the carbonate reservoir. If the warm and wet early Mars constraint is relaxed, then we find best overall agreement with present day reservoirs for initial CO2 inventories of 0.5-1.0 bar. We also find that the polar caps can have a profound effect on how the system evolves. If the initial amount of CO2 is less than some critical value, then there is not enough heating of the poles to prevent permanent caps from forming. Once formed, these caps control how the system evolves, because they set the surface pressure and, hence, the thermal environment. If the initial amount of CO2 is greater than this critical value, then caps do not form initially, but can form later on, when weathering and escape lower the surface pressure to a point at which polar heating is no longer sufficient to prevent cap formation and the collapse of the climate system. Our modeling suggests this critical initial amount of CO2 is between 1 and 2 bar, but its true value will depend on all factors affecting the polar heat budget.

Description: The flow field inside a whirling annular seal operating at a Reynolds number of 24,000 and a Taylor number of 6600 has been measured using a 3-D laser Doppler anemometer system. Two eccentricity ratios were considered, 0.10 and 0.50. The seal has a diameter of 164 mm, is 37.3 mm long, and has a clearance of 1.27 mm. The rotor was mounted eccentrically on the shaft such that the whirl ratio is 1.0 and the rotor follows a circular orbit. The mean axial velocity is not uniform around the circumference of the seal; near the inlet a region characterized by high velocity of the seal. By the exit, another region of high axial velocity is not uniform around the circumference of the seal; near the inlet a region characterized by high velocity of the seal. By the exit, another region of high axial velocity has developed, this time on the suction side of the seal. The magnitude and azimuthal distance of the migration increased with increasing whirl amplitude (eccentricity). Throughout the seal length, the azimuthal mean velocity varied inversely with the mean axial velocity. Increasing the whirl amplitude did not increase the magnitude of the azimuthal velocity at the seal exit.

Description: Details of the chondrule and Ca-Al-rich inclusion (CAI) formation during the earliest history of the solar system are imperfectly known. Because CAI's are more 'refractory' than ferromagnesian chondrules and have the lowest recorded initial Sr-87/Sr-86 ratios of any solar system materials, the expectation is that CAI's formed earlier than chondrules. But it is not known, for example, if CAI formation had stopped by the time chondrule formation began. Conventional (absolute) age-dating techniques cannot adequately resolve small age differences (less than 10(exp 6) years) between objects of such antiquity. One approach has been to look at systematic differences in the daughter products of short-lived radionuclides such as Al-26 and I-129. Unfortunately, neither system appears to be 'well-behaved.' One possible reason for this circumstance is that later secondary events have partially reset the isotopic systems, but a viable alternative continues to be large-scale (nebular) heterogeneity in initial isotopic abundances, which would of course render the systems nearly useless as chronometers. In the past two years the nature of this problem has been redefined somewhat. Examination of the Al-Mg isotopic database for all CAI's suggests that the vast majority of inclusions originally had the same initial Al-26/Al-27 abundance ratio, and that the ill-behaved isotopic systematics now observed are the results of later partial reequilibration due to thermal processing. Isotopic heterogeneities did exist in the nebula, as demonstrated by the existence of so-called FUN inclusions in CV3 chondrites and isotopically anomalous hibonite grains in CM2 chondrites, which had little or no live Al-26 at the time of their formation. But, among the population of CV3 inclusions at least, FUN inclusions appear to have been a relatively minor nebular component.

Description: A robot system for use in an enclosed environment was designed and tested. The conceptual design will be used to assist in research performed by the Controlled Ecological Life Support System (CELSS) project. Design specifications include maximum load capacity, operation at specified environmental conditions, low maintenance, and safety. The robot system must not be hazardous to the sealed environment, and be capable of stowing and deploying within a minimum area of the CELSS chamber facility. This design consists of a telescoping robot arm that slides vertically on a shaft positioned in the center of the CELSS chamber. The telescoping robot arm consists of a series of links which can be fully extended to a length equal to the radius of the working envelope of the CELSS chamber. The vertical motion of the robot arm is achieved through the use of a combination ball screw/ball spline actuator system. The robot arm rotates cylindrically about the vertical axis through use of a turntable bearing attached to a central mounting structure fitted to the actuator shaft. The shaft is installed in an overhead rail system allowing the entire structure to be stowed and deployed within the CELSS chamber. The overhead rail system is located above the chamber's upper lamps and extends to the center of the CELSS chamber. The mounting interface of the actuator shaft and rail system allows the entire actuator shaft to be detached and removed from the CELSS chamber. When the actuator shaft is deployed, it is held fixed at the bottom of the chamber by placing a square knob on the bottom of the shaft into a recessed square fitting in the bottom of the chamber floor. A support boot ensures the rigidity of the shaft. Three student teams combined into one group designed a model of the CELSS chamber robot that they could build. They investigated materials, availability, and strength in their design. After the model arm and stand were built, the class performed pre-tests on the entire system. A stability pre-test was used to determine whether the model robot arm would tip over on the stand when it was fully extended. Results showed the stand tipped when 50 Newtons were applied horizontally to the top of the vertical shaft while the arm was fully extended. This proved that it was stable. Another pre-test was the actuator slip test used to determine if there is an adequate coefficient of friction between the actuator drive wheels and drive cable to enable the actuator to fully extend and retract the arm. This pre-test revealed that the coefficient of friction was not large enough to prevent slippage. Sandpaper was glued to the drive wheel and this eliminated the slippage problem. The class preformed a fit test in the CELSS chamber to ensure that the completed robot arm is capable of reaching the entire working envelope. The robot was centered in the chamber and the arm was fully extended to the sides of the chamber. The arm was also able to retract to clear the drain pipes separating the upper and lower plant trays.

Description: We present an analysis of a cusp ion step observed between two poleward-moving events of enhanced ionospheric electron temperature. From the computed variation of the reconnection rate and the onset times of the associated ionospheric events, the distance between the satellite and the X-line can be estimated, but with a large uncertainty due to that in the determination of the low-energy cut-off of the ion velocity distribution function, f(E). Nevertheless, analysis of the time series f(t) shows the reconnection site to be on the dayside magnetopause, consistent with the pulsating cusp model, and the best estimate of the X-line location is 13 R(E) from the satellite. The ion precipitation is used to reconstruct the field-parallel part of the Cowley-D ion distribution function injected into the open low latitude boundary layer (LLBL) in the vicinity of the X-line. From this the Alfven speed, plasma density, magnetic field, parallel ion temperature, and flow velocity of the magnetosheath near the X-line can be derived.

Description: Using a spherically symmetric, self-gravitating, linear viscoelastic Earth model, we predict present-day three-dimensional surface deformation rates and baseline evolutions arising as a consequence of the late Pleistocene glacial cycles. In general, we use realistic models for the space-time geometry of the final late Pleistocene deglaciation event and incorporate a gravitationally self-consistent ocean meltwater redistribution. The predictions of horizontal velocity presented differ significantly, in both their amplitude and their spatial variation, from those presented in earlier analysis of others which adopted simplified models of both the late Pleistocene ice history and the Earth rheology. An important characteristic of our predicted velocity fields is that the melting of the Laurentide ice sheet over Canada is capable of contributing appreciably to the adjustment in Europe. The sensitivity of the predictions to variations in mantle rheology is investigated by considering a number of different Earth models, and by computing appropriate Frechet kernels. These calculations suggest that the sensitivity of the deformations to the Earth's rheology is significant and strongly dependent on the location of the site relative to the ancient ice sheet. The effects on the predictions of three-dimensional deformation rates of altering the ice history or adopting approximate models for the ocean meltwater redistribution have also been considered and found to be important (the former especially so). Finally, for a suite of Earth models we provide predictions of the velocity of a number of baselines in North America and Europe. We find that, in general, both radial and tangential motions contribute significantly to baseline length changes, and that these contributions are a strong function of the Earth model. We have, furthermore, found a set of Earth models which, together with the ICE-3G deglaciation chronology, produce predictions of baseline length changes that are consistent with very long baseline interferometry measurements of baselines within Europe.

Description: We outline a complete spectral formalism for computing high spatial resolution three-dimensional deformations arising from the surface mass loading of a spherically symmetric planet. The main advantages of the formalism are that all surface mass loads are always described using a consistent mathematical representation and that calculations of deformation fields for various spatial resolutions can be performed by simpley altering the spherical harmonic degree truncation level of the procedure. The latter may be important when incorporating improved observational constraints on a particular surface mass load, when considering potential errors in the computed field associated with mass loading having a spatial scale unresolved by the observational constraints, or when treating a number of global surface mass loads constrained with different spatial resolutions. The advantages do not extend to traditional 'Green's function' approaches which involve surface element discretizations of the global mass loads. Another advantage of the spectral formalism, over the Green's function approach, is that a posteriori analyses of the computed deformation fields are easily performed. In developing the spectral formalism, we consider specific cases where the Earth's mantle is assumed to respond as an elastic, slightly anelastic, or linear viscoelastic medium. In the case of an elastic or slightly anelastic mantle rheology the spectral response equations incorporate frequency dependent Love numbers. The formalism can therefore be used, for example, to compute the potentially resonant deformational response associated with the free core nutation and Chandler wobble eigenfunctions. For completeness, the spectral response equations include both body forces, as arise from the gravitational attraction of the Sun and the Moon, and surface mass loads. In either case, and for both elastic and anelastic mantle rheologies, we outline a pseudo-spectral technique for computing the ocean adjustment associated with the total gravitational perturbation induced by the external forcing. Three-dimensional deformations computed using the usual Love number approach are generally referenced to an origin at the center of mass of the undeformed planet. We derive a spectral technique for transforming the results to an origin located at the center of mass of the deformed planet.

Description: We use over a decade of geodetic Very Long Baseline Interferometry (VLBI) data to estimate parameters in a resonance expansion of the frequency dependence of the tidal h(sub 2) Love number within the diurnal band. The resonance is associated with the retrograde free core nutation (RFCN). We obtain a value for the real part of the resonance strength of (-0.27 +/- 0.03) x 10(exp -3); a value of -0.19 x 10(exp -3) is predicted theoretically. Uncertainties in the VLBI estimates of the body tide radial displacement amplitudes are approximately 0.5 mm (1.1 mm for the K1 frequency), but they do not yield sufficiently small Love number uncertainties for placing useful constraints on the frequency of the RFCN, given the much smaller uncertainties obtained from independent analyses using nutation or gravimetric data. We also consider the imaginary part of the tidal h(sub 2) Love number. The estimated imaginary part of the resonance strength is (0.00 +/- 0.02) x 10(exp -3). The estimated imaginary part of the nonresonant component of the Love number implies a phase angle in the diurnal tidal response of the Earth of 0.7 deg +/- 0.5 deg (lag).

Description: The flow field inside a whirling annular seal has been measured using a 3-D Laser Doppler Anemometer (LDA) system. The seal investigated has a clearance of 1.27 mm, a length of 37.3 mm and is mounted on a drive shaft with a 50 percent eccentricity ratio. This results in the rotor whirling at the same speed as the shaft rotation (whirl ratio = 1.0). The seal is operated at a Reynolds number of 12000 and a Taylor number of 6300 (3600 rpm). The 3-D LDA system is equipped with a rotary encoding system which is used to produce phase averaged measurements of the entire mean velocity vector field and Reynolds stress tensor field from 0.13 mm upstream to 0.13 mm downstream of the seal. The mean velocity field reveals a highly three-dimensional flow field with large radial velocities near the inlet of the seal as well as a recirculation zone on the rotor surface. The location of maximum mean axial velocity migrates from the pressure side of the rotor at the inlet to the suction side at the exit. Turbulence production is a maximum near the seal inlet as indicated by the rapid increase of the turbulence kinetic energy (kappa). However, turbulence production and dissipation attain equilibrium fairly quickly with kappa remaining relatively constant over the last half of the seal.

Description: Analysis of Global Positioning System (GPS) data from two sites separated by horizontal distance of only approximately 2.2 m yielded phase residuals exhibiting a systematic elevation angle dependence. One of the two GPS antennas was mounted on an approximately 1 m high concrete pillar, and the other was mounted on a standard wooden tripod. We performed elevation angle cutoff tests with these data, and established that the vertical coordinate of site position was sensitive to the minimum elevation angle (elevation cutoff) of the data analyzed. For example, the vertical coordinate of site position changed by 9.7 plus or minus 0.8 mm when the minimum elevation angle was increased from 10 to 25. We performed simulations based on a simple (ray tracing) multipath model with a single horizontal reflector, and demonstrated that the elevation angle cutoff test results and the pattern of the residual versus elevation angle could be qualitatively reproduced if the reflector were located 0.1-0.2 m beneath the antenna phase center. We therefore, hypothesized that the source of the elevation-angle-dependent error were multipath reflections and scattering and that the horizontal surface of the pillar, located a distance of approximately 0.2 m beneath the antenna phase center, was the primary reflector. We tested this hypothesis by placing microwave absorbing material between the antenna and the pillar in a number of configurations and analyzed the changes in apparent position of the antenna. The results indicate that (1) the horizontal surface of the pillar is indeed the main reflector, (2) both the concrete and the metal plate embedded in the pillar are significant reflectors, and (3) the reflection can be reduced to a great degree by the use of microwave absorbing materials. These results have significant implications for the accuracy of global GPS geodetic tracking networks which use pillar-antenna configuration identical or similar to the one used here (at the Westford WFRD GPS site).