ALBERT

Description: Submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy at the Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution June, 1980

Description: This thesis consists of three papers applying the techniques of time series analysis to geophysical data. Surface wave dispersion along the Walvis Ridge, South Atlantic Ocean, is obtained by bandpass filtering the recorded seismogram in the frequency domain. The group velocity is anomalously low in the period range of 15-50 s, and formal inversion of the data indicates both crustal thickening to 12.5 km and low shear velocity (4.25-4.35 km/s) to depths of 40-50 km. The electromagnetic induction fields at a deep ocean site northeast of Hawaii were used to determine the electrical conductivity of the earth to 400 km depth. Singular value decomposition of the data matrix indicates three degrees of freedom, suggesting source field complications and a two dimensional conductive structure. Inversion of one of the principal terms in the response function shows an abrupt rise in electrical conductivity to 0.05 mho/m near 160 km with no resolvable decrease below this. A model study suggests that moving source fields influence the induction appreciably in the other principal response tunction. A set of piston cores from the northeast Atlantic Ocean are used to construct paleomagnetic time series covering the interval 25-127 kybp. Stratigraphic control is provided by counts of planktonic toraminifera, and empirical orthogonal function analysis shows a significant decrease in sedimentation rate at the interglaciai/glacial transition. The sediments are magnetically stable and reliable relative paleointensity measurements could be obtained. Spectral analysis of the directions reveals a predominant 10 ky periodicity and no dominant looping direction.

Description: I was supported for the early parts of this work by a NSF Graduate Fellowship. The Walvis Ridge study was supported by the WHOI Education Office and the Defense Advanced Research Projects Agency. The induction study was funded by the NSF under grants OCE74-12730 and OCE77-8633, and by the WHOI Ocean Industries Program. The paleomagnetic study was supported by NSF contracts OCE77-82255 and ÖCE79-19258.

Description: Author Posting. © Elsevier B.V., 2007. This is the author's version of the work. It is posted here by permission of Elsevier B.V. for personal use, not for redistribution. The definitive version was published in Icarus 189: 72-82, doi:10.1016/j.icarus.2007.01.010.

Description: Recent observations of the south pole of Saturn’s moon Enceladus by the Cassini spacecraft have revealed an active world, powered by internal heat. In this paper, we propose that localized subsurface melting on Enceladus has produced an internal south polar sea. Evidence for this localized sea comes from the shape of Enceladus, which does not match a differentiated body at its current orbital position. We show that melting induced by the observed heat flow at the south pole produces a large enough pit to match the shape of Enceladus with a differentiated rock and ice interior. Numerical modeling of melting and ice flow shows that the sea produced beneath the south pole is stable against inflow of ductile ice from its surroundings for the duration of the heating. The shape modification due to melting also produces a negative degree-two gravity anomaly, which can reorient the spin axis of Enceladus in order to place the sea at the pole.

Description: There is a growing number of observational evidences of dynamic quasi-periodical magnetosphere response to continuously southward interplan etary magnetic field (IMF). However, traditional global MHD simulatio ns with magnetic reconnection supported by numerical dissipation and ad hoc anomalous resistivity driven by steady southward IMF often prod uce only quasi-steady configurations with almost stationary near-eart h neutral line. This discrepancy can be explained by the assumption that global MHD simulations significantly underestimate the reconnectio n rate in the magnetotail during substorm expansion phase. Indeed, co mparative studies of magnetic reconnection in small scale geometries demonstrated that traditional resistive MHD did not produce the fast r econnection rates observed in kinetic simulations. The major approxim ation of the traditional MHD approach is an isotropic fluid assumption) with zero off-diagonal pressure tensor components. The approximatio n, however, becomes invalid in the diffusion region around the reconn ection site where ions become unmagnetized and experience nongyrotropic behaviour. Deviation from gyrotropy in particle distribution functi on caused by kinetic effects manifests itself in nongyrotropic pressu re tensor with nonzero off-diagonal components. We use the global MHD code BATS-R-US and replace ad hoc parameters such as "critical curren t density" and "anomalous resistivity" with a physically motivated di ssipation model. The key element of the approach is to identify diffusion regions where the isotropic fluid MHD approximation is not applic able. We developed an algorithm that searches for locations of magnet otail reconnection sites. The algorithm takes advantage of block-based domain-decomposition technique employed by the BATS-R-US. Boundaries of the diffusion region around each reconnection site are estimated from the gyrotropic orbit threshold condition, where the ion gyroradius is equal to the distance to the reconnection site. Inside diffusion regions ions are treated as nongyrotropic fluid with nonzero off-dia gonal components of the pressure tensor. The primary kinetic mechanism controlling the dissipation in the diffusion region is incorporated into global MHD simulations in terms of spatially localized nongyrotropic corrections to the induction equation. The magnitude of the non-g yrotropic corrections to the electric field and spatial scales of the diffusion regions are calculated self-consistently at each time step of the simulation using local MHD plasma and field parameters at the reconnection site without introduction of any ad hoc parameters. We d emonstrated that magnetotail reconnection is inherently unsteady even when the solar wind is steady. Global MHD simulations with nongyrotropic corrections produce bursts of fast reconnection typically observe d in small-scale kinetic simulations. During the bursts the length of the diffusion region does not exceed 2R(sub E) approximates 12(c/ome ga * pi). The bursts of the fast reconnection last only for a few min utes. After reaching the maximum value the reconnection rate decreases while the length of the diffusion region increases. The decreased ra te, however, is still significantly larger that the steady reconnection rate characteristic for MHD simulations with reconnection supported by numerical resistivity alone. Magnetotail reconnection supported b y nongyrotropic effects results in a tailward retreat of the reconnection site with average speed of the order of 100 km/s, accompanied by magnetotail stretching and thin current sheet formation in the near-E arth plasma sheet. Overall magnetotail response to the steady low-mach-number solar wind with southward IMF exhibits quasi-periodic loading /unloading dynamics typical for frequently observed multiple substorm s.

Description: This paper presents an approach for slip prediction from a distance for wheeled ground robots using visual information as input. Large amounts of slippage which can occur on certain surfaces, such as sandy slopes, will negatively affect rover mobility. Therefore, obtaining information about slip before entering such terrain can be very useful for better planning and avoiding these areas. To address this problem, terrain appearance and geometry information about map cells are correlated to the slip measured by the rover while traversing each cell. This relationship is learned from previous experience, so slip can be predicted remotely from visual information only. The proposed method consists of terrain type recognition and nonlinear regression modeling. The method has been implemented and tested offline on several off-road terrains including: soil, sand, gravel, and woodchips. The final slip prediction error is about 20%. The system is intended for improved navigation on steep slopes and rough terrain for Mars rovers.

Description: Amino acids are considered organic molecular indicators in the search for extant and extinct life in the Solar System. Extraction of these molecules from a particulate solid matrix, such as Martian regolith, will be critical to their in situ detection and analysis. The goals of this study were to optimize a laboratory amino acid extraction protocol by quantitatively measuring the yields of extracted amino acids as a function of liquid water temperature and sample extraction time and to compare the results to the standard HCl vapor- phase hydrolysis yields for the same soil samples. Soil samples from the Yungay region of the Atacama Desert ( Martian regolith analog) were collected during a field study in the summer of 2005. The amino acids ( alanine, aspartic acid, glutamic acid, glycine, serine, and valine) chosen for analysis were present in the samples at concentrations of 1 - 70 parts- per- billion. Subcritical water extraction efficiency was examined over the temperature range of 30 - 325 degrees C, at pressures of 17.2 or 20.0 MPa, and for water- sample contact equilibration times of 0 - 30 min. None of the amino acids were extracted in detectable amounts at 30 degrees C ( at 17.2 MPa), suggesting that amino acids are too strongly bound by the soil matrix to be extracted at such a low temperature. Between 150 degrees C and 250 degrees C ( at 17.2 MPa), the extraction efficiencies of glycine, alanine, and valine were observed to increase with increasing water temperature, consistent with higher solubility at higher temperatures, perhaps due to the decreasing dielectric constant of water. Amino acids were not detected in extracts collected at 325 degrees C ( at 20.0 MPa), probably due to amino acid decomposition at this temperature. The optimal subcritical water extraction conditions for these amino acids from Atacama Desert soils were achieved at 200 degrees C, 17.2 MPa, and a water- sample contact equilibration time of 10 min.

Description: The Orbiter radiator system consists of eight individual 4.6 m x 3.2 m panels located with four on each payload bay door. Forward panels #1 and #2 are 2.3 cm thick while the aft panels #3 and #4 have a smaller overall thickness of 1.3 cm. The honeycomb radiator panels consist of 0.028 cm thick Aluminum 2024-T81 facesheets and Al5056-H39 cores. The face-sheets are topped with 0.005 in. (0.127 mm) silver-Teflon tape. The radiators are located on the inside of the shuttle payload bay doors, which are closed during ascent and reentry, limiting damage to the on-orbit portion of the mission. Post-flight inspections at the Kennedy Space Center (KSC) following the STS-115 mission revealed a large micrometeoroid/orbital debris (MMOD) impact near the hinge line on the #4 starboard payload bay door radiator panel. The features of this impact make it the largest ever recorded on an orbiter payload bay door radiator. The general location of the damage site and the adjacent radiator panels can be seen in Figure 2. Initial measurements of the defect indicated that the hole in the facesheet was 0.108 in. (2.74 mm) in diameter. Figure 3 shows an image of the front side damage. Subsequent observations revealed exit damage on the rear facesheet. Impact damage features on the rear facesheet included a 0.03 in. diameter hole (0.76 mm), a approx.0.05 in. tall bulge (approx.1.3 mm), and a larger approx.0.2 in. tall bulge (approx.5.1 mm) that exhibited a crack over 0.27 in. (6.8 mm) long. A large approx.1 in. (25 mm) diameter region of the honeycomb core was also damaged. Refer to Figure 4 for an image of the backside damage to the panel. No damage was found on thermal blankets or payload bay door structure under the radiator panel. Figure 5 shows the front facesheet with the thermal tape removed. Ultrasound examination indicated a maximum facesheet debond extent of approximately 1 in. (25 mm) from the entry hole. X-ray examinations revealed damage to an estimated 31 honeycomb cells with an extent of 0.85 in. x 1.1 in. (21.6 x 27.9 mm). Pieces of the radiator at and surrounding the impact site were recovered during the repair procedures at KSC. They included the thermal tape, front facesheet, honeycomb core, and rear facesheet. These articles were examined at JSC using a scanning electron microscope (SEM) with an energy dispersive x-ray spectrometer (EDS). Figure 6 shows SEM images of the entry hole in the facesheet. The asymmetric height of the lip may be attributed to projectile shape and impact angle. Numerous instances of a glass-fiber organic matrix composite were observed in the facesheet tape sample. The fibers were approximately 10 micrometers in diameter and variable lengths. EDS analysis indicated a composition of Mg, Ca, Al, Si, and O. Figures 7 and 8 present images of the fiber bundles, which were believed to be circuit board material based on similarity in fiber diameter, orientation, consistency, and composition. A test program was initiated in an attempt to simulate the observed damage to the radiator facesheet and honeycomb. Twelve test shots were performed using projectiles cut from a 1.6 mm thick fiberglass circuit board substrate panel. Results from test HITF07017, shown in figures 9 and 10, correlates with the observed impact features reasonably well. The test was performed at 4.14 km/sec with an impact angle of 45 degrees using a cylindrical projectile with a diameter and length of 1.25 mm. The fiberglass circuit board material had a density of 1.65 g/cu cm, giving a projectile mass of 2.53 mg. An analysis was performed using the Bumper code to estimate the probability of impact to the shuttle from a 1.25 mm diameter particle. Table 1 shows a 1.6% chance (impact odds = 1 in 62) of a 1.25 mm or larger MMOD impact on the radiators of the vehicle during a typical ISS mission. There is a 0.4% chance (impact odds = 1 in 260) that a 1.25 mm or larger MMOD particle would impact the RCC wing leading edge and nose cap during a typical miion. Figure 11 illustrates the vulnerable areas of the wing leading edge reinforced carbon-carbon (RCC), an area of the vehicle that is very sensitive to impact damage. The highlighted red, orange, yellow, and light green areas would be expected to experience critical damage if impacted by an OD particle such as the one that hit the RH4 radiator panel on STS-115.

