Polarimetric NIR reflectance measurements of regolith simulants at zero phase angle

Polarimetric reflectance measurements have been made at a wavelength of for a suite of predominantly lunar regolith simulants in support of development efforts for the BepiColombo Laser Altimeter (BELA). Measurements were made using an instrument at the University of Bern, Switzerland, that has been modified to accommodate photometric characterizations of laboratory samples down to 0^° phase angle (g) with linearly polarized illumination and a linearly polarized receiver. The data reveal opposition surges that are dominated by polarization state-conserving reflectance terms for all samples. Relative strengths of the trans-state reflectance terms are lowest for the darkest samples, suggesting that multiple scattering is responsible for state conversions. Normal albedo measurements of the lunar simulants range from 0.13 to 0.22 at zero illumination angle (i=0^°). The total reflectance of the regolith simulants at g=0^° were found to decrease with increasing i, which is inconsistent with predictions of reflectance models for Mercury. However, the g=0^° reflectance remains higher at i>0^° than would be expected for a gray Lambert surface that is scaled to the g=i=0^° reflectance value. Polarization ratios for reflectance under polarized illumination but unconstrained emission show that the samples scatter light in the in-plane polarization state more efficiently than in the transverse state at g=0^° and i>0^°. The opposite is true for g>0^° polarization ratios, which indicate that transverse polarized illumination scatters more efficiently at high g. The polarization effect appears to correlate with the sample's characteristic grain size, but the statistical basis of this trend is weak. The implications of these measurements upon the performance of the BELA instrument are discussed.     

The dark matter environment of the Abell 901/902 supercluster: a weak lensing analysis of the HST STAGES survey

We present a high-resolution dark matter reconstruction of the   z = 0.165  Abell 901/902 supercluster from a weak lensing analysis of the Hubble Space Telescope STAGES survey. We detect the four main structures of the supercluster at high significance, resolving substructure within and between the clusters. We find that the distribution of dark matter is well traced by the cluster galaxies, with the brightest cluster galaxies marking out the strongest peaks in the dark matter distribution. We also find a significant extension of the dark matter distribution of Abell 901a in the direction of an infalling X-ray group Abell 901α. We present mass, mass-to-light and mass-to-stellar mass ratio measurements of the structures and substructures that we detect. We find no evidence for variation of the mass-to-light and mass-to-stellar mass ratio between the different clusters. We compare our space-based lensing analysis with an earlier ground-based lensing analysis of the supercluster to demonstrate the importance of space-based imaging for future weak lensing dark matter 'observations'.     

Molecular line mapping of the giant molecular cloud associated with RCW 106 – II. Column density and dynamical state of the clumps

We present a fully sampled C¹⁸O (1–0) map towards the southern giant molecular cloud (GMC) associated with the H  ii region RCW 106, and use it in combination with previous ¹³CO (1–0) mapping to estimate the gas column density as a function of position and velocity. We find localized regions of significant ¹³CO optical depth in the northern part of the cloud, with several of the high-opacity clouds in this region likely associated with a limb-brightened shell around the H  ii region G333.6−0.2. Optical depth corrections broaden the distribution of column densities in the cloud, yielding a lognormal distribution as predicted by simulations of turbulence. Decomposing the ¹³CO and C¹⁸O data cubes into clumps, we find relatively weak correlations between size and linewidth, and a more sensitive dependence of luminosity on size than would be predicted by a constant average column density. The clump mass spectrum has a slope near −1.7, consistent with previous studies. The most massive clumps appear to have gravitational binding energies well in excess of virial equilibrium; we discuss possible explanations, which include magnetic support and neglect of time-varying surface terms in the virial theorem. Unlike molecular clouds as a whole, the clumps within the RCW 106 GMC, while elongated, appear to show random orientations with respect to the Galactic plane.     

The first galaxies: assembly, cooling and the onset of turbulence

We investigate the properties of the first galaxies at   z ≳ 10  with highly resolved numerical simulations, starting from cosmological initial conditions and taking into account all relevant primordial chemistry and cooling. A first galaxy is characterized by the onset of atomic hydrogen cooling, once the virial temperature exceeds  ≃10⁴ K  , and its ability to retain photoheated gas. We follow the complex accretion and star formation history of a  ≃5 × 10⁷ M_⊙  system by means of a detailed merger tree and derive an upper limit on the number of Population III (Pop III) stars formed prior to its assembly. We investigate the thermal and chemical evolution of infalling gas and find that partial ionization at temperatures  ≳10⁴ K  catalyses the formation of  H₂  and hydrogen deuteride, allowing the gas to cool to the temperature of the cosmic microwave background. Depending on the strength of radiative and chemical feedback, primordial star formation might be dominated by intermediate-mass Pop III stars formed during the assembly of the first galaxies. Accretion on to the nascent galaxy begins with hot accretion, where gas is accreted directly from the intergalactic medium and shock heated to the virial temperature, but is quickly accompanied by a phase of cold accretion, where the gas cools in filaments before flowing into the parent halo with high velocities. The latter drives supersonic turbulence at the centre of the galaxy and could lead to very efficient chemical mixing. The onset of turbulence in the first galaxies thus likely marks the transition to Pop II star formation.     

Radio observations of NGC 2808 and other globular clusters: constraints on intermediate-mass black holes

We present the results of a deep radio observation of the globular cluster NGC 2808. We show that there are no sources detected within the core of the cluster, placing constraints on both the pulsar population of the cluster and the mass of a possible intermediate-mass black hole in NGC 2808. We compare the results for this cluster with other constraints on intermediate-mass black holes derived from accretion measures. With the exception of G1 in M 31 which has previously shown radio emission, even with considerably more conservative assumptions, only the clusters with the poorest of observational constraints are consistent with falling on the   M _BH–σ  relation. This result is interpreted in terms of the fundamental differences between galaxies and globular clusters.     

Nonlinear Force-Free Modeling of Coronal Magnetic Fields. II. Modeling a Filament Arcade and Simulated Chromospheric and Photospheric Vector Fields

We compare a variety of nonlinear force-free field (NLFFF) extrapolation algorithms, including optimization, magneto-frictional, and Grad – Rubin-like codes, applied to a solar-like reference model. The model used to test the algorithms includes realistic photospheric Lorentz forces and a complex field including a weakly twisted, right helical flux bundle. The codes were applied to both forced “photospheric” and more force-free “chromospheric” vector magnetic field boundary data derived from the model. When applied to the chromospheric boundary data, the codes are able to recover the presence of the flux bundle and the field’s free energy, though some details of the field connectivity are lost. When the codes are applied to the forced photospheric boundary data, the reference model field is not well recovered, indicating that the combination of Lorentz forces and small spatial scale structure at the photosphere severely impact the extrapolation of the field. Preprocessing of the forced photospheric boundary does improve the extrapolations considerably for the layers above the chromosphere, but the extrapolations are sensitive to the details of the numerical codes and neither the field connectivity nor the free magnetic energy in the full volume are well recovered. The magnetic virial theorem gives a rapid measure of the total magnetic energy without extrapolation though, like the NLFFF codes, it is sensitive to the Lorentz forces in the coronal volume. Both the magnetic virial theorem and the Wiegelmann extrapolation, when applied to the preprocessed photospheric boundary, give a magnetic energy which is nearly equivalent to the value derived from the chromospheric boundary, but both underestimate the free energy above the photosphere by at least a factor of two. We discuss the interpretation of the preprocessed field in this context. When applying the NLFFF codes to solar data, the problems associated with Lorentz forces present in the low solar atmosphere must be recognized: the various codes will not necessarily converge to the correct, or even the same, solution. On 07/07/2007, the NLFFF team was saddened by the news that Tom Metcalf had died as the result of an accident. We remain grateful for having had the opportunity to benefit from his unwavering dedication to the problems encountered in attempting to understand the Sun’s magnetic field; Tom had completed this paper several months before his death, leading the team through the many steps described above.     

XUV Photometer System (XPS): Improved Solar Irradiance Algorithm Using CHIANTI Spectral Models

Solar soft X-ray (XUV) radiation is highly variable on all time scales and strongly affects Earth’s ionosphere and upper atmosphere; consequently, the solar XUV irradiance is important for atmospheric studies and for space weather applications. Although there have been several recent measurements of the solar XUV irradiance, detailed understanding of the solar XUV irradiance, especially its variability during flares, has been hampered by the broad bands measured in the XUV range. In particular, the simple conversion of the XUV photometer signal into irradiance, in which a static solar spectrum is assumed, overestimates the flare variations by more than a factor of two as compared to the atmospheric response to the flares. To address this deficiency in the simple conversion, an improved algorithm using CHIANTI spectral models has been developed to process the XUV Photometer System (XPS) measurements with its broadband photometers. Model spectra representative of quiet Sun, active region, and flares are combined to match the signals from the XPS and produce spectra from 0.1 to 40 nm in 0.1-nm intervals for the XPS Level 4 data product. The two XPS instruments are aboard NASA’s Solar Radiation and Climate Experiment (SORCE) and Thermosphere, Ionosphere, Mesosphere, Energetics, and Dynamics (TIMED) satellites. In addition, the XPS responsivities have been updated for the latest XPS data processing version. The new XPS results are consistent with daily variations from the previous simple conversion technique used for XPS and are also consistent with spectral measurements made at wavelengths longer than 27 nm. Most importantly, the XPS flare variations are reduced by factors of 2 – 4 at wavelengths shorter than 14 nm and are more consistent, for the first time, with atmospheric response to solar flares. Along with the details of the new XPS algorithm, several comparisons to dayglow and photoelectron measurements and model results are also presented to help verify the accuracy of the new XUV irradiance spectra.     

Spatially distributed modelling of soil erosion and sediment yield at regional scales in Spain

Initiated by the need to quantify erosion rates and the impacts of global changes on erosion, several attempts have been made to apply erosion models at regional scales. However, these models have often been directed towards on-site soil erosion estimates, emphasising sheet and rill erosion processes, and disregarding gully erosion, channel erosion and sediment transport. These models are therefore of limited use for the assessment of sediment yield, off-site impacts of erosion, and for the development of environmental management to control these impacts at regional scale. This study analyses and compares three spatially distributed models for the prediction of soil erosion and/or sediment yield at regional scales: the WATEM-SEDEM model that is based on the empirical Revised Universal Soil Loss Equation (RUSLE) in combination with a sediment transport equation, the physics-based Pan European Soil Erosion Risk Assessment model (PESERA), and a newly developed Spatially Distributed Scoring model (SPADS). The three models were applied to 61 Spanish drainage basins and model predictions were evaluated against data on measured reservoir sedimentation rates. Global data sets on land use, climate, elevation and soil characteristics were used as model input for WATEM-SEDEM and SPADS, whereas published soil erosion estimates of PESERA at 1 km² resolution were used directly. SPADS and WATEM-SEDEM provided best results after separate calibration for basins with a Sediment Delivery Ratio (SDR) higher than 5% and those with an SDR lower than 5%. In this way, SPADS explained 67% of variation in sediment yield, while WATEM-SEDEM explained 48% of the variation. PESERA represents a promising alternative to the use of empirical models at the regional scale as it can be applied to very diverse environments with little calibration. However, PESERA provides soil erosion rates and not sediment yield estimates. For most basins PESERA soil erosion rates vary between fifty and close to zero percent of total sediment yield. Two major factors may explain this discrepancy between modelled soil erosion rates and measured sediment yield. First, it may be that PESERA underestimates soil erosion under Mediterranean conditions, although PESERA soil erosion rates are of the same order of magnitude as erosion rates measured in erosion plot studies. Second, gully-, river channel erosion and sediment transport processes may be much more important than sheet- and rill erosion for regional scale sediment yield in these environments. These issues therefore require further attention in future model development. Although spatially lumped models provide better predictions of sediment yield at the basin scale, and while validation of the predicted spatial patterns of sources and sinks of sediment requires further research, spatially distributed models are expected to be of value to support management decisions regarding the assessment of on-site and off-site impacts of erosion at the regional scale.     

Climate, weather, and north polar observations from the Mars Reconnaissance Orbiter Mars Color Imager

The Mars Reconnaissance Orbiter observes Mars from a nearly circular, polar orbit. From this vantage point, the Mars Color Imager extends the ∼5 Mars years record of Mars Global Surveyor global, visible-wavelength multi-color observations of meteorological events and adds measurements at three additional visible and two ultraviolet wavelengths. Observations of the global distribution of ozone (which anti-correlates with water vapor) and water ice and dust clouds allow tracking of atmospheric circulation. Regional and local observations emphasize smaller scale atmospheric dynamics, especially those related to dust lifting and subsequent motion. Polar observations detail variations related to the polar heat budget, including changes in polar frosts and ices, and storms generated at high thermal contrast boundaries.     

The librations, shape, and icy shell of Europa

Europa, the smallest of the Galilean satellites, has a young icy surface and most likely contains an internal ocean. The primary objective of possible future missions to Europa is the unambiguous detection and characterization of a subsurface ocean. The thickness of the overlying icy shell provides important information on the thermal evolution of the satellite and on the interaction between the ocean and the surface, the latter being fundamental for astrobiology. However, the thickness is not well known, and estimates range from several hundred of meters to some ten of kilometers. Here, we investigate the use of libration (rotation variation) observations to study the interior structure of Europa and in particular its icy shell. A dynamical libration model is developed, which includes gravitational coupling between the icy shell and the heavy solid interior. The amplitude of the main libration signal at 3.55 days (the orbital period) is shown to depend on Europa's shape and structure. Models of the interior structure of Europa are constructed and the equatorial flattening of the internal layers, which are key parameters for the libration, are calculated by assuming that Europa is in hydrostatic equilibrium. Europa's flattened shape is assumed to be due to rotation and permanent tides, and we extend the classical Radau equation for rotationally flattened bodies to include also tidal deformation. We show that the presence of an ocean increases the amplitude of libration by about 10%, depending mainly on the thickness of the icy shell. Therefore, libration observations offer possibility of detection of a subsurface ocean in Europa and estimation of the thickness of its overlying icy shell.