Decoupling error for the atmospheric correction in ocean color remote sensing algorithms

This paper studies the decoupling error associated with the atmospheric correction procedures in the ocean color remote sensing algorithms. The decoupling error is caused by the lack of proper consideration of multiple scattering between the atmospheric and ocean components. In other words, the atmosphere and ocean are not coupled properly. A vector radiative transfer model for the coupled atmosphere and ocean (CAO) system based on the successive order of scattering (SOS) method is used to study the error. The inherent optical properties (IOPs) of the ocean are provided by the most updated bio-optical models. Two wavelengths are used in the study, 412 and 555 nm. For a detector located just above the ocean interface, the decoupling errors range from 0.3% to 7% at 412 nm; and from 0.3% to 3 % at 555 nm for zenith viewing angles smaller than 70°. The decoupling errors are significantly larger for larger zenith viewing angles for this detector. For a detector at the top of the atmosphere (TOA), it is hard to separate the decoupling error from the error introduced by the diffuse transmittance. If we assume the upwelling radiance is uniform just below the ocean surface when estimating the diffuse transmittance, the decoupling errors are from ?4% to 8% for zenith viewing angles smaller than 70°; and negative decoupling errors show up at mainly large zenith viewing angles.     

Development of a GPU-based high-performance radiative transfer model for the Infrared Atmospheric Sounding Interferometer (IASI)

Satellite-observed radiance is a nonlinear functional of surface properties and atmospheric temperature and absorbing gas profiles as described by the radiative transfer equation (RTE). In the era of hyperspectral sounders with thousands of high-resolution channels, the computation of the radiative transfer model becomes more time-consuming. The radiative transfer model performance in operational numerical weather prediction systems still limits the number of channels we can use in hyperspectral sounders to only a few hundreds. To take the full advantage of such high-resolution infrared observations, a computationally efficient radiative transfer model is needed to facilitate satellite data assimilation. In recent years the programmable commodity graphics processing unit (GPU) has evolved into a highly parallel, multi-threaded, many-core processor with tremendous computational speed and very high memory bandwidth. The radiative transfer model is very suitable for the GPU implementation to take advantage of the hardware’s efficiency and parallelism where radiances of many channels can be calculated in parallel in GPUs.In this paper, we develop a GPU-based high-performance radiative transfer model for the Infrared Atmospheric Sounding Interferometer (IASI) launched in 2006 onboard the first European meteorological polar-orbiting satellites, METOP-A. Each IASI spectrum has 8461 spectral channels. The IASI radiative transfer model consists of three modules. The first module for computing the regression predictors takes less than 0.004% of CPU time, while the second module for transmittance computation and the third module for radiance computation take approximately 92.5% and 7.5%, respectively. Our GPU-based IASI radiative transfer model is developed to run on a low-cost personal supercomputer with four GPUs with total 960 compute cores, delivering near 4 TFlops theoretical peak performance. By massively parallelizing the second and third modules, we reached 364× speedup for 1 GPU and 1455× speedup for all 4 GPUs, both with respect to the original CPU-based single-threaded Fortran code with the –O₂ compiling optimization. The significant 1455× speedup using a computer with four GPUs means that the proposed GPU-based high-performance forward model is able to compute one day’s amount of 1,296,000 IASI spectra within nearly 10 min, whereas the original single CPU-based version will impractically take more than 10 days. This model runs over 80% of the theoretical memory bandwidth with asynchronous data transfer. A novel CPU–GPU pipeline implementation of the IASI radiative transfer model is proposed. The GPU-based high-performance IASI radiative transfer model is suitable for the assimilation of the IASI radiance observations into the operational numerical weather forecast model.     

Propagation effects of a fractal rough ocean surface on the vertical electric field generated by lightning return strokes

Based on an improved two-dimension (2D) fractal model of rough ocean surface, the propagation effects of the rough ocean surface on the vertical electric fields generated by lightning return strokes are analyzed. The results show that the rough ocean surface has much effect on the electric field derivatives, but has no or little effect on the field peaks. The frequency above 10 MHz is attenuated significantly by the rough ocean surface, and the rapid attenuation of frequency above 10 MHz in the experimentally obtained spectrum may be taken into account the errors introduced by the roughness of the ocean surface.     

Wavelength selection of bidirectional laser transmission based on Monte Carlo method

The laser detection technology in uncertain and dynamic environments is of utmost importance in many fields. A model of transient radiative transfer of bidirectional path laser based on Monte Carlo method is developed to investigate the optimum wavelength of active detector at complex atmospheric conditions. The radiative parameters of atmosphere are calculated by HITRAN database and Mie theory at several typical atmospheric conditions including the standard atmosphere, urban aerosol, and radiation fog. Transmission characteristics for five spectral bands at the above atmospheric conditions are calculated by this model. The optimal transmission ability occurred in bands 0.2–0.5, 1.4–1.6, and 0.75–1.25 μm on the condition of standard atmosphere, urban aerosol, and radiation fog, respectively. All results provide effective reference and basic support for choosing the optimal spectral band for active detection.     

A fast linearized pseudo-spherical two orders of scattering model to account for polarization in vertically inhomogeneous scattering-absorbing media

We calculate the reflection matrix for the first two orders of scattering in a vertically inhomogeneous, scattering-absorbing medium. We take full account of polarization and perform a complete linearization (analytic differentiation) of the reflection matrix with respect to both the inherent optical properties of the medium and the surface reflection condition. Further, we compute a scalar-vector correction to the total intensity due to the effect of polarization; this correction is also fully linearized. The solar beam attenuation has been computed for a pseudo-spherical atmosphere.Results from the two orders of scattering (2OS) model have been tested against scalar intensities for an inhomogeneous atmosphere, and against Stokes vector results for a homogeneous atmosphere. We have also performed backscatter simulations of reflected sunlight in the O₂A band for a variety of geometries, and compared our results with those from a full vector multiple scattering code. Our results are exact in the center of strong lines and most inaccurate in the continuum, where polarization is least significant. The s- and p-polarized radiances are always computed very accurately. The effect of gas absorption optical depth, solar zenith angle, viewing geometry, surface albedo and wind speed (in the case of ocean glint) on the intensity, polarization and corresponding weighting functions have been investigated. It is shown that the 2OS model provides fast and reliably accurate polarization corrections to the scalar-model radiance and weighting function fields. The model can be implemented in operational retrieval algorithms as an adjunct radiative transfer code to deal with polarization effects.     

A k-distribution-based radiation code and its computational optimization for an atmospheric general circulation model

The gas absorption process scheme in the broadband radiative transfer code “mstrn8”, which is used to calculate atmospheric radiative transfer efficiently in a general circulation model, is improved. Three major improvements are made. The first is an update of the database of line absorption parameters and the continuum absorption model. The second is a change to the definition of the selection rule for gas absorption used to choose which absorption bands to include. The last is an upgrade of the optimization method used to decrease the number of quadrature points used for numerical integration in the correlated k-distribution approach, thereby realizing higher computational efficiency without losing accuracy. The new radiation package termed “mstrnX” computes radiation fluxes and heating rates with errors less than 0.6 W/m² and 0.3 K/day, respectively, through the troposphere and the lower stratosphere for any standard AFGL atmospheres. A serious cold bias problem of an atmospheric general circulation model using the ancestor code “mstrn8” is almost solved by the upgrade to “mstrnX”.     

2.7 μm fluorescence radiative dynamics and energy transfer between Er3+ and Tm3+ ions in fluoride glass under 800 nm and 980 nm excitation

A detailed study of the fluorescence radiative dynamics and energy transfer processes between Er and Tm ions in the Er³⁺/Tm³⁺ doped fluoride glass is reported. The fluorescence properties of 2.7 μm emission, other infrared and visible emissions are investigated under different selective laser excitations. Three Judd–Ofelt intensity parameters, energy transfer microparameters and efficiency have been determined and discussed. It is found that present Er³⁺/Tm³⁺ doped fluoride glass possesses large calculated emission cross section (8.98×10^–21 cm²) around 2.7 μm. The more suitable pumping scheme for laser applications at 2.7 μm laser is 980 nm excitation for Er³⁺/Tm³⁺ doped fluoride glass.     

Mid-infrared luminescence and energy transfer of Dy3+/Tm3+ doped fluorophosphate glass

This work reports the observation of emissions at 2.9 μm, 1.8 μm and 1.47 μm from Dy³⁺/Tm³⁺ codoped fluorophosphate glass upon excitation of a conventional 800 nm laser diode. Judd–Ofelt intensity parameters and radiative properties of Dy³⁺ ions in present glasses were calculated using the Judd–Ofelt theory. The mechanism and microparameters of energy transfer processes were investigated based on photoluminescence performance and lifetime measurements. The Dy³⁺/Tm³⁺ codoped fluorophosphate glass possessing advantageous spectroscopic characteristics as well as excellent thermal stability is a promising candidate for an efficient 2.9 μm laser.     

Discrete-ordinate method with matrix exponential for a pseudo-spherical atmosphere: Scalar case

We present a discrete-ordinate algorithm using the matrix-exponential solution for pseudo-spherical radiative transfer. Following the finite-element technique we introduce the concept of layer equation and formulate the discrete radiative transfer problem in terms of the level values of the radiance. The layer quantities are expressed by means of matrix exponentials, which are computed by using the matrix eigenvalue method and the Padé approximation. These solution methods lead to a compact and versatile formulation of the radiative transfer. Simulated nadir and limb radiances for an aerosol-loaded atmosphere and a cloudy atmosphere are presented along with a discussion of the model intercomparisons and timings.     

Broadband 2.84 μm luminescence properties and Judd–Ofelt analysis in Dy3+ doped ZrF4–BaF2–LaF3–AlF3–YF3 glass

2.84 μm luminescence with a bandwidth of 213 nm is obtained in Dy³⁺ doped (ZrF₄–BaF₂–LaF₃–AlF₃–YF₃) ZBLAY glass. Three intensity parameters and radiative properties have been determined from the absorption spectrum based on the Judd–Ofelt theory. The 2.84 μm emission characteristics and energy transfer mechanism upon excitation of a conventional 808 nm laser diode are investigated. The prepared Dy³⁺ doped ZBLAY glass possessing high predicted spontaneous transition probability (45.92 s^?1) along with large calculated emission cross section (1.17×10^?20 cm²) has potential applications in 2.8 μm laser.     

Low loss propagation in slow light photonic crystal waveguides at group indices up to 60

We have designed slow light photonic crystal waveguides operating in a low loss and constant dispersion window of Δλ = 2 nm around λ = 1565 nm with a group index of n_g = 60. We experimentally demonstrate a relatively low propagation loss, of 130 dB/cm, for waveguides up to 800 μm in length. This result is particularly remarkable given that the waveguides were written on an electron-beam lithography tool with a writefield of 100 μm that exhibits stitching errors of typically 10–50 nm. We reduced the impact of these stitching errors by introducing “slow–fast–slow” mode conversion interfaces and show that these interfaces reduce the loss from 320 dB/cm to 130 dB/cm at n_g = 60. This significant improvement highlights the importance of the slow–fast–slow method and shows that high performance slow light waveguides can be realised with lengths much longer than the writing field of a given e-beam lithography tool.     

New probing techniques of radiative shocks

Radiative shock waves propagating in xenon at a low pressure have been produced using 60 joules of iodine laser (λ = 1.315 μm) at PALS center. The shocks have been probed by XUV imaging using a Zn X-raylaser (λ = 21 nm) generated with a 20-ns delay after the shock creating pulse. Auxiliary high-speed silicon diodes allowed performing space- and time-resolved measurement of plasma self-emission in the visible and XUV. The results show the generation of a shock wave propagating at 60 km/s preceded by a radiative precursor. This demonstrates the feasibility of radiative shock generation using high power infrared lasers and the use of XRL backlighting as a suitable diagnostic for shock imaging.     

A new approach to the retrieval of surface properties from earthshine measurements

Instruments such as the MODIS and MISR radiometers on EOS AM-1, and POLDER on ADEOS have been deployed for the remote sensing retrieval of surface properties. Typically, retrieval algorithms use linear combinations of semi-empirical bidirectional reflectance distribution function (BRDF) kernels to model surface reflectance. The retrieval proceeds in two steps; first, an atmospheric correction relates surface BRDF to top-of-atmosphere (TOA) reflectances, then regression is used to establish the linear coefficients used in the kernel combination. BRDF kernels may also depend on a number of physical or empirical non-linear parameters (e.g. ocean wind speed for a specular BRDF); such parameters are usually assumed known. A major source of error in this retrieval comes from lack of knowledge of planetary boundary layer (PBL) aerosol properties.In this paper, we present a different approach to surface property retrieval. For the radiative transfer simulations, we use the discrete ordinate LIDORT model, which has the capability to generate simultaneous fields of radiances and weighting functions in a multiply scattering multi-layer atmosphere. Surface-atmosphere coupling due to multiple scattering and reflection effects is treated in full; the use of an atmospheric correction is not required. Further, it is shown that sensitivities of TOA reflectances to both linear and non-linear surface BRDF parameters may be established directly by explicit analytic differentiation of the discrete ordinate radiative transfer equations. Surface properties may thus be retrieved directly and conveniently from satellite measurements using standard non-linear fitting methods. In the fitting for BRDF parameters, lower-boundary aerosol properties can either be retrieved as auxiliary parameters, or they can be regarded as forward model parameter errors. We present examples of simulated radiances and surface/aerosol weighting functions for combinations of multi-angle measurements at several different wavelengths, and we perform some examples of self-consistent non-linear fitting to demonstrate feasibility for this kind of surface property retrieval.     

Picosecond molecular switch based on the influence of photogenerated electric fields on optical charge transfer transitions

A concept for a picosecond molecular switch is demonstrated using a photoinduced electron transfer reaction in a covalently linked, fixed distance donor–acceptor molecule D–A linked to a perylene-3,4-dicarboximide chromophore, C. The chromophore C possesses a strong charge transfer transition in its optical spectrum. Selective excitation of C within D–A–C using 530 nm, 130 fs laser pulses produces^1 * C, which undergoes singlet–singlet energy transfer to produce^1 * D, which in turn transfers an electron to A. If the D–A–C system is selectively excited with 416 nm, 130 fs laser pulses to produce D ⁺ – A ⁻ –C prior to excitation of C with 530 nm, 130 fs laser pulses, a 25% lower yield of^1 * C is generated. The intense local electric field produced by D ⁺ – A ⁻ causes a 15 nm electrochromic red shift of the charge transfer absorption of C. Thus, the absorption of C at 530 nm is significantly diminished by the presence of D ⁺ – A ⁻ . The need to use two laser pulses with different wavelengths to observe these effects, and the resulting picosecond time response makes it possible to consider applications of this concept in the design of molecular switches.     

海洋气溶胶多角度偏振辐射特性研究

海洋气溶胶是对流层气溶胶的重要组成部分,对全球的辐射收支平衡及气侯变化均有重要的影响。评价气溶胶的直接与间接辐射效应需要对气溶胶的性质进行深入的研究。多角度偏振为气溶胶光学物理性质研究提供了新的方法。在对可见光550 nm和近红外860 nm波段处海洋气溶胶的光学性质的研究基础上,采用矢量辐射传输模型模拟了TOA（top of atmosphere）反射率和偏振反射率与下垫面性质、观测方位角、气溶胶光学厚度之间的敏感性。模拟结果表明,海洋型气溶胶的多角度偏振信息可以有效地体现气溶胶的光学性质,可以利用多角度偏振遥感信息进行海洋气溶胶的反演,为利用多角度偏振遥感数据进行海洋气溶胶提供了理论基础。     

A single-scattering approximation for infrared radiative transfer in limb geometry in the Martian atmosphere