Description: A second-order unstructured-grid code, developed and used primarily for steady aerodynamic simulations, is applied to the synthetic jet in a cross flow. The code, FUN3D, is a vertex-centered finite-volume method originally developed by Anderson[1, 2], and is currently supported by members of the Fast Adaptive Aerospace Tools team at NASA Langley. Used primarily for design[3] and analysis[4] of steady aerodynamic configurations, FUN3D incorporates a discrete adjoint capability, and supports parallel computations using MPI. A detailed description of the FUN3D code can be found in the references given above. The code is under continuous development and contains a variety of flux splitting algorithms for the inviscid terms, two methods for computing gradients, several turbulence models, and several solution methodologies; all in varying states of development. Only the most robust and reliable components, based on experiences with steady aerodynamic simulations, were employed in this work. As applied in this work, FUN3D solves the Reynolds averaged Navier-Stokes equations using the one equation turbulence model of Spalart and Allmaras[5]. The spatial discretization is formed on unstructured meshes using a vertex-centered approach. The inviscid terms are evaluated by a flux-difference splitting formulation using least-squares reconstruction and Roe-type approximate Riemann fluxes. Green-Gauss gradient evaluations are used for viscous and turbulence modeling terms. The discrete spatial operator is combined with a backward time operator which is then solved iteratively using point or line Gauss-Seidel and local time stepping in a pseudo time. For steady flows, the physical time step is set to infinity and the pseudo time step is ramped up with the iteration count. A second-order backward in time operator is used for time accurate flows with 20 to 50 steps in the pseudo time applied at each physical time step. For this effort, FUN3D was modified to support spatially varying boundary and initial conditions, and unsteady boundary conditions. Also, a specialized in/out flow boundary condition was implemented to model the action of the diaphragm. This boundary condition is described below in more detail. The grids were generated using the internally developed codes GridEX[6] for meshing the surfaces and inviscid regions of the domain, and for CAD access; and MesherX[7] for meshing the viscous regions. Grid spacing in on the surfaces and in the inviscid regions are indirectly controlled by specifying sources. The viscous layers are generated using an advancing layer technique. MeshersX allows the user to control the spatial variation of the first step off the surface, growth rates, and the termination criterion by providing small problem dependent subroutines.

Description: carbonaceous chondrites are of the major interest since they contain one of the most primitive organic matters. However, aqueous alteration has more or less overprinted their original features in a way that needed to be assessed. That was done in the present study by comparing the mineralogy of the most altered CR1 chondrite, GRO 95577, to a less altered CR2, Renazzo. Their modal analyses were achieved thanks to a new method, based on X-ray elemental maps acquired on electron microprobe, and on IDL image treatment. It allowed the collection of new data on the composition of Renazzo and confirmed the classification of GRO 95577 as a CR1. New alteration products for CRs, vermiculite and clinochlore, were observed. The homogeneity of the Fe-poor clays in the CR1 and the distinctive matrix composition in the two chondrites suggest a wide-range of aqueous alteration on CRs. The preservation of the outlines of the chondrules in GRO 95577 and the elemental transfers of Al, Fe and Ca throughout the chondrule and of Fe and S from the matrix to the chondrule favor the idea of an asteroidal location of the aqueous alteration. From their mineralogical descriptions and modal abundances, the element repartitions in Renazzo and GRO 95577 were computed. It indicates a possible relationship between these two chondrites via an isochemical alteration process. Knowing the chemical reactions that occurred during the alteration, it was thus possible to decipher the mineralogical modal abundances in the unaltered CR body.

Description: Rubidium-strontium and samarium-neodymium isotopes of lunar meteorite LaPaz Icefield (LAP) 02205 are consistent with derivation of the parent magma from a source region similar to that which produced the Apollo 12 low-Ti olivine basalts followed by mixing of the magma with small amounts (1 to 2 wt%) of trace element-enriched material similar to lunar KREEP-rich sample SaU 169. The crystallization age of LAP 02205 is most precisely dated by an internal Rb-Sr isochron of 2991+/-14 Ma, with an initial Sr-87/Sr-88 at the time of crystallization of 0.699836+/-0.000010. Leachable REE-rich phosphate phases of LAP 02205 do not plot on a Sm-Nd mineral isochron, indicating contamination or open system behavior of the phosphates. Excluding anomalous phases from the calculation of a Sm-Nd isochron yields a crystallization age of 2992+/-85 (initial Epsilon Nd-143 = +2.9+/-0.8) that is within error of the Rb-Sr age, and in agreement with other independent age determinations for LAP 02205 from Ar-Ar and U-Pb methods. The calculated Sm-147/Nd-144 source ratios for LAP 02205, various Apollo 12 and 15 basalts, and samples with strong affinities to KREEP (SaU 169, NWA 773, 15386) are uncorrelated with their crystallization ages. This finding does not support the involvement of a common KREEP component as a heat source for lunar melting events that occurred after crystallization of the lunar magma ocean.

Description: The Os-187/Os-188 for twenty-two ureilite whole rock samples, including monomict, augite-bearing, and polymict lithologies, were examined in order to constrain the provenance and subsequent magmatic processing of the ureilite parent body (or bodies). The Re/Os ratios of most ureilites show evidence for a recent disturbance, probably related to Re mobility during weathering, and no meaningful chronological information can be extracted from the present data set. The ureilite Os-187/Os-188 ratios span a range from 0.11739 to 0.13018, with an average of 0.1258+/-0.0023 (1(sigma)), similar to typical carbonaceous chondrites, and distinct from ordinary or enstatite chondrites. The similar mean of Os-187/Os-188 measured for the ureilites and carbonaceous chondrites suggests that the ureilite parent body probably formed within the same region of the solar nebula as carbonaceous chondrites. From the narrow range of the 187Os/188Os distribution in ureilite meteorites it is further concluded that Re was not significantly fractionated from Os during planetary differentiation and was not lost along with the missing ureilitic melt component. The lack of large Re/Os fractionations requires that Re/Os partitioning was controlled by a metal phase, and thus metal had to be stable throughout the interval of magmatic processing on the ureilite parent body.

Description: This viewgraph presentation reviews the Orion Crew Exploration vehicle (CEV) and its usage in the exploration of the moon and subsequent travel to Mars. Schedules for development and testing of the CEV are shown. Also displayed are various high level design views of the CEV, the launch abort system, the Atlas Docking adapter, and the service module.

Description: Observations from the Polar spacecraft, taken during a period of northward interplanetary magnetic field (IMF) show magnetosheath ions within the magnetosphere with velocity distributions resulting from multiple merging sites along the same field line. The observations from the TIDE instrument show two separate ion energy-time dispersions that are attributed to two widely separated (-20Re) merging sites. Estimates of the initial merging times show that they occurred nearly simultaneously (within 5 minutes.) Along with these populations, cold, ionospheric ions were observed counterstreaming along the field lines. The presence of such ions is evidence that these field lines are connected to the ionosphere on both ends. These results are consistent with the hypothesis that double merging can produce closed field lines populated by solar wind plasma. While the merging sites cannot be unambiguously located, the observations and analyses favor one site poleward of the northern cusp and a second site at low latitudes.

Description: Extensive melt on the Greenland Ice Sheet has been documented by a variety of ground and satellite measurements in recent years. If the well-documented warming continues in the Arctic, melting of the Greenland Ice Sheet will likely accelerate, contributing to sea-level rise. Modeling studies indicate that an annual or summer temperature rise of 1 C on the ice sheet will increase melt by 20-50% therefore, surface temperature is one of the most important ice-sheet parameters to study for analysis of changes in the mass balance of the ice-sheet. The Greenland Ice Sheet contains enough water to produce a rise in eustatic sea level of up to 7.0 m if the ice were to melt completely. However, even small changes (centimeters) in sea level would cause important economic and societal consequences in the world's major coastal cities thus it is extremely important to monitor changes in the ice-sheet surface temperature and to ultimately quantify these changes in terms of amount of sea-level rise. We have compiled a high-resolution, daily time series of surface temperature of the Greenland Ice Sheet, using the I-km resolution, clear-sky land-surface temperature (LST) standard product from the Moderate-Resolution Imaging Spectroradiometer (MODIS), from 2000 - 2006. We also use Gravity Recovery and Climate Experiment (GRACE) data, averaged over 10-day periods, to measure change in mass of the ice sheet as it melt and snow accumulates. Surface temperature can be used to determine frequency of surface melt, timing of the start and the end of the melt season, and duration of melt. In conjunction with GRACE data, it can also be used to analyze timing of ice-sheet mass loss and gain.

Description: The Space Technology 5 (ST-5) mission successfully placed three micro-satellites in a 300 x 4500 km dawn-dusk orbit on 22 March 2006. Each spacecraft carried a boom-mounted vector fluxgate magnetometer that returned highly sensitive and accurate measurements of the geomagnetic field. These data allow, for the first time, the separation of temporal and spatial variations in field-aligned current (FAC) perturbations measured in low-Earth orbit on time scales of approximately 10 sec to 10 min. The constellation measurements are used to directly determine field-aligned current sheet motion, thickness and current density. In doing so, we demonstrate two multi-point methods for the inference of FAC current density that have not previously been possible in low-Earth orbit; 1) the "standard method," based upon s/c velocity, but corrected for FAC current sheet motion, and 2) the "gradiometer method" which uses simultaneous magnetic field measurements at two points with known separation. Future studies will apply these methods to the entire ST-5 data set and expand to include geomagnetic field gradient analyses as well as field-aligned and ionospheric currents.

Description: Man-made molecules called chlorofluorcarbons (CFCs) are broken apart in the stratosphere by high energy light, and the reactive chlorine gases that come from them cause the ozone hole. Since the ozone layer stops high energy light from reaching low altitudes, CFCs must be transported to high altitudes to be broken apart. The number of molecules per volume (the density) is much smaller at high altitudes than near the surface, and CFC molecules have a very small chance of reaching that altitude in any particular year. Many tons of CFCs were put into the atmosphere during the end of the last century, and it will take many years for all of them to be destroyed. Each CFC has an atmospheric lifetime that depends on the amount of energy required to break them apart. Two of the gases that were made the most are CFC13 and CF2C12. It takes more energy to break apart CF2C12 than CFC13, and its lifetime is about 100 years, nearly twice as long as the lifetime for CFC13. It is hard to figure out the lifetimes from surface measurements because we don't know exactly how much was released into the air each year. Atmospheric models are used to predict what will happen to ozone and other gases as the CFCs decrease and other gases like C02 continue to increase during the next century. CFC lifetimes are used to predict future concentrations and all assessment models use the predicted future concentrations. The models have different circulations and the amount of CFC lost according to the model may not match the loss that is expected according to the lifetime. In models the amount destroyed per year depends on how fast the model pushes air into the stratosphere and how much goes to high altitudes each year. This paper looks at the way the model circulation changes the lifetimes, and looks at measurements that tell us which model is more realistic. Some models do a good job reproducing the age-of-air, which tells us that these models are circulating the stratospheric air at the right speed. These same models also do a good job reproducing the amount of CFCs in the lower atmosphere where they were measured by instruments on NASA's ER-2, a research plane that flies in the lower stratosphere. The lifetime for CFC13 that is calculated using the models that do the best job matching the data is about 25% longer than most people thought. This paper shows that using these measurements to decide which models are more realistic helps us understand why their predictions are different from each other and also to decide which predictions are more likely.

Description: During summer 2006, the NASA African Monsoon Multidisciplinary Analysis (NAMMA) organized a field campaign in Africa called Special Observation Period (SOP-3), in which scientists in the field were involved in a number of surface network and aircraft measurements. One of the scientific goals of the campaign was to understand the nature and causes for tropical cyclogenesis originating out of African Easterly Waves (AEWs, westward propagating atmospheric disturbances sometimes associated with precursors of hurricanes), and the role that the Saharan Air Layer (SAL, a hot and dry air layer advecting large amounts of dust) can play in the formation or suppression of tropical cyclones. During the NAMMA campaign a high-resolution global model, the NASA GEOS-5, was operationally run by the NASA Global Modeling and Assimilation Office (GMAO) in support to the mission. The daily GEOS-5 forecasts were found to be very useful by decision-making scientists in the field as an aid to discriminate between developing and non-developing AEWs and plan the flight tracks. In the post-event analyses which were performed mostly by the Goddard Laboratory for Atmospheres, two events were highlighted: a non-developing AEW which appeared to have been suppressed by Saharan air, compared to a developing AEW which was the precursor of hurricane Helene. Both events were successfully predicted by the GEOS-5 during the real-time forecasts provided in support to the mission. In this work it is found that very steep moisture gradients and a strong thermal dipole, with relatively warm air in the mid-troposphere and cool air below, are associated with SAL in both the GEOS-5 forecasts and the NCEP analyses, even at -great distance- from the Sahara. The presence of these unusual thermodynamic features over the Atlantic Ocean, at several thousands of kilometers from the African coastline, is suggestive that SAL mixing is very minimal and that the model's capability of retaining the different properties of air masses during transport are important to represent effectively the role of dry air intrusions in the tropical circulation.

Description: The Geoscience Laser Altimeter System (GLAS), launched in January 2003, is a laser altimeter and lidar for the Earth Observing System's (EOS) ICESat mission. GLAS accommodates three, sequentially operated, diode-pumped, solid-state, Nd:YAG laser transmitters. The laser transmitter requirements, design and qualification test results for this space-based remote sensing instrument is summarized and presented

