Reduced-order modeling of light transport in tissue for real-time monitoring of brain hemodynamics using diffuse optical tomography

This paper proposes a new reconstruction method for diffuse optical tomography using reduced-order models of light transport in tissue. The models, which directly map optical tissue parameters to optical flux measurements at the detector locations, are derived based on data generated by numerical simulation of a reference model. The reconstruction algorithm based on the reduced-order models is a few orders of magnitude faster than the one based on a finite element approximation on a fine mesh incorporating a priori anatomical information acquired by magnetic resonance imaging. We demonstrate the accuracy and speed of the approach using a phantom experiment and through numerical simulation of brain activation in a rat's head. The applicability of the approach for real-time monitoring of brain hemodynamics is demonstrated through a hypercapnic experiment. We show that our results agree with the expected physiological changes and with results of a similar experimental study. However, by using our approach, a three-dimensional tomographic reconstruction can be performed in ∼3  s per time point instead of the 1 to 2 h it takes when using the conventional finite element modeling approach

Evaluation of seven European aerosol optical depth retrieval algorithms for climate analysis

Satellite data are increasingly used to provide observation-based estimates of the effects of aerosols on climate. The Aerosol-cci project, part of the European Space Agency's Climate Change Initiative (CCI), was designed to provide essential climate variables for aerosols from satellite data. Eight algorithms, developed for the retrieval of aerosol properties using data from AATSR (4), MERIS (3) and POLDER, were evaluated to determine their suitability for climate studies. The primary result from each of these algorithms is the aerosol optical depth (AOD) at several wavelengths, together with the Ångström exponent (AE) which describes the spectral variation of the AOD for a given wavelength pair. Other aerosol parameters which are possibly retrieved from satellite observations are not considered in this paper. The AOD and AE (AE only for Level 2) were evaluated against independent collocated observations from the ground-based AERONET sun photometer network and against “reference” satellite data provided by MODIS and MISR. Tools used for the evaluation were developed for daily products as produced by the retrieval with a spatial resolution of 10 × 10 km2 (Level 2) and daily or monthly aggregates (Level 3). These tools include statistics for L2 and L3 products compared with AERONET, as well as scoring based on spatial and temporal correlations. In this paper we describe their use in a round robin (RR) evaluation of four months of data, one month for each season in 2008. The amount of data was restricted to only four months because of the large effort made to improve the algorithms, and to evaluate the improvement and current status, before larger data sets will be processed. Evaluation criteria are discussed. Results presented show the current status of the European aerosol algorithms in comparison to both AERONET and MODIS and MISR data. The comparison leads to a preliminary conclusion that the scores are similar, including those for the references, but the coverage of AATSR needs to be enhanced and further improvements are possible for most algorithms. None of the algorithms, including the references, outperforms all others everywhere. AATSR data can be used for the retrieval of AOD and AE over land and ocean. PARASOL and one of the MERIS algorithms have been evaluated over ocean only and both algorithms provide good results.

Improvements in Aerosol Optical Depth Estimation Using Multiangle CHRIS/PROBA Images

A method has been developed to estimate aerosol optical depth (AOD) over land surfaces using high spatial resolution, hyperspectral, and multiangle Compact High Resolution Imaging Spectrometer (CHRIS)/Project for On Board Autonomy (PROBA) images. The CHRIS instrument is mounted aboard the PROBA satellite and provides up to 62 bands. The PROBA satellite allows pointing to obtain imagery from five different view angles within a short time interval. The method uses inversion of a coupled surface/atmosphere radiative transfer model and includes a general physical model of angular surface reflectance. An iterative process is used to determine the optimum value providing the best fit of the corrected reflectance values for a number of view angles and wavelengths with those provided by the physical model. This method has previously been demonstrated on data from the Advanced Along-Track Scanning Radiometer and is extended here to the spectral and angular sampling of CHRIS/PROBA. The values obtained from these observations are validated using ground-based sun-photometer measurements. Results from 22 image sets show an rms error of 0.11 in AOD at 550 nm, which is reduced to 0.06 after an automatic screening procedure.

Note: Active optical detection of cloud from a balloon platform

A disposable backscatter instrument is described for optical detection of cloud in the atmosphere from a balloon-carried platform. It uses an ultra-bright light emitting diode (LED) illumination source with a photodiode detector. Scattering of the LED light by cloud droplets generates a small optical signal which is separated from background light fluctuations using a lock-in technique. The signal to noise obtained permits cloud detection using the scattered LED light, even in daytime. The response is interpreted in terms of the equivalent visual range within the cloud. The device is lightweight (150 g) and low power (∼30 mA), for use alongside a conventional meteorological radiosonde.

