LU503782B1 - Bias Correction Method and System of TanSat Satellite XCO2 Retrieval Data - Google Patents
Bias Correction Method and System of TanSat Satellite XCO2 Retrieval Data Download PDFInfo
- Publication number
- LU503782B1 LU503782B1 LU503782A LU503782A LU503782B1 LU 503782 B1 LU503782 B1 LU 503782B1 LU 503782 A LU503782 A LU 503782A LU 503782 A LU503782 A LU 503782A LU 503782 B1 LU503782 B1 LU 503782B1
- Authority
- LU
- Luxembourg
- Prior art keywords
- xco
- retrieval
- value
- albedo
- der
- Prior art date
Links
- 238000012937 correction Methods 0.000 title claims abstract description 125
- 238000000034 method Methods 0.000 title claims abstract description 34
- MYMOFIZGZYHOMD-UHFFFAOYSA-N Dioxygen Chemical compound O=O MYMOFIZGZYHOMD-UHFFFAOYSA-N 0.000 claims abstract description 158
- 229910001882 dioxygen Inorganic materials 0.000 claims abstract description 158
- 238000001228 spectrum Methods 0.000 claims abstract description 83
- 239000000443 aerosol Substances 0.000 claims description 34
- 230000008030 elimination Effects 0.000 claims description 29
- 238000003379 elimination reaction Methods 0.000 claims description 29
- 238000004364 calculation method Methods 0.000 claims description 20
- 238000012546 transfer Methods 0.000 claims description 15
- 238000012417 linear regression Methods 0.000 claims description 11
- 238000012935 Averaging Methods 0.000 claims description 6
- 230000007246 mechanism Effects 0.000 claims description 2
- 101100116570 Caenorhabditis elegans cup-2 gene Proteins 0.000 claims 1
- 101100322245 Caenorhabditis elegans des-2 gene Proteins 0.000 claims 1
- 101100116572 Drosophila melanogaster Der-1 gene Proteins 0.000 claims 1
- 230000000875 corresponding effect Effects 0.000 description 75
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 description 44
- 229910002092 carbon dioxide Inorganic materials 0.000 description 22
- 239000001569 carbon dioxide Substances 0.000 description 22
- 238000010586 diagram Methods 0.000 description 20
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 7
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 7
- 229910052799 carbon Inorganic materials 0.000 description 7
- 229910052760 oxygen Inorganic materials 0.000 description 7
- 239000001301 oxygen Substances 0.000 description 7
- 230000003287 optical effect Effects 0.000 description 6
- 230000003595 spectral effect Effects 0.000 description 6
- 230000008859 change Effects 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 230000009897 systematic effect Effects 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 238000012544 monitoring process Methods 0.000 description 4
- 238000013075 data extraction Methods 0.000 description 3
- 230000004907 flux Effects 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 230000005855 radiation Effects 0.000 description 3
- 238000004497 NIR spectroscopy Methods 0.000 description 2
- 230000000630 rising effect Effects 0.000 description 2
- 238000012216 screening Methods 0.000 description 2
- GMBQZIIUCVWOCD-WWASVFFGSA-N Sarsapogenine Chemical compound O([C@@H]1[C@@H]([C@]2(CC[C@@H]3[C@@]4(C)CC[C@H](O)C[C@H]4CC[C@H]3[C@@H]2C1)C)[C@@H]1C)[C@]11CC[C@H](C)CO1 GMBQZIIUCVWOCD-WWASVFFGSA-N 0.000 description 1
- 241000669244 Unaspis euonymi Species 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 238000000711 polarimetry Methods 0.000 description 1
- 238000001556 precipitation Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000000611 regression analysis Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000005316 response function Methods 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 238000010257 thawing Methods 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
- 238000010792 warming Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/35—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
- G01N21/3504—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light for analysing gases, e.g. multi-gas analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N2021/1793—Remote sensing
- G01N2021/1795—Atmospheric mapping of gases
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/35—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
- G01N21/3504—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light for analysing gases, e.g. multi-gas analysis
- G01N2021/3531—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light for analysing gases, e.g. multi-gas analysis without instrumental source, i.e. radiometric
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/35—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
- G01N21/359—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light using near infrared light
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2201/00—Features of devices classified in G01N21/00
- G01N2201/02—Mechanical
- G01N2201/021—Special mounting in general
- G01N2201/0214—Airborne
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2201/00—Features of devices classified in G01N21/00
- G01N2201/12—Circuits of general importance; Signal processing
- G01N2201/121—Correction signals
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Spectrometry And Color Measurement (AREA)
Abstract
The present application provides a bias correction method and system of TanSat satellite XCO2 retrieval data, the method comprises the following steps: obtaining an oxygen gas albedo error of each XCO2 retrieval value according to a full-solar-spectrum of oxygen gas band and an observation spectrum and solar zenith angle of the oxygen gas band corresponding to each XCO2 retrieval value in the TanSat satellite XCO2 retrieval data; calculating an observation angle error of each XCO2 retrieval value according to a satellite observation angle and the solar zenith angle corresponding to each XCO2 retrieval value in the N of the XCO2 retrieval values; obtaining a regression coefficient of the XCO2 retrieval data according to N of the XCO2 retrieval values and a pre-obtained XCO2 foundation observation value corresponding to the XCO2 retrieval value; performing a bias correction on each XCO2 retrieval value according to each XCO2 retrieval value, the regression coefficient, and the oxygen gas albedo error, the observation angle error of each XCO2 retrieval value to improve the accuracy of retrieval data.
Description
BL-5647
Bias Correction Method and System of TanSat Satellite XCO; Retrieval Data LUS03782
The present invention belongs to the technical field of remote sensing information, and particularly relates to a bias correction method and system of Tansat satellite XCOz retrieval data.
Background Information
Since the Industrial Revolution, with the continuous progress of various technologies, the human society has developed rapidly, and the neglect of environmental effects caused by global climate change has now become the primary issue of human interests. Global climate change is mainly manifested in precipitation redistribution, frozen soil thawing and sea level rising caused by climate warming. Severe climate change has begun to endanger ecosystem balance and gradually threaten human living environment. Wherein, rising atmospheric carbon dioxide concentration is considered to be the main cause of global climate change. In order to improve the understanding of natural and artificial surface carbon sources and carbon sinks, the atmospheric carbon dioxide monitoring technology has been developed rapidly in recent years.
The atmospheric carbon dioxide column concentration (XCO-) can be effectively monitored by detecting the near-infrared spectroscopy through remote sensing. The Total Carbon Column
Observing Network (TCCON for short) based on Fourier Transform Spectrometer (FTS for short) has high observation accuracy and can provide long time series continuous observation data, but the global distribution of observation stations is not uniform enough and the spatial scale of data is too small to reflect the spatial and temporal distribution of global atmospheric carbon dioxide concentration. The development of satellite-based remote sensing technology has effectively solved this problem, and the global XCO; spatial distribution and change trend can be obtained through satellite observation. However, the space-based distribution is not as accurate as the ground measurement. The bias of XCO; is caused by many factors, such as aerosol optical depth, thin cirrus cloud, ground pressure retrieval error and so on. These biases will lead to large errors in the estimation of regional CO; flux by retrieval analysis. Therefore, the satellite XCO» products must use the high-precision ground measurement data for bias correction to meet higher precision requirements, so as to provide reference for the international community to formulate greenhouse 1
BL-5647
LU503782 gas policies and carbon tariff mechanisms.
In the bias correction of the satellite observation data, some related studies used the linear regression to perform the bias correction of the satellite observation data. Wunch etc. correlated the difference between the XCO; retrieval result of GOSAT (a remote sensing satellite for atmospheric carbon dioxide monitoring launched in Japan in 2009) and the earth's surface albedo, the difference between the retrieval result and prior earth's surface pressure, the difference between the air mass and Oz-A band spectral radiance, and the GOSAT data from Southern
Hemisphere was used as the reference value for linear regression, and bias correction was performed based on TCCON data. After Wunch etc., Cogan etc. used analog data calculated based on the atmospheric chemical transport model GEOS-Chem to perform the bias correction on the
GOSAT XCO; data obtained from the UoL-FP retrieval algorithm. Guerlet etc. performed bias correction on the GOSAT XCO; data retrieved from the Dutch Space Research Institute/Karlsruhe
Institute of Technology based on measurements from 12 TCCON sites as a reference.
The Tansat is a global carbon dioxide scientific monitoring satellite (or carbon satellite for short), which is the first satellite in China to measure atmospheric carbon dioxide concentration from space. Emission in December 2016, satellites equipped with Atmospheric Carbon-dioxide
Grating Spectroradiometer (ACGS) and Cloud and Aerosol Polarimetry Imager (CAPI) to detect atmospheric carbon dioxide from near-infrared spectroscopy and acquired clouds and aerosol data simultaneously. The results of the XCO; retrieval and validation study for TanSat show that the bias of the uncorrected TanSat XCO; retrieval result from the TCCON data is -2.11ppm. So far, there were no bias correction XCO; data of TanSat has been published yet.
Accordingly, there is a need to provide an improved solution to the above-mentioned deficiencies of the prior art.
It is an object of the present application to provide a bias correction method and system of
Tansat satellite XCO; retrieval data to solve or alleviate the above-mentioned problems in the prior art.
In order to achieve the above object, the present invention provides the following technical solutions: 2
BL-5647
LU503782
The present application provides a bias correction method of TanSat satellite XCO» retrieval data (TanSat satellite XCO» retrieved data), comprising the following steps: obtaining an oxygen gas albedo error of each XCO» retrieval value (XCO» retrieved value) according to a pre-obtained full-solar-spectrum of an oxygen gas band and a pre-obtained observation spectrum and solar zenith angle of the oxygen gas band corresponding to each XCO; retrieval value in the TanSat satellite XCO; retrieval data; wherein the XCO» retrieval data contains N of the XCO- retrieval values, and N is a positive integer; calculating an observation angle error of each XCO; retrieval value according to a satellite observation angle and the solar zenith angle corresponding to each
XCO; retrieval value in the N of the XCO; retrieval values; obtaining a regression coefficient of the XCO; retrieval data according to N of the XCO; retrieval values and a pre-obtamed XCO, foundation observation value corresponding to the XCO; retrieval value; performing a bias correction on each XCO; retrieval value according to each XCO; retrieval value, the regression coefficient of the XCO; retrieval value, and the oxygen gas albedo error and the observation angle error of each XCO; retrieval value to obtain a bias correction result of each XCO; retrieval value.
Alternatively, in any example of the present application, the step of obtaining an oxygen gas albedo error of each XCO; retrieval value according to a pre-obtained full-solar-spectrum of an oxygen gas band and a pre-obtained observation spectrum and solar zenith angle of the oxygen gas band corresponding to each XCO; retrieval value in the TanSat satellite XCO; retrieval data comprises the following steps: calculating an oxygen gas band non-absorbing region albedo corresponding to each XCO» retrieval value according to the pre-obtained full-solar-spectrum of the oxygen gas band and a pre-obtained observation spectrum and solar zenith angle of the oxygen gas band corresponding to each XCO; retrieval value in the N of the XCO; retrieval values; averaging the oxygen gas band non-absorbing region albedo corresponding to N of the XCO, retrieval values to obtain an oxygen gas albedo average value of N of the XCO; retrieval values; calculating the oxygen gas albedo error of each XCO; retrieval value according to the oxygen gas band non-absorbing region albedo corresponding to each XCO» retrieval value, and the oxygen gas albedo average value.
Alternatively, in any example of the present application, the step of calculating an observation angle error of each XCO» retrieval value according to a satellite observation angle and the solar zenith angle corresponding to each XCO; retrieval value in the N of the XCO; retrieval values 3
BL-5647
LU503782 comprises the following steps: calculating an air mass factor of each XCO» retrieval value according to the satellite observation angle and the solar zenith angle corresponding to each XCO, retrieval value in the N of the XCO; retrieval values; averaging the air mass factor of N of the
X CO; retrieval values to obtain an air mass factor average value of N of the XCO; retrieval values; calculating the observation angle error of each XCO» retrieval value according to the air mass factor corresponding to each XCO; retrieval value, and the air mass factor average value.
Alternatively, in any example of the present application, the step of obtaining a regression coefficient of the XCO» retrieval data according to N of the XCO; retrieval values and a pre-obtained XCO; foundation observation value corresponding to the XCO; retrieval value is specifically: performing linear regression on N of the XCO; retrieval values and a pre-obtained
XCO; foundation observation value corresponding to the XCO» retrieval value to obtain the regression coefficient of the XCO; retrieval data.
Alternatively, in any example of the present application, the step of performing a bias correction on each XCO; retrieval value according to each XCO; retrieval value, the regression coefficient of the XCO» retrieval value, and the oxygen gas albedo error and the observation angle error of each
X CO; retrieval value to obtain a bias correction result of each XCO; retrieval value comprises the following steps: performing an albedo error correction and an observation angle error correction on each XCO; retrieval value respectively according to each XCO» retrieval value, and the oxygen gas albedo error and the observation angle error of each XCO; retrieval value; performing a quotient operation on each XCO» retrieval value after the albedo error correction and the observation angle error correction and the regression coefficient to obtain the bias correction result of each XCO; retrieval value.
Alternatively, in any example of the present application, before the step of obtaining an oxygen gas albedo error of each XCO» retrieval value according to a pre-obtained full-solar-spectrum of an oxygen gas band and a pre-obtained observation spectrum and solar zenith angle of the oxygen gas band corresponding to each XCO» retrieval value in the TanSat satellite XCO» retrieval data, further comprising the following steps: extracting a satellite parameter and an observation parameter of the TanSat satellite; performing cloud region data elimination on cloud data, and performing aerosol high value region data elimination on aerosol data; obtaining an analog spectrum based on an atmospheric radiative transfer model of a preset retrieval parameter, 4
BL-5647
LU503782 according to the satellite parameter, the observation parameter, a cloud data performed cloud region data elimination, an aerosol data performed aerosol high value region data elimination and an atmospheric profile parameter; solving the analog spectrum and the observation spectrum in the satellite parameter to obtain a solution result based on an optimal estimation method; judging whether the solution result converges, and 1f not, updating the preset retrieval parameter in the atmospheric radiative transfer model until the solution result converges; obtaining the TanSat satellite XCO» retrieval data according to the solution result.
The present application also provides a bias correction system of TanSat satellite XCO» retrieval data, comprising: an albedo error unit configured to obtain an oxygen gas albedo error of each XCO; retrieval value according to a pre-obtained full-solar-spectrum of an oxygen gas band and a pre-obtained observation spectrum and solar zenith angle of the oxygen gas band corresponding to each XCO; retrieval value in the TanSat satellite XCO» retrieval data; wherein the XCO» retrieval data contains N of the XCO; retrieval values, and N is a positive integer; an observation angle error unit configured to calculate an observation angle error of each XCO, retrieval value according to a satellite observation angle and the solar zenith angle corresponding to each XCO; retrieval value in the N of the XCO; retrieval values; a regression coefficient unit configured to obtain a regression coefficient of the XCO; retrieval data according to N of the
XCO; retrieval values and a pre-obtained XCO; foundation observation value corresponding to the
XCO- retrieval value; a bias correction unit configured to perform a bias correction on each XCO, retrieval value according to each XCO; retrieval value, the regression coefficient of the XCO, retrieval value, and the oxygen gas albedo error and the observation angle error of each XCO, retrieval value to obtain a bias correction result of each XCO; retrieval value.
Alternatively, in any example of the present application, the albedo error unit comprises: an albedo calculation subunit configured to calculate an oxygen gas band non-absorbing region albedo corresponding to each XCO; retrieval value according to the pre-obtained full-solar-spectrum of the oxygen gas band and a pre-obtained observation spectrum and solar zenith angle of the oxygen gas band corresponding to each XCO; retrieval value in the N of the
XCO; retrieval values; an albedo average calculation subunit configured to average the oxygen gas band non-absorbing region albedo corresponding to N of the XCO; retrieval values to obtain an oxygen gas albedo average value of N of the XCO; retrieval values; an albedo error calculation
BL-5647
LU503782 subunit configured to calculate the oxygen gas albedo error of each XCO; retrieval value according to the oxygen gas band non-absorbing region albedo corresponding to each XCO», retrieval value, and the oxygen gas albedo average value.
Alternatively, in any example of the present application, the observation angle error unit comprises: an air mass factor calculation subunit configured to calculate an air mass factor of each
XCO; retrieval value according to the satellite observation angle and the solar zenith angle corresponding to each XCO; retrieval value in the N of the XCO retrieval values; an air mass factor average calculation subunit configured to average the air mass factor of N of the XCO, retrieval values to obtain an air mass factor average value of N of the XCO; retrieval values; an observation angle error calculation subunit configured to calculate the observation angle error of each XCO» retrieval value according to the air mass factor corresponding to each XCO-, retrieval value, and the air mass factor average value.
Alternatively, in any example of the present application, the regression coefficient unit is further configured to perform linear regression on N of the XCO» retrieval values and a pre-obtained
XCO; foundation observation value corresponding to the XCO» retrieval value to obtain the regression coefficient of the XCO; retrieval data.
Alternatively, in any example of the present application, the bias correction unit comprises: a first correction subunit configured to perform an albedo error correction and an observation angle error correction on each XCO-, retrieval value respectively according to each XCO; retrieval value, and the oxygen gas albedo error and the observation angle error of each XCO; retrieval value; a second correction subunit configured to perform a quotient operation on each XCO; retrieval value after the albedo error correction and the observation angle error correction and the regression coefficient to obtain the bias correction result of each XCO; retrieval value.
Alternatively, in any of the example of the present application, further comprises: a data extraction unit configured to extract a satellite parameter and an observation parameter of the
TanSat satellite; a data elimination unit configured to perform cloud region data elimination on cloud data, and performing aerosol high value region data elimination on aerosol data; an analog spectrum unit configured to obtain an analog spectrum based on an atmospheric radiative transfer model of a preset retrieval parameter, according to the satellite parameter, the observation parameter, a cloud data performed cloud region data elimination, an aerosol data performed 6
BL-5647
LU503782 aerosol high value region data elimination and an atmospheric profile parameter; a spectrum solving unit configured to solve the analog spectrum and the observation spectrum in the satellite parameter to obtain a solution result based on an optimal estimation method; a result update unit configured to judge whether the solution result converges, and if not, updating the preset retrieval parameter in the atmospheric radiative transfer model until the solution result converges; a retrieval data unit configured to obtain the TanSat satellite XCO; retrieval data according to the solution result.
Compared with the closest prior art, the technical solutions of the examples of the present application have the following beneficial effects:
In the technical solution provided in the example of the present application, obtaining an oxygen gas albedo error of each XCO; retrieval value according to a pre-obtained full-solar-spectrum of an oxygen gas band and a pre-obtained observation spectrum and solar zenith angle of the oxygen gas band corresponding to each XCO; retrieval value in the TanSat satellite XCO- retrieval data; calculating an observation angle error of each XCO; retrieval value according to a satellite observation angle and the solar zenith angle corresponding to each XCO, retrieval value; obtaining a regression coefficient of the XCO; retrieval data according to N of the
XCO; retrieval values and a pre-obtained XCO; foundation observation value corresponding to the
XCO- retrieval value; then, performing a bias correction on each XCO» retrieval value according to each XCO; retrieval value, the regression coefficient of each XCO; retrieval value, and the oxygen gas albedo error, the observation angle error of each XCO; retrieval value to obtain a bias correction result of each XCO; retrieval value. Therefore, the oxygen gas albedo error, the observation angle error and the system error of TanSat satellite observation in XCO; retrieval value are eliminated, and the XCO> retrieval data correction is finally realized, which effectively improves the accuracy of TanSat satellite observation.
The accompanying drawings, which are incorporated in and constitute a part of this application, are included to provide a further understanding of the application, and the description of the exemplary examples and illustrations of the application are intended to explain the application and are not intended to limit the application. Wherein: 7
BL-5647
LU503782
Fig. 1 is a schematic flow diagram of the bias correction method of TanSat satellite XCO» retrieval data provided according to some examples of the present application;
Fig. 2 1s a schematic flow diagram of step S101 in the bias correction method of TanSat satellite
XCO; retrieval data provided according to some examples of the present application;
Fig. 3 1s a schematic flow diagram of step S102 the bias correction method of TanSat satellite
XCO; retrieval data provided according to some examples of the present application;
Fig. 4 is a schematic flow diagram of step S104 in the bias correction method of TanSat satellite
XCO; retrieval data provided according to some examples of the present application;
Fig. 5 is a schematic flow diagram for obtaining TanSat satellite XCO> retrieval data provided according to some examples of the present application;
Fig. 6 is a schematic structure diagram of the bias correction system of TanSat satellite XCOz retrieval data provided according to some examples of the present application;
Fig. 7 is a schematic structure diagram of the albedo error unit according to some examples of the present application;
Fig. 8 is a schematic structure diagram of the observation angle error unit according to some examples of the present application;
Fig. 9 is a schematic structure diagram of the bias correction unit according to some examples of the present application;
Fig. 10 is a schematic structure diagram for obtaining TanSat satellite XCOz retrieval data provided according to some examples of the present application.
The present application will be described in detail below with reference to the accompanying drawings and in conjunction with the examples. The various examples are provided by way of interpretation of the application and not limiting the application. Indeed it will be apparent to those skilled in the art that modifications and variations may be made in the present application without departing from the scope or spirit of the application. For example features shown or described as part of one example may be used in another example to produce yet another example. It is therefore desirable that the application encompass such modifications and variations falling within the scope of the appended claims and their equivalents. 8
BL-5647
LU503782
Exemplary method
Fig. 1 is a schematic flow diagram of the bias correction method of TanSat satellite XCO» retrieval data provided according to some examples of the present application; as shown in Fig. 1, the method comprises the following steps:
Step S101, obtaining an oxygen gas albedo error of each XCO; retrieval value according to a pre-obtained full-solar-spectrum of an oxygen gas band and a pre-obtained observation spectrum and solar zenith angle of the oxygen gas band corresponding to each XCO; retrieval value in the
TanSat satellite XCO» retrieval data; wherein the XCO; retrieval data contains N of the XCO; retrieval values, and N is a positive integer;
In the example of the present application, the full-solar-spectrum of oxygen gas band can be selected according to needs, for example, selecting the full-solar-spectrum data published by Dr.
Robert L. Kurucz laboratory; the observation spectrum of the oxygen gas band corresponding to
XCO» retrieval value is observed by the Atmospheric Carbon-dioxide Grating Spectroradiometer (ACGS for short) mounted in the TanSat satellite at the observation point. The solar zenith angle can be obtained from the LIB data (carbon satellite primary calibration data) of TanSat satellite. It is to be understood that the above description is intended to be illustrative, and that the present examples are not restrictive.
Fig. 2 is a schematic flow diagram of step S101 in the bias correction method of TanSat satellite
XCO- retrieval data provided according to some examples of the present application; as shown in
Fig. 2, the step of obtaining an oxygen gas albedo error of each XCO; retrieval value according to a pre-obtained full-solar-spectrum of an oxygen gas band and a pre-obtained observation spectrum and solar zenith angle of the oxygen gas band corresponding to each XCO; retrieval value in the
TanSat satellite XCO-, retrieval data comprises the following steps:
Step S111, calculating an oxygen gas band non-absorbing region albedo corresponding to each
XCO» retrieval value according to the pre-obtained full-solar-spectrum of the oxygen gas band and a pre-obtained observation spectrum and solar zenith angle of the oxygen gas band corresponding to each XCO; retrieval value in the N of the XCO; retrieval values;
In the example of the present application, the albedo of the oxygen gas band non-absorbing region is the ratio of the irradiance obtained without absorption of other molecules in the oxygen gas band and the irradiance corresponding solar spectrum in the process of sunlight being reflected 9
BL-5647
LU503782 back to the sensor mounted on the TanSat satellite through atmospheric radiation and through the earth's surface. It is to be understood that the above description is intended to be illustrative, and that the present examples are not restrictive.
In the example of the present application, when calculating an oxygen gas band non-absorbing region albedo corresponding to each XCO retrieval value according to the pre-obtained full-solar-spectrum of the oxygen gas band and a pre-obtained observation spectrum and solar zenith angle of the oxygen gas band corresponding to each XCO; retrieval value, the oxygen gas band non-absorbing region albedo corresponding to each XCO; retrieval value is calculated according to the irradiance of the full-solar-spectrum of the oxygen gas band, the radiance of the observation spectrum of the oxygen gas band corresponding to each XCO; retrieval value and the solar zenith angle. Specifically, the oxygen gas band non-absorbing region albedo corresponding to XCO» retrieval value is calculated by the following formula (1). The formula (1) is as follows: albedo_0, = Tr RT D
Wherein, the albedo_0, is the oxygen gas band non-absorbing region albedo corresponding to the XCO; retrieval value, the radiance_O2 reasured is the radiance of the observation spectrum of the oxygen gas band corresponding to the XCO; retrieval value, and indicates the radiant flux of per unit of spectral wavelength, per unit of projected area, per unit of solid angle in w/m°/sr/nm; the irradiance O,_ is the irradiance of the full-solar-spectrum of the oxygen gas band, and indicates the radiant flux received per unit area at unit spectral wavelength in w/m?*/nm; the 6, is the solar zenith angle corresponding to XCO» retrieval value; the m is the circumference ratio (the value is 3.14). The oxygen gas albedo may be a positive number, e.g. may range from (0, 1).
It is to be understood that the above description is intended to be illustrative, and that the present examples are not restrictive.
Step S121, averaging the oxygen gas band non-absorbing region albedo corresponding to N of the XCO; retrieval values to obtain an oxygen gas albedo average value of N of the XCO; retrieval values;
In the example of the present application, the XCO; retrieval data contains N of the XCO» retrieval values, and each oxygen gas band non-absorbing region albedo is calculated corresponding to each XCO; retrieval value, and thus N of the oxygen gas band non-absorbing
BL-5647
LU503782 region albedos are obtained. By performing a weighted average calculation on the obtained N of the oxygen gas band non-absorbing region albedos, the oxygen gas albedo average value of N of the XCO; retrieval values can be obtained. Wherein, the albedo_O, indicates the oxygen gas albedo average value of N of the XCO; retrieval values. It is to be understood that the above description is intended to be illustrative, and that the present examples are not restrictive.
Step S131, calculating the oxygen gas albedo error of each XCO; retrieval value according to the oxygen gas band non-absorbing region albedo corresponding to each XCO; retrieval value, and the oxygen gas albedo average value;
In the example of the present application, when calculating the oxygen gas albedo error of each
XCO; retrieval value according to the oxygen gas band non-absorbing region albedo corresponding to each XCO; retrieval value, and the oxygen gas albedo average value, the oxygen gas albedo error of each XCO» retrieval value is calculated according to the oxygen gas band non-absorbing region albedo and oxygen gas albedo average value corresponding to each XCO, retrieval value based on a preset oxygen gas albedo correction coefficient. It is to be understood that the above description is intended to be illustrative, and that the present examples are not restrictive.
In the example of the present application, the C , indicates the oxygen gas albedo correction coefficient. The product of the difference between oxygen gas band non-absorbing region albedo and the oxygen gas albedo average value corresponding to each XCO» retrieval value and the preset oxygen gas albedo correction coefficient is the oxygen gas albedo error of each XCO; retrieval value, i.e. the oxygen gas albedo error of each XCO» retrieval value is C 2(albedo_0, — albedo_0,), wherein, the C , indicates the oxygen gas albedo correction coefficient. It is to be understood that the above description is intended to be illustrative, and that the present examples are not restrictive.
In the example of the present application, the albedo error of each XCO; retrieval value can be calculated according to the oxygen gas band non-absorbing region albedo corresponding to each
XCO- retrieval value, and the oxygen gas albedo average value, and then the oxygen gas albedo correction coefficient can be obtained by performing the regression analysis on N of the oxygen gas band non-absorbing region albedo and the albedo bias. It is to be understood that the above description is intended to be illustrative, and that the present examples are not restrictive. 11
BL-5647
LU503782
Step S102, calculating an observation angle error of each XCO; retrieval value according to a satellite observation angle and the solar zenith angle corresponding to each XCO; retrieval value in the N of the XCO; retrieval values;
In the example of the present application, in the process of forward atmospheric radiation transmission, the XCO; retrieval is an estimation of the optical path length, and in this process, influenced by clouds and aerosols and other factors, the optical path length would be wrongly estimated, and the size of the satellite observation angle and the solar zenith angle would further affect the bias of the XCO; retrieval result due to the optical path length; therefore, it is necessary to correct the XCO; retrieval result by counting the air mass factor of XCO; retrieval value calculated by the satellite observation angle and the solar zenith angle, namely, the XCO; retrieval value is required to correct by the observation angle error of XCO; retrieval value. It is to be understood that the above description is intended to be illustrative, and that the present examples are not restrictive.
Fig. 3 is a schematic flow diagram of step S102 the bias correction method of TanSat satellite
XCO- retrieval data provided according to some examples of the present application; as shown in
Fig. 3, the step of calculating an observation angle error of each XCO; retrieval value according to a satellite observation angle and the solar zenith angle corresponding to each XCO; retrieval value in the N of the XCO; retrieval values comprises the following steps:
Step S112, calculating an air mass factor of each XCO; retrieval value according to the satellite observation angle and the solar zenith angle corresponding to each XCO; retrieval value in the N of the XCO; retrieval values;
In the example of the present application, the air mass factor of XCO; retrieval value is calculated by the following formula (2). The formula (2) is as follows: airmass = om Tagen (2)
Wherein, the airmass is the air mass factor of the XCO» retrieval value, the 67 is the solar zenith angle corresponding to XCO; retrieval value, and the 8 is the satellite observation angle of the TanSat satellite. It 1s to be understood that the above description is intended to be illustrative, and that the present examples are not restrictive.
Step S122, averaging the air mass factor of N of the XCO; retrieval values to obtain an air mass 12
BL-5647
LU503782 factor average value of N of the XCO» retrieval values;
In the example of the present application, the XCO; retrieval data contains N of the XCO» retrieval value, and each air mass factor corresponding to each XCO» retrieval value is calculated, and thus N of the air mass factors are obtained. By performing a weighted average calculation on the obtained N of the air mass factors, the air mass factor average value of N of the XCO; retrieval values can be obtained. Wherein, the airmass indicates the air mass factor average value of N of the XCO; retrieval values. It is to be understood that the above description is intended to be illustrative, and that the present examples are not restrictive.
Step S132, calculating the observation angle error of each XCO; retrieval value according to the air mass factor corresponding to each XCO; retrieval value, and the air mass factor average value;
In the present example, when calculating the observation angle error of each XCO; retrieval value according to the air mass factor corresponding to each XCO; retrieval value, and the air mass factor average value, the observation angle error for each XCO» retrieval value is calculated according to the air mass factor and the air mass factor average for each XCO; retrieval value based on the preset air mass factor correction coefficient. It is to be understood that the above description is intended to be illustrative, and that the present examples are not restrictive.
In the example of the present application, the C ;indicates the air mass factor correction coefficient, which refers to the atmospheric carbon dioxide column error resulting in C , parts per million (ppm) concentration unit per caused by 1 atmospheric optical mass (airmass) unit difference. It is to be understood that the above description is intended to be illustrative, and that the present examples are not restrictive.
In the example of the present application, the product of the difference between the air mass factor corresponding to each XCO; retrieval value and the air mass factor average and the preset air mass factor correction coefficient is the observation angle error of each XCO; retrieval value, i.e. the observation angle error of each XCO; retrieval value is C , (airmass — arrmass). It is to be understood that the above description is intended to be illustrative, and that the present examples are not restrictive.
In the example of the present application, the air mass factor correction coefficient (C 4) is obtained according to actual satellite observation angle and error statistics, and the value of C ; is less than or equal to 5 ppm/airmass. For example, C ; is empirically set to 2 ppm/airmass. It is to 13
BL-5647
LU503782 be understood that the above description 1s intended to be illustrative, and that the present examples are not restrictive.
Step S103, obtaining a regression coefficient of the XCO; retrieval data according to N of the
XCO; retrieval values and a pre-obtained XCO; foundation observation value corresponding to the
X CO; retrieval value;
In the example of the present application, the XCO; foundation observation value is obtained by observing the observation point with the Fourier Transform Spectrometer, which corresponds to
TanSat satellite XCO; retrieval value in one-to-one correspondence. Wherein, the XCO, foundation observation value can be defined as ground-based observation data, and TanSat satellite XCO; retrieval value can be defined as space-based observation data. It is to be understood that the above description is intended to be illustrative, and that the present examples are not restrictive.
In the example of the present application, compared with the observation of an observation point using the Fourier Transform Spectrometer, when using a satellite to perform XCO, observation, due to the limitation of satellite-based remote sensing technology and the uncertainty of observation, the ground-based observation data 1s closer to the "real atmosphere”, and there 1s still some systematic bias in the space-based observation data. Therefore, the ground-based observation data can be used to correct the space-based observation data, and a more realistic
XCO- retrieval data can be obtained. It is to be understood that the above description is intended to be illustrative, and that the present examples are not restrictive.
In some alternative examples, when obtaining a regression coefficient of the XCO; retrieval data according to N of the XCO; retrieval values and a pre-obtained XCO; foundation observation value corresponding to the XCO; retrieval value, the linear regression is performed on N of the
XCO; retrieval values and a pre-obtained XCO; foundation observation value corresponding to the
XCO» retrieval value to obtain the regression coefficient of the XCO; retrieval data. It is to be understood that the above description is intended to be illustrative, and that the present examples are not restrictive.
In the example of the present application, the ground-based observation data is closer to the "real atmosphere”, and there are still some systematic bias in the space-based observation data.
Therefore, using the pre-obtained XCO; foundation observation value (namely, ground-based 14
BL-5647
LU503782 observation data) and N of the XCO; retrieval value of TanSat satellite (namely, space-based observation data) to perform the linear regression, the regression coefficient of space-based observation data relatrve to ground-based observation data is obtained, so as to realize the bias correction of space-based observation data by the ground-based observation data and obtain the more realistic XCO» retrieval data. It is to be understood that the above description is intended to be illustrative, and that the present examples are not restrictive.
In the example of the present application, the C 9 indicates the regression coefficient. The unary linear regression equation of N of the XCO-, retrieval values of TanSat satellite with respect to the
XCO; foundation observation value is obtained according to the pre-obtained XCO; foundation observation value and N of the XCO; retrieval values. Assuming that the foundation observation of the TanSat satellite and the Fourier Transform Spectrometer are synchronized at "zero point",
Le. when the foundation observation of the Fourier Transform Spectrometer is zero, the observation of the Atmospheric Carbon-dioxide Grating Spectroradiometer mounted on the TanSat satellite 1s also zero, so that the intercept of the obtained unary linear regression equation is zero, the regression coefficient C og can be obtained. It is to be understood that the above description is intended to be illustrative, and that the present examples are not restrictive.
Step S104, performing a bias correction on each XCO» retrieval value according to each XCO, retrieval value, the regression coefficient of the XCO; retrieval value, and the oxygen gas albedo error and the observation angle error of each XCO; retrieval value to obtain a bias correction result of each XCO; retrieval value.
In the example of the present application, when obtaining the XCO- retrieval value, the retrieval parameters are used, wherein the retrieval parameters comprise the observation mode, the scattering mode and the trace gas profile, and may specifically include the earth's surface pressure, the albedo, the carbon dioxide and oxygen gas profile, the observation angle, etc. Because the retrieval parameters will lead to uncertainty in the retrieval result, including the earth's surface pressure, the albedo, the air mass parameter and the aerosol, etc. In the retrieval process, the earth's surface pressure is not retrieved synchronously and aerosol parameter is only used as the threshold for low concentration region screening. Therefore, when correcting the XCO; retrieval data, the oxygen gas albedo error and the observation angle error need to be corrected. It is to be understood that the above description is intended to be illustrative, and that the present examples
BL-5647 …. LU503782 are not restrictive.
In the example of the present application, the correction of XCO; retrieval data is achieved by eliminating the albedo error, the observation angle error, and the systematic bias at the TanSat satellite observation in the XCO; retrieval data. It is to be understood that the above description is intended to be illustrative, and that the present examples are not restrictive.
In the example of the present application, the oxygen gas albedo error in the XCO; retrieval data is eliminated by the oxygen gas albedo correction coefficient; the observation angle error in the XCO; retrieval data is eliminated by the air mass factor correction coefficient; the systematic bias at the TanSat satellite observation is eliminated by the regression coefficient, and the XCO; retrieval data correction is finally achieved. It is to be understood that the above description is intended to be illustrative, and that the present examples are not restrictive.
Fig. 4 is a schematic flow diagram of step S104 in the bias correction method of TanSat satellite
XCO- retrieval data provided according to some examples of the present application; as shown in
Fig. 4, the step of performing a bias correction on each XCO; retrieval value according to each
XCO» retrieval value, the regression coefficient of the XCO; retrieval value, and the oxygen gas albedo error and the observation angle error of each XCO; retrieval value to obtain a bias correction result of each XCO» retrieval value comprises the following steps:
Step S114, performing an albedo error correction and an observation angle error correction on each XCO; retrieval value respectively according to each XCO; retrieval value, and the oxygen gas albedo error and the observation angle error of each XCO» retrieval value;
In the example of the present application, when performing an albedo error correction and an observation angle error correction on each XCO; retrieval value respectively according to each
XCO» retrieval value, and the oxygen gas albedo error and the observation angle error of each
XCO» retrieval value, Performing the difference operation on each XCO; retrieval value and the oxygen gas albedo error and the observation angle error of each XCO» retrieval value can complete the albedo error correction and the observation angle error correction of each XCO; retrieval value. It is to be understood that the above description is intended to be illustrative, and that the present examples are not restrictive.
In the example of the present application, the Xrgirieved indicates the XCO-, retrieval value,
Xrgirieved — C ,(airmass — airmass) — C 2(albedo_0, — albedo_0,) indicates that the 16
BL-5647
LU503782 oxygen gas albedo error and the observation angle error of the XCO» retrieval value is eliminated and the albedo error correction and the observation angle error correction of the XCO; retrieval value is completed. It is to be understood that the above description is intended to be illustrative, and that the present examples are not restrictive.
Step S124, performing a quotient operation on each XCO; retrieval value after the albedo error correction and the observation angle error correction and the regression coefficient to obtain the bias correction result of each XCO; retrieval value;
In the present example, the bias correction result of each XCO; retrieval value is calculated by the following formula (3). The formula (3) is as follows: xegreeted — K662 TC x (érmass atrmass) ~C y(atbedo 0, albedo 0) (q)
Co
Wherein, the Xggrrected is the bias correction result of XCOz retrieval value; the Xrgtrievedis the XCO; retrieval value; the C , is the air mass factor correction coefficient, and the value is 2 ppm/airmass; the airmass is the air mass factor of the XCO; retrieval value; the aırmass is the air mass factor average value of N of the XCO; retrieval values; the C , (airmass — airmass) is the observation angle error of the XCO; retrieval value; the C „is the oxygen gas albedo correction coefficient, the value is 10ppm/unit of albedo (column concentration error per unit albedo); the albedo_0, is the oxygen gas band non-absorbing region albedo corresponding to the XCO retrieval value; the albedo_0, is the oxygen gas albedo average value of N of the
XCO» retrieval values; the C 2(albedo_0, — albedo_0,) is the oxygen gas albedo error of the
XCO» retrieval value; the Co is the regression coefficient. It is to be understood that the above description is intended to be illustrative, and that the present examples are not restrictive.
In the example of the present application, the correction of XCO; retrieval data is achieved by climinating the oxygen gas albedo error, the observation angle error and the system error at TanSat satellite observation in XCO; retrieval value. Calculating the oxygen gas albedo of each XCO, retrieval value according to a pre-obtained full-solar-spectrum of an oxygen gas band and a pre-obtained observation spectrum and solar zenith angle of the oxygen gas band corresponding to each XCO» retrieval value in the TanSat satellite XCO; retrieval data, and the oxygen gas albedo error in the XCO; retrieval value is eliminated by the oxygen gas albedo correction coefficient; calculating an air mass factor of each XCO; retrieval value and an air mass factor average value of 17
BL-5647
LU503782
N of the XCO; retrieval values according to the satellite observation angle and the solar zenith angle corresponding to each XCO-, retrieval value in the N of the XCO; retrieval values, and the observation angle error in the XCO» retrieval value is eliminated by the air mass factor correction coefficient (namely, the air mass factor correction); obtaining a regression coefficient of the XCO, retrieval data according to N of the XCO; retrieval values and a pre-obtained XCO» foundation observation value corresponding to the XCO; retrieval value, and the system bias at the TanSat satellite observation is eliminated by the regression coefficient (namely, the scale bias correction), and the correction of the XCO» retrieval data is finally realized, which effectively improves the accuracy of TanSat satellite observation. It is to be understood that the above description is intended to be illustrative, and that the present examples are not restrictive.
In an application scenario, the Atmospheric Carbon-dioxide Grating Spectroradiometer (including 9 sensors) is installed on the TanSat satellite. Each observation of 9 sensors will simultaneously record the spectrum and measure the reflected sunlight from the earth's surface.
Different sensors have different spectral radiation responses. Therefore, in the XCO; retrieval value processing, the different sensors are performed spectral calibration separately. At the retrieved, a separate instrument response function is used to retrieve the XCO; data of TanSat satellite, and finally the XCO; retrieval data of TanSat satellite is obtained. It is to be understood that the above description is intended to be illustrative, and that the present examples are not restrictive.
Fig. 5 is a schematic flow diagram for obtaining TanSat satellite XCO> retrieval data provided according to some examples of the present application; as shown in Fig. 5, before the step of obtaining an oxygen gas albedo error of each XCO; retrieval value according to a pre-obtained full-solar-spectrum of an oxygen gas band and a pre-obtained observation spectrum and solar zenith angle of the oxygen gas band corresponding to each XCO; retrieval value in the TanSat satellite XCO» retrieval data, further comprising the following steps:
Step 501, extracting a satellite parameter and an observation parameter of the TanSat satellite;
In the example of the present application, the satellite parameter of the TanSat satellite comprises the observation angle, the observation time, the observation instrument, the spectral signal-to-noise ratio and the observation spectrum; the observation parameter comprises the longitude and latitude, the earth's surface albedo, the earth's surface elevation, the observation 18
BL-5647
LU503782 quality mark and the land and water area mark of the observation point. It is to be understood that the above description is intended to be illustrative, and that the present examples are not restrictive.
In the example of the present application, the XCO; data is retrieved using LIB data (carbon satellite primary calibration data) of the satellite, wherein the L1B data is in a hdf (Hierarchical
Data Format) 5 format, and thus the field corresponding to the satellite parameter and the observation parameter can be directly extracted from the L1B data as the corresponding satellite parameter and observation parameter, respectively. Thereby, large-scale monitoring of observation point can be realized and the data extraction can be performed quickly. It is to be understood that the above description is intended to be illustrative, and that the present examples are not restrictive.
Step 502, performing cloud region data elimination on cloud data, and performing aerosol high value region data elimination on aerosol data;
In the example of the present application, the cloud data is derived from FY-4A data (the observation data of Fengyun No. 4 satellite) matched in time and space, and the aerosol data is derived from FY-3C data (the observation data of Fengyun No. 3 satellite). It is to be understood that the above description is intended to be illustrative, and that the present examples are not restrictive.
In the example of the present application, when performing cloud region data elimination on cloud data, the cloud filtering and cloud screening are performed using a cloud detection product.
When performing aerosol high value region data elimination on aerosol data, the data with
Aerosol Optical Depth (AOD for short) value greater than 0.3 in the aerosol data is eliminated. It is to be understood that the above description is intended to be illustrative, and that the present examples are not restrictive.
Step S503, obtaining an analog spectrum based on an atmospheric radiative transfer model of a preset retrieval parameter, according to the satellite parameter, the observation parameter, a cloud data performed cloud region data elimination, an aerosol data performed acrosol high value region data elimination and an atmospheric profile parameter;
In an example of the present application, the retrieval parameters comprise the observation mode, the scattering mode, and may specifically include the earth's surface pressure, the albedo, 19
BL-5647
LU503782 the air mass parameter, the observation angle, etc., the atmospheric profile parameters comprise the carbon dioxide and oxygen profile, the temperature profiles, the potential profiles the humidity profile, and the pressure profile. It is to be understood that the above description is intended to be illustrative, and that the present examples are not restrictive.
In the example of the present application, the observation spectrum is observed by the
Atmospheric Carbon-dioxide Grating Spectroradiometer, and the analog spectrum is obtained based on atmospheric radiative transfer model, whereby the observation spectrum can be compared using the analog spectrum. It is to be understood that the above description is intended to be illustrative, and that the present examples are not restrictive.
Step S504, solving the analog spectrum and the observation spectrum in the satellite parameter to obtain a solution result based on an optimal estimation method;
In the example of the present application, the forward analog is performed based on the atmospheric radiative transfer model to obtain the analog spectrum, the analog value of analog spectrum is compared with the observation spectrum to generate the cost function, and an optimal solution of the retrieved atmospheric carbon dioxide and oxygen profile is obtained by minimizing the cost function, namely, the solution result (the number of carbon dioxide molecules and the number of oxygen molecules in the atmosphere). It is to be understood that the above description 1s intended to be illustrative, and that the present examples are not restrictive.
Step S505, judging whether the solution result converges, and if not, updating the preset retrieval parameter in the atmospheric radiative transfer model until the solution result converges;
In the example of the present application, the solution result does not converge that indicates the error of the retrieved atmospheric carbon dioxide and oxygen profile is also relatively large.
Therefore, the preset retrieval parameter of the atmospheric radiative transfer model needs to be updated to converge the solution result, that is, the cost function of the analog value of the analog spectrum and the retrieval value of the observation spectrum is minimal. It is to be understood that the above description is intended to be illustrative, and that the present examples are not restrictive.
Step S506, obtaining the TanSat satellite XCO: retrieval data according to the solution result.
In the example of the present application, the solution result converges (i.e. the cost function between the simulated value of analog spectrum and the inverted value of observation spectrum is
BL-5647
LU503782 minimal) indicates that atmospheric carbon dioxide and the oxygen profiles are obtamed by retrieved, and then the XCO; retrieval value can be calculated. It is to be understood that the above description is intended to be illustrative, and that the present examples are not restrictive.
In the example of the present application, when obtaining the TanSat satellite XCO; retrieval data according to the solution result, the TanSat satellite XCO» retrieval value is obtained according to the number of carbon dioxide molecules, the number of oxygen molecules and the percentage of oxygen gas in the atmosphere. Specifically, XCO; retrieval value can be obtained by the following formula (4).
Xrgirieved = I) / on
Wherein, the Xrgirieved is TanSat satellite XCOz retrieval value, the COS°! is the number of the carbon dioxide molecules (i.e. the solution result), the 05°" is the number of oxygen gas molecules,the 0” is the percentage of oxygen in the atmosphere (the value is 20.95%). It is to be understood that the above description is intended to be illustrative, and that the present examples are not restrictive.
In the example of the present application, the oxygen gas albedo error in the XCO; retrieval data is eliminated by the oxygen gas albedo correction coefficient; the observation angle error in the XCO:; retrieval data is eliminated by the air mass factor correction coefficient (namely, air mass factor correction); the systematic error (i.e. scale programming correction) at the TanSat satellite observation is eliminated by the regression coefficient, and the XCO; retrieval data correction is finally achieved; the accuracy of TanSat satellite observation data is improved effectively. It is to be understood that the above description is intended to be illustrative, and that the present examples are not restrictive.
Exemplary system
Fig. 6 1s a schematic structure diagram of the bias correction system of TanSat satellite XCO» retrieval data provided according to some examples of the present application; as shown in Fig. 6, the system comprises: an albedo error unit configured to obtain an oxygen gas albedo error of each XCO; retrieval value according to a pre-obtained full-solar-spectrum of an oxygen gas band and a pre-obtained observation spectrum and solar zenith angle of the oxygen gas band corresponding to each XCO, 21
BL-5647
LU503782 retrieval value in the TanSat satellite XCO;z retrieval data; wherein the XCO; retrieval data contains N of the XCO» retrieval values, and N is a positive integer; an observation angle error unit configured to calculate an observation angle error of each XCO» retrieval value according to a satellite observation angle and the solar zenith angle corresponding to each XCO» retrieval value in the N of the XCO; retrieval values; a regression coefficient unit configured to obtain a regression coefficient of the XCO; retrieval data according to N of the XCO; retrieval values and a pre-obtained XCO; foundation observation value corresponding to the XCO; retrieval value; a bias correction unit configured to perform a bias correction on each XCO; retrieval value according to each XCO; retrieval value, the regression coefficient of the XCO; retrieval value, and the oxygen gas albedo error and the observation angle error of each XCO; retrieval value to obtain a bias correction result of each XCO; retrieval value.
It is to be understood that the above description is intended to be illustrative, and that the present examples are not restrictive.
Fig. 7 is a schematic structure diagram of the albedo error unit according to some examples of the present application; as shown in Fig. 7, the albedo error unit comprises: an albedo calculation subunit configured to calculate an oxygen gas band non-absorbing region albedo corresponding to each XCO; retrieval value according to the pre-obtained full-solar-spectrum of the oxygen gas band and a pre-obtained observation spectrum and solar zenith angle of the oxygen gas band corresponding to each XCO; retrieval value in the N of the
X CO; retrieval values; an albedo average calculation subunit configured to average the oxygen gas band non-absorbing region albedo corresponding to N of the XCO; retrieval values to obtain an oxygen gas albedo average value of N of the XCO; retrieval values; an albedo error calculation subunit configured to calculate the oxygen gas albedo error of each
XCO; retrieval value according to the oxygen gas band non-absorbing region albedo corresponding to each XCO; retrieval value, and the oxygen gas albedo average value.
It is to be understood that the above description is intended to be illustrative, and that the present examples are not restrictive.
Fig. 8 is a schematic structure diagram of the observation angle error unit according to some 22
BL-5647
LU503782 examples of the present application; as shown in Fig. 8, the observation angle error unit comprises: an air mass factor calculation subunit configured to calculate an air mass factor of each XCO, retrieval value according to the satellite observation angle and the solar zenith angle corresponding to each XCO-, retrieval value in the N of the XCO; retrieval values; an air mass factor average calculation subunit configured to average the air mass factor of N of the XCO; retrieval values to obtain an air mass factor average value of N of the XCO; retrieval values; an observation angle error calculation subunit configured to calculate the observation angle error of each XCO; retrieval value according to the air mass factor corresponding to each XCO, retrieval value, and the air mass factor average value.
It is to be understood that the above description is intended to be illustrative, and that the present examples are not restrictive.
In some alternative examples, the regression coefficient unit is further configured to perform linear regression on N of the XCO; retrieval values and a pre-obtained XCO; foundation observation value corresponding to the XCO- retrieval value to obtain the regression coefficient of the XCO; retrieval data. It is to be understood that the above description is intended to be illustrative, and that the present examples are not restrictive.
Fig. 9 is a schematic structure diagram of the observation angle error unit according to some examples of the present application; as shown in Fig. 9, the bias correction unit comprises: a first correction subunit configured to perform an albedo error correction and an observation angle error correction on each XCO-, retrieval value respectively according to each XCO; retrieval value, and the oxygen gas albedo error and the observation angle error of each XCO» retrieval value; a second correction subunit configured to perform a quotient operation on each XCO; retrieval value after the albedo error correction and the observation angle error correction and the regression coefficient to obtain the bias correction result of each XCO; retrieval value.
It is to be understood that the above description is intended to be illustrative, and that the present examples are not restrictive.
Fig. 10 is a schematic structure diagram for obtaining TanSat satellite XCO; retrieval data 23
BL-5647
LU503782 provided according to some examples of the present application; as shown in Fig. 10, the bias correction system of TanSat satellite XCO; retrieval data further comprises: a data extraction unit 1001 configured to extract a satellite parameter and an observation parameter of the TanSat satellite; a data elimination unit 1002 configured to perform cloud region data elimination on cloud data, and performing aerosol high value region data elimination on aerosol data; an analog spectrum unit 1003 configured to obtain an analog spectrum based on an atmospheric radiative transfer model of a preset retrieval parameter, according to the satellite parameter, the observation parameter, a cloud data performed cloud region data elimination, an aerosol data performed aerosol high value region data elimination and an atmospheric profile parameter; a spectrum solving unit 1004 configured to solve the analog spectrum and the observation spectrum in the satellite parameter to obtain a solution result based on an optimal estimation method; a result update unit 1005 configured to judge whether the solution result converges, and if not, updating the preset retrieval parameter in the atmospheric radiative transfer model until the solution result converges; a retrieval data unit 1006 configured to obtain the TanSat satellite XCO» retrieval data according to the solution result.
It is to be understood that the above description is intended to be illustrative, and that the present examples are not restrictive.
The bias correction system of TanSat satellite XCO» retrieval data provided by the example of the present application can realize the various processes in the above-mentioned bias correction method examples of the TanSat satellite XCO» retrieval data and achieve the same functions and effects, and will not be described in detail herein. It is to be understood that the above description 1s intended to be illustrative, and that the present examples are not restrictive.
The foregoing is merely a preferred example of the present invention and is not intended to limit the present invention so that various modifications and changes may be made by those skilled in the art. All changes, equivalents, improvements, etc. that come within the spirit and scope of the invention are intended to be contained therein. 24
Claims (10)
1. a bias correction method of TanSat satellite XCO» retrieval data, characterized in that, comprising the following steps: obtaining an oxygen gas albedo error of each XCO» retrieval value according to a pre-obtained full-solar-spectrum of an oxygen gas band and a pre-obtained observation spectrum and solar zenith angle of the oxygen gas band corresponding to each XCO; retrieval value in the TanSat satellite XCO: retrieval data; wherein the XCO» retrieval data contains N of the XCO- retrieval values, and N is a positive integer; calculating an observation angle error of each XCO; retrieval value according to a satellite observation angle and the solar zenith angle corresponding to each XCO; retrieval value in the N of the XCO; retrieval values; obtaining a regression coefficient of the XCO; retrieval data according to N of the XCO, retrieval values and a pre-obtained XCO» foundation observation value corresponding to the XCO, retrieval value; performing a bias correction on each XCO; retrieval value according to each XCO; retrieval value, the regression coefficient of the XCO; retrieval value, and the oxygen gas albedo error and the observation angle error of each XCO; retrieval value to obtain a bias correction result of each X CO; retrieval value, wherein, the step of performing a bias correction on each XCO; retrieval value according to each XCO; retrieval value, the regression coefficient of the XCO; retrieval value, and the oxygen gas albedo error and the observation angle error of each XCO; retrieval value to obtain a bias correction result of each XCO» retrieval value comprises the following steps: performing an albedo error correction and an observation angle error correction on each XCO» retrieval value respectively according to each XCO» retrieval value, and the oxygen gas albedo error and the observation angle error of each XCO- retrieval value; performing a quotient operation on each XCO; retrieval value after the albedo error correction and the observation angle error correction and the regression coefficient to obtain the bias correction result of each XCO-, retrieval value.
2. The method according to claim 1, characterized in that, the step of obtaining an oxygen gas
BL-5647 LU503782 albedo error of each XCO; retrieval value according to a pre-obtained full-solar-spectrum of an oxygen gas band and a pre-obtained observation spectrum and solar zenith angle of the oxygen gas band corresponding to each XCO; retrieval value in the TanSat satellite XCO; retrieval data comprises the following steps: calculating an oxygen gas band non-absorbing region albedo corresponding to each XCO», retrieval value according to the pre-obtained full-solar-spectrum of the oxygen gas band and a pre-obtained observation spectrum and solar zenith angle of the oxygen gas band corresponding to each XCO- retrieval value in the N of the XCO; retrieval values; averaging the oxygen gas band non-absorbing region albedo corresponding to N of the XCO» retrieval values to obtain an oxygen gas albedo average value of N of the XCO; retrieval values; calculating the oxygen gas albedo error of each XCO; retrieval value according to the oxygen gas band non-absorbing region albedo corresponding to each XCO» retrieval value, and the oxygen gas albedo average value.
3. The method according to claim 1, characterized in that, the step of calculating an observation angle error of each XCO» retrieval value according to a satellite observation angle and the solar zenith angle corresponding to each XCO; retrieval value in the N of the XCO; retrieval values comprises the following steps: calculating an air mass factor of each XCO; retrieval value according to the satellite observation angle and the solar zenith angle corresponding to each XCO» retrieval value in the N of the XCO» retrieval values; averaging the air mass factor of N of the XCO; retrieval values to obtain an air mass factor average value of N of the XCO- retrieval values; calculating the observation angle error of each XCO; retrieval value according to the air mass factor corresponding to each XCO; retrieval value, and the air mass factor average value.
4. The method according to claim 1, characterized in that, the step of obtaining a regression coefficient of the XCO» retrieval data according to N of the XCO; retrieval values and a pre-obtained XCO; foundation observation value corresponding to the XCO; retrieval value is specifically: performing linear regression on N of the XCO; retrieval values and a pre-obtained 26
BL-5647 LU503782 XCO; foundation observation value corresponding to the XCO» retrieval value to obtain the regression coefficient of the XCO; retrieval data.
5. The method according to any one of claims 1-4, characterized in that, before the step of obtaining an oxygen gas albedo error of each XCO; retrieval value according to a pre-obtained full-solar-spectrum of an oxygen gas band and a pre-obtained observation spectrum and solar zenith angle of the oxygen gas band corresponding to each XCO; retrieval value in the TanSat satellite XCO» retrieval data, further comprising the following steps: extracting a satellite parameter and an observation parameter of the TanSat satellite; performing cloud region data elimination on cloud data, and performing aerosol high value region data elimination on aerosol data; obtaining an analog spectrum based on an atmospheric radiative transfer model of a preset retrieval parameter, according to the satellite parameter, the observation parameter, a cloud data performed cloud region data elimination, an aerosol data performed aerosol high value region data elimination and an atmospheric profile parameter; solving the analog spectrum and the observation spectrum in the satellite parameter to obtain a solution result based on an optimal estimation method; judging whether the solution result converges, and if not, updating the preset retrieval parameter in the atmospheric radiative transfer model until the solution result converges; obtaining the TanSat satellite XCO; retrieval data according to the solution result.
6. a bias correction system of TanSat satellite XCO; retrieval data, characterized in that, comprising: an albedo error unit configured to obtain an oxygen gas albedo error of each XCO; retrieval value according to a pre-obtained full-solar-spectrum of an oxygen gas band and a pre-obtained observation spectrum and solar zenith angle of the oxygen gas band corresponding to each XCO, retrieval value in the TanSat satellite XCO;z retrieval data; wherein the XCO; retrieval data contains N of the XCO» retrieval values, and N is a positive integer; an observation angle error unit configured to calculate an observation angle error of each XCO» retrieval value according to a satellite observation angle and the solar zenith angle corresponding 27
BL-5647 LU503782 to each XCO» retrieval value in the N of the XCO; retrieval values; a regression coefficient unit configured to obtain a regression coefficient of the XCO; retrieval data according to N of the XCO; retrieval values and a pre-obtained XCO; foundation observation value corresponding to the XCO; retrieval value; a bias correction unit configured to perform a bias correction on each XCO; retrieval value according to each XCO; retrieval value, the regression coefficient of the XCO; retrieval value, and the oxygen gas albedo error and the observation angle error of each XCO; retrieval value to obtain a bias correction result of each XCO- retrieval value, wherein, the bias correction unit comprises: a first correction subunit configured to perform an albedo error correction and an observation angle error correction on each XCO; retrieval value respectively according to each XCO; retrieval value, and the oxygen gas albedo error and the observation angle error of each XCO; retrieval value; a second correction subunit configured to perform a quotient operation on each XCO; retrieval value after the albedo error correction and the observation angle error correction and the regression coefficient to obtain the bias correction result of each XCO; retrieval value.
7. The system according to claim 6, characterized in that, the albedo error unit comprises: an albedo calculation subunit configured to calculate an oxygen gas band non-absorbing region albedo corresponding to each XCO; retrieval value according to the pre-obtained full-solar-spectrum of the oxygen gas band and a pre-obtained observation spectrum and solar zenith angle of the oxygen gas band corresponding to each XCO; retrieval value in the N of the X CO; retrieval values; an albedo average calculation subunit configured to average the oxygen gas band non-absorbing region albedo corresponding to N of the XCO; retrieval values to obtain an oxygen gas albedo average value of N of the XCO- retrieval values; an albedo error calculation subunit configured to calculate the oxygen gas albedo error of each XCO; retrieval value according to the oxygen gas band non-absorbing region albedo corresponding to each XCO; retrieval value, and the oxygen gas albedo average value. 28
BL-5647 LU503782
8. The system according to claim 6, characterized in that, the observation angle error unit comprises: an air mass factor calculation subunit configured to calculate an air mass factor of each XCO, retrieval value according to the satellite observation angle and the solar zenith angle corresponding to each XCO; retrieval value in the N of the XCO; retrieval values; an air mass factor average calculation subunit configured to average the air mass factor of N of the XCO; retrieval values to obtain an air mass factor average value of N of the XCO; retrieval values; an observation angle error calculation subunit configured to calculate the observation angle error of each XCO; retrieval value according to the air mass factor corresponding to each XCO, retrieval value, and the air mass factor average value.
9. The system according to claim 6, characterized in that, the regression coefficient unit is further configured to perform linear regression on N of the XCO; retrieval values and a pre-obtained XCO; foundation observation value corresponding to the XCO; retrieval value to obtain the regression coefficient of the XCO; retrieval data.
10. The system according to claim 6, characterized in that, the bias correction system of TanSat satellite XCO» retrieval data further comprises: a data extraction unit configured to extract a satellite parameter and an observation parameter of the TanSat satellite; a data elimination unit configured to perform cloud region data elimination on cloud data, and performing aerosol high value region data elimination on aerosol data; an analog spectrum unit configured to obtain an analog spectrum based on an atmospheric radiative transfer model of a preset retrieval parameter, according to the satellite parameter, the observation parameter, a cloud data performed cloud region data elimination, an aerosol data performed aerosol high value region data elimination and an atmospheric profile parameter; a spectrum solving unit configured to solve the analog spectrum and the observation spectrum in the satellite parameter to obtain a solution result based on an optimal estimation method; a result update unit configured to judge whether the solution result converges, and if not, 29
BL-5647 LU503782 updating the preset retrieval parameter in the atmospheric radiative transfer model until the solution result converges; a retrieval data unit configured to obtain the TanSat satellite XCO» retrieval data according to the solution result.
BL-5647 ANSPRÜCHE 7509768
1. Verfahren zur Abweichungskorrektur von XCO;-Inversionsdaten für TanSat-Satelliten, dadurch gekennzeichnet, dass es umfasst: Erhalten eines Sauerstoff-Albedo-Fehlers für jeden der XCO,-Inversionswerte auf der Grundlage eines zuvor erhaltenen Sonnenspektrums des Sauerstoffbandes und eines Beobachtungsspektums des Sauerstoffbandes, das jedem XCOz-Inversionswert in den zuvor erhaltenen XCOz-Inversionsdaten für TanSat-Satelliten entspricht, und eines solaren Zenitwinkels, wobei die XCOz-Inversionsdaten N XCOz-Inversionswerte enthalten, wobei N eine positive ganze Zahl ist; Berechnen eines Beobachtungswinkelfehlers für jeden der XCOz-Inversionswerte auf der Grundlage des Satellitenbeobachtungswinkels, der jedem der XCOz-Inversionswerte von N XCO;-Inversionswerten entspricht, und des Sonnenzenitwinkels; Erhalten der Regressionskoeffizient der XCOz-Inversionsdaten auf der Grundlage von N XCO--Inversionswerten und zuvor erhaltenen XCOz-Bodenbeobachtungen, die den XCO--Inversionswerten entsprechen; Durchführen der Abweichungskorrektur für jeden XCO»-Inversionswert auf der Grundlage jedes XCO--Inversionswertes, des Regressionskoeffizienten der XCO--Inversionsdaten und des Sauerstoff-Albedofehlers und des Beobachtungswinkelfehlers für jeden = XCO--Inversionswert, um ein Abweichungskorrekturergebnis für jeden XCOz-Inversionswert zu erhalten; das Durchführen der Abweichungskorrektur für jeden XCO-,-Inversionswert auf der Grundlage jedes XCO--Inversionswertes, des Regressionskoeffizienten der XCO--Inversionsdaten und des Sauerstoff-Albedofehlers und des Beobachtungswinkelfehlers für jeden = XCO--Inversionswert, um ein Abweichungskorrekturergebnis für jeden XCOz-Inversionswert zu erhalten, umfasst: Durchführen einer Albedo-Fehlerkorrektur und einer Beobachtungswinkelfehlerkorrektur auf jeden XCOz-Inversionswert auf der 1
BL-5647 Grundlage jedes XCOz-Inversionswertes und des Sauerstoff-Albedo-Fehlers bzw. des 7505792 Beobachtungswinkelfehlers für jeden XCO;-Inversionswert; Quotientbilden jedes XCO,-Inversionswertes nach der Albedo-Fehlerkorrektur und der Beobachtungswinkelfehlerkorrektur mit dem Regressionskoeffizienten, um ein Abweichungskorrekturergebnis für jeden XCOz-Inversionswert zu erhalten.
2. Verfahren nach Anspruch 1, dadurch gekennzeichnet, dass das Erhalten eines Sauerstoff-Albedo-Fehlers fiir jeden XCO-,-Inversionswert auf der Grundlage eines zuvor erhaltenen Sonnenspektrums des Sauerstoffbandes und eines Beobachtungsspektums des Sauerstoffbandes, das jedem XCOz-Inversionswert in den zuvor erhaltenen XCOz-Inversionsdaten für TanSat-Satelliten entspricht, und eines solaren Zenitwinkels umfasst: Berechnen der Albedo des nicht-absorbierenden Bereichs des Sauerstoffbandes, das jedem XCO--Inversionswert entspricht, auf der Grundlage des zuvor erhaltenen Sonnenspektrums des Sauerstoffbandes und des Beobachtungsspektums des Sauerstoffbandes, das jedem XCO-,-Inversionswert aus den zuvor N zuvor erhaltenen XCO>-Inversionswerten entspricht, und des solaren Zenitwinkels; Mitteln der Albedo des nicht-absorbierenden Bereichs des Sauerstoffbandes, der den N XCOz-Inversionen entspricht, um einen Mittelwert der Sauerstoff-Albedos der N XCO--Inversionen zu erhalten; Berechnen eines Sauerstoff-Albedo-Fehlers für jeden XCO;-Inversionswert auf der Grundlage der Albedo des nicht absorbierenden Bereichs des Sauerstoffbands, die jedem XCOz-Inversionswert entspricht, und des Sauerstoff-Albedo-Mittelwertes.
3. Verfahren nach Anspruch 1, dadurch gekennzeichnet, dass das Berechnen eines Beobachtungswinkelfehlers für jeden der XCOz-Inversionswerte auf der Grundlage des Satellitenbeobachtungswinkels, der jedem der XCO--Inversionswerte von N XCO>-Inversionswerten entspricht, und des Sonnenzenitwinkels umfasst; Berechnen eines Luftmassenfaktors für jeden XCOz-Inversionswert auf der Grundlage des Satellitenbeobachtungswinkels, der jedem XCOz-Inversionswert entspricht, und des Sonnenzenitwinkels; Mitteln der Luftmassenfaktoren für N XCOz-Inversionswerte, um einen Mittelwert 2
BL-5647 der Luftmassenfaktoren für N der XCOz-Inversionswerte zu erhalten; 0508782 Berechnen eines beobachteten Winkelfehlers für jeden XCOz-Inversionswert auf der Grundlage des Luftmassenfaktors jedes XCO»-Inversionswertes und des Luftmassenfaktor-Mittelwertes.
4. Verfahren nach Anspruch 1, dadurch gekennzeichnet, dass das Erhalten der Regressionskoeffizient der XCOz-Inversionsdaten auf der Grundlage von N XCOz-Inversionswerten und zuvor erhaltenen XCOz-Bodenbeobachtungen, die den XCOz-Inversionswerten entsprechen, derart ausgeführt ist, dass eine lineare Regression von N XCO--Inversionswerten und zuvor erhaltenen XCO»-Bodenbeobachtungen, die den XCO:-Inversionswerten entsprechen, durchgeführt wird um einen Regressionskoeffizienten für die XCOz-Inversionsdaten zu erhalten.
5. Verfahren nach einem der Ansprüche 1 bis 4, dadurch gekennzeichnet, dass vor dem Erhalten eines Sauerstoff-Albedo-Fehlers für jeden XCOz-Inversionswert auf der Grundlage eines zuvor erhaltenen Sonnenspektrums des Sauerstoffbandes und eines Beobachtungsspektums des Sauerstoffbandes, das jedem XCOz-Inversionswert in den zuvor erhaltenen XCOs-Inversionsdaten für TanSat-Satelliten entspricht, und des solaren Zenitwinkels das Verfahren ferner umfasst: Extrahieren von Satellitenparametern und Beobachtungsparametern des TanSat-Satelliten; Datenunterdrücken von bewölkten Bereichen für Wolkendaten und Datenunterdrücken von Aerosolbereichen mit hohen Aerosolwerten für Aerosoldaten; Erhalten des Analogspektrums auf der Grundlage eines atmosphärischen Strahlungstransfermodells mit vorbestimmten Inversionsparametern und auf der Basis der Satellitenparameter, der Beobachtungsparameter, der Wolkendaten mit Datenunterdrückung von bewölkten Bereichen, der Aerosoldaten mit Datenunterdrückung von Bereichen mit hohen Aerosolwerten und der atmosphärischen Profilparameter; Lösen des Analogspektrums und des Beobachtungsspektums aus den Satellitenparametern auf der Grundlage eines optimalen Schätzverfahrens, um ein 3
BL-5647 Lôsungsergebnis zu erhalten; 505782 Bestimmen, ob die Lösungsergebnisse konvergieren, und wenn nicht, Aktualisieren der vorbestimmten Inversionsparameter in dem atmosphärischen Strahlungstransportmodell, bis die Lösungsergebnisse konvergieren; Erhalten der XCO;-Inversionsdaten für TanSat-Satelliten basierend auf den Lösungsergebnissen.
6. System zur Abweichungskorrektur von XCO-,-Inversionsdaten für TanSat-Satelliten, dadurch gekennzeichnet, umfassend eine Albedo-Fehlereinheit, die so eingerichtet ist, einen Sauerstoff-Albedo-Fehler für jeden der XCOs-Inversionswerte auf der Grundlage eines zuvor erhaltenen Sonnenspektrums des Sauerstoffbandes und eines zuvor erhaltenen Beobachtungsspektums des Sauerstoffbandes, das jedem XCOz-Inversionswert in den XCOz-Inversionsdaten für TanSat-Satelliten entspricht, und eines solaren Zenitwinkels zu erhalten; wobei die XCOz-Inversionsdaten N XCOz-Inversionswerte umfassen, wobei N eine positive ganze Zahl ist; eine Beobachtungswinkelfehlereinheit, die so eingerichtet ist, einen Beobachtungswinkelfehler für jeden XCOz-Inversionswert auf der Grundlage des Satellitenbeobachtungswinkels, der jedem der XCOz-Inversionswerte von N XCO--Inversionswerten entspricht, und des solaren Zenitwinkels zu berechnen; eine Regressionskoeffizienteneinheit, die so eingerichtet ist, einen Regressionskoeffizienten für die XCOz-Inversionsdaten auf der Grundlage von N XCO--Inversionswerten und zuvor erhaltenen bodengestützten XCO;-Beobachtungen, die den XCOz-Inversionswerten entsprechen, zu erhalten; eine Abweichungskorrektureinheit, die so eingerichtet ist, eine Abweichnungskorrektur für jeden XCOz-Inversionswert auf der Grundlage jedes XCOz-Inversionswertes, eines Regressionskoeffizienten der XCOs-Inversionsdaten und eines Sauerstoff-Albedofehlers und eines Beobachtungswinkelfehlers für jeden XCOz-Inversionswert durchzuführen, um ein Abweichungskorrekturergebnis für jeden XCOz-Inversionswert zu erhalten; wobei die Abweichungskorrektureinheit vorzugsweise umfasst: 4
BL-5647 eine erste Korrekturuntereinheit, die so eingerichtet ist, jeden XCOz-Inversionswert 7505792 um den Albedofehler und den Beobachtungswinkelfehler auf der Grundlage jedes XCO>-Inversionswertes und des Sauerstoff-Albedofehlers und des Beobachtungswinkelfehlers für jeden XCOz-Inversionswert zu korrigieren; eine zweite Korrekturuntereinheit, die so eingerichtet ist, die um den Albedofehler und den Beobachtungswinkelfehler korrigierten XCO»-Inversionen mit den Regressionskoeffizienten zu quotieren, um ein Abweichungskorrekturergebnis für jede XCOs-Inversion zu erhalten.
7. System nach Anspruch 6, dadurch gekennzeichnet, dass die Albedo-Fehlereinheit umfasst: eine Albedo-Berechnungsuntereinheit, die so eingerichtet ist, die Albedo des nicht-absorbierenden Bereichs des Sauerstoffbandes, das jedem XCO--Inversionswert entspricht, auf der Grundlage des zuvor erhaltenen Sonnenspektrums des Sauerstoffbande und des Beobachtungsspektums des Sauerstoffbandes, das jedem XCO--Inversionswert aus den N zuvor erhaltenen XCOz-Inversionswerten entspricht, und des Sonnenzenitwinkels zu berechnen; eine Albedo-Mittelungsberechnungsuntereinheit, die so eingerichtet ist, die Albedo des nicht absorbierenden Bereichs des Sauerstoffbandes, das den N XCOz-Inversionswerten entspricht, zu mitteln, um einen Mittelwert der Sauerstoff-Albedos der N XCOz-Inversionswerte zu erhalten; eine Albedo-Fehlerberechnungs-Untereinheit, die so eingerichtet ist, einen Sauerstoff-Albedo-Fehler fiir jeden XCO--Inversionswert auf der Grundlage der Albedo des nicht-absorbierenden Bereichs des Sauerstoffbands, die jedem XCOz-Inversionswert entspricht, und des Sauerstoff-Albedo-Mittelwertes zu berechnen.
8. System nach Anspruch 6 dadurch gekennzeichnet, dass die Beobachtungswinkelfehlereinheit umfasst: eine Luftmassenfaktor-Berechnungsuntereinheit, die so eingerichtet ist, einen Luftmassenfaktor für jeden XCOz-Inversionswert, der jedem XCOz-Inversionswert entspricht, auf der Grundlage des Satellitenbeobachtungswinkels und des solaren
BL-5647 Zenitwinkels, der jedem XCOz-Inversionswert von N XCOz-Inversionswerten 7505792 entspricht, zu berechnen; eine Luftmassenfaktor-Mittelwertberechnungsuntereinheit, die so eingerichtet ist, die Luftmassenfaktoren der N XCO;z-Inversionswerte zu mitteln, um einen Mittelwert der Luftmassenfaktoren der N XCO,-Inversionswerte zu erhalten; eine Untereinheit zur Berechnung des Beobachtungswinkelfehlers, die so eingerichtet ist, einen Beobachtungswinkelfehler für jeden XCO--Inversionswert auf der Grundlage des Gasmassenfaktors jedes XCOs-Inversionswertes und des Gasmassenfaktordurchschnitts zu berechnen.
9. System nach Anspruch 6, dadurch gekennzeichnet, dass die Regressionskoeffizienten-Finheit ferner so eingerichtet ist, dass eine lineare Regression von N XCO--Inversionswerten und zuvor erhaltenen bodengestützten XCO»-Beobachtungen, die den XCOz-Inversionswerten entsprechen, durchgeführt wird, um einen Regressionskoeffizienten für die XCO»-Inversionsdaten zu erhalten.
10. System nach Anspruch 6, dadurch gekennzeichnet, dass das System zur Abweichungskorrektur von XCOz-Inversionsdaten für TanSat-Satelliten umfasst: eine Datenextraktionseinheit, die so eingerichtet ist, Satellitenparameter und Beobachtungsparameter des TanSat-Satelliten zu extrahieren; eine Datenunterdrückungseinheit, die so eingerichtet ist, eine Datenunterdrückung von bewôlkten Bereichen für Wolkendaten und eine Datenunterdrückung von Aerosolbereichen mit hohen Aerosolwerten für Aerosoldaten durchzuführen: eine Analogspektrumeinheit, die so eingerichtet ist, ein Analogspektrum auf der Grundlage eines atmosphärischen Strahlungstransfermodells mit vorbestimmten Inversionsparametern und auf der Basis der Satellitenparameter, der Beobachtungsparameter, der Wolkendaten mit Datenunterdrückung von bewôlkten Bereichen, der Aerosoldaten mit Datenunterdrückung von Bereichen mit hohen Aerosolwerten und der atmosphärischen Profilparameter zu erhalten; eine Spektrallôsungseinheit, die so eingerichtet ist, das Analogspektrum und das Beobachtungsspektum in den Satellitenparametern auf der Grundlage eines optimalen Schätzverfahrens zu 16sen, um ein Lôsungsergebnis zu erhalten; 6
BL-5647 eine Ergebnisaktualisierungseinheit, die so eingerichtet ist, zu bestimmen, ob die 503782 Lösungsergebnisse konvergieren, und wenn nicht, so die vorbestimmten Inversionsparameter in dem atmosphärischen Strahlungstransportmodell aktualisiert werden, bis die Lösungsergebnisse konvergieren; eine Inversionsdateneinheit, die so eingerichtet ist, die XCO»-Inversionsdaten für TanSat-Satelliten basierend auf den Lôsungsergebnissen zu erhalten.
7
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| LU503782A LU503782B1 (fr) | 2023-03-30 | 2023-03-30 | Bias Correction Method and System of TanSat Satellite XCO2 Retrieval Data |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| LU503782A LU503782B1 (fr) | 2023-03-30 | 2023-03-30 | Bias Correction Method and System of TanSat Satellite XCO2 Retrieval Data |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| LU503782B1 true LU503782B1 (fr) | 2023-10-02 |
Family
ID=88244287
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| LU503782A LU503782B1 (fr) | 2023-03-30 | 2023-03-30 | Bias Correction Method and System of TanSat Satellite XCO2 Retrieval Data |
Country Status (1)
| Country | Link |
|---|---|
| LU (1) | LU503782B1 (fr) |
-
2023
- 2023-03-30 LU LU503782A patent/LU503782B1/fr active IP Right Grant
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Leroy et al. | Retrieval of atmospheric properties and surface bidirectional reflectances over land from POLDER/ADEOS | |
| Morino et al. | Preliminary validation of column-averaged volume mixing ratios of carbon dioxide and methane retrieved from GOSAT short-wavelength infrared spectra | |
| Verhoelst et al. | Metrology of ground-based satellite validation: co-location mismatch and smoothing issues of total ozone comparisons | |
| Zhang et al. | An operational approach for generating the global land surface downward shortwave radiation product from MODIS data | |
| Scholze et al. | Reviews and syntheses: Systematic Earth observations for use in terrestrial carbon cycle data assimilation systems | |
| Li et al. | MAX-DOAS measurements in southern China: retrieval of aerosol extinctions and validation using ground-based in-situ data | |
| CN111579504B (zh) | 基于光学遥感的大气污染成分垂直分布反演方法 | |
| CN105092055B (zh) | 基于冷云目标的气象卫星太阳反射波段辐射定标方法 | |
| Schneider et al. | Accomplishments of the MUSICA project to provide accurate, long-term, global and high-resolution observations of tropospheric {H 2 O, δD} pairs–a review | |
| Berkoff et al. | Nocturnal aerosol optical depth measurements with a small-aperture automated photometer using the moon as a light source | |
| Xie et al. | Calculating NDVI for Landsat7-ETM data after atmospheric correction using 6S model: A case study in Zhangye city, China | |
| Jiménez et al. | Exploring the merging of the global land evaporation WACMOS-ET products based on local tower measurements | |
| CN114817825B (zh) | 一种基于超光谱遥感的排放源co2快速成像方法 | |
| Chang et al. | Calibration and evaluation of precipitable water vapor from MODIS infrared observations at night | |
| Mélin | Uncertainties in ocean colour remote sensing | |
| Rokotyan et al. | A posteriori calculation of δ 18 O and δD in atmospheric water vapour from ground-based near-infrared FTIR retrievals of H 2 16 O, H 2 18 O, and HD 16 O | |
| Hu et al. | Upscaling GOME-2 SIF from clear-sky instantaneous observations to all-sky sums leading to an improved SIF–GPP correlation | |
| Yao et al. | Retrieval of solar-induced chlorophyll fluorescence (SIF) from satellite measurements: comparison of SIF between TanSat and OCO-2 | |
| Palmer et al. | Potential improvements in global carbon flux estimates from a network of laser heterodyne radiometer measurements of column carbon dioxide | |
| LU503782B1 (fr) | Bias Correction Method and System of TanSat Satellite XCO2 Retrieval Data | |
| CN116167202A (zh) | 一种基于多源数据的卫星遥感二氧化碳排放反演方法 | |
| Wang et al. | Investigation of the consistency of atmospheric CO 2 retrievals from different space-based sensors: Intercomparison and spatiotemporal analysis | |
| Sun et al. | Construction of a spatially gridded heat flux map based on airborne flux Measurements using remote sensing and machine learning methods | |
| Reinhardt | On the retrieval of circumsolar radiation from satellite observations and weather model output | |
| CN112161943B (zh) | TanSat卫星XCO2反演数据的偏差校正方法和系统 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| FG | Patent granted |
Effective date: 20231002 |