Description: "Qualification" of fiber optic components holds a very different meaning than it did ten years ago. In the past, qualification meant extensive prolonged testing and screening that led to a programmatic method of reliability assurance. For space flight programs today, the combination of using higher performance commercial technology, with shorter development schedules and tighter mission budgets makes long term testing and reliability characterization unfeasible. In many cases space flight missions will be using technology within years of its development and an example of this is fiber laser technology. Although the technology itself is not a new product the components that comprise a fiber laser system change frequently as processes and packaging changes occur. Once a process or the materials for manufacturing a component change, even the data that existed on its predecessor can no longer provide assurance on the newer version. In order to assure reliability during a space flight mission, the component engineer must understand the requirements of the space flight environment as well as the physics of failure of the components themselves. This can be incorporated into an efficient and effective testing plan that "qualifies" a component to specific criteria defined by the program given the mission requirements and the component limitations. This requires interaction at the very initial stages of design between the system design engineer, mechanical engineer, subsystem engineer and the component hardware engineer. Although this is the desired interaction what typically occurs is that the subsystem engineer asks the components or development engineers to meet difficult requirements without knowledge of the current industry situation or the lack of qualification data. This is then passed on to the vendor who can provide little help with such a harsh set of requirements due to high cost of testing for space flight environments. This presentation is designed to guide the engineers of design, development and components, and vendors of commercial components with how to make an efficient and effective qualification test plan with some basic generic information about many space flight requirements. Issues related to the ~ physics of failure, acceptance criteria and lessons learned will also be discussed to assist with understanding how to approach a space flight mission in an ever changing commercial photonics industry.

Description: The ST5 payload, part of NASA s New Millennium Program headquartered at JPL, consisted of three micro satellites (approx. 30 kg each) deployed into orbit from the Pegasus XL launch. ST5 was a technology demonstration mission, intended to test new technologies for potential use for future missions. In order to meet the launch date schedule of ST 5, a different approach was required rather than the standard I&T approach used for single, room-sized satellites. The I&T phase was planned for spacecraft #1 to undergo integration and test first, followed by spacecraft #2 and #3 in tandem. A team of engineers and technicians planned and executed the integration of all three spacecraft emphasizing versatility and commonality. They increased their knowledge and efficiency through spacecraft #1 integration and testing and utilized their experience and knowledge to safely execute I&T for spacecraft #2 and #3. Each integration team member could perform many different roles and functions and thus better support activities on any of the three spacecraft. The I&T campaign was completed with STS s successful launch on March 22,2006

Description: The primary payload onboard the Cloud-Aerosol Lidar Infrared Pathfinder Satellite Observations (CALIPSO) satellite is a dual-wavelength backscatter lidar designed to provide vertical profiling of clouds and aerosols. Launched in April 2006, the first data from this new satellite was obtained in June 2006. As with any new satellite measurement capability, an immediate post-launch requirement is to verify that the data being acquired is correct lest scientific conclusions begin to be drawn based on flawed data. A standard approach to verifying satellite data is to take a similar, or validation, instrument and fly it onboard a research aircraft. Using an aircraft allows the validation instrument to get directly under the satellite so that both the satellite instrument and the aircraft instrument are sensing the same region of the atmosphere. Although there are almost always some differences in the sampling capabilities of the two instruments, it is nevertheless possible to directly compare the measurements. To validate the measurements from the CALIPSO lidar, a similar instrument, the Cloud Physics Lidar, was flown onboard the NASA high-altitude ER-2 aircraft during July- August 2006. This paper presents results to demonstrate that the CALIPSO lidar is properly calibrated and the CALIPSO Level 1 data products are correct. The importance of the results is to demonstrate to the research community that CALIPSO Level 1 data can be confidently used for scientific research.

Description: The prediction of separation bubbles on NACA 65-213 and NACA 0012 using a modified Chen-Thyson transition model is presented. The contents include: 1) Background; 2) Analysis of NACA 65-213 separation bubble using cebeci's viscous-inviscid interaction method; 3) Analysis of NACA 0012 separation bubble using navier-stokes method; and 4) Comparison with experiment.

Description: Experiments on boundary layer transition with flat, concave and convex walls and various levels of free-stream disturbance and with zero and strong streamwise acceleration have been conducted. Measurements of both fluid mechanics and heat transfer processes were taken. Examples are profiles of mean velocity and temperature; Reynolds normal and shear stresses; turbulent streamwise and cross-stream heat fluxed; turbulent Prandtl number; and streamwise variations of wall skin friction and heat transfer coefficient values. Free-stream turbulence levels were varied over the range from about 0.3 percent to about 8 percent. The effects of curvature on the onset of transition under low disturbance conditions are clear; concave curvature leads to an earlier and more rapid transition and the opposite is true for convex curvature This was previously known but little documentation of the transport processes in the flow was available

Description: Measurements on transition under different levels of adverse pressure gradient and free-stream turbulence level are described. This extensive series of investigations, which was predicated on intermittency measurement techniques, has resulted in correlations for transition length and turbulent spot formation rate. These correlations rae intended to be used in conjunction with boundary layer prediction methods and examples are given of such predictions. More effective predictions of the transition region, especially under conditions of variable pressure gradient, are dependent on a more comprehensive understanding of transition and spot behavior. It is expected that this will result in improved transition modeling.

Description: Experimental work with leading edge separation bubbles is presented to clarify the issues regarding transition in separated regions. Hot-wire measurements, in the form of oscilloscope traces, turbulence intermittency and conditionally sampled velocity distributions are given. The resulting points of transition onset and spot production rates are compared to existing correlations.

Description: A new concept and technique has been developed to directly control boundary-layer transition and turbulence. Near-wall vertical motions are directly suppressed through the application of Lorentz force. Current (j) and magnetic (b) fields are applied parallel to the boundary and normal to each other to produce a Lorentz force (j x B) normal to the boundary. This approach is called magnetic turbulence control (MTC). Experiments have been performed on flat-plate transitional and turbulent boundary layers in water seeded with a weak electrolyte.

Description: An experimental investigation of boundary layer transition in a multi-stage turbine has been completed using surface-mounted hot-film sensors. Tests were carried out using the two-stage Low Speed Research Turbine of the Aerodynamics Research Laboratory of GE Aircraft Engines. Blading in this facility models current, state-of-the-art low pressure turbine configurations. The instrumentation technique involved arrays of densely-packed hot-film sensors on the surfaces of second stage rotor and nozzle blades. The arrays were located at mid-span on both the suction and pressure surfaces. Boundary layer measurements were acquired over a complete range of relevant Reynolds numbers. Data acquisition capabilities provided means for detailed data interrogation in both time and frequency domains. Data indicate that significant regions of laminar and transitional boundary layer flow exist on the rotor and nozzle suction surfaces. Evidence of relaminarization both near the leading edge of the suction surface and along much of the pressure surface was observed. Measurements also reveal the nature of the turbulent bursts occuring within and between the wake segments convecting through the blade row. The complex character of boundary layer transition resulting from flow unsteadiness due to nozzle/nozzle, rotor/nozzle, and nozzle/rotor wake interactions are elucidated using these data. These measurements underscore the need to provide turbomachinery designers with models of boundary layer transition to facilitate accurate prediction of aerodynamic loss and heat transfer.

Description: The end-stage phase of boundary layer transition is characterized by the development of hairpin-like vortices which evolve rapidly into patches of turbulent behavior. In general, the characteristics of the evolution form this hairpin stage to the turbulent stage is poorly understood, which has prompted the present experimental examination of hairpin vortex development and growth processes. Two topics of particular relevance to the workshop focus will be covered: 1) the growth of turbulent spots through the generatio and amalgamation of hairpin-like vortices, and 2) the development of hairpin vortices during transition in an end-wall junction flow. Brief summaries of these studies are described below. Using controlled generation of hairpin vortices by surface injection in a critical laminar boundary layer, detailed flow visualization studies have been done of the phases of growth of single hairpin vortices, from the initial hairgin generation, through the systematic generation of secondary hairpin-like flow structures, culminating in the evolution to a turbulent spot. The key to the growth process is strong vortex-surface interactions, which give rise to strong eruptive events adjacent to the surface, which results in the generation of subsequent hairpin vortex structures due to inviscid-viscuous interactions between the eruptive events and the free steam fluid. The general process of vortex-surface fluid interaction, coupled with subsequent interactions and amalgamation of the generated multiple hairpin-type vortices, is demonstrated as a physical mechanism for the growth and development of turbulent spots. When a boundary layer flow along a surface encounters a bluff body obstruction extending from the surface (such as cylinder or wing), the strong adverse pressure gradients generated by these types of flows result in the concentration of the impinging vorticity into a system of discrete vortices near the end-wall juncture of the obstruction, with the extensions of the vortices engirdling the obstruction to form "necklace" or "horseshoe" vortices. Recent hydrogen bubble and particle image visualization have shown that as Reynolds number is increased for a laminar approach flow, the flow will become critical, and a destabilization of the necklace vortices results in the development of an azimuthal waviness, or "kinks", in the vortices. These vortex kinks are accentuated by Biot-Savart effects, causing portions of a distorted necklace vortex to make a rapid approach to the surface, precipitating processes of localized, three-dimensional surface interactions. These interactions result in the rapid generation, focussing, and ejection of thin tongues of surface fluid, which rapidly roll-over and appear as hairpin vortices in the junction region. Subsequent amalgamation of these hairpin vortices with the necklace vortices produces a complex transitional-type flow. A presentation of key results from both these studies will be done, emphasizing both the ubiquity of such hairpin-type flow structures in manifold transitional-type flows, and the importance of vortex-surface interactions n the development of hairpin vortices.

Description: Our research involves study of the behavior of k-epsilon turbulence models for simulation of bypass-level transition over flat surfaces and turbine blades. One facet of the research has been to assess the performance of a multitude of k-epsilon models in what we call "natural transition", i.e. no modifications to the k-e models. The study has been to ascertain what features in the dynamics of the model affect the start and end of the transition. Some of the findings are in keeping with those reported by others (e.g. ERCOFTAC). A second facet of the research has been to develop and benchmark a new multi-time scale k-epsilon model (MTS) for use in simulating bypass-level transition. This model has certain features of the published MTS models by Hanjalic, Launder, and Schiestel, and by Kim and his coworkers. The major new feature of our MTS model is that it can be used to compute wall shear flows as a low-turbulence Reynolds number type of model, i.e. there is no required partition with patching a one-equation k model in the near-wall region to a two-equation k-epsilon model in the outer part of the flow. Our MTS model has been studied extensively to understand its dynamics in predicting the onset of transition and the end-stage of the transition. Results to date indicate that it far superior to the standard unmodified k-epsilon models. The effects of protracted pressure gradients on the model behavior are currently being investigated.

Description: The transition process which takes place in the attachment-line boundary layer in the presence of gross contamination is an issue of considerable interest to wing designers. It is well known that this flow is very sensitive to the presence of isolated roughness and that transition can be initiated at a very low value of the local medium thickness Reynolds number.Moreover, once the attachment line is turbulent, the flow over the whole wing chords, top and bottom surface, will be turbulent and this has major implications for wind drag.

Description: This viewgraph presentation gives a general overview of the X-43A program. The contents include: 1) X-43A Program Overview; 2) Vehicle Description; 3) Flight 1, MIB & Return to Flight; 4) Flight 2 and Results; and 5) Flight 3 and Results.

Description: The similarity among turbulent spots observed in various transition experiments, and the rate in which they contaminate the surrounding laminar boundary layer is only cursory. The shape of the spot depends on the Reynolds number of the surrounding boundary layer and on the pressure gradient to which it and the surrounding laminar flow are exposed. The propagation speeds of the spot boundaries depend, in addition, on the location from which the spot originated and do not simply scale with the local free stream velocity. The understanding of the manner in which the turbulent manner in which the turbulent spot destabilizes the surrounding, vortical fluid is a key to the understanding of the transition process. We therefore turned to detailed observations near the spot boundaries in general and near the spanwise tip of the spot in particular.

Description: The transition from laminar to turbulent flow in a boundary layer is a complex phenomenon that may take different routes, each involving distinct stages governed by different, often not-yet unraveled dynamical principles. There are, surprisingly, questions concerning virtually every stage in the process, beginning with receptivity to external disturbances, the linear stability of spatially developing flows, different possible nonlinear end games, the formation and propagation of turbulent spots and the emergence of fully developed turbulent flow. There seems no doubt that the flow has to be seen as a forced, nonlinear spatio-temporal system, but the system is so complex that to extract simple insights is still very difficult.

Description: Experiment are being carried out to study the process by which th almost periodic disturbance waves generated naturally by the freestream evolve into turbulence. The boundary layer on a flat plate has been used for this study. The novelty of the approach is in the form of artificial excitation that is used. In this work the flow is excited artificially by deterministic white noise. The weak T-S wave created develops down stream, becomes nonlinear and blows up locally onto a highly distorted flow. These large local distortions of the mean flow allow very high frequency disturbances to grow and form into small turbulent spots. The spots arise from the excitation, and if the same noise sequence is repeated a spot will form at the same position and time instant relative to the excitation.

