Crustal S-Wave Velocity from Apparent Incidence Angles: A Case Study in Preparation for InSight

Retrieval of crustal structure and thickness of Mars is among the main goals of InSight. Here we investigate which constraints on the crust at the landing site can be provided by apparent P-wave incidence angles derived from P-receiver functions. We consider receiver functions for six different Mars models, calculated from synthetic seismograms generated via Instaseis from the Green’s function databases of the Marsquake Service, in detail. To allow for a larger range of crustal thicknesses and structures, we additionally analyze data from five broad-band stations across Central Europe. We find that the likely usable epicentral distance range for P-wave receiver functions on Mars lies between \(35^{\circ}\) and the core shadow, and can be extended to more than \(150^{\circ}\) by also using the PP-phase. Comparison to models for the spatial distribution of Martian seismicity indicates that sufficient seismicity should occur within the P-wave distance range around InSight within the nominal mission duration to allow for the application of our method. Apparent P-wave incidence angles are derived from the amplitudes of vertical and radial receiver functions at the P-wave onset within a range of period bands, up to 120 s. The apparent incidence angles are directly related to apparent S-wave velocities, which are inverted for the subsurface S-wave velocity structure via a grid search. The veracity of the forward calculated receiver functions and apparent S-wave velocities is ensured by benchmarking various algorithms against the Instaseis synthetics. Results indicate that apparent S-wave velocity curves provide valuable constraints on crustal thickness and structure, even without any additional constraints, and considering the location uncertainty and limited data quantity of InSight. S-wave velocities in the upper half of the crust are constrained best, but if reliable measurements at long periods are available, the curves also provide constraints down to the uppermost mantle. Besides, it is demonstrated that the apparent velocity curves can differentiate between crustal velocity models that are indistinguishable by other methods.     

The Solar Origin of Corotating Interaction Regions and Their Formation in the Inner Heliosphere

Corotating Interaction Regions (CIRs) form as a consequence of the compression of the solar wind at the interface between fast speed streams and slow streams. Dynamic interaction of solar wind streams is a general feature of the heliospheric medium; when the sources of the solar wind streams are relatively stable, the interaction regions form a pattern which corotates with the Sun. The regions of origin of the high speed solar wind streams have been clearly identified as the coronal holes with their open magnetic field structures. The origin of the slow speed solar wind is less clear; slow streams may well originate from a range of coronal configurations adjacent to, or above magnetically closed structures. This article addresses the coronal origin of the stable pattern of solar wind streams which leads to the formation of CIRs. In particular, coronal models based on photospheric measurements are reviewed; we also examine the observations of kinematic and compositional solar wind features at 1 AU, their appearance in the stream interfaces (SIs) of CIRs, and their relationship to the structure of the solar surface and the inner corona; finally we summarise the Helios observations in the inner heliosphere of CIRs and their precursors to give a link between the optical observations on their solar origin and the in-situ plasma observations at 1 AU after their formation. The most important question that remains to be answered concerning the solar origin of CIRs is related to the origin and morphology of the slow solar wind. This revised version was published online in August 2006 with corrections to the Cover Date.     

Corotating Interaction Regions at High Latitudes

Ulysses observed a stable strong CIR from early 1992 through 1994 during its first journey into the southern hemisphere. After the rapid latitude scan in early 1995, Ulysses observed a weaker CIR from early 1996 to mid-1997 in the northern hemisphere as it traveled back to the ecliptic at the orbit of Jupiter. These two CIRs are the observational basis of the investigation into the latitudinal structure of CIRs. The first CIR was caused by an extension of the northern coronal hole into the southern hemisphere during declining solar activity, whereas the second CIR near solar minimum activity was caused by small warps in the streamer belt. The latitudinal structure is described through the presentation of three 26-day periods during the southern CIR. The first at ∼24°S shows the full plasma interaction region including fast and slow wind streams, the compressed shocked flows with embedded stream interface and heliospheric current sheet (HCS), and the forward and reverse shocks with associated accelerated ions and electrons. The second at 40°S exhibits only the reverse shock, accelerated particles, and the 26-day modulation of cosmic rays. The third at 60°S shows only the accelerated particles and modulated cosmic rays. The possible mechanisms for the access of the accelerated particles and the CIR-modulated cosmic rays to high latitudes above the plasma interaction region are presented. They include direct magnetic field connection across latitude due to stochastic field line weaving or to systematic weaving caused by solar differential rotation combined with non-radial expansion of the fast wind. Another possible mechanism is particle diffusion across the average magnetic field, which includes stochastic field line weaving. A constraint on connection to a distant portion of the CIR is energy loss in the solar wind, which is substantial for the relatively slow-moving accelerated ions. Finally, the weaker northern CIR is compared with the southern CIR. It is weak because the inclination of the streamer belt and HCS decreased as Ulysses traveled to lower latitudes so that the spacecraft remained at about the maximum latitudinal extent of the HCS. This revised version was published online in August 2006 with corrections to the Cover Date.     

Comparison of various periodograms for sinusoid detection andfrequency estimation

With the advent of the fast Fourier transform (FFT) algorithm, the periodogram and its variants such as the Bartlett's procedure and Welch method, have become very popular for spectral analysis. However, there has not been a thorough comparison of the detection and estimation performances of these methods. Different forms of the periodogram are studied here for single real tone detection and frequency estimation in the presence of white Gaussian noise. The threshold effect in frequency estimation, that is, when the estimation errors become several orders of magnitude greater than the Cramer-Rao lower bound (CRLB), is also investigated. It is shown that the standard periodogram gives the optimum detection performance for a pure tone while the Welch method is the best detector when there is phase instability in the sinusoid. As expected, since the conventional periodogram is a maximum likelihood estimator of frequency, it generally provides the minimum mean square frequency estimation errors     

Dynamically adaptable m-best 2-D assignment algorithm andmultilevel parallelization

In recent years, there has been considerable interest within the tracking community in an approach to data association based on the m-best two-dimensional (2D) assignment algorithm. Much of the interest has been spurred by its ability to provide various efficient data association solutions, including joint probabilistic data association (JPDA) and multiple hypothesis tracking (MHT). The focus of this work is to describe several recent improvements to the m-best 2D assignment algorithm. One improvement is to utilize a nonintrusive 2D assignment algorithm switching mechanism, based on a problem sparsity threshold. Dynamic switching between two different 2D assignment algorithms, highly suited for sparse and dense problems, respectively, enables more efficient solutions to the numerous 2D assignment problems generated in the m-best 2D assignment framework. Another improvement is to utilize a multilevel parallelization enabling many independent and highly parallelizable tasks to be executed concurrently, including 1) solving the multiple 2D assignment problems via a parallelization of the m-best partitioning task, and 2) calculating the numerous gating tests, state estimates, covariance calculations, and likelihood function evaluations (used as cost coefficients in the 2D assignment problem) via a parallelization of the data association interface task. Using both simulated data and an air traffic surveillance (ATS) problem based on data from two Federal Aviation Administration (FAA) air traffic control radars, we demonstrate that efficient solutions to the data association problem are obtainable using our improvements in the m-best 2D assignment algorithm     

Evaluation of gravity-dependent membrane potential shift in Paramecium

It is still debated whether or not gravity can stimulate unicellular organisms. This question may be settled by revealing changes in the membrane potential in a manner depending on the gravitational forces imposed on the cell. We estimated the gravity-dependent membrane potential shift to be about 1 mV G⁻¹ for Paramecium showing gravikinesis at 1–5 G, on the basis of measurements of gravity-induced changes in active propulsion and those of propulsive velocity in solutions, in which the membrane potential has been measured electrophysiologically. The shift in membrane potential to this extent may occur from mechanoreceptive changes in K⁺ or Ca²⁺ conductance by about 1% and might be at the limit of electrophysiological measurement using membrane potential-sensitive dyes. Our measurements of propulsive velocity vs membrane potential also suggested that the reported propulsive force of Paramecium measured in a solution of graded densities with the aid of a video centrifuge microscope at 350 G was 11 times as large as that for −29 mV, i.e., the resting membrane potential at [K⁺]_o = 1 mM and [Ca²⁺]_o = 1 mM, and, by extrapolation, that Paramecium was hyperpolarized to −60 mV by gravity stimulation of 100- G equivalent, the value corrected by considering the reduction of density difference between the interior and exterior of the cell in the graded density solution. The estimated shift of the membrane potential from −29 mV to −60 mV by 100- G equivalent stimulation, i.e., 0.3 mV G⁻¹, could reach the magnitude entirely feasible to be measured more directly.     

Modelling solar irradiance variations: Comparison with observations,including line-ratio variations

Solar irradiance variations show a strong temporal and spectral dependence. The progression of the Sun through its activity cycle as well as solar rotation are mirrored in the irradiance variations. The spectral dependence is such that the variations are several magnitudes larger in the EUV than in the visible or infrared.We present a simple 3-component model that is based on the assumption that changes in the solar flux are exclusively due to changes in spot and facular coverage. We compare our model to observations of the spectral solar irradiance variations.Despite its simplicity, we find that the agreement between our model and the observations is surprisingly good. We also explore the reliability and the limitations of our approach by comparing observations of the solar facular contrast and of the changes in spectral line depths with our calculations.     

Drilling systems for extraterrestrial subsurface exploration

Drilling consists of 2 processes: breaking the formation with a bit and removing the drilled cuttings. In rotary drilling, rotational speed and weight on bit are used to control drilling, and the optimization of these parameters can markedly improve drilling performance. Although fluids are used for cuttings removal in terrestrial drilling, most planetary drilling systems conduct dry drilling with an auger. Chip removal via water-ice sublimation (when excavating water-ice-bound formations at pressure below the triple point of water) and pneumatic systems are also possible. Pneumatic systems use the gas or vaporization products of a high-density liquid brought from Earth, gas provided by an in situ compressor, or combustion products of a monopropellant. Drill bits can be divided into coring bits, which excavate an annular shaped hole, and full-faced bits. While cylindrical cores are generally superior as scientific samples, and coring drills have better performance characteristics, full-faced bits are simpler systems because the handling of a core requires a very complex robotic mechanism. The greatest constraints to extraterrestrial drilling are (1) the extreme environmental conditions, such as temperature, dust, and pressure; (2) the light-time communications delay, which necessitates highly autonomous systems; and (3) the mission and science constraints, such as mass and power budgets and the types of drilled samples needed for scientific analysis. A classification scheme based on drilling depth is proposed. Each of the 4 depth categories (surface drills, 1-meter class drills, 10-meter class drills, and deep drills) has distinct technological profiles and scientific ramifications.     

Experiment LEND of the NASA Lunar Reconnaissance Orbiter for high-resolution mapping of neutron emission of the Moon

The scientific objectives of neutron mapping of the Moon are presented as 3 investigation tasks of NASA's Lunar Reconnaissance Orbiter mission. Two tasks focus on mapping hydrogen content over the entire Moon and on testing the presence of water-ice deposits at the bottom of permanently shadowed craters at the lunar poles. The third task corresponds to the determination of neutron contribution to the total radiation dose at an altitude of 50 km above the Moon. We show that the Lunar Exploration Neutron Detector (LEND) will be capable of carrying out all 3 investigations. The design concept of LEND is presented together with results of numerical simulations of the instrument's sensitivity for hydrogen detection. The sensitivity of LEND is shown to be characterized by a hydrogen detection limit of about 100 ppm for a polar reference area with a radius of 5 km. If the presence of ice deposits in polar "cold traps" is confirmed, a unique record of many millions of years of lunar history would be obtained, by which the history of lunar impacts could be discerned from the layers of water ice and dust. Future applications of a LEND-type instrument for Mars orbital observations are also discussed.     

Spectral signatures of photosynthesis. II. Coevolution with other stars and the atmosphere on extrasolar worlds

As photosynthesis on Earth produces the primary signatures of life that can be detected astronomically at the global scale, a strong focus of the search for extrasolar life will be photosynthesis, particularly photosynthesis that has evolved with a different parent star. We take previously simulated planetary atmospheric compositions for Earth-like planets around observed F2V and K2V, modeled M1V and M5V stars, and around the active M4.5V star AD Leo; our scenarios use Earth's atmospheric composition as well as very low O2 content in case anoxygenic photosynthesis dominates. With a line-by-line radiative transfer model, we calculate the incident spectral photon flux densities at the surface of the planet and under water. We identify bands of available photosynthetically relevant radiation and find that photosynthetic pigments on planets around F2V stars may peak in absorbance in the blue, K2V in the red-orange, and M stars in the near-infrared, in bands at 0.93-1.1 microm, 1.1-1.4 microm, 1.5-1.8 microm, and 1.8-2.5 microm. However, underwater organisms will be restricted to wavelengths shorter than 1.4 microm and more likely below 1.1 microm. M star planets without oxygenic photosynthesis will have photon fluxes above 1.6 microm curtailed by methane. Longer-wavelength, multi-photo-system series would reduce the quantum yield but could allow for oxygenic photosystems at longer wavelengths. A wavelength of 1.1 microm is a possible upper cutoff for electronic transitions versus only vibrational energy; however, this cutoff is not strict, since such energetics depend on molecular configuration. M star planets could be a half to a tenth as productive as Earth in the visible, but exceed Earth if useful photons extend to 1.1 microm for anoxygenic photosynthesis. Under water, organisms would still be able to survive ultraviolet flares from young M stars and acquire adequate light for growth.