Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Four-Dimensional Parameter Estimation for Mixed Far-Field and Near-Field Target Localization Using Bistatic MIMO Arrays and Higher-Order Singular Value Decomposition
Previous Article in Journal
The Reconstruction of FY-4A and FY-4B Cloudless Top-of-Atmosphere Radiation and Full-Coverage Particulate Matter Products Reveals the Influence of Meteorological Factors in Pollution Events
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On-Orbit Wavelength Calibration Error Analysis of the Spaceborne Hyperspectral Greenhouse Gas Monitoring Instrument Using the Solar Fraunhofer Lines

by
Yulong Guo
1,2,
Cailan Gong
1,*,
Yong Hu
1,
Fuqiang Zheng
1 and
Yunmeng Liu
1
1
Key Laboratory of Infrared System Detection and Imaging Technologies, Shanghai Institute of Technical Physics, Chinese Academy of Sciences, Shanghai 200083, China
2
University of Chinese Academy of Sciences, Beijing 100049, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2024, 16(18), 3367; https://doi.org/10.3390/rs16183367
Submission received: 8 June 2024 / Revised: 6 September 2024 / Accepted: 8 September 2024 / Published: 10 September 2024

Abstract

:
Accurate on-orbit wavelength calibration of the spaceborne hyperspectral payload is the key to the quantitative analysis and application of observational data. Due to the high spectral resolution of general spaceborne hyperspectral greenhouse gas (GHG) detection instruments, the common Fraunhofer lines in the solar atmosphere can be used as a reference for on-orbit wavelength calibration. Based on the performances of a GHG detection instrument under development, this study simulated the instrument’s solar-viewing measurement spectra and analyzed the main sources of errors in the on-orbit wavelength calibration method of the instrument using the solar Fraunhofer lines, including the Doppler shift correction error, the instrumental measurement error, and the peak-seek algorithm error. The calibration accuracy was independently calculated for 65 Fraunhofer lines within the spectral range of the instrument. The results show that the wavelength calibration accuracy is mainly affected by the asymmetry of the Fraunhofer lines and the random error associated with instrument measurement, and it can cause calibration errors of more than 1/10 of the spectral resolution at maximum. A total of 49 Fraunhofer lines that meet the requirements for calibration accuracy were screened based on the design parameters of the instrument. Due to the uncertainty of simulation, the results in this study have inherent limitations, but provide valuable insights for quantitatively analyzing the errors of the on-orbit wavelength calibration method using the Fraunhofer lines, evaluating the influence of instrumental parameters on the calibration accuracy, and enhancing the accuracy of on-orbit wavelength calibration for similar GHG detection payloads.

Graphical Abstract

1. Introduction

In recent years, the problem of global climate change caused by excessive emissions of CO2 and other greenhouse gases has become more and more serious and has developed into a global issue of concern to all countries around the world [1,2,3,4]. Detecting the distribution of atmospheric CO2 concentrations is important for understanding global greenhouse gas emissions and uptake and improving our understanding of climate change. Spaceborne hyperspectral GHG detection instruments can acquire fine spectral information on surface targets or atmospheric molecules and detect the narrow-band absorption lines of gas molecules in the atmosphere. The GHG concentrations can be inverted based on the depth of the gas absorption spectra. Hyperspectral remote sensing has unrivaled advantages in the inversion of atmospheric composition [5,6,7]. In 2019, the IPCC explicitly proposed expanding the observation network by adding satellite remote sensing observations and supporting and validating the emission inventories through “top-down” flux calculations centered on assimilation models [8]. SCIAMACHY on ESA’s Envisat is the world’s first atmospheric tropospheric and stratospheric GHG space probe payload, which was launched in March 2002 as a joint project of Germany, the Netherlands, and Belgium [9]. The main satellites dedicated to GHG detection currently in orbit include NASA’s OCO-2, JAXA’s GOSAT series, GHGSat constellation of GHGSat Inc. (Montreal, QC, Canada), and TanSat supported by the Ministry of Science and Technology of China (MOST), the Chinese Academy of Sciences (CAS), and the China Meteorological Administration (CMA) [10,11,12,13,14]. In addition, several Integrated Earth Observation Satellites also carry GHG detection payloads, such as ESA’s Sentinel-5P and China’s FY-3D and GF-5 satellites [15,16,17].
For atmospheric hyperspectral remote sensing payloads, the radiometric and spectral calibration accuracy of the measured data is an important factor affecting the inversion of atmospheric composition [18]. The spectral calibration accuracy mainly includes the spectral resolution calibration accuracy and the wavelength calibration accuracy, which not only directly affects the positioning accuracy of atmospheric absorption spectral lines but also impacts the radiometric calibration and correction, thus affecting the accuracy of measured radiance. The Jet Propulsion Laboratory (LaKan Yada and Pasadena, CA, USA) has reported that spectral calibration errors of 5% cause significant errors in measured radiance in the solar reflectance band [19]. For hyperspectral GHG detection instruments, deviations in the spectral position, particularly at the absorption peaks and valleys, can result in decreased detection precision. Relevant studies have shown that the degree of spectral calibration center wavelength shift significantly affects the precision of spectral radiance retrieval conducted by the instrument, and a minor 1% shift in the center wavelength results in an average relative error of 0.1% in spectral radiance [20]. Therefore, minimizing center wavelength shifts can effectively reduce the error of the acquired data and enhance the retrieval accuracy of the GHG concentration. The main hyperspectral GHG detection satellites in orbit, such as OCO-2 and TanSat, require on-orbit wavelength calibration errors of <10% spectral resolution or even less to measure the atmospheric column-averaged dry air mole fractions of carbon dioxide (XCO2) with a precision of <1% (~4 ppm) [13,21,22].
To improve the reliability of quantitative spectral data measured by the hyperspectral GHG detection instruments, we need to perform accurate spectral calibrations in the laboratory to derive the relevant calibration coefficients. However, due to vibration, changes in environmental factors, the aging of instrument components, and other factors, changes in detector center wavelength position and spectral resolution will occur after the instrument is operated in orbit with the satellite, so it is necessary to carry out on-orbit spectral calibration of the instruments to calibrate the results of laboratory calibrations [23,24]. On-orbit spectral calibration of hyperspectral remote sensing loads is mainly carried out using characteristic spectral lines with known wavelengths, which are generally generated by the internal calibration light source of the payload or by an external light source. At present, one way to perform on-orbit spectral calibration is to assemble standard lamps with standard emission lines and determine the wavelength shift of the instrument according to the change in the emission lines’ positions. Another method is to add a standard filter with a known absorption line to the optical path and measure the change in absorption line position to achieve calibration. Both of the above methods have certain limitations. The first method requires power consumption to ignite the lamps, and in order to realize the calibration of different bands of the instrument, it is necessary to use a variety of standard lamps covering different wavelength ranges and a rotating mechanism to convert these lamps. The second method also requires a rotating mechanism to move filters in and out of the optical path [25,26,27,28,29]. Both methods increase the demand for satellite resources and require high precision of the rotating mechanism, which reduces the reliability of the calibration.
Spaceborne hyperspectral GHG detection instruments have a high spectral resolution and are capable of resolving the Fraunhofer lines in the solar spectrum, which are caused by the presence of elements in the Sun. The Fraunhofer lines are fixed in the solar spectrum and are not affected by absorption in the Earth’s atmosphere. They are sharper and the wavelength positions of the absorption peaks are more precise than the Earth observation spectra. This method utilizes the Sun’s intrinsic characteristic spectral lines as the standard, which reduces the influence of the internal mechanical structure of the payload and reduces the demand on satellite resources. Therefore, in recent years, hyperspectral payloads have often utilized solar Fraunhofer lines for on-orbit wavelength calibration [30,31]. The Medium Resolution Imaging Spectrometer (MERIS) on Envisat was configured both for Earth and diffuser observations and acquired data for Fraunhofer spectral calibration [32]. A solar diffuser was mounted on OCO-2, which is used to acquire routine observations of the Sun. The spectroscopic performance has been routinely monitored in orbit using measurements of the positions and shapes [21]. The variations in the center wavelength of the atmospheric carbon dioxide grating spectrometer (ACGS) carried on TanSat were calculated at the locations of the Fraunhofer lines by comparing solar-viewing measurements and a high-resolution solar reference spectrum [22].
The solar Fraunhofer lines vary greatly in strength and distribution, especially in the visible and near-infrared bands [33,34]. At the same time, the Fraunhofer lines observed by the instrument are limited by design parameters such as the spectral resolution and sampling rate. When the spectral resolution and sampling rate are high, the instruments can accurately resolve the stronger Fraunhofer lines, and on-orbit wavelength calibration based on these lines is feasible; however, when the spectral resolution and sampling rate are limited, the instruments’ observations of Fraunhofer lines may be affected by the asymmetry of the lines themselves and the adjacent absorption lines, which in turn affects the accuracy of wavelength calibration. In addition, due to noise signals and uncertainties in the radiometric calibration, some Fraunhofer lines may not be accurately resolved or the peak-seek accuracy may be reduced, which needs to be discriminated and screened. Finally, the Doppler shift correction of measured spectra and the peak-seek algorithm can also cause errors in the process of wavelength calibration. Therefore, it is important to analyze the accuracy of the on-orbit wavelength calibration method using the solar Fraunhofer lines with respect to the design parameters of the hyperspectral GHG detection instrument and to screen the solar Fraunhofer lines in the spectral range of the instrument to meet the calibration accuracy requirements.
In this study, the design parameters and performance of a GHG detection instrument under development are used as references. A generic spaceborne hyperspectral GHG monitoring instrument (GHG Monitor) is defined, whose specific design parameters are shown in Section 2.1, and the instrument’s solar observation spectrum is simulated. Seven typical Fraunhofer lines are selected as examples for each of the three bands. The main sources of errors in the wavelength calibration method based on the solar Fraunhofer lines, including the Doppler shift correction error, the instrument measurement systematic error, random error, and the peak-seek algorithm error, are highlighted through the peak-seek algorithm by fitting the simulated measured Fraunhofer lines. The effects of instrumental measurement errors due to spectral resolution, SNR, and other design parameters on the calibration accuracy are also analyzed. Finally, the results of the analysis are used to distinguish the 65 Fraunhofer lines within the design spectral range of GHG Monitor, which offers a basis and reference for the selection of Fraunhofer lines used for on-orbit high-precision wavelength calibration of general hyperspectral GHG detection instruments. This study analyzes the effect of instrument performance on the accuracy of on-orbit wavelength calibration, which is informative for the design of general GHG detection instruments.