Description: A program sponsored by the National Aeronautics and Space Administration (NASA) for the investigation of the heat transfer in the transition region of turbine vanes and blades with the object of improving the capability for predicting heat transfer is described,. The accurate prediction of gas-side heat transfer is important to the determination of turbine longevity, engine performance and developmental costs. The need for accurate predictions will become greater as the operating temperatures and stage loading levels of advanced turbine engines increase. The present methods for predicting transition shear stress and heat transfer on turbine blades are based on incomplete knowledge and are largely empirical. To meet the objectives of the NASA program, a team approach consisting of researchers from government, universities, a research institute, and a small business is presented. The research is divided into areas of experimentation, direct numerical simulation (DNS) and turbulence modeling. A summary of the results to date is given for the above research areas in a high-disturbance environment (bypass transition) with a discussion of the model development necessary for use in numerical codes.

Description: In order to understand the end-stages of boundary layer transition in low as well as high disturbance environments it is desirable to establish a unified view of the sequences of physico-mathematical phenomena that lead from laminar flow to self-sustained "bursting" in wall turbulence. The dominant driving disturbances: oncoming free turbulence, unsteady pressure fields, inhomogeneous density fields, inhomogeneities in wall geometry, all force disturbed motions within the boundary layer via multiple competitive receptivity mechanisms. For small disturbances, a sequence of instabilities then leads to sporadic local bursting very near the wall which can sustain turbulence. The local seeds of turbulence then somehow propagate to engulf quite rapidly the surrounding disturbed but still laminar regions.

Description: Quantitative observations of transitional boundary layers in regions of strong flow deceleration on an axial compressor stator blade are reported. Measurements are obtained at a fixed chordwise position, and the blade incidence was varied by changing the compressor throughflow so as to move the transition region relative to the stationary probe. It was thus possible to observe typical boundary layer behavior at various stages of transition in the turbomachine environment. The range of observations covers separating laminar flow at transition onset, and reattachment of intermittently turbulent periodically separated shear layers.

Description: Experimental work at the University of Oxford Osney Lab has demonstrated characteristics of the late-stage transition process by the use of thin-film heat transfer gauges. The development of turbulent spots has been observed in a range of environments, including flat plates, turbine blade cascade tests and wake-passing experiments. These results were taken at Mach/Reynolds numbers and gas-to-wall temperature ratios representative of gas turbines. Analyses of the spot characteristics are consistent with measurements taken in low speed experiments, and support the Schubauer and Klebanoff type of turbulent spots. The addition of simulated wakes from upstream stages has been observed to be primarily superpositional for these tests.

Description: A spatially developing direct numerical simulation has been performed for flow over a flat plate that is subjected to a one-time fluid injection through an elongated slit in the wall. The flow parameters have been chosen to closely approximate the experimental conditions of Haidari, Taylow, and Smith (AIAA-89-0964). A hairpin vortex quickly develops near the upstream end of the slit, and a pair of necklace vortices form around the slow-moving injection fluid. As seen in the experiments and reported in Haidari and Smith (in review, JFM), the hairpin vortex spawns both in-line and sidelobe secondary vortices. However, no subdsidiary vortices (those formed by the inviscid deformation of a vortex-line bundle) are observed. At later times, a set of three different types of vortices are identified: hairpin vortex structures with heads that rise away from the wall horseshoe-shaped vortices that do not rise out of the boundary layer, and quasi-streamwise vortices. These structures interact with each other and with the wall layer to generate new vortices that are similar in structure to those mentioned above, although a particular parent vortex may have an offspring that more nearly resembles another member of the set. Perturbation velocity and vertical vorticity contours reveal an arrowhead shape of the highly disturbed region that is reminiscent of a turbulent spot. Spatially averaged velocity profiles in the highly disturbed area are nonlaminar, but as yet do not show typical low-of-the-wall behavior.

Description: A series of experiments are described which examine the growth of turbulent spots on a flat plate at Reynolds and Mach numbers typical of gas-turbine blading. A short-duration piston tunnel is employed and rapid-response miniature surface-heat-transfer gauges are used to asses the state of the boundary layer. The leading- and trailing-edge velocities of spots are reported for different external pressure gradients and Mach numbers. Also, the lateral spreading angle is determined from the heat-transfer signals which demonstrate dramatically the reduction in spot growth associated with favorable pressure gradients. An associated experiment on the development of turbulent wedges is also reported where liquid-crystal heat-transfer techniques are employed in low-speed wind tunnel to visualize and measure the wedge characteristics. Finally, both liquid crystal techniques and hot-film measurements from flight tests at Mach number of 0.6 are presented.

Description: A transitional laminar boundary layer is developed on a 1m wide km long flat plate in a 0.6m deep water channel with a freestream velocity of 15-50 cm/s. A particulate dispenser under computer control ejects individual particles having diameters of 1/3 delta into the free stream. The particulates are introduced with an initial velocity of U(sub infinity) in the direction of the free stream. They have differing specific gravities of 1.03-2.7 which introduces an additional non-dimensional parameter relating the time taken to traverse the boundary layer to the convective time scale. The particulates produce a wake in the upper region of the boundary layer as they sink towards the wall. Visualization data taken over the range 5 x 10(exp 4) less than Re(sub x) less than 5 x 10(exp 5) indicate that turbulent spots are produced by the disturbances due to the wake rather than by the particulates themselves. This suggests that the spot formation process in this case may be inviscid in nature and may not be strongly influenced by the presence of the wall.

Description: Airfoils at high Reynolds numbers, in general, have small separation bubbles that are usually confined to the leading edge. Since the Reynolds number is large, the turbulence model for the transition region between the laminar and turbulent flow is not important. Furthermore, the onset of transition occurs either at separation or prior to separation and can be predicted satisfactorily by empirical correlations when the incident angle is small and can be assumed to correspond to laminar separation when the correlations do not apply, i.e., at high incidence angles.

Description: This lecture reviews current practice as well as new modeling ideas for the calculation of at least skin friction and heat transfer between the onset and end of transition.

Description: For incompressible benchmark flows, we have demonstrated the capability of the parabolized stability equations (PSE) to simulate the transition process in excellent agreement with microscopic experiments and direct Navier-Stokes simulations at modest computational cost. Encouraged by these results, we have developed the PSE methodology of three-dimensional boundary-layers in general curvilinear coordinates for the range from low to hypersonic speeds, and for both linear and nonlinear problems. For given initial and boundary conditions, the approach permits simulations from receptivity through linear and secondary instabilities into the late stages of transition where significant changes in skin friction and heat transfer coefficients occur. We have performed transition simulations for a variety of two- and three-dimensional similarity solutions and for realistic flows over swept wings at subsonic and supersonic speeds, the pressure ans suction side of turbine blades at low and medium turbulence levels, and over a blunt cone at Mach number Ma = 8. We present selected results for different transition mechanisms with emphasis on the late stage of transition and the evolution of wall-shear stress and heat transfer.

Description: The National Aeronautics and Space Administration is currently designing the Crew Exploration Vehicle (CEV) as a replacement for the Space Shuttle for manned missions to the International Space Station, as a command module for returning astronauts to the moon, and as an earth reentry vehicle for the final leg of manned missions to the moon and Mars. The CEV resembles a scaled-up version of the heritage Apollo vehicle; however, the CEV seal requirements are different than those from Apollo because of its different mission requirements. A review is presented of some of the seals used on the Apollo spacecraft for the gap between the heat shield and backshell and for penetrations through the heat shield, docking hatches, windows, and the capsule pressure hull.

Description: A viewgraph presentation describing the hypersonics program at NASA Dryden Flight Research Center is shown. The topics include: 1) X-43A Program Overview; 2) Vehicle Description; 3) Flight 1, MIB & Return to Flight; 4) Flight 2 and Results; 5) Flight 3 and Results; and 6) Concluding Remarks

Description: An overview of the NASA Glenn Research Center Drive Systems Research will be presented. The primary purpose of this research is to improve performance, reliability, and integrity of aerospace drive systems and space mechanisms. The research is conducted through a combination of in-house, academia, and through contractors. Research is conducted through computer code development and validated through component and system testing. The drive system activity currently has four major thrust areas including: thermal behavior of high speed gearing, health and usage monitoring, advanced components, and space mechanisms.

Description: The Space Technology 7 (ST7) experiment will perform an on-orbit system-level validation of two specific Disturbance Reduction System technologies: colloidal micronewton thrusters and drag-free control. The ST7 Disturbance Reduction System (DRS) is designed to maintain the spacecraft s position with respect to a free-floating test mass while limiting the residual accelerations of that test mass over the frequency range of 1 to 30 mHz. This paper presents the overall design and analysis of the spacecraft drag-free and attitude controllers, with particular attention given to its primary mission mode. These controllers close the loop between the drag-free sensors and the colloidal micronewton thrusters.

Description: Over the past two decades, risk management and risk analysis have emerged throughout the business community in the United States (US) as prominent planning and development strategies used to mitigate risk of failure and ensure a high return on investment (ROI) for business endeavors (financial and otherwise). They are generic tools that can be applied to any business regardless of the sector (i.e., government, university, private) and have been used by the Federal government in the form of institutional practices aimed at maximizing the probability of success in business activities. One US Federal agency that incorporates risk management and analysis techniques into business and/or engineering activities is the National Aeronautics and Space Administration (NASA). The present work is a discussion on mission, spacecraft and instrument design (as well as technology development) and the role of risk management, analysis and mitigation as a fundamental tool in the design process.

Description: This report documents the current status of CompactPCI(Registered TradeMark) connectors in GSFC spaceflight applications. To the extent the information is known, this report summarizes to what component quality level each NASA contractor (referred to as OEM in this report) procured the parts, and what board level and system level testing was performed. The report also provides the current status of the reliability assessment for each GSFC project based on the results of testing and FMEA (Failure Mode Effects Analysis). This report addresses how the CompactPCI(Registered TradeMark) connectors came into existence, and how these became the connector style chosen by many designers of space flight hardware. It identifies the design philosophy and the lack of robustness which has led to several known failure modes. These failure modes include fretting of connector pins during vibration, shock and thermal cycling, exposure of underplating, and increased resistance, including brief excursions to very high resistance. Each of these are signs of aging, which becomes an increasing concern for long duration orbiting space flight applications. This report addresses the mitigation strategy to replace CompactPCI(Registered TradeMark) connectors with space qualified Hypertronics 2mm cPCI connectors. The Hypertronics 2mm cPCI connectors are pin-to-pin compatible with the CompactPCI(Registered TradeMark) connectors and meet all of the same technical requirements, except the ability to hot mate, and to mate directly with a CompactPCI of the opposite gender. A detailed comparison of the CompactPCI(Registered TradeMark) connector and the Hypertronics 2mm cPCI connector is provided to describe the ruggedness of Hypertronics connector for space flight applications. Finally, this report makes recommendations for flight hardware for the future missions where the hardware is yet to be built, as well as for the hardware which has already been built with CompactPCI(Registered TradeMark) connectors.

Description: Sea level rise in the Gulf of Mexico has occurred at a rate of 1.8-2.2 mm/yr during the 20th century, or nearly the same as observed globally due to combined steric and water mass changes. Tide gauges in coastal Louisiana, however, record a substantially larger rate of rise and while a number of causal mechanisms may be responsible, their specific contribution is poorly understood. Using a realistic viscoelastic Earth model, detailed geologic parameters for south Louisiana and new GPS data, we demonstrate that Holocene sedimentary loading in the Gulf and Mississippi River delta is capable of contributing to 1-8 mm/yr of subsidence over areas of 30-0.75 x 10(exp 3) km(exp 2).

Description: A rising volcanic plume from an unknown source was observed on 9-11 August 2006 in the Vava'u Island group in the northernmost islands of Tonga [Matangi News Online, 2006]. On 12 August, the crew on board the yacht Maiken, sailing west from Vava'u to Fiji, encountered 'a vast, many miles wide, belt of densely packed pumice' floating on the water (F. Fransson personal communication, 2006). Later, the crew sailed south and discovered that the source of the pumice was a newly erupting submarine volcano near Home Reef (18.991 deg S, 174.767 deg W).

Description: The flow over the two-dimensional hump model is computed by solving the RANS equations with kappa-omega (SST) model. The governing equations, the flow equations and the turbulent equations, are solved using the 5th order accurate weighted essentially non-oscillatory (WENO) scheme for space discretization and using explicit third order total-variation-diminishing (TVD) Runge-Kutta scheme for time integration. The WENO and the TVD methods and the formulas are explained in [1] and the application of ENO method to N-S equations is given in [2]. The solution method implemented in this computation is described in detail in [3].

Description: Computational analyses have been conducted on the Wall-mounted Glauert-Goldschmied type body ("hump" model) with the Full Unstructured Navier-Stokes 2-D (FUN2D) flow solver developed at NASA LaRC. This investigation uses the time-accurate Reynolds-averaged Navier- Stokes (RANS) approach to predict aerodynamic performance of the active flow control experimental database for the hump model. The workshop is designed to assess the current capabilities of different classes of turbulent flow solution methodologies, such as RANS, to predict flow fields induced by synthetic jets and separation control geometries. The hump model being studied is geometrically similar to that previously tested both experimentally and computationally at NASA LaRC [ref. 1 and 2, respectively].

Description: NASA's MODIS instrument on board the Terra satellite is one of the premier tools to assess aerosol over land and ocean because of its high quality calibration and consistency. We analyze Terra-MODIS's seven year record of aerosol optical depth (AOD) observations to determine whether global aerosol has increased or decreased during this period. This record shows that AOD has decreased over land and increased over ocean. Only the ocean trend is statistically significant and corresponds to an increase in AOD of 0.009, or a 15% increase from background conditions. The strongest increasing trends occur over regions and seasons noted for strong biomass burning. This suggests that biomass burning aerosol dominates the increasing trend over oceans and mitigates the otherwise mostly negative trend over the continents.

Description: The interdependence of rapid ice flow, surface topography and the spatial distribution of subglacial water are examined by linking existing theories. The motivation is to investigate whether the acceleration of an ice-stream tributary contains a positive feedback that encourages the retention of subglacial water that leads to faster flow. Periodically varying surface and bed topographies are related through a linear ice-flow perturbation theory for various values of mean surface slope, perturbation amplitude and basal sliding speeds. The topographic variations lead to a periodic variation in hydraulic potential that is used to infer the tendency for subglacial water to be retained in local hydraulic potential minima. If water retention leads to enhanced basal sliding, a positive feedback loop is closed that could explain the transition from slower tributary flow to faster-streaming flow and the sustained downstream acceleration along the tributary-ice-stream system. A sensitivity study illustrates that the same range of topographic wavelengths most effectively transmitted from the bed to the surface also strongly influences the behavior of subglacial water. A lubrication index is defined to qualitatively measure the heterogeneity of the subglacial hydrologic system. Application of this index to field data shows that the transition from tributary to ice stream closely agrees with the location where subglacial water may be first stored.

Description: This study presents the first direct comparison of variations in seasonal GWS derived from GRACE TWS and simulated SM with GW-level measurements in a semiarid region. Results showed that variations in GWS and SM are the main sources controlling TWS changes over the High Plains, with negligible storage changes from surface water, snow, and biomass. Seasonal variations in GRACE TWS compare favorably with combined GWS from GW-level measurements (total 2,700 wells, average 1,050 GW-level measurements per season) and simulated SM from the Noah land surface model (R = 0.82, RMSD = 33 mm). Estimated uncertainty in seasonal GRACE-derived TWS is 8 mm, and estimated uncertainty in TWS changes is 11 mm. Estimated uncertainty in SM changes is 11 mm and combined uncertainty for TWS-SM changes is 15 mm. Seasonal TWS changes are detectable in 7 out of 9 monitored periods and maximum changes within a year (e.g. between winter and summer) are detectable in all 5 monitored periods. Grace-derived GWS calculated from TWS-SM generally agrees with estimates based on GW-level measurements (R = 0.58, RMSD = 33 mm). Seasonal TWS-SM changes are detectable in 5 out of the 9 monitored periods and maximum changes are detectable in all 5 monitored periods. Good correspondence between GRACE data and GW-level measurements from the intensively monitored High Plains aquifer validates the potential for using GRACE TWS and simulated SM to monitor GWS changes and aquifer depletion in semiarid regions subjected to intensive irrigation pumpage. This method can be used to monitor regions where large-scale aquifer depletion is ongoing, and in situ measurements are limited, such as the North China Plain or western India. This potential should be enhanced by future advances in GRACE processing, which will improve the spatial and temporal resolution of TWS changes, and will further increase applicability of GRACE data for monitoring GWS.

Description: TerraSAR-X is an advanced synthetic aperture radar satellite system for scientific and commercial applications that is realized in a public-private partnership between the German Aerospace Center (DLR) and the Astrium GmbH. TerraSAR-X was launched at June 15, 2007 on top of a Russian DNEPR-1 rocket into a 514 km sun-synchronous dusk-dawn orbit with an 11-day repeat cycle and will be operated for a period of at least 5 years during which it will provide high resolution SAR-data in the X-band. Due to the objectives of the interferometric campaigns the satellite has to comply to tight orbit control requirements, which are formulated in the form of a 250 m toroidal tube around a pre-flight determined reference trajectory (see [1] for details). The acquisition of the reference orbit was one of the main and key activities during the Launch and Early Orbit Phase (LEOP) and had to compensate for both injection errors and spacecraft safe mode attitude control thruster activities. The paper summarizes the activities of GSOC flight dynamics team during both LEOP and early Commissioning Phase, where the main tasks have been 1) the first-acquisition support via angle-tracking and GPS-based orbit determination, 2) maneuver planning for target orbit acquisition and maintenance, and 3) precise orbit and attitude determination for SAR processing support. Furthermore, a presentation on the achieved results and encountered problems will be addressed.

Description: In this paper, we present a general collision risk assessment method, which has been applied through numerical simulations to the Automated Transfer Vehicle (ATV) case. During ATV ascent towards the International Space Station, close approaches between the ATV and objects of the USSTRACOM catalog will be monitored through collision rosk assessment. Usually, collision risk assessment relies on an exclusion volume or a probability threshold method. Probability methods are more effective than exclusion volumes but require accurate covariance data. In this work, we propose to use a criterion defined by an adaptive exclusion area. This criterion does not require any probability calculation but is more effective than exclusion volume methods as demonstrated by our numerical experiments. The results of these studies, when confirmed and finalized, will be used for the ATV operations.

Description: This paper extends a previously developed method for finding spacecraft initial conditions (ICs) that minimize the drift resulting from J2 disturbances while also minimizing the fuel required to attain those ICs. It generalizes the single spacecraft optimization to a formation-wide optimization valid for an arbitrary number of vehicles. Additionally, the desired locations of the spacecraft, separate from the starting location, can be specified, either with respect to a reference orbit, or relative to the other spacecraft in the formation. The three objectives (minimize drift, minimize fuel, and maintain a geometric template) are expressed as competing costs in a linear optimization, and are traded against one another through the use of scalar weights. By carefully selecting these weights and re-initializing the formation at regular intervals, a closed-loop, formation-wide control system is created. This control system can be used to reconfigure the formations on the fly, and creates fuel-efficient plans by placing the spacecraft in semi-invariant orbits. The overall approach is demonstrated through nonlinear simulations for two formations a GEO orbit, and an elliptical orbit.

Description: Space missions to small solar system bodies must deal with multiple perturbations acting on the spacecraft. These include strong perturbations from the gravity field and solar tide, but for small bodies the most important perturbations may arise from solar radiation pressure (SRP) acting on the spacecraft. Previous research has generally investigated the effect of the gravity field, solar tide, and SRP acting on a spacecraft trajectory about an asteroid in isolation and has not considered their joint effect. In this paper a more general theoretical discussion of the joint effects of these forces is given.

Description: The Mars swing-by in the early morning of the 25th of February 2007 was one of the most critical events the Rosetta mission has experienced so far on its way to the comet Churyumov-Gerasimenko. The closest approach took place at a distance of only 250 km from the planet s surface. Missing the optimal target would have translated into considerable fuel cost. In order to achieve confidence in operating through this highly critical mission phase, a navigation analysis exercise was carried out beforehand. This paper describes the purpose and the chosen approach for this preparatory Flight Dynamics activity. It presents and discusses results of the analysis. Emphasis is put on the question of what is needed to simulate a valuable data set representative for operations. The results of the navigation analysis are compared with real data obtained during swing-by operations.

Description: Alternatives to the Tracking and Data Relay Satellite (TDRS) orbit estimation procedure were studied to develop a technique that both produces more reliable results and is more amenable to automation than the prior procedure. The Earth Observing System (EOS) Terra mission has TDRS ephemeris prediction 3(sigma) requirements of 75 meters in position and 5.5 millimeters per second in velocity over a 1.5-day prediction span. Meeting these requirements sometimes required reruns of the prior orbit determination (OD) process, with manual editing of tracking data to get an acceptable solution. After a study of the available alternatives, the Flight Dynamics Facility (FDF) began using the Real-Time Orbit Determination (RTOD(Registered TradeMark)) Kalman filter program for operational support of TDRSs in February 2007. This extended Kalman filter (EKF) is used for daily support, including within hours after most thrusting, to estimate the spacecraft position, velocity, and solar radiation coefficient of reflectivity (C(sub R)). The tracking data used are from the Bilateration Ranging Transponder System (BRTS), selected TDRS System (TDRSS) User satellite tracking data, and Telemetry, Tracking, and Command (TT&C) data. Degraded filter results right after maneuvers and some momentum unloads provided incentive for a hybrid OD technique. The results of combining EKF strengths with the Goddard Trajectory Determination System (GTDS) Differential Correction (DC) program batch-least-squares solutions, as recommended in a 2005 paper on the chain-bias technique, are also presented.

Description: For low energy spacecraft trajectories such as multi-moon orbiters for the Jupiter system, multiple gravity assists by moons could be used in conjunction with ballistic capture to drastically decrease fuel usage. In this paper, we outline a procedure to obtain a family of zero-fuel multi-moon orbiter trajectories, using a family of Keplerian maps derived by the first author previously. The maps capture well the dynamics of the full equations of motion; the phase space contains a connected chaotic zone where intersections between unstable resonant orbit manifolds provide the template for lanes of fast migration between orbits of different semimajor axes. Patched three body approach is used and the four body problem is broken down into two three-body problems, and the search space is considerably reduced by the use of properties of the Keplerian maps. We also introduce the notion of Switching Region where the perturbations due to the two perturbing moons are of comparable strength, and which separates the domains of applicability of the corresponding two Keplerian maps.

Description: The Solar Dynamics Observatory (SDO) mission is the first Space Weather Research Network mission, part of NASA s Living With a Star program.1 This program seeks to understand the changing Sun and its effects on the Solar System, life, and society. To this end, the SDO spacecraft will carry three Sun-observing instruments to geosynchronous orbit: Helioseismic and Magnetic Imager (HMI), led by Stanford University; Atmospheric Imaging Assembly (AIA), led by Lockheed Martin Space and Astrophysics Laboratory; and Extreme Ultraviolet Variability Experiment (EVE), led by the University of Colorado. Links describing the instruments in detail may be found through the SDO web site.2 The basic mission goals are to observe the Sun for a very high percentage of the 5-year mission (10-year goal) with long stretches of uninterrupted observations and with constant, high-data-rate transmission to a dedicated ground station. These goals guided the design of the spacecraft bus that will carry and service the three-instrument payload. At the time of this publication, the SDO spacecraft bus is well into the integration and testing phase at the NASA Goddard Space Flight Center (GSFC). A three-axis stabilized attitude control system (ACS) is needed both to point at the Sun accurately and to keep the roll about the Sun vector correctly positioned. The ACS has four reaction wheel modes and 2 thruster actuated modes. More details about the ACS in general and the control modes in particular can be found in Refs. [3-6]. All four of SDO s wheel-actuated control modes involve Sun-pointing controllers, as might be expected from such a mission. Science mode, during which most science data is collected, uses specialized guide telescopes to point accurately at the Sun. Inertial mode has two sub-modes, one tracks a Sun-referenced target orientation, and another maintains an absolute (star-referenced) target orientation, that both employ a Kalman filter to process data from a digital Sun sensor and two star trackers. However, this paper is concerned only with the other two modes: Safe Hold (SH) and Sun Acquisition (SA).

Description: An analytical approach for spin-stabilized spacecraft attitude prediction is presented for the influence of the residual magnetic torques and the satellite in an elliptical orbit. Assuming a quadripole model for the Earth s magnetic field, an analytical averaging method is applied to obtain the mean residual torque in every orbital period. The orbit mean anomaly is used to compute the average components of residual torque in the spacecraft body frame reference system. The theory is developed for time variations in the orbital elements, giving rise to many curvature integrals. It is observed that the residual magnetic torque does not have component along the spin axis. The inclusion of this torque on the rotational motion differential equations of a spin stabilized spacecraft yields conditions to derive an analytical solution. The solution shows that the residual torque does not affect the spin velocity magnitude, contributing only for the precession and the drift of the spin axis of the spacecraft. The theory developed has been applied to the Brazilian s spin stabilized satellites, which are quite appropriated for verification and comparison of the theory with the data generated and processed by the Satellite Control Center of Brazil National Research Institute. The results show the period that the analytical solution can be used to the attitude propagation, within the dispersion range of the attitude determination system performance of Satellite Control Center of Brazil National Research Institute.

Description: Currently two gravity field satellite missions, CHAMP and GRACE, are equipped with high sensitivity electrostatic accelerometers, measuring the non-conservative forces acting on the spacecraft in three orthogonal directions. During the gravity field recovery these measurements help to separate gravitational and non-gravitational contributions in the observed orbit perturbations. For precise orbit determination purposes all these missions have a dual-frequency GPS receiver on board. The reduced dynamic technique combines the dense and accurate GPS observations with physical models of the forces acting on the spacecraft, complemented by empirical accelerations, which are stochastic parameters adjusted in the orbit determination process. When the spacecraft carries an accelerometer, these measured accelerations can be used to replace the models of the non-conservative forces, such as air drag and solar radiation pressure. This approach is implemented in a batch least-squares estimator of the GPS High Precision Orbit Determination Software Tools (GHOST), developed at DLR/GSOC and DEOS. It is extensively tested with data of the CHAMP and GRACE satellites. As accelerometer observations typically can be affected by an unknown scale factor and bias in each measurement direction, they require calibration during processing. Therefore the estimated state vector is augmented with six parameters: a scale and bias factor for the three axes. In order to converge efficiently to a good solution, reasonable a priori values for the bias factor are necessary. These are calculated by combining the mean value of the accelerometer observations with the mean value of the non-conservative force models and empirical accelerations, estimated when using these models. When replacing the non-conservative force models with accelerometer observations and still estimating empirical accelerations, a good orbit precision is achieved. 100 days of GRACE B data processing results in a mean orbit fit of a few centimeters with respect to high-quality JPL reference orbits. This shows a slightly better consistency compared to the case when using force models. A purely dynamic orbit, without estimating empirical accelerations thus only adjusting six state parameters and the bias and scale factors, gives an orbit fit for the GRACE B test case below the decimeter level. The in orbit calibrated accelerometer observations can be used to validate the modelled accelerations and estimated empirical accelerations computed with the GHOST tools. In along track direction they show the best resemblance, with a mean correlation coefficient of 93% for the same period. In radial and normal direction the correlation is smaller. During days of high solar activity the benefit of using accelerometer observations is clearly visible. The observations during these days show fluctuations which the modelled and empirical accelerations can not follow.

Description: The Lunar Reconnaissance Orbiter (LRO) mission is the first of a series of lunar robotic spacecraft scheduled for launch in Fall 2008. LRO will spend at least one year in a low altitude polar orbit around the Moon, collecting lunar environment science and mapping data to enable future human exploration. The LRO employs a 3-axis stabilized attitude control system (ACS) whose primary control mode, the "Observing mode", provides Lunar Nadir, off-Nadir, and Inertial fine pointing for the science data collection and instrument calibration. The controller combines the capability of fine pointing with that of on-demand large angle full-sky attitude reorientation into a single ACS mode, providing simplicity of spacecraft operation as well as maximum flexibility for science data collection. A conventional suite of ACS components is employed in this mode to meet the pointing and control objectives. This paper describes the design and analysis of the primary LRO fine pointing and attitude re-orientation controller function, known as the "Observing mode" of the ACS subsystem. The control design utilizes quaternion feedback, augmented with a unique algorithm that ensures accurate Nadir tracking during large angle yaw maneuvers in the presence of high system momentum and/or maneuver rates. Results of system stability analysis and Monte Carlo simulations demonstrate that the observing mode controller can meet fine pointing and maneuver performance requirements.

Description: Europe's first polar-orbiting weather satellite, METOPA, was launched by a Soyuz launcher from Baikonur Cosmodrome on the 19th of October of 2006. The routine operations of METOP-A are conducted by EUMETSAT (European Organization for Exploitation of Meteorological Satellites) in the frame of the European Polar System mission (EPS). The METOP-A Launch and Early Orbit Phase (LEOP) operations have been performed by ESA/ESOC. The Flight Dynamics Orbit Determination and Control team (OD&C) at ESOC was in charge of correcting the S/C orbit as delivered by the launcher in such a way that EUMETSAT would be able to acquire the reference orbit with a drift-stop manoeuvre approximately two weeks after a LEOP of 3 days and Hand-Over to the EUMETSAT Control Centre (EUMETSAT-CC) in Darmstadt, Germany. The various strict constraints and the short amount of time available for ESOC operations made this task challenging. Several strategies were prepared before launch and analysed during LEOP based on the achieved injection orbit. This paper presents the different manoeuvre strategies investigated and finally applied to acquire the operational orbit, reporting as well the details of its execution and final achieved state.

Description: This paper details a High Gain Antenna (HGA) pointing algorithm which mitigates jitter during the motion of the antennas on the Solar Dynamics Observatory (SDO) spacecraft. SDO has two HGAs which point towards the Earth and send data to a ground station at a high rate. These antennas are required to track the ground station during the spacecraft Inertial and Science modes, which include periods of inertial Sunpointing as well as calibration slews. The HGAs also experience handoff seasons, where the antennas trade off between pointing at the ground station and pointing away from the Earth. The science instruments on SDO require fine Sun pointing and have a very low jitter tolerance. Analysis showed that the nominal tracking and slewing motions of the antennas cause enough jitter to exceed the HGA portion of the jitter budget. The HGA pointing control algorithm was expanded from its original form as a means to mitigate the jitter.

Description: TerraSAR-X (TSX) and TanDEM-X (TDX) are two advanced synthetic aperture radar (SAR) satellites flying in formation. SAR interferometry allows a high resolution imaging of the Earth by processing SAR images obtained from two slightly different orbits. TSX operates as a repeat-pass interferometer in the first phase of its lifetime and will be supplemented after two years by TDX in order to produce digital elevation models (DEM) with unprecedented accuracy. Such a flying formation makes indeed possible a simultaneous interferometric data acquisition characterized by highly flexible baselines with range of variations between a few hundreds meters and several kilometers [1]. TSX has been successfully launched on the 15th of June, 2007. TDX is expected to be launched on the 31st of May, 2009. A safe and robust maintenance of the formation is based on the concept of relative eccentricity/inclination (e/i) vector separation whose efficiency has already been demonstrated during the Gravity Recovery and Climate Experiment (GRACE) [2]. Here, the satellite relative motion is parameterized by mean of relative orbit elements and the key idea is to align the relative eccentricity and inclination vectors to minimize the hazard of a collision. Previous studies have already shown the pertinence of this concept and have described the way of controlling the formation using an impulsive deterministic control law [3]. Despite the completely different relative orbit control requirements, the same approach can be applied to the TSX/TDX formation. The task of TDX is to maintain the close formation configuration by actively controlling its relative motion with respect to TSX, the leader of the formation. TDX must replicate the absolute orbit keeping maneuvers executed by TSX and also compensate the natural deviation of the relative e/i vectors. In fact the relative orbital elements of the formation tend to drift because of the secular non-keplerian perturbations acting on both satellites. The goal of the ground segment is thus to regularly correct this configuration by performing small orbit correction maneuvers on TDX. The ground station contacts are limited due to the geographic position of the station and the costs for contact time. Only with a polar ground station a contact visibility is possible every orbit for LEO satellites. TSX and TDX use only the Weilheim ground station (in the southern part of Germany) during routine operations. This station allows two scheduled contact per day for the nominal orbit configuration, meaning that the satellite conditions can be checked with an interval of 12 hours. While this limitation is usually not critical for single satellite operations, the visibility constraints drive the achievable orbit control accuracy for a LEO formation if a ground based approach is chosen. Along-track position uncertainties and maneuver execution errors affect the relative motion and can be compensated only after a ground station contact.

Description: In March 2006, the Tracking and Data Relay Satellite (TDRS)-3 experienced an unexpected thrusting event, which caused significant changes to its orbit. Recovery from this anomaly was protracted, raising concerns during the Independent Review Team (IRT) investigation of the anomaly regarding the contingency response readiness. The simulations and readiness exercises discussed in this paper were part of the response to the IRT concerns. This paper explains the various levels of simulation needed to enhance the proficiency of the Flight Dynamics Facility (FDF) and supporting elements in recovery from a TDRS contingency situation. The main emergency to address is when a TDRS has experienced uncommanded, unreported, or misreported thrusting, causing a ground station to lose the ability to acquire the spacecraft, as happened in 2006. The following levels of simulation are proposed: 1) Tests that would be performed by the individual support sites to verify that internal procedures and tools are in place and up to date; 2) Tabletop simulations that would involve all of the key support sites talking through their respective operating procedures to ensure that proper notifications are made and communications links are established; and 3) Comprehensive simulations that would be infrequent, but realistic, involving data exchanges between ground sites and voice and electronic communications among the supporting elements.

Description: This paper presents the state of the art and future prospects for autonomous real-time on-orbit calibration of gyros and attitude sensors. The current practice in ground-based calibration is presented briefly to contrast it with on-orbit calibration. The technical and economic benefits of on-orbit calibration are discussed. Various algorithms for on-orbit calibration are evaluated, including some that are already operating on board spacecraft. Because Redundant Inertial Measurement Units (RIMUs, which are IMUs that have more than three sense axes) are almost ubiquitous on spacecraft, special attention will be given to calibration of RIMUs. In addition, we discuss autonomous on board calibration and how it may be implemented.

Description: Recent literature in applied estimation theory reflects growing interest in the sigma-point (also called unscented ) formulation for optimal sequential state estimation, often describing performance comparisons with extended Kalman filters as applied to specific dynamical problems [c.f. 1, 2, 3]. Favorable attributes of sigma-point filters are described as including a lower expected error for nonlinear even non-differentiable dynamical systems, and a straightforward formulation not requiring derivation or implementation of any partial derivative Jacobian matrices. These attributes are particularly attractive, e.g. in terms of enabling simplified code architecture and streamlined testing, in the formulation of estimators for nonlinear spaceflight mechanics systems, such as filter software onboard deep-space robotic spacecraft. As presented in [4], the Sigma-Point Consider Filter (SPCF) algorithm extends the sigma-point filter algorithm to the problem of consider covariance analysis. Considering parameters in a dynamical system, while estimating its state, provides an upper bound on the estimated state covariance, which is viewed as a conservative approach to designing estimators for problems of general guidance, navigation and control. This is because, whether a parameter in the system model is observable or not, error in the knowledge of the value of a non-estimated parameter will increase the actual uncertainty of the estimated state of the system beyond the level formally indicated by the covariance of an estimator that neglects errors or uncertainty in that parameter. The equations for SPCF covariance evolution are obtained in a fashion similar to the derivation approach taken with standard (i.e. linearized or extended) consider parameterized Kalman filters (c.f. [5]). While in [4] the SPCF and linear-theory consider filter (LTCF) were applied to an illustrative linear dynamics/linear measurement problem, in the present work examines the SPCF as applied to nonlinear sequential consider covariance analysis, i.e. in the presence of nonlinear dynamics and nonlinear measurements. A simple SPCF for orbit determination, exemplifying an algorithm hosted in the guidance, navigation and control (GN&C) computer processor of a hypothetical robotic spacecraft, was implemented, and compared with an identically-parameterized (standard) extended, consider-parameterized Kalman filter. The onboard filtering scenario examined is a hypothetical spacecraft orbit about a small natural body with imperfectly-known mass. The formulations, relative complexities, and performances of the filters are compared and discussed.

Description: This paper presents the work of the Navigation Working Group of the Consultative Committee for Space Data Systems (CCSDS) on development of standards addressing the transfer of orbit, attitude and tracking data for space objects. Much progress has been made since the initial presentation of the standards in 2004, including the progression of the orbit data standard to an accepted standard, and the near completion of the attitude and tracking data standards. The orbit, attitude and tracking standards attempt to address predominant parameterizations for their respective data, and create a message format that enables communication of the data across space agencies and other entities. The messages detailed in each standard are built upon a keyword = value paradigm, where a fixed list of keywords is provided in the standard where users specify information about their data, and also use keywords to encapsulate their data. The paper presents a primer on the CCSDS standardization process to put in context the state of the message standards, and the parameterizations supported in each standard, then shows examples of these standards for orbit, attitude and tracking data. Finalization of the standards is expected by the end of calendar year 2007.

Description: This paper studies the interaction between two satellites after docking. In order to maintain the docked state with uncertainty in the motion of the target vehicle, a game theoretic controller with Stackelberg strategy to minimize the interaction between the satellites is considered. The small perturbation approximation leads to LQ differential game scheme, which is validated to address the docking interactions between a service vehicle and a target vehicle. The open-loop solution are compared with Nash strategy, and it is shown that less control efforts are obtained with Stackelberg strategy.

Description: During design of the SDO Science and Inertial mode PID controllers, the decision was made to disable the integral torque whenever system stability was in question. Three different schemes were developed to determine when to disable or enable the integral torque, and a trade study was performed to determine which scheme to implement. The trade study compared complexity of the control logic, risk of not reenabling the integral gain in time to reject steady-state error, and the amount of integral torque space used. The first scheme calculated a simplified Routh criterion to determine when to disable the integral torque. The second scheme calculates the PD part of the torque and looked to see if that torque would cause actuator saturation. If so, only the PD torque is used. If not, the integral torque is added. Finally, the third scheme compares the attitude and rate errors to limits and disables the integral torque if either of the errors is greater than the limit. Based on the trade study results, the third scheme was selected. Once it was decided when to disable the integral torque, analysis was performed to determine how to disable the integral torque and whether or not to reset the integrator once the integral torque was reenabled. Three ways to disable the integral torque were investigated: zero the input into the integrator, which causes the integral part of the PID control torque to be held constant; zero the integral torque directly but allow the integrator to continue integrating; or zero the integral torque directly and reset the integrator on integral torque reactivation. The analysis looked at complexity of the control logic, slew time plus settling time between each calibration maneuver step, and ability to reject steady-state error. Based on the results of the analysis, the decision was made to zero the input into the integrator without resetting it. Throughout the analysis, a high fidelity simulation was used to test the various implementation methods.

Description: The SC "Phobos-Grunt" flight is planned to 2009 in Russia with the purpose to deliver to the Earth the soil samples of the Mars satellite Phobos. The mission will pass under the following scheme [1-4]: the SC flight from the Earth to the Mars, the SC transit on the Mars satellite orbit, the motion round the Mars on the observation orbit and on the quasi-synchronous one [5], landing on Phobos, taking of a ground and start in the direction to the Earth. The implementation of complicated dynamical operations in the Phobos vicinity is foreseen by the project. The SC will be in a disturbance sphere of gravitational fields from the Sun, the Mars and the Phobos. The SC orbit determination is carried out on a totality of trajectory measurements executed from ground tracking stations and measurements of autonomous systems onboard space vehicle relatively the Phobos. As ground measurements the radio engineering measurements of range and range rate are used. There are possible as onboard optical observations of the Phobos by a television system and ranges from the SC up to the Phobos surface by laser locator. As soon as the Phobos orbit accuracy is insufficient for a solution of a problem of landing its orbit determination will be carried out together with determination of the SC orbit. Therefore the algorithms for joint improving of initial conditions of the SC and the Phobos are necessary to determine parameters of the SC relative the Phobos motion within a single dynamical motion model. After putting on the martial satellite orbit, on the Phobos observation orbit, on the quasi-synchronous orbit in the Phobos vicinity the equipment guidance and the following process of the SC orbit determination relatively Phobos requires a priori knowledge of the Phobos orbit parameters with sufficiently high precision. These parameters should be obtained beforehand using both all modern observations and historical ones.

Description: Carrying six science instruments and three engineering payloads, the Mars Reconnaissance Orbiter (MRO) is the first mission in a low Mars orbit to characterize the surface, subsurface, and atmospheric properties with unprecedented detail. After a seven-month interplanetary cruise, MRO arrived at Mars executing a 1.0 km/s Mars Orbit Insertion (MOI) maneuver. MRO achieved a 430 km periapsis altitude with the final orbit solution indicating that only 10 km was attributable to navigation prediction error. With the last interplanetary maneuver performed four months before MOI, this was a significant accomplishment. This paper describes the navigation analyses and results during the 210-day interplanetary cruise. As of August 2007 MRO has returned more than 18 Terabits of scientific data in support of the objectives set by the Mars Exploration Program (MEP). The robust and exceptional interplanetary navigation performance paved the way for a successful MRO mission.

Description: The TanDEM-X formation employs two separate spacecraft to collect interferometric Synthetic Aperture Radar (SAR) measurements over baselines of about 1 km. These will allow the generation ofa global Digital Elevation Model (DEM) with an relative vertical accuracy of 2-4 m and a 10 m ground resolution. As part of the ground processing, the separation of the SAR antennas at the time of each data take must be reconstructed with a 1 mm accuracy using measurements from two geodetic grade GPS receivers. The paper discusses the TanDEM-X mission as well as the methods employed for determining the interferometric baseline with utmost precision. Measurements collected during the close fly-by of the two GRACE satellites serve as a reference case to illustrate the processing concept, expected accuracy and quality control strategies.

Description: The cloud of cataloged debris produced in low earth orbit by the fragmentation of the Fengyun-1C spacecraft was propagated for 15 years, taking into account all relevant perturbations. Unfortunately, the cloud resulted to be very stable, not suffering substantial debris decay during the time span considered. The only significant short term evolution was the differential spreading of the orbital planes of the fragments, leading to the formation of a debris shell around the earth approximately 7-8 months after the breakup, and the perigee precession of the elliptical orbits. Both effects will render the shell more "isotropic" in the coming years. The immediate consequence of the Chinese anti-satellite test, carried out in an orbital regime populated by many important operational satellites, was to increase significantly the probability of collision with man-made debris. For the two Italian spacecraft launched in the first half of 2007, the collision probability with cataloged objects increased by 12% for AGILE, in equatorial orbit, and by 38% for COSMO-SkyMed 1, in sun-synchronous orbit.

Description: During 2006, three ESA interplanetary spacecraft, Rosetta, Mars Express (MEX) and Venus Express (VEX), passed through superior solar conjunction. For all three spacecraft, the noise in the post-fit range-rate residuals from the orbit determination was analysed. At small Sun-Earth-Probe (SEP) angles the level was almost two orders of magnitude higher than normal. The main objective was to characterize the Doppler (rangerate) noise as a function of SEP angle. At least then the range-rate data can be appropriately weighted within the orbit determination so that the solution uncertainties are realistic. For VEX, some intervals of particularly noisy Doppler data could be correlated with unusual solar activity. For Rosetta, the biases in the range data residuals were analysed with the aim of improving the model used for calibrating the signal delay due to the solar plasma. The model, which originally had fixed coefficients, was adjusted to achieve better fits to the data. Even the relatively small Doppler biases were well represented. Using the improved model, the electron density at 20 solar radii was compared with earlier results obtained by radio science studies using Voyager 2 and Ulysses radiometric data. There is some evidence for a dependency of the density on the phase within the 11 years solar cycle.

Description: This paper describes the operations to control the Moon impact of the 3-axis stabilized spacecraft SMART-1 in September 2006. SMART-1 was launched on 27/09/2003. It was the first ESA mission to use an Electric Propulsion (EP) engine as the main motor to spiral out of the Earth gravity field and reach a scientific moon orbit [1]. During September 2005 the last EP maneuvers were performed using the remaining Xenon, in order to compensate for the 3rd body perturbations of the Sun and Earth. These operations extended the mission for an additional year. Afterwards the EP performance became unpredictable and low, so that no meaningful operation for the moon impact could be done. To move the predicted impact point on the 16/8/2006 into visibility from Earth an alternative Delta-V strategy was designed. Due to their alignment, the attitude thrusters could not be used directly to generate the Delta-V, so this strategy was based on controlled angular momentum biasing. Firing along the velocity vector around apolune, the remaining Hydrazine left from the attitude control budget was used, to shift the impact to the required coordinates.

Description: In this paper, optimum trajectories in Earth Transfer Orbit (ETO) for a lunar transportation system are proposed. This paper aims at improving the payload ratio of the reusable orbital transfer vehicle (OTV), which transports the payload from Low Earth Orbit (LEO) to Lunar Low Orbit (LLO) and returns to LEO. In ETO, we discuss ballistic flight using chemical propulsion, multi-impulse flight using electrical propulsion, and aero-assisted flight using aero-brake. The feasibility of the OTV is considered.

Description: The Lunar CRater Observation and Sensing Satellite (LCROSS) was competitively selected by the National Aeronautical and Space Administration (NASA) Exploration Systems Mission Directorate (ESMD) as a low-cost (〈 $80M) 1000 kg secondary payload to be launched with the Lunar Reconnaissance Orbiter (LRO) in October of 2008. LCROSS is a lunar impactor mission that will investigate the presence or absence of water in a permanently shadowed crater. Following launch, trans-lunar injection (TLI) and separation from LRO, LCROSS will remain attached to the launch vehicle's approximately 2300 kg spent Earth Departure Upper Stage (EDUS) and will guide it toward an impact of a permanently shadowed crater at the lunar South Pole. Hours prior to impact, LCROSS will separate from the EDUS and perform a braking maneuver that will allow the spacecraft to take measurements of the resulting EDUS impact ejecta cloud for several minutes, before impacting the crater as well. As a cost-capped secondary mission that must accommodate specific LRO launch dates, LCROSS faces unique challenges and constraints that must be carefully reconciled in order to satisfy an ambitious set of science observation requirements. This paper examines driving mission requirements and constraints and describes the trajectory design and navigation strategy that shape the LCROSS mission.

Description: The Mars Reconnaissance Orbiter reached Mars on March 10, 2006 and performed a Mars orbit insertion maneuver of 1 km/s to enter into a large elliptical orbit. Three weeks later, aerobraking operations began and lasted about five months. Aerobraking utilized the atmospheric drag to reduce the large elliptical orbit into a smaller, near circular orbit. At the time of MRO aerobraking, there were three other operational spacecraft orbiting Mars and the navigation team had to minimize the possibility of a collision. This paper describes the daily operations of the MRO navigation team during this time as well as the collision avoidance strategy development and implementation.

Description: The biggest and most advanced Earth Observation Satellite in-orbit, developed by the European Space Agency (ESA) and its member states, is Envisat. It was launched on March 1, 2002 by an Ariane V from French Guyana and holds a total of 10 multi-disciplinary Earth observation instruments, among which an Advanced Synthetic Aperture Radar (ASAR). The ASAR user community requested the Flight Dynamics division of the European Space Operations Centre (ESOC) to investigate how the orbit control maintenance strategy for Envisat could be changed to optimize ASAR interferometry opportunities overall and in addition support the International Polar Year 2007/2008 initiative. The Polar Regions play a pivotal role in understanding our planet and our impact on it as they are recognized as sensitive barometers of environmental change. One of the main themes of the International Polar Year 2007/2008 is therefore the study of Earth s changing ice and snow, and its impact on our planet and our lives. Naturally, ESA would like to support this very important initiative. This paper presents the investigations that have been conducted to support these requests in the best possible way. It discusses the orbit maintenance strategy that has been in place since its launch, ensuring the actual orbit to be within 1 km of a so-called reference orbit, and presents the new orbit maintenance strategy that is aimed at improving/increasing the opportunities for Envisat ASAR interferometry, while preserving the fuel on board the spacecraft. The hydrazine on-board Envisat happens to be a precious resource as only approximately 300 kg of it was available at launch, like ERS-2. The difference being however that the mass of Envisat is approximately 3.2 times that of ERS-2.

Description: First this paper describes a method how the trajectory of the launcher can be modelled for the contingency analysis without having much information about the launch vehicle itself. From a dense sequence of state vectors a velocity profile is derived which is sufficiently accurate to enable the Flight Dynamics Team to integrate parts of the launcher trajectory on its own and to simulate contingency cases by modifying the velocity profile. Then the paper focuses on the thorough visibility analysis which has to follow the contingency case or burn performance simulations. In the ideal case it is possible to identify a ground station which is able to acquire the satellite independent from the burn performance. The correlations between the burn performance and the pointing at subsequent ground stations are derived with the aim of establishing simple guidelines which can be applied quickly and which significantly improve the chance of acquisition at subsequent ground stations. In the paper the method is applied to the Soyuz/Fregat launch with the MetOp satellite. Overall the paper shows that the launcher trajectory modelling with the simulation of contingency cases in connection with a ground station visibility analysis leads to a proper selection of ground stations and acquisition methods. In the MetOp case this ensured successful contact of all ground stations during the first hour after separation without having to rely on any early orbit determination result or state vector update.

Description: The ESA-funded "Cross-Scale Technology Reference Study has been carried out with the primary aim to identify and analyse a mission concept for the investigation of fundamental space plasma processes that involve dynamical non-linear coupling across multiple length scales. To fulfill this scientific mission goal, a constellation of spacecraft is required, flying in loose formations around the Earth and sampling three characteristic plasma scale distances simultaneously, with at least two satellites per scale: electron kinetic (~10 km), ion kinetic (~100-2000 km), magnetospheric fluid (~3000-15000 km). The key Cross-Scale mission drivers identified are the number of S/C, the space segment configuration, the reference orbit design, the transfer and deployment strategy, the inter-satellite localization and synchronization process and the mission operations. This paper presents a comprehensive overview of the mission design and analysis for the Cross-Scale concept and outlines a technically feasible mission architecture for a multi-dimensional investigation of space plasma phenomena. The main effort has been devoted to apply a thorough mission-level trade-off approach and to accomplish an exhaustive analysis, so as to allow the characterization of a wide range of mission requirements and design solutions.