Optical properties of Saharan dust aerosol and contribution from the coarse mode as measured during the Fennec 2011 aircraft campaign

New in-situ aircraft measurements of Saharan dust originating from Mali, Mauritania and Algeria taken during the Fennec 2011 aircraft campaign over a remote part of the Sahara Desert are presented. Size distributions extending to 300 μm are shown, representing measurements extending further into the coarse mode than previously published for airborne Saharan dust. A significant coarse mode was present in the size distribution measurements with effective diameter (deff) from 2.3 to 19.4 μm and coarse mode volume median diameter (dvc) from 5.8 to 45.3 μm. The mean size distribution had a larger relative proportion of coarse mode particles than previous aircraft measurements. The largest particles (with deff >12 μm, or dvc >25 μm) were only encountered within 1 km of the ground. Number concentration, mass loading and extinction coefficient showed inverse relationships to dust age since uplift. Dust particle size showed a weak exponential relationship to dust age. Two cases of freshly uplifted dust showed quite different characteristics of size distribution and number concentration. Single Scattering Albed (SSA) values at 550 nm calculated from the measured size distributions revealed high absorption ranging from 0.70 to 0.97 depending on the refractive index. SSA was found to be strongly related to deff. New instrumentation revealed that direct measurements, behind Rosemount inlets, overestimate SSA by up to 0.11 when deff is greater than 2 μm. This is caused by aircraft inlet inefficiencies and sampling losses. Previous measurements of SSA from aircraft measurements may also have been overestimates for this reason. Radiative transfer calculations indicate that the range of SSAs during Fennec 2011 can lead to underestimates in shortwave atmospheric heating rates by 2.0 to 3.0 times if the coarse mode is neglected. This will have an impact on Saharan atmospheric dynamics and circulation,which should be taken into account by numerical weather prediction and climate models.

A study on the sensitivities of simulated aerosol optical properties to composition and size distribution using airborne measurements

We present a flexible framework to calculate the optical properties of atmospheric aerosols at a given relative humidity based on their composition and size distribution. The similarity of this framework to climate model parameterisations allows rapid and extensive sensitivity tests of the impact of uncertainties in data or of new measurements on climate relevant aerosol properties. The data collected by the FAAM BAe-146 aircraft during the EUCAARI-LONGREX and VOCALS-REx campaigns have been used in a closure study to analyse the agreement between calculated and measured aerosol optical properties for two very different aerosol types. The agreement achieved for the EUCAARI-LONGREX flights is within the measurement uncertainties for both scattering and absorption. However, there is poor agreement between the calculated and the measured scattering for the VOCALS-REx flights. The high concentration of sulphate, which is a scattering aerosol with no absorption in the visible spectrum, made the absorption measurements during VOCALS-REx unreliable, and thus no closure study was possible for the absorption. The calculated hygroscopic scattering growth factor overestimates the measured values during EUCAARI-LONGREX and VOCALS-REx by ∼30% and ∼20%, respectively. We have also tested the sensitivity of the calculated aerosol optical properties to the uncertainties in the refractive indices, the hygroscopic growth factors and the aerosol size distribution. The largest source of uncertainty in the calculated scattering is the aerosol size distribution (∼35%), followed by the assumed hygroscopic growth factor for organic aerosol (∼15%), while the predominant source of uncertainty in the calculated absorption is the refractive index of organic aerosol (28–60%), although we would expect the refractive index of black carbon to be important for aerosol with a higher black carbon fraction.

Impact of atmospheric transport on the evolution of microphysical and optical properties of Saharan dust

Saharan dust affects the climate by altering the radiation balance and by depositing minerals to the Atlantic Ocean. Both are dependent on particle size. We present aircraft measurements comprising 42 profiles of size distribution (0.1–300 µm), representing freshly uplifted dust, regional aged dust, and dust in the Saharan Air Layer (SAL) over the Canary Islands. The mean effective diameter of dust in SAL profiles is 4.5 µm smaller than that in freshly uplifted dust, while the vertical structure changes from a low shallow layer (0–1.5 km) to a well-mixed deep Saharan dust layer (0–5 km). Size distributions show a loss of 60 to 90% of particles larger than 30 µm 12 h after uplift. The single scattering albedo (SSA) increases from 0.92 to 0.94 to 0.95 between fresh, aged, and SAL profiles: this is enough to alter heating rates by 26%. Some fresh dust close to the surface shows SSA as low as 0.85

