A method for constructing a parameterization scheme for sand or dust emission based on surface albedo

Abstract

The present invention provides a method for constructing a sand or dust emission parameterization scheme based on surface albedo, which belongs to the field of sand or dust numerical simulation. The method comprises: based on wind tunnel experiments, performing stress measurement and light radiation measurement on a study area, and constructing a frictional wind speed calculation model; dividing the study area into grids, obtaining the surface albedo of each grid by acquiring satellite remote sensing data of each grid, and calculating using the frictional wind speed calculation model to obtain the frictional wind speed of each grid; constructing a sand or dust emission formula, and substituting the frictional wind speed of each grid into the sand or dust emission formula to calculate a sand or dust emission inventory for each grid, thereby completing the construction of the sand or dust emission parameterization scheme. The present invention solves the problems that simulation results and measured results of existing sand or dust emission parameterization schemes are uncertain, and are subject to real-time limitations and spatial singleness limitations.

Description

Technical Field

[0001] The present invention belongs to the field of sand and dust numerical simulation, and in particular relates to a method for constructing a sand or dust emission parameterization scheme based on surface albedo. Background Art

[0002] Dust plays a crucial role in the Earth system, not only directly impacting the atmospheric radiation properties but also indirectly affecting the Earth's radiation balance and water cycle. Dust pollution can cause respiratory and cardiovascular diseases and even death. Residents living in dust-prone areas for extended periods of time are at a significantly increased risk of developing desert pneumoconiosis. In my country's arid and semi-arid regions, sandstorms occur frequently. Dust travels near the surface, increasing particulate matter concentrations, causing air pollution and reducing visibility. When dust rises to a certain height, it can be easily transported over long distances, affecting air quality along the transmission path.

[0003] Currently, numerous sand and dust parameterization schemes have been developed, supported by field observations, wind tunnel experiments, and wind-blown sand physics. These schemes are mostly based on the silhouette coverage of the windward side of the land. Through extensive field measurements and land surveys, different land types are identified and the frictional wind speeds for each type of land are derived. However, these parameterization schemes can only simulate fixed surface types, which are spatially fixed, temporally static, and have a single surface condition. Therefore, the simulated results have a certain degree of uncertainty compared to the measured results. Summary of the Invention

[0004] In response to the above-mentioned deficiencies in the prior art, the present invention provides a method for constructing a sand or dust emission parameterization scheme based on surface albedo, which solves the problems that the simulation results and measured results of the existing sand or dust emission parameterization schemes are uncertain, and are subject to real-time limitations and spatial uniformity limitations.

[0005] In order to achieve the above objectives, the technical solution adopted by the present invention is as follows: On the one hand, the present invention provides a method for constructing a sand or dust emission parameterization scheme based on surface albedo, comprising the following steps:

[0006] S1. Based on wind tunnel experiments and light radiation experiments, stress and light radiation measurements are conducted in the study area to construct a friction wind speed calculation model;

[0007] S2. Divide the study area into grids, obtain satellite remote sensing data for each grid, obtain the surface albedo of each grid, and use the frictional wind speed calculation model to calculate the frictional wind speed of each grid;

[0008] S3. Construct a dust emission formula and substitute the frictional wind speed of each grid into the dust emission formula to calculate the dust emission list of each grid and complete the construction of a parameterized sand or dust emission scheme.

[0009] The beneficial effects of the present invention are as follows: the present invention performs stress measurement based on wind tunnel experiments and performs light radiation measurement based on light radiation experiments, constructs a response relationship between surface albedo and frictional wind speed, and utilizes the high temporal and spatial resolution surface albedo provided by satellite remote sensing data to realize the construction of a sand or dust raising parameterization scheme based on surface albedo, improves the accuracy of the simulation results and the measured results of the parameterization scheme, and reduces the influence of real-time limitations and spatial singleness limitations on the sand or dust raising parameterization scheme.

[0010] Furthermore, the S1 includes the following steps:

[0011] S101. Based on wind tunnel experiments, measure the shear stress of the ground surface and the lateral stress of each roughness element to obtain the experimental friction wind speed;