Description: A method is presented for generating and maintaining a lunar mapping orbit using continuous low-thrust hardware. Optimal control theory is used to maintain a lunar orbit that is low-altitude, near-polar, and Sun-synchronous; three typical requirements for a successful lunar mapping mission. The analysis of the optimal control problem leads to the commonly seen two-point boundary value problem, which is solved using a simple indirect shooting algorithm. Simulations are presented for a 50-day mapping duration, in which it is shown that a very tight control is achieved with thrust levels below 1 N for a 1000 kg spacecraft. A straightforward approach for using the method presented to compute missions of any duration is also discussed.

Description: The Lunar Reconnaissance Orbiter (LRO) is scheduled to launch in 2008 as the first mission under NASA's Vision for Space Exploration. Following several weeks in a quasi-frozen commissioning orbit, LRO will fly in a 50 km mean altitude lunar polar orbit. During the one year mission duration, the orbital dynamics of a low lunar orbit force LRO to perform periodic sets of stationkeeping maneuvers. This paper explores the characteristics of low lunar orbits and explains how the LRO stationkeeping plan is designed to accommodate the dynamics in such an orbit. The stationkeeping algorithm used for LRO must meet five mission constraints. These five constraints are to maintain ground station contact during maneuvers, to control the altitude variation of the orbit, to distribute periselene equally between northern and southern hemispheres, to match eccentricity at the beginning and the end of the sidereal period, and to minimize stationkeeping deltaV. This paper addresses how the maneuver plan for LRO is designed to meet all of the above constraints.

Description: Since 2000, Boeing 702 satellites have used electric propulsion for transfer to geostationary orbits. The use of the 25cm Xenon Ion Propulsion System (25cm XIPS) results in more than a tenfold increase in specific impulse with the corresponding decrease in propellant mass needed to complete the mission when compared to chemical propulsion[1]. In addition to more favorable mass properties, with the use of XIPS, the 702 has been able to achieve orbit insertions with higher accuracy than it would have been possible with the use of chemical thrusters. This paper describes the experience attained by using the 702 XIPS ascent strategy to transfer satellite to geosynchronous orbits.

Description: The Lunar Reconnaissance Orbiter (LRO) undergoes a series of thruster maneuvers to attain lunar orbit. The first of the series of lunar orbit insertion (LOI) maneuvers is crucial to the success of the mission. Therefore, it is important to characterize the disturbances acting on the spacecraft during this phase of the mission. This paper focuses on the internal disturbance force caused by fuel slosh and its impact on attitude control. During the first LOI maneuver (LOI-1), approximately 50% of the total fuel mass is used or roughly 25% of the spacecraft s wet mass, during the 38-minute burn. The forces imparted on the spacecraft from the fuel are dependent on the fill level of the two fuel tanks. During LOI-1, the fill level in both tanks varies greatly and thus so does the disturbance level caused by the fuel. It is therefore necessary to account for the time-varying mass properties of the spacecraft and the effects of the varying fuel levels during the entire 38-minute maneuver. Two simulations are developed in Mathworks s Simulink to analyze the fuel slosh effect. The first model, a baseline model, is a rigid body dynamics model where the fuel slosh is not modeled. The second is a multibody model, developed using a multibody dynamics toolbox, where each of the two fuel tanks and the remaining spacecraft body are treated as separate rigid bodies. The simulations are executed in a piece-wise fashion to account for the time-varying mass properties, and to accommodate the multibody toolbox. Disturbances caused by fuel slosh during both lunar and mission orbit insertions will be analyzed through simulation of different dynamics models. Results of the analysis will show the effects of the slosh disturbance on the spacecraft s attitude.

Description: Optimization techniques are critical when investigating Earth to Mars trajectories since they have the potential of reducing the total (delta)V of a mission. A deep space maneuver (DSM) executed during the cruise may improve a trajectory by reducing the total mission V. Nonetheless, DSMs not only may improve trajectory performance (from an energetic point of view) but also open up new families of trajectories that would satisfy very specific mission requirements not achievable with ballistic trajectories. In the following pages, various specific examples showing the potential advantages of the usage of broken plane maneuvers will be introduced. These examples correspond to possible scenarios for Earth to Mars trajectories during the next decade (2010-2020).

Description: ExoMars is ESA s next mission to planet Mars. The probe is aimed for launch either in 2013 or in 2016. The project is currently undergoing Phase B1 studies under ESA management and Thales Alenia Space Italia project leadership. In that context, DEIMOS Space is responsible for the Mission Analysis and Design for the interplanetary and the entry, descent and landing (EDL) activities. The present mission baseline is based on an Ariane 5 or Proton M launch in 2013 of a spacecraft Composite bearing a Carrier Module (CM) and a Descent Module (DM). A back-up option is proposed in 2016. This paper presents the current status of the interplanetary mission design from launch up to the start of the EDL phase.

Description: SIMBOL-X is a high energy new generation telescope covering by a single instrument a continuous energy range starting at classical X-rays and extending to hard X-rays, i.e. from 0.5 to 80 keV. It is using in this field a focalizing payload which until now was used for energy below 10 keV only, via the construction of a telescope distributed on two satellites flying in formation. SIMBOL-X permits a gain of two orders of magnitude in sensibility and spatial resolution in comparison to state of the art hard X-rays instruments. The mirror satellite will be in free flight on a high elliptical orbit and will target the object to observe very precisely, thus focusing the hard X-ray emission thanks to this mirror module. At the focal point area which is situated 20 meters behind the mirror satellite, the detector satellite maintains its position on a forced orbit thanks to a radio link with the mirror satellite and a lateral displacement sensor using a beam emitted onboard the mirror satellite. This configuration is said "formation flying". The location of the detector satellite shall be very finely tuned as it carries the focal plane of this distributed telescope. To provide science measurements, the Simbol-X orbit has been chosen High elliptic (HEO), which means elliptical orbit with a high perigee altitude. Preliminary studies where made with an orbit with an altitude of the perigee of 44000km and altitude of the apogee of 253000km. The orbit was seven days ground track repeated in order to maintain a perigee pass over the Malindi ground station to download scientific telemetry. But as studies went on, difficulties in mass budget, link budget, perigee maintenance and formation flying maintenance were raised. This was mainly due to the vicinity of the Moon and its disturbing effect on the satellites orbits. Alternative orbits have been proposed in order to demonstrate the feasibility of the mission. The problematic of bringing the two satellites from their injection orbit to their operational orbit 20 m apart from each other and then maintain this configuration is very challenging. It requires theoretical development of the relative motion between two satellites in high eccentric orbit with large differential disturbance on the two bodies. This paper will present the mission analysis for the Simbol-X satellites with the complex problematic of doing formation flying in high elliptic orbit.

Description: Within the framework of a partnership agreement, EADS ASTRIUM has worked since June 2006 for the CNES formation flying experiment on the PRISMA mission. EADS ASTRIUM is responsible for the anti-collision function. This responsibility covers the design and the development of the function as a Matlab/Simulink library, as well as its functional validation and performance assessment. PRISMA is a technology in-orbit testbed mission from the Swedish National Space Board, mainly devoted to formation flying demonstration. PRISMA is made of two micro-satellites that will be launched in 2009 on a quasi-circular SSO at about 700 km of altitude. The CNES FFIORD experiment embedded on PRISMA aims at flight validating an FFRF sensor designed for formation control, and assessing its performances, in preparation to future formation flying missions such as Simbol X; FFIORD aims as well at validating various typical autonomous rendezvous and formation guidance and control algorithms. This paper presents the principles of the collision avoidance function developed by EADS ASTRIUM for FFIORD; three kinds of maneuvers were implemented and are presented in this paper with their performances.

Description: The problem of determining the orbit of a space object from measurements based on one pass through the field of view of a radar is not a new one. Extensive research in this area has been carried out in the USA and Russia since the late 50s when these countries started the development of ballistic missile defense (BMD) and Early Warning systems. In Russia these investigations got additional stimulation in the early 60s after the decision to create a Space Surveillance System, whose primary task would be the maintenance of the satellite catalog. These problems were the focus of research interest until the middle 70s when the appropriate techniques and software were implemented for all radars. Then for more than 20 years no new research papers appeared on this subject. This produced an impression that all the problems of track determination based on one pass had been solved and there was no need for further research. In the late 90s interest in this problem arose again in relation to the following. It was estimated that there would be greater than 100,000 objects with size greater than 1-2 cm and collision of an operational spacecraft with any of these objects could have catastrophic results. Thus, for prevention of hazardous approaches and collisions with valuable spacecraft the existing satellite catalog should be extended by at least an order of magnitude This is a very difficult scientific and engineering task. One of the issues is the development of data fusion procedures and the software capable of maintaining such a huge catalog in near real time. The number of daily processed measurements (of all types, radar and optical) for such a system may constitute millions, thus increasing the number of measurements by at least an order of magnitude. Since we will have ten times more satellites and measurements the computer effort required for the correlation of measurements will be two orders of magnitude greater. This could create significant problems for processing data close to real time even for modern computers. Preliminary "compression" of data for one pass through the field of view of a sensor can significantly reduce the requirements to computers and data communication. This compression will occur when all the single measurements of the sensor are replaced by the orbit determined on their basis. The single measurement here means the radar parameters (range, azimuth, elevation, and in some cases range rate) measured by a single pulse.

Description: Precision mean element (PME) satellite theories play a key role in orbit dynamics analyses. These theories employ: nonsingular orbital elements comprehensive force models Generalized Method of Averaging Numerical interpolation concepts The Draper Semianalytical Satellite Theory (DSST) (Refs. 1 - 6), whose development was led by the author, and the independently-developed Universal Semianalytical Method (USM) (Ref. 7) are examples of such theories. These theories provide the capability to tailor the force modeling to meet the desired computational speed vs. accuracy trade-off. The flexibility of such theories is demonstrated by their ability to include complicated atmosphere density models and spacecraft models in the perturbation theory context. The value of high speed satellite theories, in this era of computational plenty, is that they allow new ways of looking at astrodynamical problems such as orbit design (Refs. 8, 9) and atmosphere density updating (Refs. 10, 11). In the mid to late-1980 s, the geodynamics community led the development of very precise geopotential models such as GEM T2 and GEM T3 (Ref. 12), and with the subsequent analysis of the TOPEX flight data, JGM-2 and JGM-3 (Ref. 13). These were high degree and order geopotentials, at least 50 x 50. In 1993, the DSST implementation in the GTDS program was extended to include the 50 x 50 geopotential models (Ref. 14). The 50 x 50 geopotential, J2000 integration coordinate system, and solid Earth tide capabilities were integrated in GTDS by Scott Carter (Ref. 15). This capability demonstrated 1 m accuracy versus the TOPEX Precise Orbit Ephemerides. Subsequently the DSST Standalone program was also extended to include high degree and order geopotential models (Ref. 5). More recently GTDS has been hosted in the Linux PC environment. However, all of these efforts have been limited to modeling the motion of an artificial Earth satellite. They did not consider the additional complexities associated with lunar, planetary, or other natural satellite orbiters. Such complexities include: additional coordinate systems (associated with the direction of the north pole of rotation and the prime meridian of the new central bodies) (Ref. 16) normalized gravity model coefficients (desirable for high degree and order fields) (Ref. 17) indirect oblateness

Description: The paper presents the attitude reorientation taking the advantage of solar radiation pressure without use of any fuel aboard. The strategy had been adopted to make Hayabusa spacecraft keep pointed toward the Sun for several months, while spinning. The paper adds the above mentioned results reported in Sedona this February showing another challenge of combining ion engines propulsion tactically balanced with the solar radiation torque with no spin motion. The operation has been performed since this March for a half year successfully. The flight results are presented with the estimated solar array panel diffusion coefficient and the ion engine's swirl torque.

Description: In the present day, orbit determination by Global Positioning System (GPS) is not unusual. Especially for low-cost small satellites, position determination by an on-board GPS receiver provides a cheap, reliable and precise method. However, the original purpose of GPS is for ground users, so the transmissions from all of the GPS satellites are directed toward the Earth s surface. Hence there are some restrictions for users above the GPS constellation to detect those signals. On the other hand, a desire for precise orbit determination for users in orbits higher than GPS constellation exists. For example, the next Japanese Very Long Baseline Interferometry (VLBI) mission "ASTRO-G" is trying to determine its orbit in an accuracy of a few centimeters at apogee. The use of GPS is essential for such ultra accurate orbit determination. This study aims to construct a method for precise orbit determination for such high orbit users, especially in High Elliptical Orbits (HEOs). There are several approaches for this objective. In this study, a hybrid method with GPS and an accelerometer is chosen. Basically, while the position cannot be determined by an on-board GPS receiver or other Range and Range Rate (RARR) method, all we can do to estimate the user satellite s position is to propagate the orbit along with the force model, which is not perfectly correct. However if it has an accelerometer (ACC), the coefficients of the air drag and the solar radiation pressure applied to the user satellite can be updated and then the propagation along with the "updated" force model can improve the fitting accuracy of the user satellite s orbit. In this study, it is assumed to use an accelerometer available in the present market. The effects by a bias error of an accelerometer will also be discussed in this paper.