Infrared optical coatings for the EarthCARE Multispectral Imager

The Earth Cloud, Aerosol and Radiation Explorer mission (EarthCARE) Multispectral Imager (MSI) is a radiometric instrument designed to provide the imaging of the atmospheric cloud cover and the cloud top surface temperature from a sun-synchronous low Earth orbit. The MSI forms part of a suite of four instruments destined to support the European Space Agency Living Planet mission on-board the EarthCARE satellite payload to be launched in 2016, whose synergy will be used to construct three-dimensional scenes, textures and temperatures of atmospheric clouds and aerosols. The MSI instrument contains seven channels: four solar channels to measure visible and short-wave infrared wavelengths, and three channels to measure infrared thermal emission. In this paper, we describe the optical layout of the infrared instrument channels, thin-film multilayer designs, the coating deposition method and the spectral system throughput for the bandpass interference filters, dichroic beam splitters, lenses and mirror coatings to discriminate wavelengths at 8.8, 10.8, & 12.0 µm. The rationale for the selection of thin-film materials, spectral measurement technique, and environmental testing performance are also presented.

Phase calibration of the EISCAT Svalbard Radar interferometer using optical satellite signatures

The link between natural ion-line enhancements in radar spectra and auroral activity has been the subject of recent studies but conclusions have been limited by the spatial and temporal resolution previously available. The next challenge is to use shorter sub-second integration times in combination with interferometric programmes to resolve spatial structure within the main radar beam, and so relate enhanced filaments to individual auroral rays. This paper presents initial studies of a technique, using optical and spectral satellite signatures, to calibrate the received phase of a signal with the position of the scattering source along the interferometric baseline of the EISCAT Svalbard Radar. It is shown that a consistent relationship can be found only if the satellite passage through the phase fringes is adjusted from the passage predicted by optical tracking. This required adjustment is interpreted as being due to the vector between the theoretical focusing points of the two antennae, i.e. the true radar baseline, differing from the baseline obtained by survey between the antenna foot points. A method to obtain a measurement of the true interferometric baseline using multiple satellite passes is outlined.

Simultaneous optical and radar signatures of poleward-moving auroral forms

Dayside poleward moving auroral forms (PMAFs) were detected between 06:30 and 07:00 UT on December 16, 1998, by the meridian scanning photometer and the all-sky camera at Ny Alesund, Svalbard. Simultaneous SuperDARN HF radar measurements permitted the study of the associated ionospheric velocity pattern. A good general agreement is observed between the location and movement of velocity enhancements (flow channels) and the PMAFs. Clear signatures of equatorward flow were detected in the vicinity of PMAFs. This flow is believed to be the signature of a return flow outside the reconnected Aux tube, as predicted by the Southwood model. The simulated signatures of this model reproduce globally the measured signatures, and differences with the experimental data can be explained by the simplifications of the model. Proposed schemes of the flow modification due to the presence of several flow channels and the modification of cusp and region 1 field-aligned currents at the time of sporadic reconnection events are shown to fit well with the observations.

Auroral bright spot sequence near 1400 MLT: Coordinated optical and ion drift observations

Optical observations of a dayside auroral brightening sequence, by means of all-sky TV cameras and meridian scanning photometers, have been combined with EISCAT ion drift observations within the same invariant latitude-MLT sector. The observations were made during a January 1989 campaign by utilizing the high F region ion densities during the maximum phase of the solar cycle. The characteristic intermittent optical events, covering ∼300 km in east-west extent, move eastward (antisunward) along the poleward boundary of the persistent background aurora at velocities of ∼1.5 km s−1 and are associated with ion flows which swing from eastward to westward, with a subsequent return to eastward, during the interval of a few minutes when there is enhanced auroral emission within the radar field of view. The breakup of discrete auroral forms occurs at the reversal (negative potential) that forms between eastward plasma flow, maximizing near the persistent arc poleward boundary, and strong transient westward flow to the south. The reported events, covering a 35 min interval around 1400 MLT, are embedded within a longer period of similar auroral activity between 0830 (1200 MLT) and 1300 UT (1600 MLT). These observations are discussed in relation to recent models of boundary layer plasma dynamics and the associated magnetosphere-ionosphere coupling. The ionospheric events may correspond to large-scale wave like motions of the low-latitude boundary layer (LLBL)/plasma sheet (PS) boundary. On the basis of this interpretation the observed spot size, speed and repetition period (∼10 min) give a wavelength (the distance between spots) of ∼900 km in the present case. The events can also be explained as ionospheric signatures of newly opened flux tubes associated with reconnection bursts at the magnetopause near 1400 MLT. We also discuss these data in relation to random, patchy reconnection (as has recently been invoked to explain the presence of the sheathlike plasma on closed field lines in the LLBL). In view of the lack of IMF data, and the existing uncertainty on the location of the open-closed field line boundary relative to the optical events, an unambiguous discrimination between the different alternatives is not easily obtained.

A coupled cloud physics–radiation parameterization of the bulk optical properties of cirrus and its impact on the Met Office unified model global atmosphere 5.0 configuration

A new coupled cloud physics–radiation parameterization of the bulk optical properties of ice clouds is presented. The parameterization is consistent with assumptions in the cloud physics scheme regarding particle size distributions (PSDs) and mass–dimensional relationships. The parameterization is based on a weighted ice crystal habit mixture model, and its bulk optical properties are parameterized as simple functions of wavelength and ice water content (IWC). This approach directly couples IWC to the bulk optical properties, negating the need for diagnosed variables, such as the ice crystal effective dimension. The parameterization is implemented into the Met Office Unified Model Global Atmosphere 5.0 (GA5) configuration. The GA5 configuration is used to simulate the annual 20-yr shortwave (SW) and longwave (LW) fluxes at the top of the atmosphere (TOA), as well as the temperature structure of the atmosphere, under various microphysical assumptions. The coupled parameterization is directly compared against the current operational radiation parameterization, while maintaining the same cloud physics assumptions. In this experiment, the impacts of the two parameterizations on the SW and LW radiative effects at TOA are also investigated and compared against observations. The 20-yr simulations are compared against the latest observations of the atmospheric temperature and radiative fluxes at TOA. The comparisons demonstrate that the choice of PSD and the assumed ice crystal shape distribution are as important as each other. Moreover, the consistent radiation parameterization removes a long-standing tropical troposphere cold temperature bias but slightly warms the southern midlatitudes by about 0.5 K.

Polymorphism and optical properties in [NH4][InSe2]

The solvothermal synthesis and characterization of two indium selenides with stoichiometry [NH4][InSe2] is described. Yellow [NH4][InSe2] (1), which exhibits a layered structure, was initially prepared in an aqueous solution of trans-1,4-diaminocyclohexane, and subsequently using a concentrated ammonia solution. A red polymorph of one-dimensional character, [NH4][InSe2] (2), was obtained using 3,5-dimethylpyridine as solvent. [NH4][InSe2] (1) crystallizes in the non-centrosymmetric space group Cc (a=11.5147(6), b=11.3242(6), c=15.9969(9) Å and β=100.354(3)°). The structural motif of the layers is the In4Se10 adamantane unit, composed of four corner-linked InSe4 tetrahedra. These units are linked by their corners, forming [InSe2]− layers which are stacked back to back along the c-direction, and interspaced by [NH4]+cations. The one-dimensional polymorph, (2), crystallizes in the tetragonal space group, I4/mcm (a=8.2519(16), c=6.9059 (14) Å). This structure contains infinite chains of edge-sharing InSe4 tetrahedra separated by [NH4]+ cations.

Determination of the embedded thermo-optical expansion coefficients of PbTe and ZnSe thin film infrared multilayers

This paper reports the first derived thermo-optical properties for vacuum deposited infrared thin films embedded in multilayers. These properties were extracted from the temperature-dependence of manufactured narrow bandpass filters across the 4-17 µm mid-infrared wavelength region. Using a repository of spaceflight multi-cavity bandpass filters, the thermo-optical expansion coefficients of PbTe and ZnSe were determined across an elevated temperature range 20-160 ºC. Embedded ZnSe films showed thermo-optical properties similar to reported bulk values, whilst the embedded PbTe films of lower optical density, deviate from reference literature sources. Detailed knowledge of derived coefficients is essential to the multilayer design of temperature-invariant narrow bandpass filters for use in non-cooled infrared detection systems. We further present manufacture of the first reported temperature-invariant multi-cavity narrow bandpass filter utilizing PbS chalcogenide layer material.