Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Impact of solar background radiation on the accuracy of wind observations of spaceborne Doppler wind lidars based on their orbits and optical parameters

Open Access Open Access

Abstract

Due to the quantum properties of light, solar background radiation (SBR) is the main source of noise in daytime wind observations of spaceborne Doppler wind lidars (DWLs). In previous works, the impact of SBR on the observation accuracy of spaceborne lidars was assessed mainly using the default or worst-case scenarios. We assessed the impact of SBR on the observations of spaceborne DWLs using the global distributions of SBR in summer and winter, which were obtained based on their orbit parameters, view geometry and optical parameters. Three experiments illustrate that the uncertainty in wind observations increases with an increase in the quantiles of SBR. The uncertainties of the whole profiles of wind are greater than 2 m s−1 in the troposphere and 3 m s−1 in the stratosphere when the quantile of the SBR is greater than 85% in summer and 95% in winter, which do not satisfy the accuracy expectations of the European Space Agency (ESA) for spaceborne DWLs. The facts indicate that the impact of SBR cannot be negligible for the observations of spaceborne DWLs. Based on the orbit parameters, view geometry, and optical parameters of new spaceborne DWLs, engineers can assess the impact of SBR on the accuracy of wind observations from a global perspective using the method proposed in this paper.

© 2019 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

1. Introduction

For a long time, global wind field was mainly measured by the network of radiosondes, wind profiler radars, et al. However, observations of wind over large portions of the Tropics and major oceans are quite incomplete owing to the limits imposed by costs and the geography [1]. The first spaceborne Doppler wind lidar (DWL) mission ADM-Aeolus (Atmospheric Dynamics Mission, ADM) operated by the European Space Agency (ESA) was launched successfully on August 22, 2018 to close the gaps in the global wind observations, and it is expected to improve the accuracy of numerical weather prediction system (NWP) forecasts. The study of Weissmann et al. shows that the forecasts errors of geopotential height, wind, and humidity in the troposphere can be reduced with the high accuracy and spatial resolution wind profiles provided by Aeolus [2]. In the Tropics, only the observed wind can reconstruct the proper atmospheric mass/wind balance [3]. Outside the Tropics, the accurate wind observations appear to be more beneficial than temperature observations in the improvement of NWP forecasts quality [4]. The studies have confirmed the significance of spaceborne DWLs.

However, the observations of spaceborne DWLs do not always improve the NWP forecasts. Experiments in meteorological data assimilation shows that the poor quality observations often fail to have a significant beneficial impact on NWP. According to the accuracy requirement of the ESA for spaceborne DWLs, the uncertainty of the wind observations in the boundary layer, troposphere, and stratosphere should be less than 1, 2, and 3 m s−1, respectively [1], which is the criterion used in this paper to assess the observation accuracy of spaceborne DWLs. The solar background radiation (SBR) is the main source of noise in daytime observations and, thus, have a negative effect on the spaceborne direct detection wind lidars observations. Firstly, the SBR increases the minimum detectability of aerosol/cloud layers for a given atmospheric scene. Secondly, it lowers the signal-to-noise ratio (SNR) and the accuracy of the observed wind speed. The simulation of McGill et al. shows that the uncertainty of the observed wind speed during daytime can be 10%~20% higher than during nighttime [5]. The ESA have tested the effect of the solar background noise (SBN) on the Aeolus wind observations using three SBR scenarios: the nighttime condition, default cases, and worst cases which equals to 154 mW m−2 sr−1 nm−1. The results show that the largest uncertainty of Rayleigh wind without the impact of SBR is 5 m s−1, whereas it is up to 11 m s−1 in the worst cases on the Rayleigh channel. The experiments show that the SBR has a considerable influence on the accuracy of daytime spaceborne DWLs observations on the Rayleigh channel. In addition, the differences of uncertainty in the Mie channel between nighttime condition and the worst cases are small, which indicates that SBR has a lower impact on the Mie channel [6]. To minimize the impact of SBR on the Rayleigh channel of Aeolus, the dawn-dusk orbit has been selected [7].

Furthermore, the SBR has an impact on other functions of spaceborne lidars, such as the measurement of aerosol, methane, and CO2. For example, differences in the measurements taken with the cloud–aerosol lidar with orthogonal polarization (CALIOP) and the ground based lidar show that the uncertainty of the attenuated backscatter coefficient in daytime is 14% larger than that in nighttime. Further, the comparisons between CALIOP and airborne high spectral resolution lidar measurements show that the uncertainty of attenuated backscatter coefficients in daytime and nighttime are approximately 20% and 6.3%, respectively, which shows the negative influence of SBR [8]. For spaceborne methane lidars, the SBR at 1645.6 nm is small at a sun-synchronous dawn-dusk polar orbit. However, it would also decrease the SNR of returned signal from methane gas [9,10]. Assuming the spaceborne CO2 lidar operates on a sun-synchronous orbit with an equator-crossing time of approximately 1:30 a.m./p.m., the uncertainty of the measured CO2 density caused by SBN is approximately 0.9 ppmv [11].

As is mentioned above, SBR would affect the observations of many kinds of spaceborne lidars. In this study, we assess the impacts of SBR on spaceborne DWLs based on the global distributions of SBR, unlike previous works which assessed it based on the default or worst cases. In the assessment, the orbit parameters and view geometry of the spaceborne DWLs are set according to that of Aeolus. In the future, spaceborne DWLs with different orbit parameters and view geometries will be launched, which, by implementing the method proposed in this paper, will allow engineers to assess the impact of SBR more accurately. The remainder of this paper is organized as follows: in section 2, the methodology to assess the impact of SBR on the wind observations is described; section 3 presents the databases used to derive the global distributions of SBR; section 4 illustrates the global distributions of SBR at 355 nm in summer and winter; in section 5 we describe the implementation of the proposed method with three experiments; and in section 6 we summarize the results.

2. Assessment methodology

2.1. Poisson noise

In spaceborne lidars, photodetectors are used to detect the intensity of the signal backscattered from atmosphere. However, even if the backscatter signal has constant energy, the number of the photoelectrons detected by the photodetectors is inherently uncertain due to the quantum nature of light. Theoretical studies the number of photoelectrons emitted on the detectors follow a Poisson distribution whose variance is equal to the mean, as follows:

Δnp2¯=n¯p
where n¯p denotes the mean number of photoelectrons emitted and Δnp is the standard deviation [12]. The noise of the emitted photoelectrons caused by the quantum nature of light is also called the Poisson noise.

The laser signal backscattered from atmosphere would direct toward the telescopes of spaceborne DWLs along with the signal arisen from the SBR. And the SBR would affect the uncertainty of the wind speed observations through Poisson noise [13]. In operation mode, we subtract the photon count caused by SBR (also called SBN in this paper) from the total signal received on the photodetectors [14,15]. In addition, the SBN signal can be obtained by averaging of the signal derived from high altitudes of the atmosphere (generally greater than 65 km). However, due to Poisson noise, the SBR would also affect the accuracy of observed wind speed after the subtraction. Assuming the energy of SBR can excite n¯p photon counts, the detector of the lidar may receive np photon counts with uncertainty of Δnp. Through the averaged photon counts detected from high altitudes of the atmosphere, we subtract n¯p photon counts from the total photon counts received by detector. Therefore, the SBR adds uncertainty of Δnp to the total backscatter signal from atmosphere, which leads to random error in wind observations.

2.2. Theoretical uncertainty in wind observations on the Rayleigh channel

Several techniques can be used to retrieve the Doppler shift for Rayleigh channel, such as the edge technique and double-edge techniques [16,17]. The latter was used in this paper, similar to the Aeolus [18]. According to the study of Tan et al. [19], the uncertainty of wind observations on Rayleigh channel can be estimated as

σHLOS=(vHLOSRATMσRATM)2+(vHLOSTσT)2+(vHLOSPσP)2+(vHLOSρσρ)2
where vHLOS is the horizontal line-of-sight (HLOS) wind component, vHLOS/RATM, νHLOS/T, νHLOS/P, and νHLOS/ρ are the local sensitivities of the νHLOS wind observations on Rayleigh response function RATM, temperature T, pressure P, and scattering ratio ρ, respectively. The Rayleigh response function RATM can be written as
RATM=NANBNA+NB
where NA and NB denote the photon counts received by Rayleigh channels A and B (also called useful signal in this paper).

The SBR affects the uncertainty in wind speed observations mainly through the response function. In the operation mode, the values of temperature and pressure are obtained from the European Centre for Medium-Range Weather Forecasts (ECWMF) data assimilation system. The scattering ratio can be derived from the intensity of the signal measured by the Rayleigh and Mie channels [20]. In this paper, we assumed that the values of temperature, pressure, and scattering ratio are accurate (no uncertainty). The uncertainty of wind speed observations for Rayleigh channel was calculated as

σHLOS=vHLOSRATMσRATM

In Eq. (4), the partial derivate of the νHLOS with respect to the Rayleigh response RATM is a function of temperature and pressure, as is shown in Fig. 1.

 figure: Fig. 1

Fig. 1 Lookup table of the partial derivative of HLOS wind component to Rayleigh response, which is determined by temperature and pressure.

Download Full Size | PDF

The uncertainty in the response RATM can be obtained by

σRATM=2NA+NBNB2σA2+NA2σB2,
where σA and σB are the uncertainties in NA and NB. Considering the impact of SBR and the dark current of Accumulation Charge Coupled Device (ACCD) detector on the accuracy of wind observations, the uncertainty in NA and NB can be calculated as
σA2=NA+NS,A+σRay2,σB2=NB+NS,B+σRay2
where the NS,A and NS,B denote the SBN received by Rayleigh channel A and B, σRay denotes the dark current in the ACCD on Rayleigh channel.

Equations (2)-(6) require the useful signals and SBN received by the two Rayleigh channels to assess the impact of SBR on the accuracy of the wind observations. In sections 2.3 and 2.4, the methods used to achieve this are presented.

2.3. Simulation of spaceborne DWLs

To obtain the useful signals NA and NB, a simulation system of spaceborne DWLs was developed according to the optical structure and instrument parameters of Aeolus, whose optical structure is shown in Fig. 2 and main instrument parameters are shown in Table 1. The instrument parameters in Table 1 mainly obtained from the ADM-Aeolus level 1B algorithm theoretical basis document (ATBD) which was released in 2006 [21]. In the subsequent development process, some parameters may have been improved, so the output of our simulation model may not be completely consistent with the observations of Aeolus. The workflow of our simulation model is described as follows, 1) the laser pulses are sent out via transmitted optics; 2) the signal backscattered from the atmosphere is received by the receiver optics and is directed toward Mie channel through interference filter and polarizing beam splitter. The Fizeau interference of the Mie channel can transmit the central part of the spectrum of the backscattered signals, including the narrow-band backscattered signals from aerosols and part of the wide-band backscattered signals from molecules; 3) the other signal is reflected towards the Rayleigh channel. The Rayleigh channel is composed of two Fabry–Pérot (F-P) interferometers, which acts as spectral filters for the backscattered signals from molecules. The photons of the beam are directed to the F-P A and, subsequently, to the ACCD detection unit. NA photons are detected by the ACCD. The other photons are reflected to the F-P B and the subsequent ACCD detection unit, and NB photons are detected.

 figure: Fig. 2

Fig. 2 Optical structure of spaceborne DWLs in the simulation system.

Download Full Size | PDF

Tables Icon

Table 1. Instrument parameters used for simulations.

The input parameters of the simulation system include the profiles of temperature, pressure, wind vector, and the extinction and backscatter coefficients of aerosols and clouds. The following steps were implemented to derive the useful signals NA and NB: 1) the intensity of the signal SX(z) received by the telescope of the spaceborne DWL was obtained using the lidar equation as is shown in Eq. (7); 2) the spectrum of the signal received by the telescope was simulated according to the Rayleigh spectrum, Mie spectrum, and Doppler shift; 3) the spectrum of transmittance and reflectivity of the Mie channel were simulated, and the spectrum reflected towards the Rayleigh channel was derived; 4) the spectra of transmittance and reflectivity of F-P A, and the spectrum of transmittance of F-P B were simulated. The intensity of signal received by the ACCD on the Rayleigh channel was obtained and converted to photon counts. Through the simulation system, the useful signal in Rayleigh channels A and B can be derived. Further details of the simulations can be found in the ADM-Aeolus level 1B ATBD and the thesis of Praffrath [21–23].

SX(z)=EλhcπAr24βX(z)R(z)2τ(z)2ΔZLOSTRXTTX
where the subscript X denotes the type of scattering, ‘a’ for aerosols and ‘m’ for molecules. E denotes the shoot energy of single laser, λ denotes the wavelength of the laser, Ardenotes the diameter of the telescope, βX(z) denotes the backscatter coefficients of the scattering particles in the measured bin. τ(z) denotes one-way atmospheric transmission, R(z) denotes the range from the measured bins to the satellite. The LOS length of the bin is ΔZLOS=ΔZ/cos(ϕ˜), where ΔZ denotes the geometric depth of the vertical range bins, ϕ˜ denote the view zenith angle of spaceborne DWLs. TRX and TTX denote the transmit and receive optics transmission. h denotes the Plank’s constant, c denotes the velocity of light.

2.4. Conversion from SBR to SBN

According to subsection 2.2, to assess the impact of SBR on the accuracy of wind observations on the Rayleigh channel, NS,A and NS,B, which are caused by SBR should be derived. To obtain NS,A and NS,B, we should firstly convert the SBR to SBN using Eq. (8).

Ns=nEQEOLsφRAr2π4ΔλΔthcλ,
where SBN NS denotes the photon counts which are caused by SBR Ls. n denotes the number of laser shots accumulated for a measurement, EQdenotes the quantum efficiency of ACCD detector, and EO denotes the total optical efficiency of the lidar transmitter. For the receiver system, φR denotes the field of view, Δλ denotes the bandwidth of the interference filter of the receiver, Δt=2ΔZ/c denotes the laser detection time [13]. The vertical range bins used in this paper were obtained from Fig. 5(a) in the paper of Marseille et al [24].

For spaceborne DWLs operating at 355 nm, the SBR originates mainly from the following contributions: 1) radiation originating from the sun, after propagating through the atmosphere, and scattered into the telescope without reaching the Earth; 2) direct radiation transmitted through the atmosphere without being scattered and reflected by the earth; 3) radiation scattered by the atmosphere onto the Earth and reflected back towards the telescope [25]. Accordingly, SBR is determined by the relative position between the sun and the satellite, the atmospheric optics, and the surface reflectance. The relative position between sun and satellite can be derived through the parameters of the orbit and view geometry of spaceborne DWLs. The global distributions of SBR with a spatial resolution of 1° × 1° were also derived based on the global databases of atmospheric molecules, aerosol optical properties, and surface albedo.

Once the SBN is obtained, its spectrum can be assumed to follow a uniform distribution with Δλ bandwidth, and its integral energy is equal to the energy of NS photons. Then, the NS,A and NS,B can be obtained using the simulation system.

3. Databases and processing method

To assess the impact of SBR on the accuracy of wind observations from a global perspective, we derived the distributions of SBR received by spaceborne DWLs in summer and winter. In this section, the databases used to obtain the seasonal global distributions of SBR and their corresponding processing method are present.

3.1. Methodology to derive the global distributions of SBR

The global distributions of SBR under clear skies for spaceborne DWLs in summer and winter at 355 nm was derived using radiative transfer model (RTM) libRadtran [26,27]. The detailed process is illustrated in Fig. 3. First, the longitudes and latitudes of the off-nadir points and their corresponding solar geometries can be derived using satellite orbit simulation software based on the orbit parameters and the view geometry of spaceborne DWLs. Next, the optical properties of aerosols, atmospheric molecules, and surface albedo of the off-nadir points can be obtained through the databases of the lidar climatology of vertical aerosol structure for space-based lidar simulation studies (LIVAS) [28], Global Ozone Monitoring Experiments (GOME) [29], and Ozone Monitoring Instrument (OMI) ozone [30]. Then, the TOA radiance of each off-nadir point can be derived using libRadtran based on its atmospheric and surface optical properties, view geometry, and solar geometry. Finally, we divided the earth into a 1° × 1° grid, each grid would include several off-nadir points, and we picked the maximum TOA radiance as the SBR in one grid.

 figure: Fig. 3

Fig. 3 Process to obtain the global distribution of SBR for spaceborne DWL.

Download Full Size | PDF

3.2. LIVAS database

The global optical properties of aerosols were obtained from the LIVAS database [28]. The database provides averaged profiles of aerosol optical properties in a uniform grid of 1° × 1° with high vertical resolution of CALIOP products, including the extinction coefficient (EXT), backscatter coefficient, and aerosol types. The database divides the aerosol into six types according to CALIOP, namely polluted continental, smoke, dust, polluted dust, clean marine, and clean continental, respectively. In each bin of the vertical profiles, the occurrence frequency of each aerosol type is also provided. To calculate the TOA radiance using libRadtran, the profiles of EXT, single scattering albedo (SSA), and asymmetry parameter (ASY) are required. However, the profiles of SSA and ASY were not provided in LIVAS. We obtained the two parameters in the following steps: firstly, the SSA and ASY of each aerosol types at 355 nm were derived based on the bi-lognormal number-based particle size distribution of each aerosol types [31] and their complex refractive indices. The size distribution parameters and SSA and ASY values of each aerosol type were provided in Table 2. Secondly, the weighted vertical profiles of SSA and ASY were obtained based on the vertical profile of aerosol type occurrence frequency using Eqs. (9) and (10).

ϖ(z)=i=16ηi(z)ϖi,
g(z)=i=16ηi(z)gi,
where ϖ(z) and g(z) denotes the weighted mean SSA and ASY at altitude z, respectively, ϖi and gi denotes the SSA and ASY for the ith aerosol type, respectively. ηi(z)denotes the occurrence frequency of the ith aerosol type.

Tables Icon

Table 2. Size distribution parameters and SSA, ASY coefficient of six aerosol types at 355 nm.

3.3. GOME database

We obtained the global distributions of surface albedo with a spatial resolution of 1° × 1° from the lambert-equivalent reflectivity (LER) of GOME database [29]. The database includes a 1-nm-wide wavelength bins centered at 335 nm, which is used as the surface albedo at 355 nm here. Study shows that the average difference of LER between 335 nm and 355 nm is negligible with a standard deviation within 0.013. The database is made for each month of the year, and the data sets in July and January were defined as the surface albedo in summer and winter, respectively.

3.4. OMI ozone database

As is shown in Fig. 3, the profiles of the atmospheric molecules are needed to derive the SBR for spaceborne DWLs. Considering that the ozone shows absorption bands in the ultraviolet region, the global database of ozone was taken into account in deriving SBR. The database was obtained from the OMI ozone database, which consists of temperature, pressure, and ozone profiles by binning ozonesonde and satellite data over 12 months in 18 latitude (each 10° wide) bins [30]. The molecular density profiles can be derived from temperature and pressure.

4. Global distributions of SBR at 355 nm

To assess the impact of SBR on the accuracy of wind observations for spaceborne DWLs of which the orbit parameters and view geometry were set according to Table 1, the global distributions of SBR at 355 nm under clear skies in summer and winter were derived using the method mentioned in subsection 3.1 and are shown in Fig. 4.

 figure: Fig. 4

Fig. 4 Global distributions of SBR at 355 nm received by the spaceborne DWLs of which the orbit and view geometry were set according to Table 1. The isolines of SBR are related to the percentiles of 50%, 60%, 70%, 80%, 90%, and 100%, respectively. (a) summer, (b) winter.

Download Full Size | PDF

Figure 4 illustrates that: 1) the SBR of most regions in the southern and northern hemispheres are equal to zero in summer and winter, respectively, which indicates that approximately half of the total wind observations is free from SBN; 2) in summer, the maximum SBR, up to 169.03 mW m−2 sr−1 nm−1, can be found in the northern pole. In winter, the maximum SBR, up to 90.98 mW m−2 sr−1 nm−1, can be found in the southern pole; 3) the SBR in the Pole regions is greater than that in other surface types, mainly due to the large surface albedo in the polar regions. According to Koelemeijer et al, the surface albedo on snow/ice land is much larger than other land cover types [29]. The differences of SBR on other regions are also small, probably due to approximately the same amount of surface albedo on other land cover types, which illustrates that the surface albedo has a great impact on the SBR under clear skies.

We also derived the frequency of SBR, as is shown in Fig. 5. The frequency of zero SBR values is nearly 50% (not shown in Figs. 5(a) and 5(b)). As is illustrated in Fig. 5, the SBR in summer is larger than that in winter. The mean SBR is 46.85 mW m−2 sr−1 nm−1 with a standard deviation of 26.09 mW m−2 sr−1 nm−1 in the northern hemisphere in summer, and the mean SBR is 33.65 mW m−2 sr−1 nm−1 with a standard deviation of 21.84 mW m−2 sr−1 nm−1 in the southern hemisphere in winter. The quantile statistics of SBR is presented in Table 3. The SBR is equal to zero when the quantiles are lower than 50%. The isolines in Fig. 4 and Table 3 illustrate that the values SBR range from ~0 to 60 mW m−2 sr−1 nm−1 in greater than 90% regions of Earth. And the SBR is extremely large only in the polar regions.

 figure: Fig. 5

Fig. 5 Frequency of SBR values. The frequency of zero SBR was not shown which is nearly 50%. (a) summer, (b) winter.

Download Full Size | PDF

Tables Icon

Table 3. Quantile SBR statistics from the Aeolus telescope.

5. Impact of SBR on the accuracy of wind observations

The accuracy of wind observations on the Rayleigh channel is mainly affected by atmospheric constituents, parameters of optical structure, and the SBR. In this section, we designed three experiments to assess the impact of SBR on the accuracy of wind speed observations on the Rayleigh channel for spaceborne DWLs.

5.1. Sensitivity of the uncertainty in the wind observations on atmospheric constituents

We assessed the performance of wind observations on six atmospheric climatological constituents, which were obtained from the Air Force Geophysical Laboratory (AFGL) atmospheric constituents profiles, including midlatitude summer, midlatitude winter, subarctic summer, subarctic winter, tropical, and the US76 model [32]. The aerosol model was set to continental clean, which was derived from the optical properties of aerosols and clouds (OPAC) data sets [33]. Two levels of SBR were used, including no SBR and the worst case SBR, which is equal to 169.03 mW m−2 sr−1 nm−1, as is illustrated in Fig. 6. The red bold lines indicates the accuracy expectations for the wind observations of spaceborne DWLs from the ESA. The main findings derived from Fig. 6 are listed as follows:

 figure: Fig. 6

Fig. 6 The profiles of wind observation uncertainty under various atmospheric constitutes. The red solid lines illustrate the accuracy requirement of the ESA for spaceborne DWLs. (a) No SBR; (b) The value of SBR equals to 169.03 mW m−2 sr−1 nm−1.

Download Full Size | PDF

  • (1) The comparisons between the results shown in Fig. 6 and the study of ESA, which was mentioned in Introduction section, were made. The ESA’s results show that: 1) The uncertainty in wind observations is ~2-5 m s−1 at the altitude between 16 and 26 km; ~1-3 m s−1 at the altitudes between 2 and 16 km; and greater than 3 m s−1 at the altitude lower than 2 km, without the impact of SBR, which is consistent with the results shown in Fig. 6(a). 2) When the value of SBR is 154 mW m−2 sr−2 nm−1, the uncertainty in the wind speed observations is ~3-12 m s−1 at the altitudes between 16 km and 26 km; ~2-5 m s−1 at the altitudes between 2 km and 16 km; and greater than 3 m s−1 at altitude lower than 2 km, which is consistent with our results shows in Fig. 6(b), where the SBR equals 169.03 mW m−2 sr−1 nm−1 [6]. The comparisons demonstrate the accuracy of our assessment model.
  • (2) The uncertainties in wind speed observations below 2 km are all higher than 2 m s−1 in the Rayleigh channel. In the operation mode, in most cases the wind observations in the boundary layer are derived from the Mie channel.
  • (3) Without consideration of SBR, the differences of the uncertainties in wind speed observations caused by the differences in atmospheric constituents are small. According to Fig. 6(a), the differences range between 0.02 m s−1 and 0.29 m s−1 in the troposphere and between 0.17 m s−1 and 0.22 m s−1 in the stratosphere.
  • (4) The SBR influences the accuracy of wind speed observations. Without the impact of SBR, the wind speed observations in most bins of the troposphere and stratosphere can satisfy the accuracy requirement for spaceborne DWLs of the ESA. When the SBR approaches 169.03 mW m−2 sr−1 nm−1, the uncertainty in the wind observations becomes large. Considering the worst case SBR, the differences of the uncertainties of wind observations caused by the differences in atmospheric constituents become large, which ranging from 0.06 to 0.93 m s−1 in the troposphere, and from 0.78 to 1.69 m s−1 in the stratosphere, respectively.
  • (5) The comparison between Figs. 6(a) and 6(b) illustrates that the uncertainty in the wind speed observations can be influenced by both atmospheric constituents and SBR and the latter can amplify these differences.

5.2. Sensitivity of the uncertainty in wind observations on SBR

In this subsection, the sensitivity of the uncertainty in the wind observations on SBR was tested. In the experiment, the atmospheric constituent was set to the US76 model, and the aerosol model was set as continental clean. The values of SBR ranged from 0 to 100 mW m−2 sr−1 nm−1, with 20 mW m−2 sr−1 nm−1 steps. The profiles of the uncertainty in wind observations under different levels of SBR are shown in Fig. 7. According to the global distributions of SBR illustrated in Fig. 4, the corresponding quantiles of the SBR are labelled in the legend of Fig. 7. The mean uncertainty in the wind observations on the Rayleigh channel in the troposphere and the stratosphere under different levels of SBR is given in Table 4. We can draw the following conclusions from Fig. 7 and Table 4:

 figure: Fig. 7

Fig. 7 Sensitivity of the uncertainty in wind observations on SBR. The label on the brackets denotes the corresponding quantiles of the global distribution of SBR on summer shown in Fig. 4(a).

Download Full Size | PDF

Tables Icon

Table 4. Mean of uncertainty in wind observations on Rayleigh channel under different levels of SBR.

  • (1) In the troposphere, taking the uncertainty in the wind observations on the Rayleigh channel as an example, with a SBR between 0 and 20 mW m−2 sr−1 nm−1, the mean difference of the uncertainties is 0.22 m s−1 and, while with a SBR between 20 and 40 mW m−2 sr−1 nm−1, it is 0.19 m s−1. With the SBR increasing from 40, 60, 80, to 100 mW m−2 sr−1 nm−1, the uncertainty increases by 0.18, 0.17, and 0.16 m s−1, respectively. Therefore, the mean uncertainty in the wind observations increases by 0.18 and 0.69 m s−1 when 20 mW m−2 sr−1 nm−1 SBR increases in the troposphere and the stratosphere, respectively. The increase of the uncertainty in the wind observations is large compared to the accuracy requirement of the ESA.
  • (2) In the troposphere, the uncertainty in wind observations on the Rayleigh channel at each 20 mW m−2 sr−1 nm−1 step increases by 0.22, 0.20, 0.18, 0.17, and 0.16 m s−1 when the SBR increases from 0 to 100 mW m−2 sr−1 nm−1. The result shows that the increase of the uncertainty becomes smaller as the SBR increases, which indicates that the uncertainty in wind observations on the Rayleigh channel is more sensitive to changes in SBR at lower SBR.
  • (3) As indicated by some bins in the vertical profile, the uncertainty in wind observations is beyond the accuracy requirements of the ESA even no SBR is considered, especially between 22.25 km and 26.25 km in the stratosphere and between 11.25 km and 16.50 km in the troposphere. The large uncertainty in the former layer are mainly caused by the relative small density of molecules which lead to the weak backscatter signals. And the large uncertainty in the latter layer are due to that the depth of vertical bins is changed from 2 km to 1 km at the altitude of 16.5 km. As is shown in Eq. (8), lower ΔZLOS will lead to lower intensity of backscatter signals, which also illustrated that the appropriate configurations of vertical sampling scenarios can improve the accuracy of spaceborne DWL measurements [34].

From Fig. 7 and Table 4, we can conclude that the uncertainty in the wind observations for spaceborne DWLs is sensitive to the level of SBR. Therefore, the impact of SBR cannot be neglected in the design of the spaceborne DWLs.

5.3. Wind observations uncertainty profiles under various amounts of SBR

In section 4, we obtained the global distributions of SBR and derived its quantiles statistics. In this subsection, the wind observations uncertainty profiles under different amounts of SBR are calculated, setting the atmospheric constituent and aerosol model were set as in subsection 5.2. Because the value of SBR is equal to zero when the percentiles are lower than 50%, the uncertainty in wind speed observations corresponding to the percentiles of SBR ranging from 50% to 100% with 5 percentile increments were derived and are illustrated in Fig. 8. The following findings can be obtained:

 figure: Fig. 8

Fig. 8 Profiles uncertainty in wind observations on the Rayleigh channel under different quantiles of SBR ranges from 50% to 100%, here the 100% percentile of SBR means the maximum value of SBR, and it’s related values of SBR are labelled in the top of the figures. The white contour illustrates the accuracy requirement for spaceborne DWLs of the ESA. (a) summer, (b)winter.

Download Full Size | PDF

  • (1) In general, the uncertainty caused by SBR in the troposphere is lower than that in the stratosphere in a specific profile of wind observations. That is mainly due to the larger molecular density in the troposphere than stratosphere, which leads to a larger intensity of Rayleigh backscatter signal.
  • (2) Compared with accuracy requirement of the ESA, the uncertainty of wind observations is less accurate when the percentile of SBR reaches a threshold value. The threshold of each bin can be estimated through the isolines illustrated in Fig. 8. In summer, the uncertainty of the whole profiles of wind observations in the troposphere and stratosphere are greater than 2 m s−1 and 3 m s−1 when the percentiles of SBR are greater than 65% and 85%, respectively, under the atmospheric and surface scenario used in this study. In winter, less accurate wind observations are made in the troposphere and stratosphere when the percentiles of SBR are greater than 75% and 95%, respectively. According to the isolines in Fig. 4, the less accurate wind observations mainly occur in regions close to the northern pole and the southern pole where the related percentiles of SBR is greater than 90%.
  • (3) Comparing the Figs. 8(a) and 8(b), the impact of SBR on the uncertainty in the wind speed observations for the spaceborne DWLs in winter is less than that in summer, which is due to the SBR in winter is smaller than that in summer.

6. Summary and conclusions

The SBR is the major source of noise on wind observations of spaceborne DWLs during daytime. To eliminate its negative effects, engineers detect and subtract the SBN from the total signal. However, the SBR still affects the accuracy of the observations due to Poisson noise. In this study, we assessed the impact of SBR on the accuracy of wind observations for spaceborne DWLs simulated based on the orbit parameters and optical structure of Aeolus. According to previous studies, the SBR has little influence on the observations in the Mie channel for Aeolus; hence, only the observations in the Rayleigh channel were considered in this study.

To assess the impact of SBR on the accuracy of wind observations from a global perspective, a simulation system of spaceborne DWLs was developed and the SBR was converted into SBN according to the instrument parameters. The global distributions of SBR at 355 nm in summer and winter were derived using the LIVAS, GOME, and OMI ozone databases. According to the distributions, we assessed the uncertainty in the wind speed observations under various amounts of SBR. The main findings are listed as following:

  • (1) Tests show that the sensitivity of the uncertainty in the wind observations on atmospheric constituents is small, which indicates that the uncertainty in the wind observation in the Rayleigh channel is less sensitive to atmospheric constituents.
  • (2) The mean increasement of uncertainty in wind observations is approximately 0.18 m s−1 in the troposphere and 0.69 m s−1 in the stratosphere with an increasement of 20 mW m−2 sr−1 nm−1 SBR value. The result indicates that the impact of SBR cannot be negligible for the measurements of spaceborne DWLs.
  • (3) The global distributions of SBR at 355 nm in summer and winter illustrate that the off-nadir points in one hemisphere are always free from SBR and the signal received from the off-nadir points in the other hemisphere are contaminated by SBN, which is attributed to the sun-synchronous, dawn-dusk orbit. The worst cases occur the in polar regions, which are mainly caused by the large surface albedo in the polar regions. Furthermore, the values of SBR in summer are larger than that in winter
  • (4) The profiles of the uncertainty in the wind observations under various amounts of SBR were derived based on global distributions. By comparing each bin of a specific profile, it was found that the uncertainty in the wind observations in the troposphere is lower than that in the stratosphere and the wind observations in some bins of a profile are less accurate compared with the accuracy requirement of the ESA, even when the SBN is not considered. The comparisons among the uncertainty of the bins at same height under different percentiles of SBR further illustrate that: 1) the uncertainty in the wind speed observations increases with increasing percentiles of SBR; 2) the uncertainty does not meet the accuracy requirement of the ESA when the percentiles of SBR are larger than a threshold value, which is different for different bins. Moreover, the impact of SBR on the wind observations is lower in winter than in summer.

In summary, the SBR has a considerable impact on the accuracy of spaceborne DWLs, especially when the values of SBR are larger than a threshold value, which indicates that the impact of SBR cannot be negligible while designing spaceborne DWLs mission. When the spaceborne DWLs operate on the sun-synchronous, dawn-dusk orbit, nearly the observations in one hemisphere are always free of SBR. For the other hemisphere, only the values of SBR in the regions close to the Pole are large, where the SBR will have significant impact on the uncertainty in wind speed observations compared with the accuracy requirement of the ESA. In the future, several spaceborne DWLs will be launched in different orbits for different purposes. For example, wind has clear diurnal variation characteristics and spaceborne DWLs that operates on orbits with different local times of ascending node may be needed. Furthermore, the corresponding global distributions of SBR would increase compared with the sun-synchronous, dawn-dusk orbit. We can assess the impact of SBR on different orbit parameters using the method mentioned in this paper, and improve the instrument parameters to compensate for the increased uncertainty in the wind speed observations caused by the increase of SBR.

Appendix

The descriptions and units of the variables used in the equations of this manuscript were listed in Table 5.

Tables Icon

Table 5. The descriptions and units of the variables used in this paper.

Funding

National Natural Science Foundation of China (NSFC) (41575020).

Acknowledgments

Thanks for the helpful discussions provided by Dr. Karsten Schmidt from DLR, Dr. Gert-Jan Marseille and Dr. Ad Stoffelen from Royal Netherlands Meteorological Institute. Thanks for the suggestions provided by Dr. Claudia Emde in running the libRadtran. The LIVAS database was provided by Dr. V. Amiridis from Institude for space applications and remote sensing, National observatory of Athens. The GOME databases was provided by R. B. A. Koelemeijer from Royal Netherlands Meteorological Institute. We also thank the anonymous reviewers for the great help in the improvement of this paper.

References

1. A. Stoffelen, J. Pailleux, E. Kallen, J. M. Vaughan, L. Isaksen, P. Flamant, W. Wergen, E. Andersson, H. Schyberg, A. Culoma, R. Meynart, M. Endemann, and P. Ingmann, “The atmospheric dynamics mission for global wind field measurement,” Bull. Am. Meteorol. Soc. 86(1), 73–88 (2005). [CrossRef]  

2. M. Weissmann and C. Cardinali, “Impact of airborne Doppler lidar observations on ECMWF forecasts,” Q. J. R. Meteorol. Soc. 133, 107–116 (2007). [CrossRef]  

3. N. Zagar, N. Gustafsson, and E. Kallen, “Dynamical response of equatorial waves in four-dimensional variational data assimilation,” Tellus, Ser. A, Dyn. Meterol. Oceanogr. 56(1), 29–46 (2004). [CrossRef]  

4. A. Cress and W. Wergen, “Impact of profile observations on the German Weather Service’s NWP system,” Meteorol. Z. (Berl.) 10(2), 91–101 (2001). [CrossRef]  

5. M. J. McGill, W. D. Hart, J. A. McKay, and J. D. Spinhirne, “Modeling the performance of direct-detection Doppler lidar systems including cloud and solar background variability,” Appl. Opt. 38(30), 6388–6397 (1999). [CrossRef]   [PubMed]  

6. M. Rennie, CCN6 Results: Further Chain-of-Processors Testing of L2B Results and Testing of CCN6 L2B Processor Algorithm Updates (European Centre for Medium-Range Weather Forecasts, 2017).

7. A. Heliere, J. L. Bezy, P. Bensi, and P. Ingmann, “System definition of the ESA Earth Explorer WALES mission,” Proc. SPIE 4881, 24–32 (2003).

8. R. R. Rogers, C. A. Hostetler, J. W. Hair, R. A. Ferrare, Z. Liu, M. D. Obland, D. B. Harper, A. L. Cook, K. A. Powell, M. A. Vaughan, and D. M. Winker, “Assessment of the CALIPSO Lidar 532 nm attenuated backscatter calibration using the NASA LaRC airborne High Spectral Resolution Lidar,” Atmos. Chem. Phys. 11(3), 1295–1311 (2011). [CrossRef]  

9. C. Kiemle, M. Quatrevalet, G. Ehret, A. Amediek, A. Fix, and M. Wirth, “Sensitivity studies for a space-based methane lidar mission,” Atmos. Meas. Tech. 4(10), 2195–2211 (2011). [CrossRef]  

10. C. Kiemle, S. R. Kawa, M. Quatrevalet, and E. V. Browell, “Performance Simulations for a Spaceborne Methane Lidar Mission,” J. Geophys. Res. Atmos. 119(7), 4365–4379 (2014). [CrossRef]  

11. S. R. Kawa, J. Mao, J. B. Abshire, G. J. Collatz, X. Sun, and C. J. Weaver, “Simulation studies for a space-based CO2 lidar mission,” Tellus B Chem. Phys. Meterol. 62, 759–769 (2010). [CrossRef]  

12. Z. Liu, W. Hunt, M. Vaughan, C. Hostetler, M. McGill, K. Powell, D. Winker, and Y. Hu, “Estimating random errors due to shot noise in backscatter lidar observations,” Appl. Opt. 45(18), 4437–4447 (2006). [CrossRef]   [PubMed]  

13. T. Y. Nakajima, T. Imai, O. Uchino, and T. Nagai, “Influence of daylight and noise current on cloud and aerosol observations by spaceborne elastic scattering lidar,” Appl. Opt. 38(24), 5218–5228 (1999). [CrossRef]   [PubMed]  

14. W. H. Hunt, D. M. Winker, M. A. Vaughan, K. A. Powell, P. L. Lucker, and C. Weimer, “CALIPSO Lidar Description and Performance Assessment,” J. Atmos. Ocean. Technol. 26(7), 1214–1228 (2009). [CrossRef]  

15. C. Zhang, X. Sun, R. Zhang, and Y. Liu, “Simulation and assessment of solar background noise for spaceborne lidar,” Appl. Opt. 57(31), 9471–9479 (2018). [CrossRef]   [PubMed]  

16. C. Flesia and C. L. Korb, “Theory of the double-edge molecular technique for Doppler lidar wind measurement,” Appl. Opt. 38(3), 432–440 (1999). [CrossRef]   [PubMed]  

17. J. A. McKay, “Modeling of direct detection Doppler wind lidar. I. The edge technique,” Appl. Opt. 37(27), 6480–6486 (1998). [CrossRef]   [PubMed]  

18. G. J. Marseille and A. Stoffelen, “Simulation of wind profiles from a space-borne Doppler wind lidar,” Q. J. R. Meteorol. Soc. 129, 3079–3098 (2003). [CrossRef]  

19. D. G. H. Tan, E. Anderson, J. De Kloe, G. Marseille, A. Stoffelen, P. Poli, M. Denneulin, A. Dabas, D. Huber, O. Reitebuch, P. Flamant, O. Le Rille, and H. Nett, “The ADM-Aeolus wind retrieval algorithms,” Tellus, Ser. A, Dyn. Meterol. Oceanogr. 60(2), 191–205 (2008). [CrossRef]  

20. P. Flamant, J. Cuesta, M. Denneulin, A. Dabas, and D. Huber, “ADM-Aeolus retrieval algorithms for aerosol and cloud products,” Tellus, Ser. A, Dyn. Meterol. Oceanogr. 60(2), 273–286 (2008). [CrossRef]  

21. O. Reitebuch, U. Paffrath, and I. Leike, ATBD: ADM-Aeolus Level 1B Product (European Space Research and Technology Centre, 2006).

22. U. Paffrath, Performance Assessment of the Aeolus Doppler Wind Lidar Prototype (Ludwig-Maximilians-Universität München, 2006).

23. R. W. Zhang, X. J. Sun, W. Yan, J. Zhao, L. Liu, Y. Li, C. L. Zhang, and J. H. Zhou, “Simulation of frequency discrimination for spaceborne Doppler wind lidar (II): Study on the retrieval of atmospheric wind speed for Rayleigh channel based on Fabry-Perot interferometer,” Wuli Xuebao 63(14), 147–156 (2014).

24. G. J. Marseille, A. Stoffelen, J. De Kloe, K. Houchi, H. Kornich, H. Schyberg, A. Straume, and O. Lerille, “ADM-Aeolus vertical sampling scenarios,” in Proceedings of the 8th International Sysposium on Tropospheric Profiling (2009).

25. A. Mousivand, W. Verhoef, M. Menenti, and B. Gorte, “Modeling Top of Atmosphere Radiance over Heterogeneous Non-Lambertian Rugged Terrain,” Remote Sens. 7(6), 8019–8044 (2015). [CrossRef]  

26. B. Mayer and A. Kylling, “Technical note: The libRadtran software package for radiative transfer calculations - description and examples of use,” Atmos. Chem. Phys. 5(7), 1855–1877 (2005). [CrossRef]  

27. C. Emde, R. Burasschnell, A. Kylling, B. Mayer, J. Gasteiger, U. Hamann, J. Kylling, B. Richter, C. Pause, T. Dowling, and L. Bugliaro, “The libRadtran software package for radiative transfer calculations (Version 2.0),” Geosci. Model Dev. 8(12), 10237–10303 (2015). [CrossRef]  

28. V. Amiridis, E. Marinou, A. Tsekeri, U. Wandinger, A. Schwarz, E. Giannakaki, R. Mamouri, P. Kokkalis, I. Binietoglou, S. Solomos, T. Herekakis, S. Kazadzis, E. Gerasopoulos, E. Proestakis, M. Kottas, D. Balis, A. Papayannis, C. Kontoes, K. Kourtidis, N. Papagiannopoulos, L. Mona, G. Pappalardo, O. Le Rille, and A. Ansmann, “LIVAS: a 3-D multi-wavelength aerosol/cloud database based on CALIPSO and EARLINET,” Atmos. Chem. Phys. 15(13), 7127–7153 (2015). [CrossRef]  

29. R. Koelemeijer, J. F. de Haan, and P. Stammes, “A database of spectral surface reflectivity in the range 335-772 nm derived from 5.5 years of GOME observations,” J. Geophys. Res. 108, 171–181 (2003). [CrossRef]  

30. R. McPeters, M. Kroon, G. Labow, E. Brinksma, D. Balis, I. Petropavlovskikh, J. P. Veefkind, P. K. Bhartia, and P. F. Levelt, “Validation of the Aura Ozone Monitoring Instrument total column ozone product,” J. Geophys. Res. Atmos. 113, D15S14 (2008). [CrossRef]  

31. Y. B. Zhou, X. J. Sun, C. L. Zhang, R. W. Zhang, Y. Li, and H. R. Li, “3D aerosol climatology over East Asia derived from CALIOP observations,” Atmos. Environ. 152, 503–518 (2017). [CrossRef]  

32. G. P. Anderson, S. A. Clough, F. X. Kneizys, J. H. Chetwynd, and E. P. Shettle, AFGL (Air Force Geophysical Laboratory) atmospheric constituent profiles (0. 120km). Environmental research papers (Air Force Geophysical Laboratory, 1986).

33. M. Hess, P. Koepke, and I. Schult, “Optical properties of aerosols and clouds: The software package OPAC,” Bull. Am. Meteorol. Soc. 79(5), 831–844 (1998). [CrossRef]  

34. J. De Kloe, G. J. Marseille, and A. Stoffelen, ADM-Aeolus, Vertical Aeolus Measurement Positioning (The Royal Netherlands Meteorological Institute, 2007).

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (8)

Fig. 1
Fig. 1 Lookup table of the partial derivative of HLOS wind component to Rayleigh response, which is determined by temperature and pressure.
Fig. 2
Fig. 2 Optical structure of spaceborne DWLs in the simulation system.
Fig. 3
Fig. 3 Process to obtain the global distribution of SBR for spaceborne DWL.
Fig. 4
Fig. 4 Global distributions of SBR at 355 nm received by the spaceborne DWLs of which the orbit and view geometry were set according to Table 1. The isolines of SBR are related to the percentiles of 50%, 60%, 70%, 80%, 90%, and 100%, respectively. (a) summer, (b) winter.
Fig. 5
Fig. 5 Frequency of SBR values. The frequency of zero SBR was not shown which is nearly 50%. (a) summer, (b) winter.
Fig. 6
Fig. 6 The profiles of wind observation uncertainty under various atmospheric constitutes. The red solid lines illustrate the accuracy requirement of the ESA for spaceborne DWLs. (a) No SBR; (b) The value of SBR equals to 169.03 mW m−2 sr−1 nm−1.
Fig. 7
Fig. 7 Sensitivity of the uncertainty in wind observations on SBR. The label on the brackets denotes the corresponding quantiles of the global distribution of SBR on summer shown in Fig. 4(a).
Fig. 8
Fig. 8 Profiles uncertainty in wind observations on the Rayleigh channel under different quantiles of SBR ranges from 50% to 100%, here the 100% percentile of SBR means the maximum value of SBR, and it’s related values of SBR are labelled in the top of the figures. The white contour illustrates the accuracy requirement for spaceborne DWLs of the ESA. (a) summer, (b)winter.

Tables (5)

Tables Icon

Table 1 Instrument parameters used for simulations.

Tables Icon

Table 2 Size distribution parameters and SSA, ASY coefficient of six aerosol types at 355 nm.

Tables Icon

Table 3 Quantile SBR statistics from the Aeolus telescope.

Tables Icon

Table 4 Mean of uncertainty in wind observations on Rayleigh channel under different levels of SBR.

Tables Icon

Table 5 The descriptions and units of the variables used in this paper.

Equations (10)

Equations on this page are rendered with MathJax. Learn more.

Δ n p 2 ¯ = n ¯ p
σ HLOS = ( v HLOS R ATM σ R ATM ) 2 + ( v HLOS T σ T ) 2 + ( v HLOS P σ P ) 2 + ( v HLOS ρ σ ρ ) 2
R ATM = N A N B N A + N B
σ HLOS = v HLOS R ATM σ R ATM
σ R ATM = 2 N A + N B N B 2 σ A 2 + N A 2 σ B 2 ,
σ A 2 = N A + N S,A + σ Ray 2 , σ B 2 = N B + N S,B + σ Ray 2
S X (z)= Eλ hc π A r 2 4 β X (z) R (z) 2 τ (z) 2 Δ Z LOS T RX T TX
N s =n E Q E O L s φ R A r 2 π 4 ΔλΔt h c λ ,
ϖ(z)= i=1 6 η i (z) ϖ i ,
g(z)= i=1 6 η i (z) g i ,
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.