CN104297816A - Calibration correcting method for ultraviolet radiation meter - Google Patents
Calibration correcting method for ultraviolet radiation meter Download PDFInfo
- Publication number
- CN104297816A CN104297816A CN201410580265.3A CN201410580265A CN104297816A CN 104297816 A CN104297816 A CN 104297816A CN 201410580265 A CN201410580265 A CN 201410580265A CN 104297816 A CN104297816 A CN 104297816A
- Authority
- CN
- China
- Prior art keywords
- radiation
- attenuation rate
- time
- formula
- meter
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01W—METEOROLOGY
- G01W1/00—Meteorology
- G01W1/18—Testing or calibrating meteorological apparatus
Landscapes
- Environmental & Geological Engineering (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Atmospheric Sciences (AREA)
- Biodiversity & Conservation Biology (AREA)
- Ecology (AREA)
- Environmental Sciences (AREA)
- Photometry And Measurement Of Optical Pulse Characteristics (AREA)
Abstract
The invention discloses a calibration correcting method for an ultraviolet radiation meter. The calibration correcting method comprises the steps that the real ultraviolet irradiance reaching the ground at a point is expressed as F(lambda), the ultraviolet irradiance output by a model is expressed as F[m] (lambda), and then the formula F(lambda)=F[m] (lambda) is obtained; the ultraviolet irradiance observed by an aged ultraviolet radiation meter is expressed as F[ab] (lambda), and then attenuation rate can be expressed as the formula D=1- , and the attenuation rate D is a function related to time, namely, D=F(t), wherein D is the monthly attenuation rate, t is the number of service months of the meter; 85 fine cloudless days are selected, and the ultraviolet irradiances of 106 time points are adopted; the time series of D is obtained according to the formula, the attenuation rate of all the months are averaged, so that the monthly average attenuation rate D[mean] is obtained, fitting is conducted on the monthly average attenuation rate D[mean] and t, the attenuation rate obtained after fitting changes exponentially along with time, the relevance (R<2>) reaches 0.96, and please see the fitting equation in the equation (4); the attenuation rate of the ultraviolet radiation meter gradually becomes larger within short time, namely D=0.00116t<2>, and then a calibration correcting formula F'(t)= is obtained. The calibration correcting method is simple and economical and reduces the use cost of the meter.
Description
Technical field
The present invention relates to climatological observation instrument, be specifically related to a kind of demarcation bearing calibration of ultraviolet light radiation meter.
Background technology
UV radiation refers to that wavelength is the radiation of 100 ~ 400 nm.Generally according to the difference of wavelength, UV radiation is divided into UVA (315 ~ 400 nm), UVB (280 ~ 315 nm) and UVC (100 ~ 280 nm) (Kaufman et al. 2002).Wherein UVC is difficult to arrive ground by high-rise Ozone Absorption substantially, and UVB has very important impact to health, affects crop growth simultaneously, and UVA can make the skin darkening of people, there is material impact to animal and plant growth.UV radiation has obvious impact to health.The sunshine medium ultraviolet radiation arriving earth surface can help human body synthesis of vitamin d, K, but also can damage skin and eyes, Immunosuppression system (Penner et al. 1992).According to the erythema effect (erythema) (McKinley et al. 1987) of UV radiation to skin, Burrows and Austin etc. propose computing method (the Burrows et al. 1994 of ultraviolet index (UV index); Austin et al. 1994), to characterize the extent of injury of UV radiation to human body.UV radiation is one of photochemically reactive energy, although UV radiation proportion in solar radiation energy is less, but because its light quantum energy level is higher, especially the biological effect such as photochemical effect, photosynthesis of plant that produces of UVB very significantly (Lars 2007).
Along with UV radiation problem more and more pay close attention to by people, meteorological department has all over China carried out UV radiation monitoring and forecast (Wu 2001 successively in recent years; Shen et al. 2002).Many scientific researches both domestic and external and service unit have installed ultraviolet light radiation meter, carry out Measurement of UV irradiance research.Along with the increase of service time, Aging equipment, observation data there will be decay.Tan et al. (2008) thinks, the main cause of Aging equipment is that sensing element is aging and diffusion sheet is contaminated.At present, China there is no for the fixed standard of Solar UV radiation exposures scale, can only return the country of origin and again demarcate, otherwise, according to the requirement of national metrological regulations, can not continue to use.For the instrument having certain spectral width, can not demarcate as light source using standard lamp, then go to measure Solar UV radiation exposures degree (Angstrom et al.1962) under physical environment.Owing to lacking autonomous demarcation means, ultraviolet light radiation meter is sent back to genuine and demarcate the certain expense of needs, and consuming time longer.
Summary of the invention
The object of the invention is to, for the above-mentioned defect of prior art, provide a kind of demarcation bearing calibration of ultraviolet light radiation meter, the maintenance cost of ultraviolet light radiation meter can be reduced and upgrade cost.
A kind of demarcation bearing calibration of ultraviolet light radiation meter is provided; It corrects in following condition: ground receiver to UV radiation comprise direct ultraviolet radiation and scattering UV radiation; Direct ultraviolet radiation is made up of molecular scattering, gas absorption delustring and aerosol extinction three part, and scattering UV radiation mainly comprises molecular scattering, atmospheric aerosol scattering and cloud and mist scattering; Its bearing calibration is as follows:
If the true uv radiation intensity on the arrival ground of certain time is
, the uv radiation intensity that pattern exports is
, then obtain following formula:
(1)
With the uv radiation intensity of aging ultraviolet light radiation meter observation be
, then attenuation rate can be expressed as:
(2)
Attenuation rate D is the function with time correlation, that is:
(3)
Wherein D is month by month attenuation rate, and t is that instrument uses month;
Select ceiling unlimited day totally 85 D, adopt the uv radiation intensity (UVA) of 106 time points.Obtain the uv radiation intensity of pattern output by TUV mode computation, by the sample mean in every month, then carry out matching with the time; Obtain the time series of D according to formula (2), attenuation rate is monthly done on average, obtain monthly attenuation rate D
mean, then carry out matching with t, the attenuation rate after matching is exponentially property change in time, correlativity (R
2) reaching 0.96, fit equation is shown in formula (4); Ultraviolet light radiation meter is within shorter service time, and attenuation rate can be approximately linear attenuation. and time long when deployed, rate of decay becomes large gradually;
(4)
According to fitting result, obtaining demarcation updating formula is:
(5)
Wherein
for revised radiation value.
Tool of the present invention has the following advantages: the uv radiation intensity on the arrival ground utilizing the TUV pattern simulation in this method to calculate is a kind of simple, economic bearing calibration to correct UVS-AB-T ultraviolet light radiation meter; Reduce the use cost of instrument.
Accompanying drawing explanation
Fig. 1 is that TUV pattern of the present invention orders the UV radiation loose some schematic diagram calculated and observe;
Fig. 2 is the fit correlation schematic diagram that attenuation rate D of the present invention and instrument use moon number t;
Fig. 3 is before the present invention corrects and corrects rear Table A and table B falls apart point data schematic diagram.
Embodiment
Shown in ginseng Fig. 1,2,3; The T-shaped ultraviolet light radiation meter of UVS-AB – two that are arranged on the basic weather station of country of Guangdong Province being to the production of Kipp & Zonen company carries out correction demarcation, the set-up time of two ultraviolet light radiation meters is respectively in Dec, 2007 (title Table A), with in May, 2010 (claiming table B), the commencement of commercial operation time is respectively in January, 2008 and in July, 2010.And by U.S. Center for Atmospheric Research (NCAR, http://ncar.ucar.edu) based on the employing troposphere UV radiation of Mie-scattering lidar exploitation and the radiative transmission mode TUV pattern (Madroniich 1993) of some visible light radiation, a kind of method of simple economy is proposed, the data of UV radiation comparing ultraviolet light radiation meter are exported by radiative transmission mode, analyze ultraviolet light radiation meter data decay characteristics, further proposition correction method, and with the new ultraviolet light radiation meter installed, revised data are checked.The uv radiation intensity that the uv radiation intensity calculated by fitting mode and ultraviolet light radiation meter are recorded, finds that attenuation rate and time are that power exponent relation decays.According to fitting result, obtain the calibration formula of ultraviolet light radiation meter:
, wherein
for revised radiation value.Utilize calibration formula, correct the data of old ultraviolet light radiation meter, further with the ultraviolet light radiation meter Data Comparison of a same model of newly installing, correct formula to check.Research finds, the data dependence of two instruments is very high, and the use value that can be improved instrument by the method corrected is described.Before correcting, the attenuation rate of Table A reaches 23.2%, and after correcting, the irradiance average of old table is 99.7% of table B, respond well.
The observation of uv radiation intensity is extremely important for the calculating of health, photochmeical reaction speed, but due to Aging equipment, makes instrument monitoring data occur decay.Utilize the uv radiation intensity that ultraviolet radiation transmission pattern (TUV) analog computation obtains, monitor with ultraviolet light radiation meter the uv radiation intensity data fitting obtained, set up the linear attenuation formula with time correlation, the UV radiation data of decay are corrected.
Pattern and data:
Two ultraviolet light radiation meters that this method embodiment adopts are the T-shaped ultraviolet light radiation meter of UVS-AB – that Kipp & Zonen company produces, be arranged on the basic station (station number 59289) of Dongguan City plate ridge country, infield is separated by about 2 m.Country's basic stance is in botanical garden, and it has an area of about 1 km for green plants, mainly fruit tree, ornamental plant and careless face, away from urban district and road, affects less by Local anesthesia source.The observation spectrum district of TVS-AB-T type ultraviolet light radiation meter is 315 ~ 400nm (UVA) and 280 ~ 315nm (UVB).The set-up time of two ultraviolet light radiation meters is respectively in Dec, 2007 (title Table A) and in May, 2010 (claiming table B), and the commencement of commercial operation time is respectively in January, 2008 and in July, 2010.
TUV pattern (Madroniich 1993) is U.S. Center for Atmospheric Research (NCAR, http://ncar.ucar.edu) be based on Mie-scattering lidar exploitation the troposphere UV radiation of employing and the radiative transmission mode of some visible light radiation, pattern wavelength 21 ~ 735 nm., resolution is 0.05nm.TUV pattern have employed ozone profile and Elterman. (1968) the gasoloid profile of United States standard atmosphere.TUV pattern main input parameter comprises: longitude and latitude, date, large gas column ozone concentration, surface albedo, 550 nm aerosol optical depths, gasoloid single scattering albedo, wave length exponent α and surface pressure etc.Main Output rusults is downward irradiance (w/ (the um m of any wave band
2)).The error of the uv radiation intensity that TUV pattern exports is within 5%.
In the parameter of TUV pattern input, 550 nm aerosol optical depths are from Hong Kong University of Science and Thchnology. and large gas column ozone concentration setting value is 320 DU (Tang et al. 1990).On the south the Qinling Mountains utilizing MODIS satellite Retrieval according to Wang et al (2004), surface albedo annual mean is between 0.1 ~ 0.14, and setting surface albedo is the general surface albedo 0.1 in careless face and forest land.Gasoloid single scattering albedo is ground, the Pearl River Delta plot mean 0.92 of Xia et al. (2013) inverting.According to the fitting result of Tan et al. (2009), South China's Wavelength Indices
value scope be 1.2 ~ 1.6, mean value is 1.4, thus setting Wavelength Indices be
be 1.4.
Ground receiver to UV radiation comprise direct ultraviolet radiation and scattering UV radiation.Direct ultraviolet radiation is made up of molecular scattering (Rayleigh), gas absorption delustring (as ozone, oxygen, steam etc.) and aerosol extinction 3 part, and scattering UV radiation mainly comprises molecular scattering (Rayleigh), atmospheric aerosol scattering and cloud and mist scattering.Cloud and mist radiation is very large on the impact of UV radiation, and for removing the impact of cloud and mist radiation, the UV radiation data picking ceiling unlimited day are herein corrected.
The true uv radiation intensity on the arrival ground of certain time is
, the uv radiation intensity that pattern exports is
, then can think:
(1)
With the uv radiation intensity of aging ultraviolet light radiation meter observation be
, then attenuation rate can be expressed as:
(2)
Attenuation rate D is the function with time correlation, that is:
(3)
Wherein d is month by month attenuation rate, and t is that instrument uses month.
Select ceiling unlimited day totally 85 d in January, 2008 to 2013 4 months, because the AOD data of MODIS inverting is limited, so altogether adopt the uv radiation intensity (UVA) of 106 time points.Obtain the uv radiation intensity of pattern output by TUV mode computation, contrast with the uv radiation intensity of observation, as shown in Figure 1.Both correlativity (R
2) be 0.592, there is the good degree of correlation, but also have certain deviation.Occur that the main cause departed from may be main from two aspects.First be that the Characteristics of The Aerosol parameter that TUV pattern exports has error.The gasoloid single scattering albedo of pattern input, large gas column ozone concentration and Wavelength Indices are fixed value, and in fact these parameters are not changeless.But owing to lacking inverting means at present, these parameters data accurately cannot be obtained.Next is the impact cannot removing micro-cloud and mist according to naked-eye observation ceiling unlimited weather completely, causes the uv radiation intensity of ground observation to occur error.These two kinds of errors are all random generations, so the positive and negative attribute of error is also random.
In order to reduce the impact that error is brought, first by the sample mean in every month, then carry out matching with the time.The time series of D will be obtained according to formula (2), and attenuation rate monthly be done on average, obtains monthly attenuation rate D
mean, then carry out matching with t, result as shown in Figure 2.Attenuation rate exponentially property change is in time found, correlativity (R after matching
2) reaching 0.96, fit equation is shown in formula (4).These are different from the analysis result of Tan et al. (2008).Ultraviolet light radiation meter is within shorter service time, and attenuation rate can be approximately linear attenuation. and time long when deployed, rate of decay becomes large gradually.
(4)
According to fitting result, obtaining calibration formula is:
(5)
Wherein
for revised radiation value.
This method is tested:
Consult following table: utilize formula (4) errata A data, by table B data, it is tested.In order to avoid other factors affect, select totally 21 d of the ceiling unlimited day after in July, 2010.During this period of time, the attenuation rate of new table is lower.Every day from 5 h to 20 h totally 16 time time ultraviolet irradiation degrees of data.Fig. 3 represent correct before and correct the scatter diagram of rear Table A and table B data.Fig. 3 shows, and the data dependence of two tables is very high, R
2up to 0.949, the aging comparatively even of instrument is described, the use value of instrument can be improved by the method corrected.Because Table A in use also can be decayed, so also use formula (5) to correct it.Irradiance average before and after Table A and table B being corrected is listed in table 1. table 1 and is shown, and the irradiance average before and after table B corrects is more or less the same, this is because table B is shorter for service time, decays less; Compared with table B, correct 76.8% of the value average out to table B of the uv radiation intensity recorded of front Table A, attenuation rate reaches 23.2%; After correcting, the measured value of Table A is about 99.7% of table B, differs minimum.
Radiation mean value (w/ (um m before and after table 1 Table A and table B correct
2))
| Before Table A is corrected | After Table A is corrected | Before table B corrects | After table B corrects |
| 12.47 | 16.19 | 16.16 | 16.24 |
The calibration formula of ultraviolet light radiation meter:
, wherein
for revised radiation value.
Claims (1)
1. the demarcation bearing calibration of a ultraviolet light radiation meter; It corrects in following condition: ground receiver to UV radiation comprise direct ultraviolet radiation and scattering UV radiation; Direct ultraviolet radiation is made up of molecular scattering, gas absorption delustring and aerosol extinction three part, and scattering UV radiation mainly comprises molecular scattering, atmospheric aerosol scattering and cloud and mist scattering; Its bearing calibration is as follows:
If the true uv radiation intensity on the arrival ground of certain time is
, the uv radiation intensity that pattern exports is
, then obtain following formula:
(1)
With the uv radiation intensity of aging ultraviolet light radiation meter observation be
, then attenuation rate can be expressed as:
(2)
Attenuation rate D is the function with time correlation, that is:
(3)
Wherein D is month by month attenuation rate, and t is that instrument uses month;
Select ceiling unlimited day totally 85 D, adopt the uv radiation intensity of 106 time points; Obtain the uv radiation intensity of pattern output by TUV mode computation, by the sample mean in every month, then carry out matching with the time; Obtain the time series of D according to formula (2), attenuation rate is monthly done on average, obtain monthly attenuation rate D
mean, then carry out matching with t, the attenuation rate after matching is exponentially property change in time, correlativity (R
2) reaching 0.96, fit equation is shown in formula (4); Ultraviolet light radiation meter is within shorter service time, and attenuation rate can be approximately linear attenuation. and time long when deployed, rate of decay becomes large gradually;
(4)
According to fitting result, obtaining demarcation updating formula is:
(5)
Wherein
for revised radiation value.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410580265.3A CN104297816B (en) | 2014-10-27 | 2014-10-27 | The demarcation bearing calibration of ultraviolet light radiation meter |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410580265.3A CN104297816B (en) | 2014-10-27 | 2014-10-27 | The demarcation bearing calibration of ultraviolet light radiation meter |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104297816A true CN104297816A (en) | 2015-01-21 |
| CN104297816B CN104297816B (en) | 2017-06-13 |
Family
ID=52317618
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410580265.3A Expired - Fee Related CN104297816B (en) | 2014-10-27 | 2014-10-27 | The demarcation bearing calibration of ultraviolet light radiation meter |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104297816B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106485002A (en) * | 2016-10-13 | 2017-03-08 | 云南省农业科学院甘蔗研究所 | Estimate solar radiation and the method for Caulis Sacchari sinensis potential production in complicated landform climatic province |
| CN110411982A (en) * | 2019-07-04 | 2019-11-05 | 云南省水利水电勘测设计研究院 | A method of surface albedo is estimated using surface observing data |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4627284A (en) * | 1985-07-22 | 1986-12-09 | Spectral Sciences, Inc. | Ultraviolet absorption hygrometer |
| JP2005147953A (en) * | 2003-11-18 | 2005-06-09 | Asaza Kikin | Environmental data instrumentation system, method, and program, summarizing server and sensor terminal used for the environmental data instrumentation system |
| CN100406912C (en) * | 2004-04-23 | 2008-07-30 | 富士施乐株式会社 | Method of measuring ultraviolet radiation and ultraviolet measuring device |
| CN100582692C (en) * | 2007-02-22 | 2010-01-20 | 富士通株式会社 | Apparatus having an input device and a display, method of controlling apparatus and recording medium |
-
2014
- 2014-10-27 CN CN201410580265.3A patent/CN104297816B/en not_active Expired - Fee Related
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4627284A (en) * | 1985-07-22 | 1986-12-09 | Spectral Sciences, Inc. | Ultraviolet absorption hygrometer |
| JP2005147953A (en) * | 2003-11-18 | 2005-06-09 | Asaza Kikin | Environmental data instrumentation system, method, and program, summarizing server and sensor terminal used for the environmental data instrumentation system |
| CN100406912C (en) * | 2004-04-23 | 2008-07-30 | 富士施乐株式会社 | Method of measuring ultraviolet radiation and ultraviolet measuring device |
| CN100582692C (en) * | 2007-02-22 | 2010-01-20 | 富士通株式会社 | Apparatus having an input device and a display, method of controlling apparatus and recording medium |
Non-Patent Citations (1)
| Title |
|---|
| 王炳忠等: "宽带紫外辐射表校准方法的新探索", 《太阳能学报》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106485002A (en) * | 2016-10-13 | 2017-03-08 | 云南省农业科学院甘蔗研究所 | Estimate solar radiation and the method for Caulis Sacchari sinensis potential production in complicated landform climatic province |
| CN106485002B (en) * | 2016-10-13 | 2019-04-19 | 云南省农业科学院甘蔗研究所 | In the method for complicated landform climatic province estimation sugarcane potential production |
| CN110411982A (en) * | 2019-07-04 | 2019-11-05 | 云南省水利水电勘测设计研究院 | A method of surface albedo is estimated using surface observing data |
| CN110411982B (en) * | 2019-07-04 | 2021-09-24 | 云南省水利水电勘测设计研究院 | Method for estimating earth surface albedo by using ground meteorological observation data |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104297816B (en) | 2017-06-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Gueymard | Interdisciplinary applications of a versatile spectral solar irradiance model: A review | |
| Pei et al. | Evolution of light use efficiency models: Improvement, uncertainties, and implications | |
| Gui et al. | Evaluation of radiosonde, MODIS-NIR-Clear, and AERONET precipitable water vapor using IGS ground-based GPS measurements over China | |
| Liu et al. | Ultraviolet radiation over China: Spatial distribution and trends | |
| Trnka et al. | Global solar radiation in Central European lowlands estimated by various empirical formulae | |
| Kaskaoutis et al. | Investigation into the wavelength dependence of the aerosol optical depth in the Athens area | |
| Wohlfahrt et al. | Quantifying nighttime ecosystem respiration of a meadow using eddy covariance, chambers and modelling | |
| Li et al. | Monthly ratios of PAR to global solar radiation measured at northern Tibetan Plateau, China | |
| Finch et al. | Photosynthetically active radiation regimes in a southern African savanna environment | |
| Zakey et al. | Atmospheric turbidity over Egypt | |
| Ide et al. | Assessment of canopy photosynthetic capacity and estimation of GPP by using spectral vegetation indices and the light–response function in a larch forest | |
| Sano et al. | Carbon dioxide exchange of a larch forest after a typhoon disturbance | |
| CN115730176A (en) | Ultraviolet hyperspectral atmospheric ozone vertical column concentration inversion method | |
| Hu et al. | Variation characteristics of ultraviolet radiation derived from measurement and reconstruction in Beijing, China | |
| Hu et al. | Ultraviolet radiation spatio-temporal characteristics derived from the ground-based measurements taken in China | |
| CN104297816A (en) | Calibration correcting method for ultraviolet radiation meter | |
| Serrano et al. | Daily and annual variations of erythemal ultraviolet radiation in Southwestern Spain | |
| García et al. | Comparison of observed and modeled cloud-free longwave downward radiation (2010–2016) at the high mountain BSRN Izaña station | |
| Malinovic-Milicevic et al. | Reconstruction of the erythemal UV radiation data in Novi Sad (Serbia) using the NEOPLANTA parametric model | |
| Bosch et al. | Global and direct photosynthetically active radiation parameterizations for clear-sky conditions | |
| 原薗芳信 et al. | Applications of MODIS-visible bands index, greenery ratio to estimate CO2 budget of a rice paddy in Japan | |
| Emetere | Gaussian algorithm for retrieving and projecting aerosols optical depth: A case study of Monrovia-Liberia | |
| Hu et al. | Variation characteristics of ultraviolet radiation over the North China Plain | |
| Pinker et al. | Aerosol radiative properties in the semiarid Western United States | |
| Gao et al. | Comparative analyses of the ultraviolet-B flux over the continental United State based on the NASA total ozone mapping spectrometer data and USDA ground-based measurements |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170613 Termination date: 20201027 |