[0012] S102. Based on the light radiation experiment, measure the light radiation of the study area to obtain the experimental surface albedo;

[0013] S103. Based on the parameterization scheme of silhouette cover, the silhouette cover is replaced by the surface albedo, and a frictional wind speed calculation model is constructed according to the experimental frictional wind speed and the experimental surface albedo.

[0014] Furthermore, the S103 includes the following steps:

[0015] S1031. Based on the parameterization scheme of silhouette coverage, extract the fixed nonlinear relationship between frictional wind speed and silhouette coverage;

[0016] S1032, performing normalization processing and inversion processing on the experimental surface albedo to obtain a processed surface albedo;

[0017] S1033. Construct a frictional wind speed calculation model based on the experimental frictional wind speed, the processed surface albedo, and the fixed nonlinear relationship between the frictional wind speed and the silhouette coverage.

[0018] Furthermore, the expression of the friction wind speed calculation model is as follows:

[0019] ;

[0020] in, represents the frictional wind speed, is the free wind speed at high altitude, represents the normalized surface albedo, A 、 B 、 C 、D Both represent model parameters.

[0021] The beneficial effects of the above further scheme are as follows: the present invention replaces the silhouette cover with the surface albedo to establish the relationship between the surface albedo and the frictional wind speed, providing a basis for the subsequent direct calculation of the frictional wind speed using the surface albedo; and simplifies the acquisition of the frictional wind speed, thereby improving the accuracy and real-time performance of the frictional wind speed.

[0022] Furthermore, the S2 includes the following steps:

[0023] S201, performing splicing, integration, and interpolation processing on multiple satellite tile data occupied by the study area to obtain satellite remote sensing data covering the entire study area;

[0024] S202, dividing the study area into grids and obtaining the surface reflectance of each grid;

[0025] S203 , calculating the frictional wind speed of each grid using a frictional wind speed calculation model based on the surface albedo of each grid, to obtain the frictional wind speed of each grid.

[0026] The beneficial effects of the above-mentioned further scheme are as follows: the present invention divides the study area into grids, provides data processing for satellite remote sensing data, and then calculates the friction wind speed of each grid, thereby realizing the construction of a sand or dust parameterization scheme for a larger geographical area, reducing the difficulty of data acquisition, and improving the accuracy, real-time and universality of the construction of the sand or dust parameterization scheme.

[0027] Furthermore, the S3 includes the following steps:

[0028] S301. Obtain the critical friction wind speed based on the soil particle size in each grid in the study area;

[0029] S302. Construct a dust emission formula based on the critical friction wind speed, the air density in each grid in the study area, the acceleration of gravity, the proportion of frozen soil, the snow cover, the proportion of dust particles of a specific size emitted into the air, and the proportion of clay.

[0030] S303: Substitute the frictional wind speed of each grid into the dust emission formula to calculate the dust emission inventory of each grid, and complete the construction of the sand or dust emission parameterization scheme.

[0031] Furthermore, the dust emission formula is expressed as follows:

[0032] ;

[0033] ;

[0034] in, Indicates the vertical emission of dust, Represents a constant, represents the air density in each grid in the study area, represents the acceleration due to gravity, represents the critical friction wind speed, represents the soil particle size in each grid in the study area, Indicates soil moisture content, represents the horizontal emission of dust, represents the proportion of frozen soil, represents the snow coverage, Indicates the proportion of dust particles of a specific size emitted into the air. Indicates the proportion of clay.

[0035] In order to achieve the above object, according to a second aspect of the present invention, there is provided a computer device, characterized in that it includes a memory for storing a computer program;

[0036] The processor is configured to implement the steps of the method for constructing a parameterized sand or dust emission scheme based on surface albedo as described above when executing the computer program.

[0037] In order to achieve the above-mentioned purpose, according to the third aspect of the present invention, a computer-readable storage medium is provided, characterized in that a computer program is stored on the computer-readable storage medium, and when the computer program is executed by a processor, the steps of the method for constructing a sand or dust emission parameterization scheme based on surface albedo as described above are implemented.

[0038] The beneficial effects of the above further scheme are as follows: the present invention constructs a dust emission formula and obtains a dust emission inventory based on surface albedo through the frictional wind speed of each grid; compared with the previous emission inventory based on silhouette cover, the present invention improves the accuracy and reliability of the simulation results. BRIEF DESCRIPTION OF THE DRAWINGS

[0039] Figure 1 Flow chart of the method of the present invention. DETAILED DESCRIPTION

[0040] The specific embodiments of the present invention are described below to facilitate understanding of the present invention by those skilled in the art. However, it should be clear that the present invention is not limited to the scope of the specific embodiments. For those skilled in the art, as long as various changes are within the spirit and scope of the present invention as defined and determined by the appended claims, these changes are obvious, and all inventions and creations utilizing the concepts of the present invention are protected.

[0041] Example 1

[0042] like Figure 1As shown, the present invention provides a method for constructing a sand or dust emission parameterization scheme based on surface albedo, and its implementation method is as follows:

[0043] S1. Based on wind tunnel experiments and light radiation experiments, stress and light radiation measurements are performed on the study area to build a friction wind speed calculation model. The specific steps are as follows:

[0044] S101. Based on wind tunnel experiments, measure the shear stress of the ground surface and the lateral stress of each roughness element to obtain the experimental friction wind speed;

[0045] S102. Based on the light radiation experiment, measure the light radiation of the study area to obtain the experimental surface albedo;

[0046] S103. Based on the parameterization scheme of silhouette cover, the silhouette cover is replaced by the surface albedo. According to the experimental frictional wind speed and the experimental surface albedo, a frictional wind speed calculation model is constructed. The specific steps are as follows:

[0047] S1031. Based on the parameterization scheme of silhouette coverage, extract the fixed nonlinear relationship between frictional wind speed and silhouette coverage;

[0048] S1032, performing normalization processing and inversion processing on the experimental surface albedo to obtain a processed surface albedo;

[0049] S1033. Construct a frictional wind speed calculation model based on the experimental frictional wind speed, the processed surface albedo, and the fixed nonlinear relationship between the frictional wind speed and the silhouette coverage.

[0050] In this embodiment, experiments are conducted in a wind tunnel to measure surface stress, including surface shear stress and lateral stress of each roughness element, to obtain experimental frictional wind speed. Furthermore, for different surface environments in the study area, light radiation measurements are performed to obtain experimental surface albedo.

[0051] Based on the parameterization scheme of silhouette coverage, the fixed nonlinear relationship between frictional wind speed and silhouette coverage is extracted; the expression of the parameterization scheme of silhouette coverage is as follows:

[0052] ;

[0053] ;

[0054] in, represents the magnitude of the surface shear stress, 、 、 as well as All represent constants, represents the air density, represents the frictional wind speed, is the free wind speed at high altitude, Indicates the silhouette coverage, that is, the light from At angular incidence, the ratio of the projection of the surface obstacle on the ground to the entire surface area can be calculated from the above formula: and silhouette coverage There is a fixed nonlinear relationship between them, and the fixed nonlinear relationship between friction wind speed and silhouette coverage is extracted;

[0055] The experimental surface albedo is normalized and inverted to obtain a processed surface albedo; the inversion process, i.e., obtaining the non-reflective portion of the surface, is used to process the fixed nonlinear relationship between the frictional wind speed and the silhouette cover by using the experimental processed surface albedo (dark sky surface albedo), and a frictional wind speed calculation model is constructed based on the experimental frictional wind speed;

[0056] The expression of the friction wind speed calculation model is as follows:

[0057] ;

[0058] in, Represents silhouette coverage, using normalized surface albedo Replace silhouette coverage , and organize the original expression to obtain the following formula:

[0059] ;

[0060] in, represents the normalized surface albedo, A 、 B 、 C 、 D All represent model parameters, which are obtained by parameter estimation based on experimental data; in this embodiment, their values are: , , , .

[0061] S2. Divide the study area into grids. Obtain satellite remote sensing data for each grid to obtain the surface albedo of each grid. Calculate the frictional wind speed using the frictional wind speed calculation model to obtain the frictional wind speed of each grid. The specific steps are as follows:

[0062] S201, performing splicing, integration, and interpolation processing on multiple satellite tile data occupied by the study area to obtain satellite remote sensing data covering the entire study area;

[0063] S202, dividing the study area into grids and obtaining the surface reflectance of each grid;

[0064] S203 , calculating the frictional wind speed of each grid using a frictional wind speed calculation model based on the surface albedo of each grid, to obtain the frictional wind speed of each grid.

[0065] In this embodiment, multiple satellite tile data covering the study area are spliced, integrated, and interpolated to obtain satellite remote sensing data covering the entire study area. The study area is then gridded, and the satellite remote sensing data of each grid is selected to obtain the surface albedo of each grid. The frictional wind speed calculation model is then used to calculate the frictional wind speed of each grid.

[0066] Since the surface albedo can be directly obtained from satellite remote sensing data, the method of the present invention is simpler to implement and has real-time and universal applicability.

[0067] S3. Construct a dust emission formula and incorporate the frictional wind speed of each grid into the dust emission formula to calculate the dust emission inventory for each grid and complete the construction of a parameterized sand or dust emission scheme. The specific steps are as follows:

[0068] S301. Obtain the critical friction wind speed based on the soil particle size in each grid in the study area;

[0069] S302. Construct a dust emission formula based on the critical friction wind speed, the air density in each grid in the study area, the acceleration of gravity, the proportion of frozen soil, the snow cover, the proportion of dust particles of a specific size emitted into the air, and the proportion of clay.

[0070] S303: Substitute the frictional wind speed of each grid into the dust emission formula to calculate the dust emission inventory of each grid, and complete the construction of the sand or dust emission parameterization scheme.

[0071] In this embodiment, the critical friction wind speed is obtained based on the soil particle size in each grid in the study area. When studying PM10, 10 Based on the MB (Marticorena-Bergametti) dust emission scheme, a dust emission formula is constructed according to the critical friction wind speed, the air density in each grid in the study area, and the gravitational acceleration. The specific expression of the dust emission formula is as follows:

[0072] ;

[0073] ;

[0074] in, Indicates the vertical emission of dust, Represents a constant, represents the air density in each grid in the study area, represents the acceleration due to gravity, represents the critical friction wind speed, represents the soil particle size in each grid in the study area, Indicates soil moisture content, Indicates the horizontal emission of dust (unit ), represents the proportion of frozen soil, represents the snow coverage, Indicates the proportion of dust particles of a specific size emitted into the air. Indicates the proportion of clay;

[0075] Substitute the frictional wind speed of each grid into the dust emission formula to calculate the vertical and horizontal emissions of dust for each grid and obtain the dust emission list.

[0076] This paper intends to use surface albedo instead of silhouette cover to improve the accuracy of frictional wind speed. Through wind tunnel experiments and light radiation experiments, the response relationship between surface albedo and frictional wind speed is constructed. The high temporal and spatial resolution surface albedo provided by satellite remote sensing data is used to establish a spatially heterogeneous and temporally dynamic frictional wind speed dataset, and complete the construction of a sand or dust emission parameterization scheme based on surface albedo.

[0077] Example 2

[0078] In this embodiment, the frictional wind speed is obtained based on the surface albedo, and the atmospheric particulate matter concentration is further simulated.

[0079] The PM10 concentration in Tianjin for each month of 2023 was simulated using the method of the present invention. The comparison between the predicted results and the measured values is shown in Table 1.

[0080] Table 1

[0081]

[0082] According to the data in Table 1, in practice, the maximum error of the method of the present invention does not exceed 30%, and it is accurate and reliable. Moreover, due to its simple implementation, it does not require a large amount of field measurement work. With sufficient computing power, it can also realize the real-time prediction of atmospheric particulate matter concentration.

[0083] Example 3

[0084] In this embodiment, a computer device is provided, characterized in that it includes a memory for storing a computer program;

[0085] A processor is configured to implement the steps of the method for constructing a parameterized sand or dust emission scheme based on surface albedo as described in Example 1 when executing the computer program.

[0086] In this embodiment, a computer-readable storage medium is also provided, characterized in that a computer program is stored on the computer-readable storage medium, and when the computer program is executed by a processor, the steps of the method for constructing a sand or dust emission parameterization scheme based on surface albedo as described in Example 1 are implemented.

Claims

1. A method for constructing a parameterized sand or dust emission scheme based on surface albedo, characterized in that: The following steps are involved: S1. Based on wind tunnel experiments and light radiation experiments, stress and light radiation measurements are conducted in the study area to construct a friction wind speed calculation model; S2. Divide the study area into grids, obtain satellite remote sensing data for each grid, obtain the surface albedo of each grid, and use the frictional wind speed calculation model to calculate the frictional wind speed of each grid; S3. Construct a dust emission formula and incorporate the frictional wind speed of each grid into the dust emission formula to calculate the dust emission inventory for each grid, completing the construction of a parameterized sand or dust emission scheme. The expression of the friction wind speed calculation model is as follows: in, represents the frictional wind speed, is the free wind speed at high altitude, represents the normalized surface albedo, A 、 B 、 C 、 D All represent model parameters; The dust emission formula is as follows: in, Indicates the vertical emission of dust, Represents a constant, represents the air density in each grid in the study area, represents the acceleration due to gravity, represents the critical friction wind speed, represents the soil particle size in each grid in the study area, Indicates soil moisture content, represents the horizontal emission of dust, represents the proportion of frozen soil, represents the snow coverage, Indicates the proportion of dust particles of a specific size emitted into the air. Indicates the proportion of clay.

2. The method for constructing a parameterized sand or dust emission scheme based on surface albedo according to claim 1, characterized in that: Said S1 comprises the following steps: S101. Based on wind tunnel experiments, measure the shear stress of the ground surface and the lateral stress of each roughness element to obtain the experimental friction wind speed; S102. Based on the light radiation experiment, measure the light radiation of the study area to obtain the experimental surface albedo; S103. Based on the parameterization scheme of silhouette cover, the silhouette cover is replaced by the surface albedo, and a frictional wind speed calculation model is constructed according to the experimental frictional wind speed and the experimental surface albedo.

3. The method for constructing a parameterized sand or dust emission scheme based on surface albedo according to claim 2, characterized in that: The S103 includes the following steps: S1031. Based on the parameterization scheme of silhouette coverage, extract the fixed nonlinear relationship between frictional wind speed and silhouette coverage; S1032, performing normalization processing and inversion processing on the experimental surface albedo to obtain a processed surface albedo; S1033. Construct a frictional wind speed calculation model based on the experimental frictional wind speed, the processed surface albedo, and the fixed nonlinear relationship between the frictional wind speed and the silhouette coverage.

4. The method for constructing a parameterized sand or dust emission scheme based on surface albedo according to claim 1, wherein: The S2 comprises the following steps: S201, performing splicing, integration, and interpolation processing on multiple satellite tile data occupied by the study area to obtain satellite remote sensing data covering the entire study area; S202, dividing the study area into grids and obtaining the surface albedo of each grid; S203 , calculating the frictional wind speed of each grid using a frictional wind speed calculation model based on the surface albedo of each grid, to obtain the frictional wind speed of each grid.

5. The method for constructing a parameterized sand or dust emission scheme based on surface albedo according to claim 1, wherein: The S3 includes the following steps: S301. Obtain the critical friction wind speed based on the soil particle size in each grid in the study area; S302. Construct a dust emission formula based on the critical friction wind speed, the air density in each grid in the study area, the acceleration of gravity, the proportion of frozen soil, the snow cover, the proportion of dust particles of a specific size emitted into the air, and the proportion of clay. S303: Substitute the frictional wind speed of each grid into the dust emission formula to calculate the dust emission inventory of each grid, and complete the construction of the sand or dust emission parameterization scheme.

6. A computer device, characterized in that: including a memory for storing a computer program; A processor, configured to implement the steps of the method for constructing a sand or dust emission parameterization scheme based on surface albedo as described in any one of claims 1 to 5 when executing the computer program.

7. A computer-readable storage medium, characterized in that The computer-readable storage medium stores a computer program, which, when executed by a processor, implements the steps of the method for constructing a sand or dust emission parameterization scheme based on surface albedo as described in any one of claims 1 to 5.

Images