We present a single-scattering approximation for infrared radiative transfer in limb geometry in the Martian atmosphere. It is based on the assumption that the upwelling internal radiation field is dominated by a surface with a uniform brightness temperature. It allows the calculation of the scattering source function for individual aerosol types, mixtures of aerosol types, and mixtures of gas and aerosol. The approximation can be applied in a Curtis-Godson radiative transfer code and is used for operational retrievals from Mars Climate Sounder measurements. Radiance comparisons with a multiple scattering model show good agreement in the mid- and far-infrared although the approximate model tends to underestimate the radiances in realistic conditions of the Martian atmosphere. Relative radiance differences are found to be about 2% in the lowermost atmosphere, increasing to ∼10% in the middle atmosphere of Mars. The increasing differences with altitude are mostly due to the increasing contribution to limb radiance of scattering relative to emission at the colder, higher atmospheric levels. This effect becomes smaller toward longer wavelengths at typical Martian temperatures. The relative radiance differences are expected to produce systematic errors of similar magnitude in retrieved opacity profiles.     

Discrete-ordinate radiative transfer in a stratified medium with first-order rotational Raman scattering

Rotational Raman scattering (RRS) by air molecules in the Earth's atmosphere is predominantly responsible for the Ring effect: Fraunhofer and absorption-feature filling-in observed in UV/visible backscatter spectra. Accurate determination of RRS effects requires detailed radiative transfer (RT) treatment. In this paper, we demonstrate that the discrete-ordinate RT equations may be solved analytically in a multi-layer multiple scattering atmosphere in the presence of RRS treated as a first-order perturbation. Based on this solution, we develop a generic pseudo-spherical RT model LIDORT-RRS for the determination of backscatter radiances with RRS included; the model will generate output at arbitrary viewing geometry and optical thickness. Model comparisons with measured RRS filling-in effects from OMI observations show very good agreement. We examine telluric RRS filling-in effects for satellite-view backscatter radiances in a spectral range covering the ozone Huggins absorption bands. The model is also used to investigate calcium H and K Fraunhofer filling-in through cloud layers in the atmosphere.     

A vector radiative transfer model for coupled atmosphere and ocean systems with a rough interface

We report on an exact vector (polarized) radiative transfer (VRT) model for coupled atmosphere and ocean systems. This VRT model is based on the successive order of scattering (SOS) method, which virtually takes all the multiple scattering processes into account, including atmospheric scattering, oceanic scattering, reflection and transmission through the rough ocean surface. The isotropic Cox–Munk wave model is used to derive the ref and transmission matrices for the rough ocean surface. Shadowing effects are included by the shadowing function. We validated the SOS results by comparing them with those calculated by two independent codes based on the doubling/adding and Monte Carlo methods. Two error analyses related to the ocean color remote sensing are performed in the coupled atmosphere and ocean systems. One is the scalar error caused by ignoring the polarization in the whole system. The other is the error introduced by ignoring the polarization of the light transmitted through the ocean interface. Both errors are significant for the cases studied. This code fits for the next generation of ocean color study because it converges fast for absorbing medium as, for instance, ocean.     

Dependence of cross-relaxation on temperature and concentration from the 1D2 level of Pr3+ in YPO4

We report on the role of cross-relaxation in the decay of the ¹D₂ level of trivalent Pr in YPO₄ in crystals with Pr concentrations of 0.1%, 1%, 2%, and 5%. We have found that the ¹D₂ level decay is purely radiative in the low-doped system. As the Pr concentration is increased, the ¹D₂ luminescence is quenched due to a cross-relaxation energy transfer between two Pr ions. The temporal behavior of the ¹D₂ luminescence following pulsed excitation has been monitored in each sample at temperatures between 30 K and 300 K, and all decay curves were fit to the Yokota–Tanimoto model. The decay times decrease as temperature increases, due to an increase in both the radiative rate and the energy transfer rate with temperature. There is little evidence of diffusion at any temperature, even in the more concentrated samples. We have also fit the decay curves using the LumiTrans computer simulation. A comparison of the fits to the decay curves of the two methods is presented.