Laboratory observations and simulations of phase reddening

The visible reflectance spectrum of many Solar System bodies changes with changing viewing geometry for reasons not fully understood. It is often observed to redden (increasing spectral slope) with increasing solar phase angle, an effect known as phase reddening. Only once, in an observation of the martian surface by the Viking 1 lander, was reddening observed up to a certain phase angle with bluing beyond, making the reflectance ratio as a function of phase angle shaped like an arch. However, in laboratory experiments this arch-shape is frequently encountered. To investigate this, we measured the bidirectional reflectance of particulate samples of several common rock types in the 400–1000 nm wavelength range and performed ray-tracing simulations. We confirm the occurrence of the arch for surfaces that are forward scattering, i.e. are composed of semi-transparent particles and are smooth on the scale of the particles, and for which the reflectance increases from the lower to the higher wavelength in the reflectance ratio. The arch shape is reproduced by the simulations, which assume a smooth surface. However, surface roughness on the scale of the particles, such as the Hapke and van Horn (Hapke, B., van Horn, H. [1963]. J. Geophys. Res. 68, 4545–4570) fairy castles that can spontaneously form when sprinkling a fine powder, leads to monotonic reddening. A further consequence of this form of microscopic roughness (being indistinct without the use of a microscope) is a flattening of the disk function at visible wavelengths, i.e. Lommel–Seeliger-type scattering. The experiments further reveal monotonic reddening for reflectance ratios at near-IR wavelengths. The simulations fail to reproduce this particular reddening, and we suspect that it results from roughness on the surface of the particles. Given that the regolith of atmosphereless Solar System bodies is composed of small particles, our results indicate that the prevalence of monotonic reddening and Lommel–Seeliger-type scattering for these bodies results from microscopic roughness, both in the form of structures built by the particles and roughness on the surface of the particles themselves. It follows from the singular Viking 1 observation that the surface in front of the lander was composed of semi-transparent particles, and was smooth on the scale of the particle size.

Uncalibrated Near-Light Photometric Stereo

In this work we solve the uncalibrated photometric stereo problem with lights placed near the scene. We investigate different image formation models and find the one that best fits our observations. Although the devised model is more complex than its far-light counterpart, we show that under a global linear ambiguity the reconstruction is possible up to a rotation and scaling, which can be easily fixed. We also propose a solution for reconstructing the normal map, the albedo, the light positions and the light intensities of a scene given only a sequence of near-light images. This is done in an alternating minimization framework which first estimates both the normals and the albedo, and then the light positions and intensities. We validate our method on real world experiments and show that a near-light model leads to a significant improvement in the surface reconstruction compared to the classic distant illumination case.

Spectrophotometric properties of the nucleus of comet 67P/Churyumov-Gerasimenko from the OSIRIS instrument onboard the ROSETTA spacecraft

Context. The Rosetta mission of the European Space Agency has been orbiting the comet 67P/Churyumov-Gerasimenko (67P) since August 2014 and is now in its escort phase. A large complement of scientific experiments designed to complete the most detailed study of a comet ever attempted are onboard Rosetta. Aims. We present results for the photometric and spectrophotometric properties of the nucleus of 67P derived from the OSIRIS imaging system, which consists of a Wide Angle Camera (WAC) and a Narrow Angle Camera (NAC). The observations presented here were performed during July and the beginning of August 2014, during the approach phase, when OSIRIS was mapping the surface of the comet with several filters at different phase angles (1.3 degrees-54 degrees). The resolution reached up to 2.1 m/px. Methods. The OSIRIS images were processed with the OSIRIS standard pipeline, then converted into I/F. radiance factors and corrected for the illumination conditions at each pixel using the Lommel-Seeliger disk law. Color cubes of the surface were produced by stacking registered and illumination-corrected images. Furthermore, photometric analysis was performed both on disk-averaged photometry in several filters and on disk-resolved images acquired with the NAC orange filter, centered at 649 ran, using Hapke modeling. Results. The disk-averaged phase function of the nucleus of 67P shows a strong opposition surge with a G parameter value of -0.13 +/- 0.01 in the HG system formalism and an absolute magnitude H-v(1, 1, 0) = 15.74 +/- 0.02 mag. The integrated spectrophotometry in 20 filters covering the 250-1000 nm wavelength range shows a red spectral behavior, without clear absorption bands except for a potential absorption centered at similar to 290 rim that is possibly due to SO2 ice. The nucleus shows strong phase reddening, with disk-averaged spectral slopes increasing from 11%/( 100 nm) to 16%/(100 nm) in the 1.3 degrees-54 degrees phase angle range. The geometric albedo of the comet is 6.5 +/- 0.2% at 649 nm, with local variations of up to similar to 16% in the Hapi region. From the disk-resolved images we computed the spectral slope together with local spectrophotometry and identified three distinct groups of regions (blue, moderately red, and red). The Hapi region is the brightest, the bluest in term of spectral slope, and the most active surface on the comet. Local spectrophotometry shows an enhancement of the flux in the 700-750 nm that is associated with coma emissions.

