CN115730176A - Ultraviolet hyperspectral atmospheric ozone vertical column concentration inversion method - Google Patents
Ultraviolet hyperspectral atmospheric ozone vertical column concentration inversion method Download PDFInfo
- Publication number
- CN115730176A CN115730176A CN202211090748.6A CN202211090748A CN115730176A CN 115730176 A CN115730176 A CN 115730176A CN 202211090748 A CN202211090748 A CN 202211090748A CN 115730176 A CN115730176 A CN 115730176A
- Authority
- CN
- China
- Prior art keywords
- ozone
- vertical column
- polynomial
- fitting
- atmospheric
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Abstract
The invention relates to ozone concentration inversion, in particular to an ultraviolet hyperspectral atmospheric ozone vertical column concentration inversion method, which determines a fitting equation based on Lambert-beer law; constructing a linear fitting matrix M, constructing a vector
From the linear fit matrix M and the vector
Converting the fitting equation into a matrix form; solving the concentration value of the ozone in the inclined column by a least square method; calculating an air quality factor AMF; calculating a vertical column concentration value of ozone by using the inclined column concentration value of ozone and an air quality factor AMF; the technical scheme provided by the invention can be thatThe method effectively overcomes the defects of low accuracy and complex calculation of the inversion result of the concentration of the atmospheric ozone vertical column in the prior art.
Description
Technical Field
The invention relates to ozone concentration inversion, in particular to an ultraviolet hyperspectral atmospheric ozone vertical column concentration inversion method.
Background
The atmospheric ozone is mainly concentrated in the stratosphere, only 10% of the ozone is present in the troposphere, and the ozone content of the atmosphere between the stratosphere and the troposphere is less. Ozone in the stratosphere is the main substance (220 nm-300 nm) for absorbing solar ultraviolet radiation, and prevents the solar ultraviolet radiation from penetrating the atmosphere and directly irradiating the earth surface. The ozone layer is a vital presence for all living beings on earth due to the protective effect of ozone on the earth's surface. Research shows that the total amount of ozone in the atmosphere is reduced by 1%, harmful solar ultraviolet radiation reaching the earth surface is increased by 1.5-12%, the probability of human skin cancer is increased by 3%, and the human skin cancer is also attacked by cataract, immune system deficiency, growth retardation and other diseases. Therefore, the detection research on the total amount of ozone in the stratosphere is of great significance not only for the research of the atmospheric science theory, but also for the research of the global ecological environment response.
At present, the international society has higher and higher requirements on real-time monitoring of air quality, and ozone is increasingly concerned by people as an important index for reflecting the air quality. In view of the important role of ozone, it is necessary to obtain high-quality observation data as a basis for studying the distribution and long-term variation tendency of ozone in the stratosphere and convection zone regions.
Zhang Lei analyzes and compares the error characteristics of the total amount of atmospheric ozone inverted by satellite detection in different periods according to the ozone total amount measured value of a Brewer spectrometer at a south pole Zhongshan station in 1993-2015; dou Xin and the like comparatively analyze the consistency of detection results of two ozone total amount observation instruments, namely Dobson and Brewer in 2014-2016, in Hebei Xianghe atmosphere comprehensive observation test stations of atmospheric physics research institute of Chinese academy of sciences; chen Tao evaluating the quality information of the foundation and satellite observed ozone total data by comparing 2008-2012 Lasa station foundation observed ozone total with three satellite inversion products; the Hassan Bencherif carries out the trend estimation of the total amount of ozone according to the ground data and satellite observation data of 1998-2017 of a south Africa Ailin station; liu Li et al, by observing the total amount of ozone through the ground, performed a test analysis on the total amount of ozone of a Fengyun No. three (FY-3A) meteorological satellite TOU (2009, 7 months to 2013, 12 months). Long-term continuous and reliable ozone observation is an important basis for studying ozone changes and their causes of formation.
The global observation technology means of the total amount of atmospheric ozone mainly comprises satellite remote sensing and ground observation, and the satellite remote sensing and the ground observation have advantages and disadvantages and are complementary to each other. The satellite remote sensing technology can provide global and long-term observation of atmospheric ozone based on the powerful functions of the satellite remote sensing technology, and the satellite remote sensing technology is rapidly developed since birth, and the ground-based observation also has irreplaceable advantages. Firstly, the foundation instrument is easy to maintain and calibrate, the total ozone amount observation data has higher stability and continuity, the satellite remote sensing is influenced by sensor calibration, cloud pollution and earth surface reflectivity, and the uncertainty of the inversion result is larger than that of the foundation instrument, so the foundation observation result is also commonly used for correcting the satellite-borne instrument. Secondly, although satellite observation can cover a wider space range and has higher spatial resolution, the satellite has very limited scanning times of the same place every day, so that relative to a foundation instrument, the time resolution of satellite observation is not high enough, short-term significant change of ozone content cannot be captured, and the requirement of detecting rapidly-developed atmospheric pollution events in real time cannot be met.
The ground ozone observation station is generally provided with two kinds of ozone total amount observation instruments, namely Dobson and Brewer, which are channel-type solar photometers, only have a few ozone detection wave bands, can eliminate the influence of other trace gases on ozone inversion only through the difference value of one or two wavelength pairs, and has less available information amount.
Disclosure of Invention
Technical problem to be solved
Aiming at the defects in the prior art, the invention provides an ultraviolet hyperspectral atmospheric ozone vertical column concentration inversion method, which can effectively overcome the defects of low accuracy and complex calculation of atmospheric ozone vertical column concentration inversion results in the prior art.
(II) technical scheme
In order to realize the purpose, the invention is realized by the following technical scheme:
an ultraviolet hyperspectral atmospheric ozone vertical column concentration inversion method comprises the following steps:
s1, determining a fitting equation based on Lambert-beer law:
wherein, I 0i For reference spectrum, R is the day-to-earth distance correction factor, I i For actually measured spectra,. Tau.s KNOWNi Optical thickness, σ, for atmospheric molecular Rayleigh and Raman scattering and aerosol Mie scattering ji Is the standard absorption cross section of gas, qs j Is the concentration value of the gas in the inclined column, P SMOi Fitting a polynomial, P, to the Slow absorbing Structure OFFSi Fitting a polynomial, P, to the stray light correction WLCi Fitting a polynomial to the wavelength correction;
s2, constructing a linear fitting matrix M, wherein the linear fitting matrix M is a matrix with n rows, and each row of the linear fitting matrix M comprises the following data:
wherein n is the number of wavelengths,
is a standard absorption cross section of ozone,
is a standard absorption cross-section for nitrogen dioxide,
is the standard absorption cross section of sulfur dioxide, sigma HCHO Is a standard absorption of formaldehydeCross section, λ i For each pixel point I wavelength, I i For the measured spectrum of each wavelength,
for the measured spectrum I i Is the measured spectrum I i A derivative with respect to wavelength;
s3, constructing a vector
(Vector)
Is an n x 1 column vector, vector
Each column of data of (a):
wherein the content of the first and second substances,
the final reference spectrum after the day-to-ground distance correction is obtained;
s4, according to the linear fitting matrix M and the vector
Converting fitting equations into matrix form
S5, solving the column vector through a least square method
Concentration value of middle ozone
S6, calculating an air quality factor AMF through the following formula:
wherein r is the distance from the center of the earth to the observation station, h EFF To absorb the effective height of the gas, ZA * The sun zenith angle is corrected by atmospheric refraction;
s7, calculating the concentration value of the ozone in the vertical column by the following formula
Preferably, the final reference spectrum
The calculation method of (2) comprises:
using known data of solar irradiance on top of the atmosphere as a reference spectrum I 0i Calculating a day-to-ground distance correction coefficient R according to the time lapse and the local longitude and latitude to obtain a final reference spectrum
Preferably, the sun-ground distance correction coefficient R is calculated by using the following formula:
wherein D =2 π N/365, N is the product of the year.
Preferably, said measured spectrum I i The measurement method of (2), comprising:
collecting direct solar light every 5 minutes from sunrise to sunset and generating continuous hyperspectral solar irradiation at ultraviolet bandThe intensity data is used as the measured spectrum I i 。
Preferably, the optical thickness of the atmospheric molecular Rayleigh and Raman scattering and aerosol Miss scattering τ s KNOWNi Calculated using the formula:
τs KNOWNi ＝σ Ms ·qs SCA
wherein σ Ms Is the molecular scattering absorption cross section, qs SCA Is a standard inclined column concentration value, P is the estimated atmospheric pressure of an observation station, and is calculated by the altitude of the observation station, P STAN Is at standard atmospheric pressure, AMF SCA Is a standard air quality factor.
Preferably, the slow absorbing structure fits a polynomial P SMOi The expansion formula of (c) is as follows:
wherein p is SMOn Is a polynomial coefficient, nsmo is the highest order of the polynomial fit, λ i For the center wavelength of each pixel point i.
Preferably, the stray light correction fitting polynomial P OFFSi The expansion formula of (c) is as follows:
wherein the content of the first and second substances,
for the measured spectrum I i Average value of p OFFSn Is the polynomial coefficient, noffs is the highest order of polynomial fit, λ i For the center wavelength of each pixel point i.
Preferably, the wavelength correction fitting polynomial P WLCi The expansion formula of (c) is as follows:
wherein p is WLCn For each polynomial fit coefficient value, nwlc is the highest order of the polynomial fit, λ i For the center wavelength of each pixel point i.
Preferably, the column vector
Including O to be solved 3 、NO 2 、SO 2 The concentration value of HCHO, and unknown p WLCn 、p OFFSn 、p SMOn Parameter, said column vector
Is represented as follows:
(III) advantageous effects
Compared with the prior art, the inversion method of the concentration of the ultraviolet hyperspectral atmospheric ozone vertical column provided by the invention has the following beneficial effects:
1) Aiming at continuous and narrow-band high spectral resolution data of ultraviolet-band solar direct light obtained by an ultraviolet high-spectrum solar radiation instrument, an inversion algorithm of atmospheric ozone vertical column concentration is researched, and compared with a differential absorption spectrum technology for solar scattered light measurement, the inversion algorithm does not need to carry out complicated radiation transmission calculation and is not influenced by Ring effect;
2) In the inversion process, factors which influence the accuracy of the inversion result, such as Rayleigh and Raman scattering of atmospheric molecules, mie scattering and absorption of aerosol, weakening of illumination intensity by a slow absorption structure in the atmosphere, reduction of measured optical thickness caused by stray light, deviation of the wavelength corresponding to the measured data from the original standard wavelength, and the like, are subjected to error elimination, so that the accuracy of the concentration inversion result of the atmospheric ozone vertical column is effectively ensured.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below. It is obvious that the drawings in the following description are only some embodiments of the invention, and that for a person skilled in the art, other drawings can be derived from them without inventive effort.
FIG. 1 is a schematic flow diagram of the present invention;
fig. 2 to 4 are inversion result graphs obtained by inverting the atmospheric ozone vertical column concentration of the sienna qinling observation station in different time periods by using the inversion method of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention. It is to be understood that the embodiments described are only a few embodiments of the present invention, and not all embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
An ultraviolet hyperspectral atmospheric ozone vertical column concentration inversion method is disclosed, and as shown in figure 1, S1, a fitting equation is determined based on Lambert-beer law:
wherein, I 0i For reference spectra, R is the coefficient of day-to-earth distance correction, I i For actually measured spectra,. Tau.s KNOWNi Optical thickness, σ, for atmospheric molecular Rayleigh and Raman scattering and aerosol Mie scattering ji Is the standard absorption cross section of gas, qs j Is the concentration value of the gas, P SMOi Fitting a polynomial, P, to the slow-absorbing structure OFFSi Fitting a polynomial, P, to the stray light correction WLCi Fitting a polynomial to the wavelength correction;
s2, constructing a linear fitting matrix M, wherein the linear fitting matrix M is a matrix with n rows, and each row of the linear fitting matrix M is as follows:
wherein n is the number of wavelengths,
is a standard absorption cross section of ozone,
is a standard absorption cross-section for nitrogen dioxide,
is the standard absorption cross section of sulfur dioxide, sigma HCHO Is the standard absorption cross section of formaldehyde, lambda i For each pixel point I wavelength, I i For the measured spectrum of each wavelength,
for the measured spectrum I i Is the measured spectrum I i A derivative with respect to wavelength;
s3, constructing a vector
(Vector)
Is an n x 1 column vector, vector
Each column of data of (a):
wherein, the first and the second end of the pipe are connected with each other,
the final reference spectrum after the day-to-ground distance correction is obtained;
s4, according to the linear fitting matrix M and the vector
Converting fitting equations into matrix form
S5, solving the column vector through a least square method
Concentration value of middle ozone
S6, calculating an air quality factor AMF through the following formula:
wherein r is the distance from the center of the earth to the observation station, h EFF To absorb the effective height of the gas, ZA * The sun zenith angle is corrected by atmospheric refraction;
s7, calculating the concentration value of the ozone in the vertical column by the following formula
Final reference spectrum
The calculation method of (2) comprises:
using known data of solar irradiance on top of the atmosphere as a reference spectrum I 0i Calculating a day-to-ground distance correction coefficient R according to the time lapse and the local longitude and latitude to obtain a final reference spectrum
The daily and terrestrial distance correction coefficient R is calculated by the following formula:
wherein D =2 π N/365, N is the product of the year.
Measured spectrum I i The measurement method of (2), comprising:
collecting direct solar light once every 5 minutes from sunrise to sunset, generating continuous hyperspectral solar irradiance data of an ultraviolet waveband, and taking the data as an actually measured spectrum I i 。
Optical thickness ts of atmospheric molecular Rayleigh and Raman scattering and aerosol Mie scattering KNOWNi Calculated using the formula:
τs KNOWNi ＝σ Ms ·qs SCA
wherein σ Ms Is the molecular scattering absorption cross section, qs SCA Is a standard inclined column concentration value, P is the estimated atmospheric pressure of an observation station, and is calculated by the altitude of the observation station, P STAN Is at standard atmospheric pressure, AMF SCA Is a standard air quality factor.
Fitting polynomial P of slow absorption structure SMOi The expansion formula of (c) is as follows:
wherein p is SMOn Is a polynomial coefficient, nsmo is the highest order of the polynomial fit, λ i For the center wavelength of each pixel point i.
Stray light correction fitting polynomial P OFFSi The expansion formula of (c) is as follows:
wherein the content of the first and second substances,
for the measured spectrum I i Average value of p OFFSn Is the polynomial coefficient, noffs is the highest order of polynomial fit, λ i For the center wavelength of each pixel point i.
Wavelength correction fitting polynomial P WLCi The expansion formula of (c) is as follows:
wherein p is WLCn For each polynomial fit coefficient value, nwlc is the highest order of the polynomial fit, λ i For the center wavelength of each pixel point i.
Column vector
Including O to be solved 3 、NO 2 、SO 2 The concentration value of HCHO, and unknown p WLCn 、p OFFSn 、p SMOn Parameter, column vector
Is represented as follows:
when ultraviolet radiation passes through the atmosphere, the ultraviolet radiation is attenuated in an e-exponential manner after being scattered and absorbed by atmospheric molecules, scattered and absorbed by aerosol and absorbed by ozone, namely Lambert-beer law:
wherein, λ is wavelength, I (λ) is instrument observation spectrum, I 0 (lambda) is the solar spectrum of the top of the atmospheric layer, R is the correction coefficient of the distance between the sun and the ground corresponding to the measurement time, n EX Number of extinction processes in the atmosphere,. Tau.s j (λ) is the optical thickness of the extinction process j at the wavelength λ.
Taking the logarithm of both sides of equation (1) can be transformed into the following equation:
due to the fact that in the formula (2)
Considering only the weakening part of the trace gas in the atmosphere on the illumination intensity, and not considering the Rayleigh and Raman scattering of atmospheric molecules, the Mie scattering and absorption of aerosol, and the weakening of the illumination intensity by slow absorption structures in the atmosphere, the formula (2) is rewritten to obtain the following formula:
wherein the optical thickness of atmospheric molecular Rayleigh and Raman scattering and aerosol Mie scattering τ s KNOWNi Calculated using the formula:
τs KNOWNi ＝σ Ms ·qs SCA
in the above formula, σ Ms Is the molecular scattering absorption cross section, qs SCA Is a standard inclined column concentration value, P is the estimated atmospheric pressure of an observation station, and is calculated by the altitude of the observation station, P STAN Is at standard atmospheric pressure, AMF SCA Is a standard air quality factor.
Wherein the slow absorbing structure is fitted with a polynomial P SMOi The expansion formula of (c) is as follows:
in the above formula, p SMOn Is a polynomial coefficient, nsmo is the highest order of the polynomial fit, λ i For the center wavelength of each pixel point i.
In real-world measurements, stray light is an unavoidable problem for many optical systems. Stray light in measurement noise includes stray light of the spectrometer itself and stray light entering the detector from the outside, and mainly comes from scattered light of optical devices (gratings, plane mirrors, etc.), reflected light of the wall of the spectrometer, reflected light of unused spectral bands of the wall near the focal plane, reflected light of the surface of the detector, and the like.
The presence of stray light may result in a reduction in the measured optical thickness, thereby reducing the inversion results of the gas concentration. In the actual measurement, the intensity of stray light entering the detector can be reduced by adding a filter, and meanwhile, the influence of residual stray light can be removed by adopting a polynomial fitting method in the spectrum processing process, namely, a stray light correction fitting polynomial P is added in the formula (3) OFFSi ：
Wherein, the stray light is corrected to fit a polynomial P OFFSi The expansion formula of (c) is as follows:
in the above formula, the first and second carbon atoms are,
for the measured spectrum I i Average value of p OFFSn Is the polynomial coefficient, noffs is the highest order of polynomial fit, λ i For the center wavelength of each pixel point i.
Generally, after a spectrometer is calibrated in a laboratory, the wavelength corresponding to each channel of the detector can be accurately obtained. However, during actual measurement, the spectrometer may be affected by external environmental factors such as temperature and humidity changes of the working environment, small fluctuation of the working voltage, and mechanical vibration, and the wavelength corresponding to the measurement data will deviate from the original standard wavelength. Thus, the wavelength of the measurement data can be calibrated accordingly during data processing by adding a wavelength correction fitting polynomial P to equation (4) WLCi And correcting the measured spectrum:
wherein the wavelength correction fits a polynomial P WLCi The expansion formula of (c) is as follows:
in the above formula, p WLCn For each polynomial fit coefficient value, nwlc is the highest order of the polynomial fit, λ i For the center wavelength of each pixel point i.
Since equation (5) is a non-linear equation, it can be obtained by linearizing it:
4 gases are mainly considered to be sun in ultraviolet bandInfluence of irradiance, optical thickness value of each gas
The expansion is as follows:
in the above formula, σ ji Are each O 3 、NO 2 、SO 2 Standard absorption cross section of HCHO.
For polynomial P SMOi 、P OFFSi 、P WLCi Unfolding is carried out and formula (6) is rearranged to obtain the following formula:
in the above formula, nwlc =1,noffs =0,nsmo =4.
Since the values on the left side of equation (7) are known, σ in the right side ji Since it is also known that the formula (7) can be rewritten in a matrix form
Wherein the vector
Is an n x 1 column vector, vector
Each column of data of (a):
column vector
Including O to be solved 3 、NO 2 、SO 2 The concentration value of HCHO, and unknown p WLCn 、p OFFSn 、p SMOn Parameter, column vector
Is represented as follows:
the linear fitting matrix M is a matrix of n rows, each row being O 3 、NO 2 、SO 2 The absorption cross section value of HCHO corresponding to the wavelength can be solved
And the center wavelength lambda of each pixel point i i Thus, the data in the linear fit matrix M are all known values. Each row of data of the linear fitting matrix M is as follows:
since it is continuous spectrum data, the number of wavelengths n must be larger than the column vector
The number of unknowns, so the column vector can be solved by the least square method
Concentration value of middle ozone
Because of the concentration value of ozone
Depending on the way of observation of the instrument and the various meteorological conditions at the time, it is therefore also necessary to switch to ozone independently of the way of observationConcentration value of (2)
Which represents the integrated concentration of the trace gas concentration along a vertical path through the atmosphere.
Air quality factor AMF as ozone inclined column concentration value
Concentration value of ozone on vertical column
Thus calculating the vertical column concentration value of ozone by the following formula
When direct sunlight observation is used, the air quality factor AMF is calculated by:
in the above formula, r is the distance from the center of the earth to the observation station, h EFF To absorb the effective height of the gas, ZA * The solar zenith angle is corrected by atmospheric refraction.
Through the solved air quality factor AMF and the concentration value of ozone in the inclined column
The concentration value of ozone in the vertical column can be solved
Fig. 2 to 4 are result graphs obtained by inverting the atmospheric ozone vertical column concentration of the sienna qinling observation station at different times (corresponding to 23 days 6 month, 30 days 6 month and 3 days 7 month in 2022 respectively) by using the inversion method in the invention, wherein a red curve is an inversion result of an ultraviolet hyperspectral ozone observer, and black x point is the total amount of ozone columns of an Aura satellite ozone Observer (OMI).
Since the transit time of a satellite is typically local 13: around 45, therefore OMI observations and 13: and (4) comparing the concentration inversion values of the ozone vertical column at the 45-time. 13, the three days: the relative deviations of the inversion values of the ozone vertical column concentration at time 45 and the OMI observations were 2DU (0.6%), 14DU (4.5%) and 7DU (2.2%), respectively. It can be seen that the concentration of the ozone vertical column obtained by inversion by the method is closer to the total amount of the OMI ozone column.
The above examples are only intended to illustrate the technical solution of the present invention, but not to limit it; although the present invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some technical features may be equivalently replaced; such modifications and substitutions do not depart from the spirit and scope of the corresponding technical solutions.
Claims (9)
1. An ultraviolet hyperspectral atmospheric ozone vertical column concentration inversion method is characterized by comprising the following steps: the method comprises the following steps:
s1, determining a fitting equation based on Lambert-beer law:
wherein, I 0i For reference spectra, R is the coefficient of day-to-earth distance correction, I i For actually measured spectra,. Tau.s KNOWNi Optical thickness, σ, for atmospheric molecular Rayleigh and Raman scattering and aerosol Mie scattering ji Is the standard absorption cross section of gas, qs j Is the concentration value of the gas in the inclined column, P SMOi Fitting a polynomial, P, to the Slow absorbing Structure OFFSi Fitting for stray light correctionPolynomial of formula P WLCi Fitting a polynomial to the wavelength correction;
s2, constructing a linear fitting matrix M, wherein the linear fitting matrix M is a matrix with n rows, and each row of the linear fitting matrix M is as follows:
wherein n is the number of wavelengths,
is a standard absorption cross section of ozone,
is a standard absorption cross-section for nitrogen dioxide,
standard absorption cross section, σ, for sulfur dioxide HCHO Is the standard absorption cross section of formaldehyde, lambda i For each pixel point I wavelength, I i For the measured spectrum of each wavelength,
for the measured spectrum I i Is the measured spectrum I i A derivative with respect to wavelength;
s3, constructing a vector
(Vector)
Is an n x 1 column vector, vector
Each column of data of (a):
wherein, the first and the second end of the pipe are connected with each other,
the final reference spectrum after the day-to-ground distance correction is obtained;
s4, according to the linear fitting matrix M and the vector
Converting fitting equations into matrix form
S5, solving the column vector through a least square method
Concentration value of middle ozone
S6, calculating an air quality factor AMF through the following formula:
wherein r is the distance from the center of the earth to the observation station, h EFF To absorb the effective height of the gas, ZA * The sun zenith angle is corrected by atmospheric refraction;
s7, calculating the concentration value of the ozone in the vertical column by the following formula
2. The ultraviolet hyperspectral atmospheric ozone vertical column concentration inversion method according to claim 1, characterized by comprising the following steps: the final reference spectrum
The calculation method of (2) comprises:
using known data of solar irradiance on top of the atmosphere as a reference spectrum I 0i Calculating a day-to-ground distance correction coefficient R according to the time lapse and the local longitude and latitude to obtain a final reference spectrum
3. The ultraviolet hyperspectral atmospheric ozone vertical column concentration inversion method according to claim 2, characterized by comprising the following steps: the sun-ground distance correction coefficient R is calculated by adopting the following formula:
wherein D =2 π N/365, N is the product of the year.
4. The ultraviolet hyperspectral atmospheric ozone vertical column concentration inversion method according to claim 1 is characterized in that: said measured spectrum I i The measurement method of (2), comprising:
collecting direct solar light once every 5 minutes from sunrise to sunset, generating continuous hyperspectral solar irradiance data of an ultraviolet waveband, and taking the data as an actually measured spectrum I i 。
5. The ultraviolet hyperspectral atmospheric ozone vertical column concentration inversion method according to claim 1, characterized by comprising the following steps: optical thickness of the atmospheric molecular Rayleigh and Raman scattering and aerosol Mie scattering τ s KNOWNi By usingThe following formula is calculated:
τs KNOWNi ＝σ Ms ·qs SCA
wherein σ Ms Is the molecular scattering absorption cross section, qs SCA Is a standard inclined column concentration value, P is the estimated atmospheric pressure of an observation station, and is calculated by the altitude of the observation station, P STAN Is at standard atmospheric pressure, AMF SCA Is a standard air quality factor.
6. The ultraviolet hyperspectral atmospheric ozone vertical column concentration inversion method according to claim 1, characterized by comprising the following steps: the slow absorbing structure fitting polynomial P SMOi The expansion formula of (c) is as follows:
wherein p is SMOn Is a polynomial coefficient, nsmo is the highest order of the polynomial fit, λ i For the center wavelength of each pixel point i.
7. The ultraviolet hyperspectral atmospheric ozone vertical column concentration inversion method according to claim 1, characterized by comprising the following steps: the stray light correction fitting polynomial P OFFSi The expansion formula of (c) is as follows:
wherein the content of the first and second substances,
for the measured spectrum I i Average value of p OFFSn Is the polynomial coefficient, noffs is the highest order of polynomial fit, λ i For the center wavelength of each pixel point i.
8. The ultraviolet hyperspectral atmospheric ozone vertical column concentration inversion method according to claim 1, characterized by comprising the following steps: the wavelength correction fitting polynomial P WLCi The expansion formula of (c) is as follows:
wherein p is WLCn For each polynomial fit coefficient value, nwlc is the highest order of the polynomial fit, λ i For the center wavelength of each pixel point i.
9. The ultraviolet hyperspectral atmospheric ozone vertical column concentration inversion method according to claim 1 is characterized in that: the column vector
Including O to be solved 3 、NO 2 、SO 2 The concentration value of HCHO, and unknown p WLCn 、p OFFSn 、p SMOn Parameter, said column vector
Is represented as follows:
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211090748.6A CN115730176A (en) | 2022-09-07 | 2022-09-07 | Ultraviolet hyperspectral atmospheric ozone vertical column concentration inversion method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211090748.6A CN115730176A (en) | 2022-09-07 | 2022-09-07 | Ultraviolet hyperspectral atmospheric ozone vertical column concentration inversion method |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN115730176A true CN115730176A (en) | 2023-03-03 |
Family
ID=85293145
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211090748.6A Pending CN115730176A (en) | 2022-09-07 | 2022-09-07 | Ultraviolet hyperspectral atmospheric ozone vertical column concentration inversion method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115730176A (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116297279A (en) * | 2023-05-18 | 2023-06-23 | 至芯半导体(杭州)有限公司 | Method, system, device and equipment for detecting concentration of formaldehyde gas/VOC gas |
| CN117009818A (en) * | 2023-04-12 | 2023-11-07 | 中国科学院空天信息创新研究院 | Atmospheric CO based on polynomial correction 2 Fusion method and system |
-
2022
- 2022-09-07 CN CN202211090748.6A patent/CN115730176A/en active Pending
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117009818A (en) * | 2023-04-12 | 2023-11-07 | 中国科学院空天信息创新研究院 | Atmospheric CO based on polynomial correction 2 Fusion method and system |
| CN117009818B (en) * | 2023-04-12 | 2024-01-23 | 中国科学院空天信息创新研究院 | Atmospheric CO based on polynomial correction 2 Fusion method and system |
| CN116297279A (en) * | 2023-05-18 | 2023-06-23 | 至芯半导体(杭州)有限公司 | Method, system, device and equipment for detecting concentration of formaldehyde gas/VOC gas |
| CN116297279B (en) * | 2023-05-18 | 2023-12-19 | 至芯半导体(杭州)有限公司 | Method, system, device and equipment for detecting concentration of formaldehyde gas/VOC gas |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Timofeyev et al. | Ground-based spectroscopic measurements of atmospheric gas composition near Saint Petersburg (Russia) | |
| Kaskaoutis et al. | Investigation into the wavelength dependence of the aerosol optical depth in the Athens area | |
| Schneising et al. | Three years of greenhouse gas column-averaged dry air mole fractions retrieved from satellite–Part 2: Methane | |
| Schneising et al. | Three years of greenhouse gas column-averaged dry air mole fractions retrieved from satellite–Part 1: Carbon dioxide | |
| Yoshida et al. | Retrieval algorithm for CO 2 and CH 4 column abundances from short-wavelength infrared spectral observations by the Greenhouse gases observing satellite | |
| Li et al. | MAX-DOAS measurements in southern China: retrieval of aerosol extinctions and validation using ground-based in-situ data | |
| CN115730176A (en) | Ultraviolet hyperspectral atmospheric ozone vertical column concentration inversion method | |
| Herman et al. | Comparison of ozone retrievals from the Pandora spectrometer system and Dobson spectrophotometer in Boulder, Colorado | |
| CN111579504B (en) | Atmospheric pollution component vertical distribution inversion method based on optical remote sensing | |
| Salerno et al. | Novel retrieval of volcanic SO2 abundance from ultraviolet spectra | |
| CN111707622B (en) | Method for measuring and calculating vertical distribution and transport flux of atmospheric water vapor based on foundation MAX-DOAS | |
| Sun et al. | Deriving the slit functions from OMI solar observations and its implications for ozone-profile retrieval | |
| Wang et al. | Vertical profiles of tropospheric ozone from MAX‐DOAS measurements during the CINDI‐2 campaign: Part 1—Development of a new retrieval algorithm | |
| Herman et al. | Ultraviolet radiation at the Earth's surface | |
| KR101968087B1 (en) | Apparatus for retrieving column amounts of ozone and nitrogen dioxide in the air using observation data of hyperspectral imaging sensor | |
| Chan et al. | Dispersive infrared spectroscopy measurements of atmospheric CO2 using a Fabry–Pérot interferometer sensor | |
| Yu et al. | Impact of 3D Cloud Structures on the Atmospheric Trace Gas Products from UV-VIS Sounders–Part II: impact on NO 2 retrieval and mitigation strategies | |
| Ortega et al. | The CU 2-D-MAX-DOAS instrument–Part 2: Raman scattering probability measurements and retrieval of aerosol optical properties | |
| Virolainen et al. | Ground-based measurements of total ozone content by the infrared method | |
| Gkertsi et al. | DOAS-based total column ozone retrieval from Phaethon system | |
| Meloni et al. | Comparison of ground‐based and Total Ozone Mapping Spectrometer erythemal UV doses at the island of Lampedusa in the period 1998–2003: Role of tropospheric aerosols | |
| Park et al. | Effect of temperature-dependent cross sections on O4 slant column density estimation by a space-borne UV–visible hyperspectral sensor | |
| Griggs et al. | Comparison of spectrally resolved outgoing longwave radiation over the tropical Pacific between 1970 and 2003 using IRIS, IMG, and AIRS | |
| An et al. | Tomographic retrieval algorithm of OH concentration profiles using double spatial heterodyne spectrometers | |
| Bernhard et al. | Comparison of measured and modeled spectral ultraviolet irradiance at Antarctic stations used to determine biases in total ozone data from various sources |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination |