Troposphere ozone column concentration calculation method and troposphere ozone column concentration calculation system

Abstract

The invention discloses a troposphere ozone column concentration calculation method and a troposphere ozone column concentration calculation system, and belongs to the technical field of environmental monitoring, and the method specifically comprises the following steps: obtaining ozone mass mixing ratio data of each isobaric surface and troposphere top pressure data, and unifying the resolution of the two data to obtain a combined data set; and calculating the ozone column concentration between the adjacent isobaric surfaces, and summing the column concentrations below the top height of the troposphere to obtain the ozone column concentration of the troposphere. According to the troposphere ozone column concentration calculation method and the troposphere ozone column concentration calculation system, the ozone mass mixing ratio data obtained from the fifth generation atmosphere reanalysis data set of global climate of the European mid-term weather forecast center and the troposphere top pressure data obtained from the reanalysis data set of the National Environmental Forecasting Center / National Atmosphere Research Center are adopted; the dependence on a data source is reduced, and the calculation precision of the troposphere ozone column concentration is improved.

Description

Technical Field

[0001] The present invention relates to the technical field of environmental monitoring, and in particular to a method and system for calculating tropospheric ozone column concentration. Background Art

[0002] The tropospheric ozone column concentration is one of the important indicators of environmental monitoring. The tropospheric ozone column concentration is used to evaluate air quality, climate change research and weather forecasting. Monitoring ozone concentration helps to deeply understand these chemical processes and provide data support for atmospheric chemistry research. After searching, the patent with announcement number CN109799549B discloses a method and device for converting gaseous pollutant concentration column. The four-dimensional grid data forecast results of physical quantities are obtained through model forecasting. In the four-dimensional grid data forecast results, the four-dimensional grid forecast results of preset gaseous pollutant concentrations are screened. According to the four-dimensional grid forecast results of preset gaseous pollutant concentrations, the preset layer preset gaseous pollutant concentration text file is obtained through meteorological mapping system processing. The layer-by-layer text file is used as the input field, and the model layer data is converted to column concentration according to the summation method to obtain the tropospheric ozone concentration column concentration. Through conversion, the model layer ozone concentration can be converted into the tropospheric ozone column concentration, reducing the ozone concentration from four dimensions to three dimensions, improving the forecast efficiency of tropospheric ozone concentration in atmospheric chemistry models, and enhancing the application effect of scientific research and business units in the environmental protection and meteorological industries. However, the tropospheric ozone column concentration calculated using a single data source is highly dependent on the data source, and the accuracy cannot be guaranteed. Summary of the invention

[0003] The purpose of the present invention is to provide a method and system for calculating tropospheric ozone column concentration to solve the above technical problems.

[0004] To achieve the above object, the present invention provides a method for calculating the tropospheric ozone column concentration, the specific steps are as follows:

[0005] Step S1: obtaining ozone mass mixing ratio data and tropopause pressure data of each isobaric surface, and performing interpolation data processing on the two data to unify the horizontal resolution of the tropopause pressure data and the ozone mass mixing ratio data to obtain a combined data set;

[0006] Step S2: extracting data with the same temporal frequency and the same spatial range in the combined data set;

[0007] Step S3: Calculate the ozone column concentration between adjacent isobaric surfaces based on the extracted data, and sum the ozone column concentrations between isobaric surfaces below the tropopause height to obtain the tropospheric ozone column concentration.

[0008] Preferably, in step S1, ozone mass mixing ratio data are obtained from the fifth generation atmospheric reanalysis data set of global climate of the European Centre for Medium-Range Weather Forecasts, and the time parameters in the data are converted into a YYYYMM format, and the ozone mass mixing ratio data include ozone mass mixing ratios at fixed time intervals on different isobaric surfaces;

[0009] The tropopause pressure data are obtained from the National Centers for Environmental Prediction / National Center for Atmospheric Research reanalysis dataset, and the time parameters in the data are converted into the YYYYMM format. The tropopause pressure data are the atmospheric tropopause pressure at fixed time intervals.

[0010] Preferably, the data resolution is converted by using an interpolation method so that the horizontal resolution of the tropospheric top pressure data and the ozone mass mixing ratio data remain consistent;

[0011] The tropopause pressure data were inverted to make the order of the pressure layer dimensions in the tropopause pressure data and the ozone mass mixing ratio data consistent.

[0012] Preferably, step S3 is specifically as follows:

[0013] Step S31: Convert the ozone mass mixing ratio of each isobaric surface in the combined data set to obtain a conversion factor CONV k ;

[0014] Step S32: Calculate the pressure difference ΔP between adjacent air layers below the tropopause height based on the tropopause pressure data in the combined data set k ;

[0015] Step S33: calculating the ozone column concentration of adjacent isobaric surfaces;

[0016] The calculation formula is as follows:

[0017] C k =CONV k *ΔP k *O 3k (:, i, j);

[0018] Among them, C k is the ozone column concentration between the kth adjacent isobaric surfaces, O 3k (:, i, j) is the ozone concentration in the set space, where i represents longitude and j represents latitude;

[0019] Step S34: Calculate the tropospheric ozone column concentration;

[0020] The calculation formula is as follows:

[0021]

[0022] Where DU is the tropospheric ozone column concentration, and n is the number of pressure isobaric surfaces below the tropospheric height.

[0023] A system based on the above-mentioned method for calculating tropospheric ozone column concentration includes:

[0024] A data acquisition module, the data acquisition module includes a first acquisition submodule and a second acquisition submodule for respectively acquiring ozone mass mixing ratio data and tropopause pressure data;

[0025] A data conversion module, the data conversion module includes a time parameter conversion submodule, a resolution conversion submodule and a dimension conversion submodule connected in sequence, the time parameter conversion submodule is used to convert the time parameter in the data into a YYYYMM format, the resolution conversion submodule is used to perform resolution conversion on the tropopause pressure data, and the dimension conversion submodule is used to perform dimension conversion on the tropopause pressure data;

[0026] A data calculation module, the data calculation module includes a pressure layer ozone column concentration calculation submodule and a tropospheric ozone column concentration calculation submodule, and the tropospheric ozone column concentration calculation submodule is connected to the pressure layer ozone column concentration calculation submodule;

[0027] An output module, used for outputting the tropospheric ozone column concentration, the output module is connected with the tropospheric ozone column concentration calculation submodule;

[0028] The first acquisition submodule and the second acquisition submodule are both connected to the time parameter conversion submodule, and the time parameter conversion submodule is connected to the atmospheric pressure layer ozone column concentration calculation submodule.

[0029] Therefore, the present invention adopts the above-mentioned method and system for calculating the tropospheric ozone column concentration, which has the beneficial effects of: using the ozone mass mixing ratio data obtained from the fifth-generation atmospheric reanalysis data set of the European Center for Medium-Term Weather Forecasts global climate and the tropospheric top pressure data obtained from the reanalysis data set of the National Center for Environmental Prediction / National Center for Atmospheric Research of the United States, and processing the data at the same time so that the resolutions of the two data are consistent, reducing the dependence on the data source and improving the calculation accuracy of the tropospheric ozone column concentration.

[0030] The technical solution of the present invention is further described in detail below through the accompanying drawings and embodiments. BRIEF DESCRIPTION OF THE DRAWINGS

[0031] Figure 1 This is a flow chart of a method for calculating tropospheric ozone column concentration according to the present invention;

[0032] Figure 2 This is a block diagram of the system principle of the present invention. DETAILED DESCRIPTION

[0033] In the description of the present invention, it should be noted that the terms "upper", "lower", "inside", "outside", etc. indicate positions or positional relationships based on the positions or positional relationships shown in the accompanying drawings, or the positions or positional relationships in which the invented product is usually placed when in use. They are only for the convenience of describing the present invention and simplifying the description, and do not indicate or imply that the device or element referred to must have a specific orientation, be constructed and operated in a specific orientation, and therefore cannot be understood as a limitation on the present invention. In the description of the present invention, it should also be noted that, unless otherwise clearly specified and limited, the terms "setting", "installation", and "connection" should be understood in a broad sense, for example, it can be a fixed connection, a detachable connection, or an integral connection; it can be a mechanical connection or an electrical connection; it can be a direct connection, or it can be an indirect connection through an intermediate medium, or it can be a connection between the two elements. For those of ordinary skill in the art, the specific meanings of the above terms in the present invention can be understood according to specific circumstances.

[0034] The embodiments of the present invention are described in detail below in conjunction with the accompanying drawings.

[0035] like Figure 1 As shown, a method for calculating the tropospheric ozone column concentration, the specific steps are as follows:

[0036] Step S1: Obtain ozone mass mixing ratio data and tropopause pressure data of each isobaric surface, and perform interpolation data processing on the two data to unify the horizontal resolution of the tropopause pressure data and the ozone mass mixing ratio data to obtain a combined data set.

[0037] The ozone mass mixing ratio data are obtained from the fifth-generation atmospheric reanalysis dataset of the global climate of the European Centre for Medium-Range Weather Forecasts, and the time parameters in the data are converted into the YYYYMM format. The ozone mass mixing ratio data include the ozone mass mixing ratios at fixed time intervals on different isobaric surfaces.

[0038] The tropopause pressure data is obtained from the National Center for Environmental Prediction / National Center for Atmospheric Research reanalysis data set, and the time parameters in the data are converted to YYYYMM format. The tropopause pressure data is the atmospheric tropopause pressure at fixed time intervals. The interpolation method is used to convert the resolution so that the horizontal resolution of the tropopause pressure data and the ozone mass mixing ratio data are consistent. The data from two different sources and different resolutions are interpolated to unify the resolution for the main calculation. The tropopause pressure data is reversed to keep the order of the pressure layer dimensions in the tropopause pressure data and the ozone mass mixing ratio data consistent.

[0039] Step S2: extracting data with the same temporal frequency and the same spatial range in the combined data set.

[0040] Step S3: Calculate the ozone column concentration between adjacent isobaric surfaces based on the extracted data, and sum the ozone column concentrations between isobaric surfaces below the tropopause height to obtain the tropospheric ozone column concentration.

[0041] Step S3 is as follows:

[0042] Step S31: Convert the ozone mass mixing ratio of each isobaric surface in the combined data set to obtain a conversion factor CONV k ;

[0043] Step S32: Calculate the pressure difference ΔP between adjacent air layers below the tropopause height based on the tropopause pressure data in the combined data set k ;

[0044] Step S33: calculating the ozone column concentration of adjacent isobaric surfaces;

[0045] The calculation formula is as follows:

[0046] C k =CONV k *ΔP k *O 3k (:, i, j);

[0047] Among them, C k is the ozone column concentration between the kth adjacent isobaric surfaces, O 3k (:, i, j) is the ozone concentration in the set space, where i represents longitude and j represents latitude;

[0048] Step S34: Calculate the tropospheric ozone column concentration;

[0049]

[0050] Where DU is the tropospheric ozone column concentration, and n is the number of pressure isobaric surfaces below the tropospheric height.

[0051] like Figure 2 As shown, a system based on the above-mentioned method for calculating tropospheric ozone column concentration includes:

[0052] A data acquisition module, the data acquisition module includes a first acquisition submodule and a second acquisition submodule for respectively acquiring ozone mass mixing ratio data and tropopause pressure data;

[0053] A data conversion module, the data conversion module includes a time parameter conversion submodule, a resolution conversion submodule and a dimension conversion submodule connected in sequence, the time parameter conversion submodule is used to convert the time parameter in the data into a YYYYMM format, the resolution conversion submodule is used to perform resolution conversion on the tropopause pressure data, and the dimension conversion submodule is used to perform dimension conversion on the tropopause pressure data;

[0054] A data calculation module, the data calculation module includes a pressure layer ozone column concentration calculation submodule and a tropospheric ozone column concentration calculation submodule, and the tropospheric ozone column concentration calculation submodule is connected to the pressure layer ozone column concentration calculation submodule;

[0055] An output module, used for outputting the tropospheric ozone column concentration, the output module is connected with the tropospheric ozone column concentration calculation submodule;

[0056] The first acquisition submodule and the second acquisition submodule are both connected to the time parameter conversion submodule, and the time parameter conversion submodule is connected to the atmospheric pressure layer ozone column concentration calculation submodule.

[0057] Finally, it should be noted that the above embodiments are only used to illustrate the technical solution of the present invention rather than to limit it. Although the present invention has been described in detail with reference to the preferred embodiments, those skilled in the art should understand that they can still modify or replace the technical solution of the present invention with equivalents, and these modifications or equivalent replacements cannot cause the modified technical solution to deviate from the spirit and scope of the technical solution of the present invention.

Claims

1. A method for calculating tropospheric ozone column concentration, characterized in that: The specific steps are as follows: Step S1: obtaining ozone mass mixing ratio data and tropopause pressure data of each isobaric surface, and performing interpolation data processing on the two data to unify the horizontal resolution of the tropopause pressure data and the ozone mass mixing ratio data to obtain a combined data set; Step S2: extracting data with the same temporal frequency and the same spatial range in the combined data set; Step S3: Calculate the ozone column concentration between adjacent isobaric surfaces based on the extracted data, and sum the ozone column concentrations between isobaric surfaces below the tropopause height to obtain the tropospheric ozone column concentration.

2. The method for calculating tropospheric ozone column concentration according to claim 1, characterized in that: In step S1, The ozone mass mixing ratio data are obtained from the fifth generation atmospheric reanalysis data set of the global climate of the European Center for Medium-Range Weather Forecasts, and the time parameters in the data are converted into the YYYYMM format. The ozone mass mixing ratio data include the ozone mass mixing ratios at fixed time intervals on different isobaric surfaces; The tropopause pressure data are obtained from the National Centers for Environmental Prediction / National Center for Atmospheric Research reanalysis dataset, and the time parameters in the data are converted into the YYYYMM format. The tropopause pressure data are the atmospheric tropopause pressure at fixed time intervals.

3. A method for calculating tropospheric ozone column concentration according to claim 2, characterized in that: The interpolation method is used to convert the data resolution so that the horizontal resolution of the tropospheric top pressure data and the ozone mass mixing ratio data remain consistent; The tropopause pressure data were inverted to make the order of the pressure layer dimensions in the tropopause pressure data and the ozone mass mixing ratio data consistent.

4. A method for calculating tropospheric ozone column concentration according to claim 3, characterized in that: Step S3 is as follows: Step S31: Convert the ozone mass mixing ratio of each isobaric surface in the combined data set to obtain a conversion factor CONV k ; Step S32: Calculate the pressure difference ΔP between adjacent air layers below the tropopause height based on the tropopause pressure data in the combined data set k ; Step S33: calculating the ozone column concentration of adjacent isobaric surfaces; The calculation formula is as follows: C k =CONV k *ΔP k *O 3k (:,i,j); Among them, C k is the ozone column concentration between the kth adjacent isobaric surfaces, O 3k (:, i, j) is the ozone concentration in the set space, where i represents longitude and j represents latitude; Step S34: Calculate the tropospheric ozone column concentration; The calculation formula is as follows: Where DU is the tropospheric ozone column concentration, and n is the number of pressure isobaric surfaces below the tropospheric height.

5. A system based on the method for calculating tropospheric ozone column concentration according to claim 4, characterized in that: include: A data acquisition module, the data acquisition module includes a first acquisition submodule and a second acquisition submodule for respectively acquiring ozone mass mixing ratio data and tropopause pressure data; A data conversion module, the data conversion module includes a time parameter conversion submodule, a resolution conversion submodule and a dimension conversion submodule connected in sequence, the time parameter conversion submodule is used to convert the time parameter in the data into a YYYYMM format, the resolution conversion submodule is used to perform resolution conversion on the tropopause pressure data, and the dimension conversion submodule is used to perform dimension conversion on the tropopause pressure data; A data calculation module, the data calculation module includes a pressure layer ozone column concentration calculation submodule and a tropospheric ozone column concentration calculation submodule, and the tropospheric ozone column concentration calculation submodule is connected to the pressure layer ozone column concentration calculation submodule; An output module, used for outputting the tropospheric ozone column concentration, the output module is connected with the tropospheric ozone column concentration calculation submodule; The first acquisition submodule and the second acquisition submodule are both connected to the time parameter conversion submodule, and the time parameter conversion submodule is connected to the atmospheric pressure layer ozone column concentration calculation submodule.

Images