Geomorphology and spectrophotometry of Philae's landing site on comet 67P/Churyumov-Gerasimenko

Context. On 12 November 2014 the European mission Rosetta succeeded in delivering a lander, named Philae, on the surface of one of the smallest, low-gravity and most primitive bodies of the solar system, the comet 67P/Churyumov-Gerasimenko (67P). Aims. The aim of this paper is to provide a comprehensive geomorphological and spectrophotometric analysis of Philae's landing site (Agilkia) to give an essential framework for the interpretation of its in situ measurements. Methods. OSIRIS images, coupled with gravitational slopes derived from the 3D shape model based on stereo-photogrammetry were used to interpret the geomorphology of the site. We adopted the Hapke model, using previously derived parameters, to photometrically correct the images in orange filter (649.2 nm). The best approximation to the Hapke model, given by the Akimov parameter-less function, was used to correct the reflectance for the effects of viewing and illumination conditions in the other filters. Spectral analyses on coregistered color cubes were used to retrieve spectrophotometric properties. Results. The landing site shows an average normal albedo of 6.7% in the orange filter with variations of similar to 15% and a global featureless spectrum with an average red spectral slope of 15.2%/100 nm between 480.7 nm (blue filter) and 882.1 nm (near-IR filter). The spatial analysis shows a well-established correlation between the geomorphological units and the photometric characteristics of the surface. In particular, smooth deposits have the highest reflectance a bluer spectrum than the outcropping material across the area. Conclusions. The featureless spectrum and the redness of the material are compatible with the results by other instruments that have suggested an organic composition. The observed small spectral variegation could be due to grain size effects. However, the combination of photometric and spectral variegation suggests that a compositional differentiation is more likely. This might be tentatively interpreted as the effect of the efficient dust-transport processes acting on 67P. High-activity regions might be the original sources for smooth fine-grained materials that then covered Agilkia as a consequence of airfall of residual material. More observations performed by OSIRIS as the comet approaches the Sun would help interpreting the processes that work at shaping the landing site and the overall nucleus.

Mapping spectroscopic uncertainties into prospective methane retrieval errors from Sentinel-5 and its precursor

Sentinel-5 (S5) and its precursor (S5P) are future European satellite missions aiming at global monitoring of methane (CH4) column-average dry air mole fractions (XCH4). The spectrometers to be deployed onboard the satellites record spectra of sunlight backscattered from the Earth's surface and atmosphere. In particular, they exploit CH4 absorption in the shortwave infrared spectral range around 1.65 mu m (S5 only) and 2.35 mu m (both S5 and S5P) wavelength. Given an accuracy goal of better than 2% for XCH4 to be delivered on regional scales, assessment and reduction of potential sources of systematic error such as spectroscopic uncertainties is crucial. Here, we investigate how spectroscopic errors propagate into retrieval errors on the global scale. To this end, absorption spectra of a ground-based Fourier transform spectrometer (FTS) operating at very high spectral resolution serve as estimate for the quality of the spectroscopic parameters. Feeding the FTS fitting residuals as a perturbation into a global ensemble of simulated S5- and S5P-like spectra at relatively low spectral resolution, XCH4 retrieval errors exceed 0.6% in large parts of the world and show systematic correlations on regional scales, calling for improved spectroscopic parameters.

Model-based changes in global annual mean surface temperature change (Delta T_g) and radiative forcing due to land ice albedo changes (Delta R_[LI]) over the last 5 Myr, supplementary material

It is still an open question how equilibrium warming in response to increasing radiative forcing - the specific equilibrium climate sensitivity S - depends on background climate. We here present palaeodata-based evidence on the state dependency of S, by using CO2 proxy data together with a 3-D ice-sheet-model-based reconstruction of land ice albedo over the last 5 million years (Myr). We find that the land ice albedo forcing depends non-linearly on the background climate, while any non-linearity of CO2 radiative forcing depends on the CO2 data set used. This non-linearity has not, so far, been accounted for in similar approaches due to previously more simplistic approximations, in which land ice albedo radiative forcing was a linear function of sea level change. The latitudinal dependency of ice-sheet area changes is important for the non-linearity between land ice albedo and sea level. In our set-up, in which the radiative forcing of CO2 and of the land ice albedo (LI) is combined, we find a state dependence in the calculated specific equilibrium climate sensitivity, S[CO2,LI], for most of the Pleistocene (last 2.1 Myr). During Pleistocene intermediate glaciated climates and interglacial periods, S[CO2,LI] is on average ~ 45 % larger than during Pleistocene full glacial conditions. In the Pliocene part of our analysis (2.6-5 Myr BP) the CO2 data uncertainties prevent a well-supported calculation for S[CO2,LI], but our analysis suggests that during times without a large land ice area in the Northern Hemisphere (e.g. before 2.82 Myr BP), the specific equilibrium climate sensitivity, S[CO2,LI], was smaller than during interglacials of the Pleistocene. We thus find support for a previously proposed state change in the climate system with the widespread appearance of northern hemispheric ice sheets. This study points for the first time to a so far overlooked non-linearity in the land ice albedo radiative forcing, which is important for similar palaeodata-based approaches to calculate climate sensitivity. However, the implications of this study for a suggested warming under CO2 doubling are not yet entirely clear since the details of necessary corrections for other slow feedbacks are not fully known and the uncertainties that exist in the ice-sheet simulations and global temperature reconstructions are large.

Documentation of sediment cores during SONNE Expedition SO228 from the western tropical Pacific

The core descriptions (chapter 7) summarize the most important results of the analysis of each sediment core following procedures applied during ODP/IODP expeditions. All cores were opened, described, and color-scanned. In the core descriptions the first column displays the lithological data that are based on visual analysis of the core and are supplemented by information from binocular and smear slide analyses. The sediment classification largely follows ODP/IODP convention. Lithological names consist of a principal name based on composition, degree of lithification, and/or texture as determined from visual description and microscopic observations. In the structure column the intensity of bioturbation together with individual or special features (turbidites, volcanic ash layers, plant debris, shell fragments, etc.) is shown. The hue and chroma attributes of color were determined by comparison with the Munsell soil color charts and are given in the color column in the Munsell notation. A GretagMacbethTM Spectrolino spectrophotometer was used to measure percent reflectance values of sediment color at 36 wavelength channels over the visible light range (380-730 nm) on all of the cores. The digital reflectance data of the spectrophotometer readings were routinely obtained from the surface (measured in 1 cm steps) of the split cores (archive half). The Spectrolino is equipped with a measuring aperture with folding mechanism allowing an exact positioning on the split core and is connected to a portable computer. The data are directly displayed within the software package Excel and can be controlled simultaneously. From all the color measurements, for each core the red/blue ratio (700 nm/450 nm) and the lightness are shown together with the visual core description. The reflectance of individual wavelengths is often significantly affected by the presence of minor amounts of oxyhydroxides or sulphides. To eliminate these effects, we used the red/blue ratio and lightness.

Hourly averages of corrected light scattering integrals in six wavelength investigated during the Atlantic Expedition of RV "Meteor" 1969 (Table 2)

With a 6-channel integrating nephelometer spectral scattering properties of the atmospheric aerosol have been measured during the third part of the Atlantic Expedition 1969. A meridional cross section of light scattering integrals in the wavelength range 0.475 µm to 0.924 µm was recorded reaching from 10° S to 60° N along 30° W. With a new algorithm the time series of hourly scattering spectra was inverted yielding a first meridional cross section of the median radius of the number size distribution in situ. Three air mass regimes could be distinguished in the course of the experiment, the first one being the extremely clean air of the SE-trade south of the ITC. An abrupt increase in light scattering marked the hemispheric change when the ship entered the NE-trade which was heavily loaded with Sahara dust. North of the trade region the ship sailed through maritime North Atlantic air masses with highly variable light scattering and a slow decrease in median radius with latitude.

Time series of black carbon aerosol measurements from Neumayer Station, 1999-2011

Continuous black carbon (BC) observations were conducted from 1999 through 2009 by an Aethalometer (AE10) and from 2006 through 2011 by a Multi-Angle Absorption Photometer (MAAP) at Neumayer Station (NM) under stringent contamination control. Considering the respective observation period, BC concentrations measured by the MAAP were somewhat higher (median ± standard deviation: 2.1 ± 2.0 ng/m**3) compared to the AE10 results (1.6 ± 2.1 ng/m**3). Neither for the AE10 nor for the MAAP data set a significant long-term trend could be detected. Consistently a pronounced seasonality was observed with both instruments showing a primary annual maximum between October and November and a minimum in April with a maximum/minimum ratio of 4.5/1.6 = 3.8 and 2.7/0.64 = 4.2 for the MAAP and AE10 data, respectively. Occasionally a secondary summer maximum in January/February was visible. With the aim to assess the impact of BC on optical properties of the aerosol at NM, we evaluated the BC data along with particle scattering coefficients measured by an integrating nephelometer. We found the mean single scattering albedo of w550 = 0.992 ± 0.0090 (median: 0.994) at a wavelength of 550 nm with a range of values from 0.95 to 1.0.