2. Materials and Methods

2.1. Instrument Design Parameters

GHG Monitor is a three-channel grating spectrometer that employs one near-infrared and two short-wave infrared spectral bands for the continuous monitoring of atmospheric CO2 concentration by measuring the atmospheric fine spectrum in the O2 and CO2 absorption wavebands. GHG Monitor utilizes a grating spectroscopy method that optimizes grating parameters to achieve high diffraction efficiency in narrow spectral bands, and it observes the solar spectrum through a solar diffuser to enable absolute radiometric and wavelength calibrations of the instrument in orbit.
GHG Monitor’s band settings center at the following wavelengths: 0.76 µm, O2 absorption band (O2A); 1.61 µm, weak CO2 absorption band (WCO2); 2.06 µm, strong CO2 absorption band (SCO2). The WCO2 band is sensitive to the change in CO2 concentration in the atmospheric boundary layer, and the absorption degree of the spectral line is approximately linear with the CO2 concentration, which is an ideal spectral band for the inversion of near-surface high-precision CO2 concentration. The SCO2 band is more sensitive to suspended particles in the atmosphere and can detect clouds, scattering, and changes in atmospheric pressure and humidity caused by suspended particles [18]. The O2A band can provide accurate surface pressure information, and when combined with the SCO2 band, the scattering problems caused by clouds and aerosols can be well mitigated, which can effectively improve the inversion accuracy of atmospheric CO2 concentration [35]. The spectral resolution represents the full width at half maxima (FWHM) of the instrument line shape (ILS). As the spectral resolution of the instrument increases, the sharper the absorption lines in the measured spectrum become, and the higher the sensitivity to atmospheric CO2 concentration gets. However, as the spectral resolution increases, the instrument detects lower radiant energy and puts forward very high requirements for SNR [36,37]. To meet the requirements of the retrieval accuracy of CO2 concentration, it is necessary to consider the indicators of the instrument comprehensively. The main design parameters of the GHG Monitor are shown in Table 1.
For a hyperspectral remote sensor with N channels, the output of channel i is shown in the following equation:
X i = +   E λ d 2 π β 2 A ρ ( λ ) τ ( λ ) τ o ( λ ) S s ( λ ) R d ( λ ) R e ( λ ) d λ
where E λ is the solar irradiance at wavelength λ ; β is the instantaneous field of view of the system; A is the effective optical area of the system; ρ ( λ ) is the reflectivity of the ground; τ ( λ ) is the spectral transmittance of the atmosphere; τ o ( λ ) is the optical efficiency of the system; S s ( λ ) is the transfer function of the dispersion system of the sensors; R d ( λ ) is the spectral responsivity of the detector; and R e ( λ ) is the spectral responsivity of the electronics system. The spectral response function of the i th channel of the system is expressed as follows:
S i = τ o ( λ ) S s ( λ ) R d ( λ ) R e ( λ )
It can usually be characterized by a super-Gaussian function, i.e., the relative spectral response R λ of each spectral channel can be expressed by Equation (3), and Figure 1 shows super-Gaussian functions of different orders n:
R λ = a · e λ λ c c n
where a is the fitting coefficient, λ c is the sampling center wavelength of each spectral channel, n is the super-Gaussian order, and the spectral resolution is expressed by the following equation:
F W H M = 2 c · ( ln ( 2 ) ) 1 n

2.2. Solar Reference Spectrum

In this study, the high-resolution Kurucz solar spectrum was used as the input for the calculation of the simulation spectra of the instrumental solar calibration observations, and it was also used as the reference spectrum for the wavelength calibration method based on the Fraunhofer lines. Its sampling resolution in the O2A band, WCO2 band, and SCO2 band of the GHG Monitor was 0.001 nm, 0.001 nm, and 0.0025 nm, respectively, and the spectral resolution was one order of magnitude higher than that of the GHG Monitor, meeting the requirements regarding the on-orbit wavelength calibration accuracy of the instrument [38,39].
For the calibration method based on the theoretical solar eigen-spectral line, it was first necessary to select the appropriately located Fraunhofer line as the calibration standard. Table 2 and Figure 2 provide a total of 21 typical reference Fraunhofer lines in three bands, and the exact positions of the absorption peaks were calculated using the extremum algorithm, which regards the positions of the absorption peak minimum as the exact positions of Fraunhofer lines. The Kurucz solar reference spectrum in the SCO2 band was interpolated to 0.001 nm using the seven-point Lagrange interpolation method.

2.3. Simulation of Measurement Spectrum

The simulated instrumental measurement spectrum I s i m u λ was obtained by convolving the high-resolution Kurucz solar reference spectrum I r e f λ over the full band at the specified ILS and sampling rate:
I s i m u λ = I r e f λ R λ
Differences in the spectral resolution of each probe element in the band were not considered in this study. The sampling center wavelength λ c of each spectral channel was calculated using the first-order spectral dispersion equation, i.e.,
λ c = λ 0 + c o l u m n F W H M / N s r
where λ 0 denotes the starting wavelength of the band, N s r is the sampling rate, and c o l u m n is the serial number of the spectral dimension probe element.
Considering the instrumental observation parameters, including solar zenith angle, incidence angle, Sun–Earth distance, satellite relative velocity to the Sun, diffuser BRDF, A/D quantization efficiency [40], etc., the simulated ideal DN values of the solar observation without considering the dark background signals and noises could be calculated, as shown in Figure 3b. Typical values were used for the instrumental observation parameters, and 14-bit quantization was used for the A/D conversion. By adding the instrument noise signal and radiometric calibration uncertainty, the simulated observation spectral radiance could be calculated. Figure 4 shows the whole simulation process of the instrumental solar measurement spectrum.
For the estimation of instrumental noise, the observation spectra were simulated at the system level in terms of the SNR. The SNR was calculated in terms of the electron number. The noise electron number N of the system can be considered to be composed of four components, including photon noise S generated by the input signal, scattering noise S d a r k generated by the dark current of the detector, scattering noise S b a c k generated by the thermal background, and electronic readout noise N r e a d .
N = S + S d a r k + S b a c k + N r e a d 2
where S represents the spectral irradiance measured by the detector and can be expressed by the following equation:
S = L λ , r , ϕ × τ × Q E × A a p × ω × N s r × Δ λ × t i n t
where L λ , r , ϕ is the spectral irradiance detected by the optical system at the solar zenith angle ϕ and ground albedo r , τ is the transmittance of the optical system, Q E is the quantum efficiency of the detector, A a p is the area of the entry pupil, ω is the instantaneous field-of-view stereo angle of the system, and Δ λ is the spectral bandwidth of the detector’s probe element, the value of which is equal to the spectral resolution of the instrument divided by the sampling rate N s r , while t i n t is the integration time. Equation (9) gives the instrument SNR:
S N R = S N
In this study, the SNR at different wavelengths was calculated using the typical spectral irradiance for solar observation of the instrument as input. Ten noise values were calculated for each channel signal based on the SNR combined with the A/D quantization error, and then a random array of numbers was used to add noise to the simulated measurement signals for each channel.
The inter-channel relative radiometric calibration uncertainty is the difference in the consistency of the radiometric calibration of the GHG Monitor across the spectral dimensions of the instrument. For different spectral dimension detectors of the instrument in the same spatial dimension, there is an inconsistency in the radiometric response, which needs to be corrected in the laboratory and during on-orbit operation [41]. However, due to the existence of uncertainties in the radiometric calibration, there is a certain error in the relative radiometric calibration between spectral channels, i.e., the inconsistency in the radiometric response between spectral channels of the instrument cannot be accurately corrected, which leads to uncertainties in the quantitative measurement spectra.
The effect of adding the relative radiometric calibration uncertainty between channels takes the form of adding a percentage uncertainty to the ideal simulated measurement spectra. The radiance of the simulated spectra after adding relative radiometric calibration uncertainty to each channel is given in the following equation:
X _ U i = X i × ( 1 + P u n c e r a i n t y )
where X i is the ideal measured spectral radiance of the i th spectral channel; P u n c e r a i n t y represents the relative radiometric calibration uncertainty between the channels, which is a random value in the interval σ r , σ r ; and σ r is the maximum value of uncertainty.
Figure 5 shows the comparison of the simulated spectrum of the instrumental solar observation in the SCO2 band before and after adding measurement uncertainty caused by instrumental noise and relative radiometric calibration uncertainty.

2.4. Wavelength Calibration Method

On-orbit wavelength calibration based on the solar Fraunhofer lines determines the position of the absorption peak center by extracting the observed spectra in the vicinity of the Fraunhofer line and using a peak-seek algorithm, typically the Gaussian fitting algorithm. The Gaussian fitting algorithm uses the Gaussian function
y x = A e x p ( x μ 2 2 σ 2 )
to optimally fit the Fraunhofer line of the normalized measurement spectra. It uses the least squares
a l g o r i t h m ( y x I ( x ) ) 2 = m i n
as the basis for optimization, and the central position of the fitted optimal Gaussian function μ o p t as the peak wavelength position, as shown in Figure 6. The wavelength offset is obtained by comparing the fitted peak-seek positions with the positions of the Fraunhofer absorption peaks in the reference spectrum. Figure 7 shows the flow of the on-orbit wavelength calibration method based on the Fraunhofer lines. In the simulation of the measurement spectra in this study, no additional instrumental wavelength offsets were added, so the offset can be regarded as a fixed error for this on-orbit wavelength calibration method under the specific instrumental design parameters.

3. Results and Discussions

3.1. Doppler Shift Correction Error

When the instrument is in orbit for solar observation, the relative velocity between the instrument and the Sun will produce a wavelength shift, called the Doppler shift, which was first discovered and proposed by the Austrian mathematician Doppler. This reflects the relationship between the frequency of the signal and the speed of motion, which is specifically manifested as follows: when the wave source and the receiver are close to each other, the wave is compressed, the wavelength becomes shorter, and the received frequency is greater than the emitted frequency; when the wave source and the receiver are far away from each other, it will produce the opposite effect; the larger the relative speed of motion between the wave source and the receiver, the greater the effect produced.
According to the theory of the Doppler effect, the observed solar spectral frequency increases when the satellite and the Sun are close to each other and decreases when they are far away from each other. The Doppler effect can be calculated from the relative velocities of the satellite and the Sun, and the revision of the Doppler effect satisfies the following equation:
f d = f ( 1 + V r e l c )
where f d denotes the Doppler shift, c is the speed of light, and f is the frequency of the original solar irradiance. V r e l is the relative velocity of the satellite and the Sun, and when V r e l is positive, it means that the two are close to each other.
Taking V r e l = 7.0   k m · s 1 as an example, the wavelength shifts due to the Doppler effect in the three bands are shown in Figure 8. Because of the narrow wavelength range of bands, the wavelength shift due to the Doppler effect at the central wavelength of different spectral channels is very small, which is negligible compared to the spectral resolution of the instrument, and therefore the wavelength shifts of the bands due to the Doppler effect can be analyzed as a whole. The average wavelength shifts of the three bands generated by the Doppler effect are 0.0178 nm, 0.0376 nm, and 0.0476 nm, respectively, which are about FWHM/2 of each band.
When the instrument is calibrated in orbit, the relative velocity calculation error between the satellite and the Sun due to on-orbit timing errors, etc., will cause a Doppler shift correction error. From the derivation, it can be seen that δ relative velocity calculation error will cause the same δ Doppler shift correction error, i.e., 2% relative velocity calculation error will cause about FWHM/100 error in the wavelength shift correction. Figure 9 shows the effect of ±5% relative velocity calculation error on the average wavelength shift of the three bands caused by the Doppler effect, using V r e l = 7.0   k m · s 1 as an example. The average wavelength shift of the three bands caused by a 5% relative velocity calculation error is 0.00089 nm, 0.00190 nm, and 0.00240 nm, respectively.

3.2. Systematic Error

When the instrument is loaded and adjusted, the spectral resolution, sampling rate, and ILS of the instrument and the central wavelength of each spectral channel probe are determined without taking into account the wavelength drift and ILS change. At this time, due to the spectral broadening caused by the instrumental line function, the observed spectrum obtained after convolution of the solar spectrum near the Fraunhofer lines may be affected by the adjacent absorption lines or the asymmetry of the Fraunhofer line themselves, resulting in a deviation between the peak position after convolution and the actual Fraunhofer peak position. This deviation is due to the instrument design parameters, and it is determined under the specified parameters, which can be regarded as the systematic error of this wavelength calibration method. In this case, the peak-seek deviation is mainly caused by the adjacent absorption peaks of the Fraunhofer line and its own asymmetry, and the systematic error is larger for the solar Fraunhofer lines that are more affected by this factor. Taking the positions of the Fraunhofer lines at 1612.039 nm and 1613.034 nm in the WCO2 band as an example, the convolution spectra of the two absorption lines are shown in Figure 10b. The deviations between the peak positions after applying the peak-seek algorithm using Gaussian fitting to the normalized simulated observed spectra and the real solar Fraunhofer peak positions are 0.0115 nm and 0.0004 nm, respectively, as shown in Figure 10c. It can be seen that there is a large deviation due to the effect of the adjacent absorption on the left side of the Fraunhofer line at 1612.039 nm. The relatively symmetrical Fraunhofer line at 1613.034 nm is less affected, at an almost negligible level.
The central wavelength’s position for each sampling spectral point can vary between positive and negative half-wavelength sampling intervals, with which the deviation of the peak-seek position of the observed Fraunhofer line varies. As shown in Figure 11, the systematic error fluctuates around the mean value with the shift in the center wavelength position of each sampling point, but the fluctuation range is small. This result also verifies that the wavelength calibration method based on the solar Fraunhofer lines is reasonable under the spectral resolution of the instrument. With the wavelength shift in the sampling point, the deviation of the peak-seek positions of the observed Fraunhofer lines from the positions of the absorption peaks of the solar reference Fraunhofer lines is still affected by the systematic error only, but almost not at all by the instrumental wavelength shift factor. Therefore, in the actual instrument’s on-orbit wavelength calibration, after removing the systematic errors, the deviation of the peak-seek position of the Fraunhofer peak can be regarded as the actual wavelength shift of the instrument.
Increasing the spectral resolution can improve the depiction of the observed spectra, especially the spectra at the absorption peak positions. Therefore, under other conditions, increasing the spectral resolution can enable the instrument to measure the spectra of the absorption peaks of the solar Fraunhofer lines more accurately, thus reducing the systematic error, especially for the more asymmetric lines. The variation in the peak-seek error at 1612.039 nm with the spectral resolution is shown in Figure 12. Figure 13 shows the change in the peak-seek systematic error with the increase in the spectral resolution of the instrument. It can be seen that, for the solar Fraunhofer lines with large systematic errors due to strong asymmetry, such as at 1598.513 nm, increasing the spectral resolution has a more obvious reduction effect on the systematic errors, while for the Fraunhofer lines with small systematic errors, the effect is negligible.
The super-Gaussian ILS has a gentler peak position as the order n increases. As can be seen from the results shown in Figure 14, it has little effect on the peak-seek error of the Fraunhofer lines. In addition, under the condition of a certain spectral resolution, the higher the number of sampling points, the smoother and more regular the measured spectral curve is, and the finer the absorption lines characterized by the curve. However, from the analysis, increasing the spectral sampling rate at a certain spectral resolution will not effectively improve the accuracy of the peak-seek position of the Fraunhofer line.
From the above analysis, under the condition of a certain spectral resolution, sampling rate, and ILS, due to the asymmetry of the Fraunhofer lines and the influence of the adjacent absorption lines, after the instrumental measurement, the position of the peak-seek using the Gaussian fitting algorithm will have a fixed deviation from the actual solar Fraunhofer peak position. This error can be significantly reduced with the enhancement of the spectral resolution of the instrument but with the super-Gaussian function order of ILS and sampling point of the center wavelength position shift having less impact. In practice, since the design parameters of the instrument have already been determined when the solar Fraunhofer lines are utilized for on-orbit wavelength calibrations, the Fraunhofer lines that do not satisfy the wavelength calibration accuracy due to asymmetry need to be screened and eliminated. In addition, since the systematic error fluctuates less with the shift in the central wavelength position of each sampling point, the effect of the systematic error can be reduced by removing the average deviation, so as to improve the wavelength calibration accuracy. Figure 15 shows the maximum deviation of the peak-seek positions of the measured spectra from the actual Fraunhofer peak positions after removing the influence of the average systematic error. It can be seen that the peak-seek position errors of the 21 Fraunhofer lines are significantly reduced after the average systematic errors are corrected, indicating that the peak-seek errors of each Fraunhofer line in the range of the observed bands vary less with the central wavelength position of the sampling point. In addition, the errors of 1612.039 nm and 1619.025 nm in the WCO2 band, which are larger than the wavelength calibration accuracy of FWHM/10, meet the requirements after correction. It can be seen that the screening of better symmetry of the Fraunhofer line, as well as the correction of the systematic error, is essential to further improving the wavelength calibration accuracy.

3.3. Random Error

In addition to the systematic error associated with the instrument design parameters, the determination of the wavelength position of the observed solar Fraunhofer lines is also affected by random errors due to the noise signal and the uncertainty of the radiometric calibration. The random error mainly has a greater impact on the absorption peaks with weaker intensity. Taking the SCO2 band as an example, the effects of instrumental SNR and inter-channel relative radiometric calibration uncertainty on the wavelength position seeking of the observed Fraunhofer lines are analyzed.
In order to analyze the effect of observation noise on the wavelength calibration accuracy based on the Fraunhofer lines, the Gaussian fitting of the simulated measurement spectra with the addition of the noise signals is performed to seek the peak at the location of the Fraunhofer lines. Figure 16 shows the peak-seek obtained using Gaussian fitting before and after adding noise at the 2045.0 nm Fraunhofer line, and it can be seen that random fluctuations occur in the observed spectral signals under the influence of noise, which brings uncertainty to the fitted peak-seek position. The experiment was repeated 400 times to statistically measure the effect of the noise signal on the peak-seek error. The peak-seek error is first processed to remove the systematic error before statistics as shown in Figure 17, so that the peak-seek error only characterizes the effect of random error.
Calculating the peak-seek error of each Fraunhofer line in experiments, Figure 18 shows the effect of the peak-seek error of the three Fraunhofer lines at 2045.0 nm, 2054.431 nm, and 2074.254 nm with the decrease in SNR. It can be seen that as the SNR decreases, the three Fraunhofer lines show a gradual dispersion of the random error, the maximum random error increases, and the absorption intensity of the Fraunhofer line is affected to a greater extent.
The variation in the maximum peak-seek error with SNR for each Fraunhofer line in the SCO2 band was statistically analyzed. As can be seen from the statistical curve in Figure 19, the maximum random error gradually increases with the decrease in the SNR, especially for the Fraunhofer lines at 2045.0 nm, 2054.431 nm, and 2056.96 nm, which is due to the fact that the peak absorption intensity of the three lines is weaker. The Fraunhofer line at 2054.431 nm is most affected by the noise signal, and the phenomenon of fitting peak shift and a large random error occurs when the SNR is higher, so the Fraunhofer line should be considered to be eliminated during wavelength calibration. In addition, the random errors of the peak-seek at 2045.0 nm and 2056.96 nm when the SNR is lower than 700 exceed FWHM/20, which affects the wavelength calibration accuracy to a certain extent. The other Fraunhofer lines are less affected by the noise signal, and the random error generated by the measurement noise has a negligible effect on the wavelength calibration accuracy.
The process of analyzing the random errors generated by the radiometric calibration uncertainty is consistent with that used to analyze the effect of the noise signal. Since the radiometric calibration uncertainty affects the observed spectra in a percentage multiplicative form, it has a greater impact than the superimposed form of the noise signal. As shown in Figure 20, the maximum random errors for the Fraunhofer lines at 2045.0 nm, 2054.431 nm, and 2056.96 nm are still larger. Among them, the peak-seek error of 2054.431 nm is too large when the uncertainty is more than 0.1%, and the peak-seek error of 2045.0 nm, 2056.96 nm is greater than FWHM/10 when the uncertainty is more than 0.2%. In practical applications, it is necessary to control inter-channel relative radiometric calibration accuracy at 0.2% or less to meet the wavelength calibration accuracy requirements or eliminate those Fraunhofer lines that are highly affected by the radiometric calibration’ uncertainty.
For the O2A and WCO2 bands, the peak-seek errors are less affected by random errors due to the higher SNR and the high intensity of the Fraunhofer lines. The random error of the peak-seek due to measurement noise and inter-channel relative radiometric calibration uncertainty is less than FWHM/20 for the Fraunhofer lines in both bands. Therefore, the random errors of wavelength calibration using the Fraunhofer lines are negligible for the O2A and WCO2 bands.

3.4. Peak-Seek Algorithm Error

In addition to the Gaussian fitting algorithm, the center-of-mass algorithm and the cubic spline fitting algorithm are also commonly used for peak wavelength detection in spectral analysis, as shown in Figure 21. The center-of-mass algorithm determines the center wavelength of the measured spectrum by calculating the weighted average of the signals. It is calculated by first multiplying the wavelength value of each spectral data point by the measurement value at its corresponding wavelength and then adding all the results and dividing them by the sum of the measurement values to obtain the weighted average wavelength, which serves as the peak position of the spectrum. Let λ i and I i be the wavelength value and spectral radiance at data point i , respectively, and then the peak wavelength of this spectral signal in the range of data points i ~ k is
λ c = i = 0 k λ i × I i i = 0 k I i
The cubic spline fitting algorithm approximates the measured spectral signal data in the form of the following equation:
f = p i = 0 n 1 w i ( y i f x i ) 2 + ( 1 p ) x 0 x n 1 λ ( x ) ( f n ( x ) ) 2 d x ,
and then the peak wavelength position of the signal is obtained by utilizing the extreme value peak-seek method for the fitted signal. w i is the i th element of the weights, λ ( x ) is the segmented constant function, and p is the equilibrium parameter, which takes values in the interval [ 0 ,   1 ] to make the fitted curves smooth and close to the data points.
In order to analyze the influence of peak-seek algorithms on the accuracy of the peak-seek position for the observed solar Fraunhofer lines, three algorithms were used in turn to calculate the peak position and statistically estimate the peak-seek position error. Among them, the center-of-mass algorithm is sensitive to the length of selected sampling points, so 3 × N s r sampling points were intercepted for the calculation. In addition, the length of the selected sampling points will also have a certain degree of influence on the results of the cubic spline function fitting algorithm and the Gaussian fitting algorithm, and 6 × N s r sampling points were intercepted for calculation in this study.
Comparing the results of the three peak-seek algorithms shown in Figure 22, the center-of-mass algorithm is the most sensitive to the asymmetry of the Fraunhofer lines, followed by the Gaussian fitting algorithm, and the cubic spline function fitting algorithm has a higher peak-seek accuracy for the asymmetry Fraunhofer lines. However, compared with the cubic spline fitting algorithm, the Gaussian fitting algorithm is less affected by the central wavelength shift of the sampling point, so the maximum peak-seek error after systematic error correction is smaller than the cubic spline fitting algorithm. As shown in Figure 23, the two Fraunhofer lines at 1612.039 nm and 1613.034 nm are still analyzed, and for the cubic spline fitting peak-seek algorithm, the maximum deviation of peak-seek error due to the central wavelength shift at the two points is 5.9 × 10 4 nm and 5.5 × 10 4 nm, respectively, while for the Gaussian fitting algorithm, the maximum deviation at the two points is 5.1 × 10 5 nm and 4.1 × 10 5 nm. Therefore, in the actual wavelength calibration, the Gaussian fitting algorithm has higher correction accuracy for the wavelength shift of the instrument, but it is necessary to select the absorption lines with a higher degree of asymmetry or to calculate and revise the systematic error.
In addition, for the simulated measurement spectra obtained by convolving the ILS in the form of a super-Gaussian function, the peak-seek calculations were also considered by fitting the super-Gaussian function of the corresponding order. For the simulated spectra calculated using the super-Gaussian ILS (n = 2.5 and n = 3.0), the Fraunhofer lines in the three bands were fitted with the corresponding order super-Gaussian function, and the results are shown in Figure 24. From the results, it can be seen that the peak-seek obtained using the corresponding super-Gaussian function of ILS does not significantly affect the accuracy of the peak-seek position.

3.5. Results of Screening

Using the above analytical method, the wavelength calibration error was calculated for all the solar Fraunhofer lines in the three bands. Firstly, the Fraunhofer lines with small absorption intensity or narrow wavelength width and the lines greatly affected by the adjacent absorption peaks and their asymmetry were eliminated. The remaining available solar Fraunhofer lines were as follows: 21 in the O2A band, 25 in the WCO2 band, and 19 in the SCO2 band. The 65 Fraunhofer lines were distinguished according to the wavelength calibration systematic error, and the results are shown in Figure 25. Among them, 56 lines met the requirement of FWHM/10 wavelength calibration accuracy, including 46 lines with a calibration accuracy less than FWHM/20, and 9 lines needed to be corrected for systematic error to meet the calibration accuracy requirements.
In addition, the random errors of wavelength calibration of 19 Fraunhofer lines in the SCO2 band were analyzed, as shown in Figure 26, showing the random errors of peak-seek positions of Fraunhofer lines in the SCO2 band at SNR = 700 and Uncertainty = 0.2%. It can be seen that, affected by the random error, the number of Fraunhofer lines in the SCO2 band that meet the requirements of FWMH/10 wavelength calibration accuracy is reduced to 11, and two of them at 2050.339 nm and 2070.436 nm are subject to systematic error correction.

4. Conclusions

It is feasible to utilize the solar Fraunhofer lines in the solar observation spectrum for on-orbit wavelength calibration under the design parameters of GHG Monitor, and the accuracy can meet the requirements of less than FWHM/10. The error source of this method is mainly due to the asymmetry of the Fraunhofer lines and the random error of the instrument measurement, which can reach a maximum of more than FWHM/10 and can be reduced by improving the spectral resolution, suppressing the noise of the instrument, and improving the accuracy of the radiometric calibration. The Doppler shift correction error caused by the relative velocity calculation error between the satellite and the Sun is limited, and its influence on the wavelength calibration accuracy is less than FWHM/40 under the condition that the relative velocity calculation error is less than 5%, which is almost negligible. Regarding the peak-seek algorithm, the center-of-mass method is the least applicable, and the cubic spline function fitting algorithm has a higher peak-seek accuracy for the asymmetry Fraunhofer lines compared with the Gaussian fitting, but the Gaussian fitting algorithm is less affected by the central wavelength shift of the sampling point.
Under the current parameter design, the number of Fraunhofer lines that can meet the wavelength calibration accuracy of less than FWHM/10 without considering the systematic error reaches 19 in the O2A band and 21 in the WCO2 band, which is sufficient for the on-orbit calibration. For the SCO2 band, some of the Fraunhofer lines are greatly affected by the measurement noise and uncertainty of radiometric calibrations, especially due to the relative radiometric calibration uncertainty between spectral channels. There are nine directly usable Fraunhofer lines in the SCO2 band considering measurement uncertainty, and therefore it is necessary to accurately correct the differences in the radiometric responses between different spectral channels in both the laboratory and on-orbit phases.
Although this study provides valuable insights into the sources of errors in the wavelength calibration method based on the solar Fraunhofer lines, it is essential to acknowledge its inherent limitations. Firstly, in practical on-orbit wavelength calibrations of hyperspectral GHG detection instruments using the solar Fraunhofer lines, the actual values of each instrumental parameter, including ILS, central wavelength, etc., should be used to simulate instrumental solar-viewing measurement. In the simulations in this study, the ideal ILS of each probe element was assumed, and its possible aberrations as well as subtle spectral resolution differences between channels were not considered. Secondly, on-orbit observations of instruments are affected by a number of factors, making it difficult to simulate them accurately. When studying the effect of random errors on the accuracy of wavelength calibration, this study used a random array of numbers to add noise and inter-channel relative radiometric calibration uncertainty to the simulated measurement spectra, which might bring uncertainty to the results. Thirdly, in addition to the discussions in this study, the sources of errors in the on-orbit wavelength calibration using the solar Fraunhofer lines also include uncertainties in solar reference spectra, BRDF of the solar diffuser, and so on, all of which can affect the accuracy on-orbit wavelength calibration of instruments. To more accurately assess the accuracy of on-orbit wavelength calibration using the solar Fraunhofer lines, it is necessary to further investigate each source of uncertainty. Moreover, further study of the complex relationships between instrument parameters and measurements will contribute to a comprehensive understanding of the detection process, leading to more accurate simulation of instrument observations.
In conclusion, based on the results of this study, it is necessary to select the Fraunhofer lines with a higher degree of symmetry and absorption intensity in solar-viewing measurement spectra for calibration and conduct a comprehensive analysis of various sources of calibration error in accordance with the on-orbit performance of the instrument. This study analyzed the main sources of errors in the wavelength calibration method based on the solar Fraunhofer lines and provides a quantitative assessment of the impact of instrument parameters on the accuracy of wavelength calibration. These insights not only aid in enhancing the accuracy of this on-orbit wavelength calibration method but also offer avenues for optimizing the design of similar GHG detection payloads. Furthermore, it should be noted that multi-combined on-orbit calibration modes, including solar observation and Earth observation, should be used in combination to obtain more accurate spectral correction results.

Author Contributions

Conceptualization, C.G. and Y.H.; methodology, Y.G.; formal analysis, Y.G. and F.Z.; resources, C.G. and Y.H.; data curation, Y.G. and F.Z.; writing—original draft preparation, Y.G.; writing—review and editing, Y.G. and Y.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Key R&D Program of China (Grant No. 2023YFB3907501): Laboratory calibration and data processing of ultra-wide carbon pollution co-monitor (infrared channel); the National Key R&D Program of China (Grant No. 2022YFB3902000): Very long wavelength infrared (VLWIR) hyperspectral imaging technology; and the Shanghai Pilot Program for Basic Research–Chinese Academy of Sciences, Shanghai Branch (JCYJ-SHFY-2022-004).

Data Availability Statement

The Kurucz solar spectrum used in this study is available at http://kurucz.harvard.edu/, which was accessed on 27 March 2024, and other simulation data are contained within this article.

Acknowledgments

We are thankful to all anonymous reviewers for providing their constructive comments on this study.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Friedlingstein, P.; O’Sullivan, M.; Jones, M.W.; Andrew, R.M.; Hauck, J.; Olsen, A.; Peters, G.P.; Peters, W.; Pongratz, J.; Sitch, S.; et al. Global Carbon Budget 2020. Earth Syst. Sci. Data 2020, 12, 3269–3340. [Google Scholar] [CrossRef]
  2. Feng, L.; Palmer, P.I.; Zhu, S.H.; Parker, R.J.; Liu, Y. Tropical methane emissions explain large fraction of recent changes in global atmospheric methane growth rate. Nat. Commun. 2022, 13, 1378. [Google Scholar] [CrossRef]
  3. Bruhwiler, L.; Basu, S.; Butler, J.H.; Chatterjee, A.; Dlugokencky, E.; Kenney, M.A.; McComiskey, A.; Montzka, S.A.; Stanitski, D. Observations of greenhouse gases as climate indicators. Clim. Chang. 2021, 165, 12. [Google Scholar] [CrossRef] [PubMed]
  4. Gautam, Y.K.; Sharma, K.; Tyagi, S.; Ambedkar, A.K.; Chaudhary, M.; Singh, B.P. Nanostructured metal oxide semiconductor-based sensors for greenhouse gas detection: Progress and challenges. R. Soc. Open Sci. 2021, 8, 201324. [Google Scholar] [CrossRef]
  5. Chatterjee, A.; Gierach, M.M.; Sutton, A.J.; Feely, R.A.; Crisp, D.; Eldering, A.; Gunson, M.R.; O’Dell, C.W.; Stephens, B.B.; Schimel, D.S. Influence of El Nino on atmospheric CO over the tropical Pacific Ocean: Findings from NASA’s OCO-2 mission. Science 2017, 358, eaam5776. [Google Scholar] [CrossRef] [PubMed]
  6. Liu, Y.; Wang, J.; Che, K.; Cai, Z.; Yang, D.; Wu, L. Satellite remote sensing of greenhouse gases: Progress and trends. Natl. Remote Sens. Bull. 2021, 25, 53–64. [Google Scholar] [CrossRef]
  7. Liu, J.J.; Baskaran, L.; Bowman, K.; Schimel, D.; Bloom, A.A.; Parazoo, N.C.; Oda, T.; Carroll, D.; Menemenlis, D.; Joiner, J.; et al. Carbon Monitoring System Flux Net Biosphere Exchange 2020 (CMS-Flux NBE 2020). Earth Syst. Sci. Data 2021, 13, 299–330. [Google Scholar] [CrossRef]
  8. Cai, B.; Zhu, S.; Yu, S.; Dong, H.; Zhang, C.; Wang, C.; Zhu, J.; Gao, Q.; Fang, S.; Pan, X.; et al. The interpretation of 2019 refinement to the 2006 IPCC guidelines for national greenhouse gas inventory. Environ. Eng. 2019, 37, 1–11. [Google Scholar]
  9. Buchwitz, M.; Schneising, O.; Burrows, J.P.; Bovensmann, H.; Reuter, M.; Notholt, J. First direct observation of the atmospheric CO year-to-year increase from space. Atmos. Chem. Phys. 2007, 7, 5341–5342. [Google Scholar] [CrossRef]
  10. Yokota, T.; Yoshida, Y.; Eguchi, N.; Ota, Y.; Tanaka, T.; Watanabe, H.; Maksyutov, S. Global Concentrations of CO2 and CH4 Retrieved from GOSAT: First Preliminary Results. Sola 2009, 5, 160–163. [Google Scholar] [CrossRef]
  11. Frankenberg, C.; Pollock, R.; Lee, R.A.M.; Rosenberg, R.; Blavier, J.F.; Crisp, D.; O’Dell, C.W.; Osterman, G.B.; Roehl, C.; Wennberg, P.O.; et al. The Orbiting Carbon Observatory (OCO-2): Spectrometer performance evaluation using pre-launch direct sun measurements. Atmos. Meas. Tech. 2015, 8, 301–313. [Google Scholar] [CrossRef]
  12. Suto, H.; Kataoka, F.; Kikuchi, N.; Knuteson, R.O.; Butz, A.; Haun, M.; Buijs, H.; Shiomi, K.; Imai, H.; Kuze, A. Thermal and near-infrared sensor for carbon observation Fourier transform spectrometer-2 (TANSO-FTS-2) on the Greenhouse gases Observing SATellite-2 (GOSAT-2) during its first year in orbit. Atmos. Meas. Tech. 2021, 14, 2013–2039. [Google Scholar] [CrossRef]
  13. Liu, Y.; Wang, J.; Yao, L.; Chen, X.; Cai, Z.N.; Yang, D.X.; Yin, Z.S.; Gu, S.Y.; Tian, L.F.; Lu, N.M.; et al. The TanSat mission: Preliminary global observations. Sci. Bull. 2018, 63, 1200–1207. [Google Scholar] [CrossRef] [PubMed]
  14. Jervis, D.; McKeever, J.; Durak, B.O.A.; Sloan, J.J.; Gains, D.; Varon, D.J.; Ramier, A.; Strupler, M.; Tarrant, E. The GHGSat-D imaging spectrometer. Atmos. Meas. Tech. 2021, 14, 2127–2140. [Google Scholar] [CrossRef]
  15. Veefkind, J.P.; Aben, I.; McMullan, K.; Förster, H.; de Vries, J.; Otter, G.; Claas, J.; Eskes, H.J.; de Haan, J.F.; Kleipool, Q.; et al. TROPOMI on the ESA Sentinel-5 Precursor: A GMES mission for global observations of the atmospheric composition for climate, air quality and ozone layer applications. Remote Sens. Environ. 2012, 120, 70–83. [Google Scholar] [CrossRef]
  16. Xiong, W. Greenhouse gases Monitoring Instrument (GMI) on GF-5 satellite (invited). Infrared Laser Eng. 2019, 48, 0303002. [Google Scholar] [CrossRef]
  17. Shao, C.; Gu, M.; Qi, C.; Wu, C. Detection of zero path difference position for FY-3D hyper-spectral infrared atmospheric sounder. Opt. Precis. Eng. 2020, 28, 2573–2580. [Google Scholar] [CrossRef]
  18. Crisp, D.; Fisher, B.M.; O’Dell, C.; Frankenberg, C.; Basilio, R.; Bösch, H.; Brown, L.R.; Castano, R.; Connor, B.; Deutscher, N.M.; et al. The ACOS CO2 retrieval algorithm—Part II: Global XCO2 data characterization. Atmos. Meas. Tech. 2012, 5, 687–707. [Google Scholar] [CrossRef]
  19. Green, R.O. Spectral calibration requirement for earth-looking imaging spectrometers in the solar-reflected spectrum. Appl. Opt. 1998, 37, 683–690. [Google Scholar] [CrossRef]
  20. Li, S.Z.; Cheng, L.; Yang, H.C.; Wang, Z.W.; Ding, L. Assessment of the Influence of Instrument Parameters on the Detection Accuracy of Greenhouse-Gases Absorption Spectrometer-2 (GAS-2). Atmosphere 2023, 14, 1418. [Google Scholar] [CrossRef]
  21. Crisp, D.; Pollock, H.R.; Rosenberg, R.; Chapsky, L.; Lee, R.A.M.; Oyafuso, F.A.; Frankenberg, C.; O’Dell, C.W.; Bruegge, C.J.; Doran, G.B.; et al. The on-orbit performance of the Orbiting Carbon Observatory-2 (OCO-2) instrument and its radiometrically calibrated products. Atmos. Meas. Tech. 2017, 10, 59–81. [Google Scholar] [CrossRef]
  22. Bi, Y.; Wang, Q.; Yang, Z.; Liu, C.; Lin, C.; Tian, L.; Zhang, N.; Wang, Y. TanSat ACGS On-orbit Wavelength Calibration Using the Solar Fraunhofer Lines. Chin. J. Atmos. Sci. 2022, 46, 645–652. [Google Scholar]
  23. Voors, R.; Dobber, M.; Dirksen, R.; Levelt, P. Method of calibration to correct for cloud-induced wavelength shifts in the Aura satellite’s Ozone Monitoring Instrument. Appl. Opt. 2006, 45, 3652–3658. [Google Scholar] [CrossRef]
  24. Xian, G.; Yan, C.; Wu, C.; Zhang, J. Effect of temperature on airborne imaging spectrometer optical properties. Infrared Laser Eng. 2015, 44, 1647–1653. [Google Scholar]
  25. Santantonio, N.; Bartoloni, A.; Bini, A.; Esposito, P.G.; Lorenzini, S. A Multi Parameter Retrieval Approach to HypSEO Spectral and Radiometric Calibration. In IEEE International Geoscience and Remote Sensing Symposium, Proceedings of IGARSS 2002, Toronto, ON, Canada, 24–28 June 2002; IEEE: New York, NY, USA, 2002; pp. 2735–2737. [Google Scholar]
  26. Mao, Y.L.; He, T.; Ren, F.; Huang, J. The On-orbit Calibration Method of Hyper-Spectral Fluorescence Imaging Spectrometer. In Fifth Symposium on Novel Optoelectronic Detection Technology and Application, Xi’an, China, 24–26 October 2018; SPIE: Bellingham, WA, USA, 2019; p. 11023. [Google Scholar]
  27. Li, B.; Li, M.; Xu, P. High Precision On-board Calibration for Spaceborne Atmospheric Monitoring Spectrometer. Spacecr. Recovery Remote Sens. 2018, 39, 60–69. [Google Scholar]
  28. Wang, T.X.; Yan, G.J.; Ren, H.Z.; Mu, X.H. Improved Methods for Spectral Calibration of On-Orbit Imaging Spectrometers. IEEE Trans. Geosci. Remote Sens. 2010, 48, 3924–3931. [Google Scholar] [CrossRef]
  29. Sornig, M.; Fischer, S.; Chlebek, C.; Muecke, M.; Honold, H.P.; Heider, B. The Hyperspectral Instrument Onboard EnMAP, Overview and Current Status. In International Conference on Space Optics-Icso 2018, Chiania, Greece, 9–12 October 2018; SPIE: Bellingham, WA, USA, 2018; p. 11180. [Google Scholar]
  30. Wang, Y.; Shi, E.; Wang, H. Method of in-flight wavelength calibration for wide band hyper-spectral imager. Infrared Laser Eng. 2019, 48, 303006. [Google Scholar] [CrossRef]
  31. Mori, S.; Cho, Y.C.R.; Tabata, H.; Yumoto, K.; Böttger, U.; Buder, M.; Dietz, E.; Hagelschuer, T.; Hübers, H.W.; Kameda, S.; et al. Fraunhofer line-based wavelength-calibration method without calibration targets for planetary lander instruments. Planet. Space Sci. 2024, 240, 105835. [Google Scholar] [CrossRef]
  32. Delwart, S.; Preusker, R.; Bourg, L.; Santer, R.; Ramon, D.; Fischer, J. MERIS in-flight spectral calibration. Int. J. Remote Sens. 2007, 28, 479–496. [Google Scholar] [CrossRef]
  33. Kuli-Zade, D.M.; Mamedov, S.G.; Aliyeva, Z.F. The fraunhofer lines classifications by asymmetry changes inside the profiles in the spectrum of the Sun. Astrophys. Bull. 2017, 72, 337–339. [Google Scholar] [CrossRef]
  34. Heath, D.F.; Georgiev, G. A New Approach for Spectroradiometric Calibration Consistency on the Ground and in Space. In Earth Observing Missions and Sensors: Development, Implementation, and Characterization, 2nd ed.; SPIE: Bellingham, WA, USA, 2012; p. 8528. [Google Scholar]
  35. O’Brien, D.M.; Mitchell, R.M.; English, S.A.; Da Costa, G.A. Airborne measurements of air mass from OA-band absorption spectra. J. Atmos. Ocean. Technol. 1998, 15, 1272–1286. [Google Scholar] [CrossRef]
  36. Mao, J.P.; Kawa, S.R. Sensitivity studies for space-based measurement of atmospheric total column carbon dioxide by reflected sunlight. Appl. Opt. 2004, 43, 914–927. [Google Scholar] [CrossRef] [PubMed]
  37. Kuang, Z.M.; Margolis, J.; Toon, G.; Crisp, D.; Yung, Y. Spaceborne measurements of atmospheric CO by high-resolution NIR spectrometry of reflected sunlight: An introductory study: Art. no. 1716. Geophys. Res. Lett. 2002, 29, 11. [Google Scholar] [CrossRef]
  38. Fontenla, J.; White, O.R.; Fox, P.A.; Avrett, E.H.; Kurucz, R.L. Calculation of solar irradiances. I. Synthesis of the solar spectrum. Astrophys. J. 1999, 518, 480–499. [Google Scholar] [CrossRef]
  39. Chance, K.; Kurucz, R.L. An improved high-resolution solar reference spectrum for earth’s atmosphere measurements in the ultraviolet, visible, and near infrared. J. Quant. Spectrosc. Radiat. Transf. 2010, 111, 1289–1295. [Google Scholar] [CrossRef]
  40. Li, S.; Li, C.L.; Xu, S.Y.; Zhang, H.; Zheng, Y.Q. Preflight radiometric calibration of a carbon dioxide spectrometer. Meas. Sci. Technol. 2019, 30, 055401. [Google Scholar] [CrossRef]
  41. Yang, Z.D.; Zhen, Y.Q.; Yin, Z.S.; Lin, C.; Bi, Y.M.; Liu, W.; Wang, Q.; Wang, L.; Gu, S.Y.; Tian, L.F. Prelaunch Radiometric Calibration of the TanSat Atmospheric Carbon Dioxide Grating Spectrometer. IEEE Trans. Geosci. Remote Sens. 2018, 56, 4225–4233. [Google Scholar] [CrossRef]
Figure 1. Super-Gaussian functions of different orders n and the other parameters of the function are as follows: a = 1 , λ c = 2060.0 , and F W H M = 0.1 . When n = 2 , the function is identical to a general Gaussian function.
Figure 1. Super-Gaussian functions of different orders n and the other parameters of the function are as follows: a = 1 , λ c = 2060.0 , and F W H M = 0.1 . When n = 2 , the function is identical to a general Gaussian function.
Remotesensing 16 03367 g001
Figure 2. Positions of 21 typical Fraunhofer lines in solar reference irradiance spectra of the three bands, where the positions of Fraunhofer lines are marked by a red “x”.
Figure 2. Positions of 21 typical Fraunhofer lines in solar reference irradiance spectra of the three bands, where the positions of Fraunhofer lines are marked by a red “x”.
Remotesensing 16 03367 g002
Figure 3. Simulations of the ideal instrumental solar observation values in the WCO2 band, where the ILS was simulated using a general Gaussian function with FWHM = 0.08 nm: (a) the high-resolution Kurucz solar irradiance spectrum in the WCO2 band; (b) simulated instrumental solar measurement DN values.
Figure 3. Simulations of the ideal instrumental solar observation values in the WCO2 band, where the ILS was simulated using a general Gaussian function with FWHM = 0.08 nm: (a) the high-resolution Kurucz solar irradiance spectrum in the WCO2 band; (b) simulated instrumental solar measurement DN values.
Remotesensing 16 03367 g003
Figure 4. Flow chart of the simulation to calculate instrumental solar measurement spectrum.
Figure 4. Flow chart of the simulation to calculate instrumental solar measurement spectrum.
Remotesensing 16 03367 g004
Figure 5. Simulation spectrum of the instrumental solar observation in the SCO2 band before (red line) and after (blue line) adding measurement uncertainty, where SNR = 750:[email protected] × 1019  p h o t o n · s 1 · m 2 · s r 1 · μ m 1 and σ r = 0.2%.
Figure 5. Simulation spectrum of the instrumental solar observation in the SCO2 band before (red line) and after (blue line) adding measurement uncertainty, where SNR = 750:[email protected] × 1019  p h o t o n · s 1 · m 2 · s r 1 · μ m 1 and σ r = 0.2%.
Remotesensing 16 03367 g005
Figure 6. Peak-seek using the Gaussian fitting algorithm to the inverted normalized simulated measurement Fraunhofer lines at 1596.446 nm in the WCO2 band: the blue dots represent the inverted normalized simulated measurement values; the black line is the fitted Gaussian function curve; the red “+” is the peak position of fitted Gaussian function curve; and the blue dashed line marks the actual position of the Fraunhofer absorption peaks in the reference spectrum. The peak-seek position of the measured Fraunhofer line at 1596.446 nm is 1596.4427 nm, and the offset of peak-seek is 0.0033 nm.
Figure 6. Peak-seek using the Gaussian fitting algorithm to the inverted normalized simulated measurement Fraunhofer lines at 1596.446 nm in the WCO2 band: the blue dots represent the inverted normalized simulated measurement values; the black line is the fitted Gaussian function curve; the red “+” is the peak position of fitted Gaussian function curve; and the blue dashed line marks the actual position of the Fraunhofer absorption peaks in the reference spectrum. The peak-seek position of the measured Fraunhofer line at 1596.446 nm is 1596.4427 nm, and the offset of peak-seek is 0.0033 nm.
Remotesensing 16 03367 g006
Figure 7. The flow of the on-orbit wavelength calibration method based on the solar Fraunhofer lines.
Figure 7. The flow of the on-orbit wavelength calibration method based on the solar Fraunhofer lines.
Remotesensing 16 03367 g007
Figure 8. The wavelength shifts at the central wavelength of different spectral channels in the three bands caused by the Doppler effect. The Doppler shifts of the starting and ending wavelengths of the three bands are 0.000327 nm, 0.000677 nm, and 0.000910 nm, respectively.
Figure 8. The wavelength shifts at the central wavelength of different spectral channels in the three bands caused by the Doppler effect. The Doppler shifts of the starting and ending wavelengths of the three bands are 0.000327 nm, 0.000677 nm, and 0.000910 nm, respectively.
Remotesensing 16 03367 g008
Figure 9. The average Doppler shift correction error of the three bands caused by ±5% relative velocity calculation when V r e l = 7.0   k m · s 1 (red “x”).
Figure 9. The average Doppler shift correction error of the three bands caused by ±5% relative velocity calculation when V r e l = 7.0   k m · s 1 (red “x”).
Remotesensing 16 03367 g009
Figure 10. Peak-seek algorithm applied to the simulated measurement of Fraunhofer lines at 1612.039 nm and 1613.034 nm in the WCO2 band: (a) the high-resolution Kurucz solar spectrum at 1612.039 nm and 1613.034 nm; (b) the simulated measurement spectra at the ILS with a general Gaussian function and FWHM = 0.08 nm; (c) peak-seek algorithm using Gaussian fitting to the normalized simulated measurement of Fraunhofer lines.
Figure 10. Peak-seek algorithm applied to the simulated measurement of Fraunhofer lines at 1612.039 nm and 1613.034 nm in the WCO2 band: (a) the high-resolution Kurucz solar spectrum at 1612.039 nm and 1613.034 nm; (b) the simulated measurement spectra at the ILS with a general Gaussian function and FWHM = 0.08 nm; (c) peak-seek algorithm using Gaussian fitting to the normalized simulated measurement of Fraunhofer lines.
Remotesensing 16 03367 g010
Figure 11. The peak-seek errors (blue squares) of Fraunhofer lines at 1612.039 nm and 1613.034 nm fluctuate around the mean value (red dashed line) affected by the shift in the central wavelength position of each sampling point.
Figure 11. The peak-seek errors (blue squares) of Fraunhofer lines at 1612.039 nm and 1613.034 nm fluctuate around the mean value (red dashed line) affected by the shift in the central wavelength position of each sampling point.
Remotesensing 16 03367 g011
Figure 12. The peak-seek errors at 1612.039 nm when (a) FWHM = 0.10 nm and (b) FWHM = 0.06 nm. The systematic error caused by the asymmetry of the Fraunhofer line decreases from 0.0138 nm to 0.0093 nm as the spectral resolution changes from 0.10 nm to 0.06 nm.
Figure 12. The peak-seek errors at 1612.039 nm when (a) FWHM = 0.10 nm and (b) FWHM = 0.06 nm. The systematic error caused by the asymmetry of the Fraunhofer line decreases from 0.0138 nm to 0.0093 nm as the spectral resolution changes from 0.10 nm to 0.06 nm.
Remotesensing 16 03367 g012
Figure 13. Variation in the peak-seek error of the 7 Fraunhofer lines in the WCO2 band with spectral resolution, which changes from 0.06 nm to 0.12 nm.
Figure 13. Variation in the peak-seek error of the 7 Fraunhofer lines in the WCO2 band with spectral resolution, which changes from 0.06 nm to 0.12 nm.
Remotesensing 16 03367 g013
Figure 14. Variation in the peak-seek error of the 7 Fraunhofer lines in the WCO2 band with the super-Gaussian function order of ILS, which changes from 1.5 to 3.5.
Figure 14. Variation in the peak-seek error of the 7 Fraunhofer lines in the WCO2 band with the super-Gaussian function order of ILS, which changes from 1.5 to 3.5.
Remotesensing 16 03367 g014
Figure 15. Comparison of the maximum peak-seek errors (“x”) of the 21 Fraunhofer lines (a) before and (b) after removing the average systematic error in three bands.
Figure 15. Comparison of the maximum peak-seek errors (“x”) of the 21 Fraunhofer lines (a) before and (b) after removing the average systematic error in three bands.
Remotesensing 16 03367 g015aRemotesensing 16 03367 g015b
Figure 16. The peak-seek obtained by Gaussian fitting (a) before and (b) after adding noise at the 2045.000 nm Fraunhofer line. The peak-seek error increases from 0.0008 nm to 0.0027 nm due to the influence of the noise signal (SNR = 750), and thus the random error due to noise is 0.0019 nm.
Figure 16. The peak-seek obtained by Gaussian fitting (a) before and (b) after adding noise at the 2045.000 nm Fraunhofer line. The peak-seek error increases from 0.0008 nm to 0.0027 nm due to the influence of the noise signal (SNR = 750), and thus the random error due to noise is 0.0019 nm.
Remotesensing 16 03367 g016
Figure 17. Peak-seek errors of the Fraunhofer lines at (a) 2056.96 nm and (b) 2063.533 nm before (left figure) and after (right figure) removing the systematic error. The blue data points denote the peak-seek position error values of 400 experiments.
Figure 17. Peak-seek errors of the Fraunhofer lines at (a) 2056.96 nm and (b) 2063.533 nm before (left figure) and after (right figure) removing the systematic error. The blue data points denote the peak-seek position error values of 400 experiments.
Remotesensing 16 03367 g017
Figure 18. Peak-seek errors of the Fraunhofer lines at (a) 2045.0 nm, (b) 2054.431 nm, and (c) 2074.254 nm when the SNR is 800, 650, and 500, where the red line is the boundary of FWHM/10.
Figure 18. Peak-seek errors of the Fraunhofer lines at (a) 2045.0 nm, (b) 2054.431 nm, and (c) 2074.254 nm when the SNR is 800, 650, and 500, where the red line is the boundary of FWHM/10.
Remotesensing 16 03367 g018
Figure 19. Variation in the peak-seek error of the 7 Fraunhofer lines in the SCO2 band with SNR changing from 500 to 850.
Figure 19. Variation in the peak-seek error of the 7 Fraunhofer lines in the SCO2 band with SNR changing from 500 to 850.
Remotesensing 16 03367 g019
Figure 20. Variation in the peak-seek error of the Fraunhofer lines in the SCO2 band with the inter-channel relative radiometric calibration uncertainty changing from 0.05% to 0.3%.
Figure 20. Variation in the peak-seek error of the Fraunhofer lines in the SCO2 band with the inter-channel relative radiometric calibration uncertainty changing from 0.05% to 0.3%.
Remotesensing 16 03367 g020
Figure 21. Peak-seek obtained using the (a) center-of-mass algorithm and (b) cubic spline fitting algorithm to the normalized simulated measurement Fraunhofer lines at 1596.446 nm in the WCO2 band.
Figure 21. Peak-seek obtained using the (a) center-of-mass algorithm and (b) cubic spline fitting algorithm to the normalized simulated measurement Fraunhofer lines at 1596.446 nm in the WCO2 band.
Remotesensing 16 03367 g021
Figure 22. The peak-seek errors of the Fraunhofer lines (a) before and (b) after removing the systematic error in three bands using the Gaussian fitting (red “x”), center-of-mass (blue “x”), and cubic spline fitting (green “x”) algorithms.
Figure 22. The peak-seek errors of the Fraunhofer lines (a) before and (b) after removing the systematic error in three bands using the Gaussian fitting (red “x”), center-of-mass (blue “x”), and cubic spline fitting (green “x”) algorithms.
Remotesensing 16 03367 g022
Figure 23. The peak-seek errors of Fraunhofer lines at 1612.039 nm and 1613.034 nm affected by the shift in the center wavelength position of each sampling point using the cubic spline fitting algorithm. The red dashed line is the reference line of the average value of peak-seek errors (blue squares).
Figure 23. The peak-seek errors of Fraunhofer lines at 1612.039 nm and 1613.034 nm affected by the shift in the center wavelength position of each sampling point using the cubic spline fitting algorithm. The red dashed line is the reference line of the average value of peak-seek errors (blue squares).
Remotesensing 16 03367 g023
Figure 24. The peak-seek errors of the Fraunhofer lines in three bands using the super-Gaussian function fitting algorithm with the corresponding order of ILS: (a) the order of the super-Gaussian function is 2.5; (b) the order of the super-Gaussian function is 3.0. The peak-seek errors are after removing the systematic error.
Figure 24. The peak-seek errors of the Fraunhofer lines in three bands using the super-Gaussian function fitting algorithm with the corresponding order of ILS: (a) the order of the super-Gaussian function is 2.5; (b) the order of the super-Gaussian function is 3.0. The peak-seek errors are after removing the systematic error.
Remotesensing 16 03367 g024
Figure 25. Selection of Fraunhofer lines in the three bands according to the systematic error: the black “x” marks indicate that the peak-seek systematic error of the Fraunhofer line is >FWHM/10; the blue triangle marks indicate that the systematic error is <FWHM/10 and >FWHM/20; the red five-pointed star marks indicate that the systematic error is <FWHM/20.
Figure 25. Selection of Fraunhofer lines in the three bands according to the systematic error: the black “x” marks indicate that the peak-seek systematic error of the Fraunhofer line is >FWHM/10; the blue triangle marks indicate that the systematic error is <FWHM/10 and >FWHM/20; the red five-pointed star marks indicate that the systematic error is <FWHM/20.
Remotesensing 16 03367 g025
Figure 26. The random errors of peak-seek positions of 19 Fraunhofer lines in the SCO2 band (SNR = 700, uncertainty = 0.2%): the blue dashed line is the reference line of FWHM/10, and the position of Fraunhofer line is marked by a green “x” if the random error is less than FWHM/10, while it is marked by a red “x” if the random error is greater than FWHM/10.
Figure 26. The random errors of peak-seek positions of 19 Fraunhofer lines in the SCO2 band (SNR = 700, uncertainty = 0.2%): the blue dashed line is the reference line of FWHM/10, and the position of Fraunhofer line is marked by a green “x” if the random error is less than FWHM/10, while it is marked by a red “x” if the random error is greater than FWHM/10.
Remotesensing 16 03367 g026
Table 1. Main design parameters of GHG Monitor.
Table 1. Main design parameters of GHG Monitor.
BandO2AWCO2SCO2
Detection ObjectiveAerosol, surface air pressureCO2CO2, H2O
Wavelength (nm)757.5~772.51595~16252040~2080
Spectral Resolution (nm)0.040.080.10
Sampling Rate3
SNR
( p h o t o n s 1 m 2 s r 1 μ m 1 )
350:[email protected] × 1019340:[email protected] × 1019230:[email protected] × 1018
Dynamic Range
( p h o t o n s 1 m 2 s r 1 μ m 1 )
1.1 × 1017~1.4 × 10215.7 × 1016~4.9 × 10203.1 × 1016~1.7 × 1020
Spatial Resolution (km)<3
Wavelength Calibration Accuracy<FWHM/10
Absolute Calibration Accuracy5%
Inter-channel Relative Calibration Accuracy0.2%
Table 2. The wavelength position of 21 typical Fraunhofer lines in the three bands.
Table 2. The wavelength position of 21 typical Fraunhofer lines in the three bands.
BandFraunhofer Lines/nm
O2A758.812761.909764.238765.972767.178768.239770.110
WCO21596.4461598.5131600.2121606.4421612.0391613.0341619.025
SCO22045.0002054.4312056.9602058.9832063.5332070.4362074.254
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Guo, Y.; Gong, C.; Hu, Y.; Zheng, F.; Liu, Y. On-Orbit Wavelength Calibration Error Analysis of the Spaceborne Hyperspectral Greenhouse Gas Monitoring Instrument Using the Solar Fraunhofer Lines. Remote Sens. 2024, 16, 3367. https://doi.org/10.3390/rs16183367

AMA Style

Guo Y, Gong C, Hu Y, Zheng F, Liu Y. On-Orbit Wavelength Calibration Error Analysis of the Spaceborne Hyperspectral Greenhouse Gas Monitoring Instrument Using the Solar Fraunhofer Lines. Remote Sensing. 2024; 16(18):3367. https://doi.org/10.3390/rs16183367

Chicago/Turabian Style

Guo, Yulong, Cailan Gong, Yong Hu, Fuqiang Zheng, and Yunmeng Liu. 2024. "On-Orbit Wavelength Calibration Error Analysis of the Spaceborne Hyperspectral Greenhouse Gas Monitoring Instrument Using the Solar Fraunhofer Lines" Remote Sensing 16, no. 18: 3367. https://doi.org/10.3390/rs16183367

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop