Gravity Recovery and Interior Laboratory (GRAIL): Mapping the Lunar Interior from Crust to Core

The Gravity Recovery and Interior Laboratory (GRAIL) is a spacecraft-to-spacecraft tracking mission that was developed to map the structure of the lunar interior by producing a detailed map of the gravity field. The resulting model of the interior will be used to address outstanding questions regarding the Moon’s thermal evolution, and will be applicable more generally to the evolution of all terrestrial planets. Each GRAIL orbiter contains a Lunar Gravity Ranging System instrument that conducts dual-one-way ranging measurements to measure precisely the relative motion between them, which in turn are used to develop the lunar gravity field map. Each orbiter also carries an Education/Public Outreach payload, Moon Knowledge Acquired by Middle-School Students (MoonKAM), in which middle school students target images of the Moon for subsequent classroom analysis. Subsequent to a successful launch on September 10, 2011, the twin GRAIL orbiters embarked on independent trajectories on a 3.5-month-long cruise to the Moon via the EL-1 Lagrange point. The spacecraft were inserted into polar orbits on December 31, 2011 and January 1, 2012. After a succession of 19 maneuvers the two orbiters settled into precision formation to begin science operations in March 1, 2012 with an average altitude of 55 km. The Primary Mission, which consisted of three 27.3-day mapping cycles, was successfully completed in June 2012. The extended mission will permit a second three-month mapping phase at an average altitude of 23 km. This paper provides an overview of the mission: science objectives and measurements, spacecraft and instruments, mission development and design, and data flow and data products.     

Atmosphere-ionosphere coupling from convectively generated gravity waves

Ionospheric variability impacts operational performances of a variety of technological systems, such as HF communication, Global Positioning System (GPS) navigation, and radar surveillance. The ionosphere is not only perturbed by geomagnetic inputs but is also influenced by atmospheric tides and other wave disturbances propagating from the troposphere to high altitudes. Atmospheric Gravity Waves (AGWs) excited by meteorological sources are one of the largest sources of mesoscale variability in the ionosphere. In this paper, Total Electron Content (TEC) data from networks of GPS receivers in the United States are analyzed to investigate AGWs in the ionosphere generated by convective thunderstorms. Two case studies of convectively generated gravity waves are presented. On April 4, 2014 two distinct large convective systems in Texas and Arkansas generated two sets of concentric AGWs that were observed in the ionosphere as Traveling Ionospheric Disturbances (TIDs). The period of the observed TIDs was 20.8 min, the horizontal wavelength was 182.4 km, and the horizontal phase speed was 146.4 m/s. The second case study shows TIDs generated from an extended squall line on December 23, 2015 stretching from the Gulf of Mexico to the Great Lakes in North America. Unlike the concentric wave features seen in the first case study, the extended squall line generated TIDs, which exhibited almost plane-parallel phase fronts. The TID period was 20.1 min, its horizontal wavelength was 209.6 km, and the horizontal phase speed was 180.1 m/s. The AGWs generated by both of these meteorological events have large vertical wavelength (>100 km), which are larger than the F2 layer thickness, thus allowing them to be discernible in the TEC dataset.     

A deorbiter CubeSat for active orbital debris removal

This paper introduces a mission concept for active removal of orbital debris based on the utilization of the CubeSat form factor. The CubeSat is deployed from a carrier spacecraft, known as a mothership, and is equipped with orbital and attitude control actuators to attach to the target debris, stabilize its attitude, and subsequently move the debris to a lower orbit where atmospheric drag is high enough for the bodies to burn up. The mass and orbit altitude of debris objects that are within the realms of the CubeSat’s propulsion capabilities are identified. The attitude control schemes for the detumbling and deorbiting phases of the mission are specified. The objective of the deorbiting maneuver is to decrease the semi-major axis of the debris orbit, at the fastest rate, from its initial value to a final value of about 6471?km (i.e., 100?km above Earth considering a circular orbit) via a continuous low-thrust orbital transfer. Two case studies are investigated to verify the performance of the deorbiter CubeSat during the detumbling and deorbiting phases of the mission. The baseline target debris used in the study are the decommissioned KOMPSAT-1 satellite and the Pegasus rocket body. The results show that the deorbiting times for the target debris are reduced significantly, from several decades to one or two years.     

Application of time-variable process noise in terrestrial reference frames determined from VLBI data

In recent years, Kalman filtering has emerged as a suitable technique to determine terrestrial reference frames (TRFs), a prime example being JTRF2014. The time series approach allows variations of station coordinates that are neither reduced by observational corrections nor considered in the functional model to be taken into account. These variations are primarily due to non-tidal geophysical loading effects that are not reduced according to the current IERS Conventions (2010). It is standard practice that the process noise models applied in Kalman filter TRF solutions are derived from time series of loading displacements and account for station dependent differences. So far, it has been assumed that the parameters of these process noise models are constant over time. However, due to the presence of seasonal and irregular variations, this assumption does not truly reflect reality. In this study, we derive a station coordinate process noise model allowing for such temporal variations. This process noise model and one that is a parameterized version of the former are applied in the computation of TRF solutions based on very long baseline interferometry data. In comparison with a solution based on a constant process noise model, we find that the station coordinates are affected at the millimeter level.     

Biochemical constraints in a protobiotic earth devoid of basic amino acids: the "BAA(-) world"

It has been hypothesized in this journal and elsewhere, based on surveys of published data from prebiotic synthesis experiments and carbonaceous meteorite analyses, that basic amino acids such as lysine and arginine were not abundant on prebiotic Earth. If the basic amino acids were incorporated only rarely into the first peptides formed in that environment, it is important to understand what protobiotic chemistry is possible in their absence. As an initial test of the hypothesis that basic amino acid negative [BAA(-)] proteins could have performed at least a subset of protobiotic chemistry, the current work reports on a survey of 13 archaeal and 13 bacterial genomes that has identified 61 modern gene sequences coding for known or putative proteins not containing arginine or lysine. Eleven of the sequences found code for proteins whose functions are well known and important in the biochemistry of modern microbial life: lysine biosynthesis protein LysW, arginine cluster proteins, copper ion binding proteins, bacterial flagellar proteins, and PE or PPE family proteins. These data indicate that the lack of basic amino acids does not prevent peptides or proteins from serving useful structural and biochemical functions. However, as would be predicted from fundamental physicochemical principles, we see no fossil evidence of prebiotic BAA(-) peptide sequences capable of interacting directly with nucleic acids.     

Energy dependence of the rigidity spectrum of Forbush decrease of galactic cosmic ray intensity

We show that rigidity spectrum of Forbush decrease (Fd) of galactic cosmic ray (GCR) intensity in September 9–23, 2005 clearly depends on energy. We calculated rigidity spectrum of the Fd based on the neutron monitors and Nagoya muon telescope channels’ data divided in three groups according to their cut off rigidities. We found that temporal changes of rigidity spectrum exponent γ are approximately similar for all cut off rigidity groups, but γ values are the larger the higher are cut off rigidities. We conclude that rigidity spectrum of Fd is hard for lower energy range and is soft for the higher energy range. We believe that an energy dependence of the power law rigidity spectrum of Fd is observed owing to the preferential convection–diffusion mechanism during Fd in September 9–23, 2005. It is a reflection of an influence of the temporal changes of the structure of the interplanetary magnetic field (IMF) turbulence in different range of frequency f during Fd. Particularly, a decisive role in formation of the character of the rigidity spectrum belongs to the changes of the exponent ν of the power spectral density (PSD) of the IMF turbulence (PSD ∝ f^−ν). The exponent ν is greater for high frequency region of the IMF turbulence (responsible for scattering of low rigidity particles of GCR), than for low frequency region of the IMF turbulence (being responsible for scattering of higher rigidity particles). Also, we challenge to estimate an existence of slab/2D structure of solar wind turbulence during the Fd in September 9–23, 2005 based on the distribution of average turbulence energy among the IMF’s components.     

An application of the weighted horizontal magnetic gradient to solar compact and eruptive events

We propose to apply the weighted horizontal magnetic gradient (${\mathit{WG}}_M$), introduced in Korsós et al., 2015, for analysing the pre-flare and pre-CME behaviour and evolution of Active Regions (ARs) using the SDO/HMI-Debrecen Data catalogue. To demonstrate the power of investigative capabilities of the ${\mathit{WG}}_M$ method, in terms of flare and CME eruptions, we studied two typical ARs, namely, AR 12158 and AR 12192. The choice of ARs represent canonical cases. AR 12158 produced an X1.6 flare with fast “halo” CME ($v_{\mathit{linear}}$ = 1267 $\text{km}{\text{s}}^{-1}$) while in AR 12192 there occurred a range of powerful X-class eruptions, i.e. X1.1, X1.6, X3.1, X1.0, X2.0 and X2.0-class energetic flares, interestingly, none with an accompanying CME. The value itself and temporal variation of ${\mathit{WG}}_M$ is found to possess potentially important diagnostic information about the intensity of the expected flare class. Furthermore, we have also estimated the flare onset time from the relationship of duration of converging and diverging motions of the area-weighted barycenters of two subgroups of opposite magnetic polarities. This test turns out not only to provide information about the intensity of the expected flare-class and the flare onset time but may also indicate whether a flare will occur with/without fast CME. We have also found that, in the case when the negative polarity barycenter has moved around and the positive one “remained” at the same coordinates preceding eruption, the flare occurred with fast “halo” CME. Otherwise, when both the negative and the positive polarity barycenters have moved around, the AR produced flares without CME. If these properties found for the movement of the barycenters are generic pre-cursors of CME eruption (or lack of it), identifying them may serve as an excellent pre-condition for refining the forecast of the lift-off of CMEs.     

Attitude stabilization of a satellite by magnetic coils

Stabilization problem for a satellite is considered. The only measurement is of the geomagnetic field in the satellite coordinates. The control is carried out by a magnetic moment of current coils (magnetorquers) mounted on the satellite body. The stabilizer constructed in this work solves the problems of magnetic and gravitational stabilization. Qualitative analysis and results of numerical simulation are presented. The results of simulation show that the proposed stabilization system is reliable, and has an appropriate accuracy and does not need powerful sources of energy, and therefore can be used for attitude control of small satellites.     

For aerospace programs,the crises continue

The Legislative Update column reports on President Bush's supplemental budget request for the war in Iraq, the formation of a defense caucus in the House of Representatives, stabilized funding for airport security, the Department of Homeland Security budget, NASA's budget request, funding for NASA's aeronautics program, and NASA's new authorization. NASA's FY04 budget request is $15.47 billion, less than 1% more than FY03. The aeronautics program budget request for FY04 is $959 million, a decrease of 8.5% from FY03. Congress is working to introduce a multiyear authorization for NASA during the 108th session.     

A model for the propagation of the westward traveling surge

A unified model is developed for the propagation of the Westward Traveling Surge (WTS) which can explain the diversity in the observed surge characteristics. The direction of the surge motion depends on the presence of polarization charges on the poleward boundary. This is related to the efficiency with which the poleward ionospheric currents are closed off into the magnetosphere by the field-aligned currents. Inclusion of the electron-ion recombination rate modifies the surge propagation velocity and leads to explicit expressions for the conductivity profile. Sufficient precipitation current is required to overcome electron-ion recombination in order for the surge to expand. When the precipitating current is less than this threshold the WTS retreats. Therefore, the model describes the ionospheric response to both the expansion and recovery phases of the magnetic substorm.