CN114088593A - Method and device for determining sea salt discharge flux - Google Patents

Abstract

The invention provides a method and a device for determining sea salt discharge flux, and belongs to the field of environmental science. The method comprises the following steps: receiving a sea salt discharge flux calculation request, wherein the sea salt discharge flux calculation request at least comprises calculation scheme configuration parameters; determining a target calculation scheme according to the calculation scheme configuration parameters; acquiring the required meteorological information and geographic information according to the target calculation scheme; and determining sea salt discharge flux according to the target calculation scheme, the meteorological information and the geographic information. By adopting the invention, the multi-scale and multi-scene applicability of sea salt discharge flux simulation can be improved.

Description

Method and device for determining sea salt discharge flux

Technical Field

The invention relates to the field of environmental science, in particular to a method and a device for determining sea salt discharge flux.

Background

Atmospheric aerosols, also known as atmospheric particulates, are a generic term for solid or liquid particulates suspended in the atmosphere. The aerosol contains various components such as sulfate, nitrate, ammonium salt, black carbon, organic carbon, sand dust, sea salt, etc., which are closely related to atmospheric pollution and climate change, and are currently important research subjects in the fields of atmospheric environment and global climate change. Sea salt aerosols are aerosols of natural origin present in the atmosphere, which are mainly produced by foam explosions or wave impact breakup in spray. Sea salt aerosols are an important component of atmospheric aerosols, and have a significant impact on the concentration and composition of particulate matter (PM10 and PM2.5) in coastal regions. On a global and regional scale, sea salt contributes significantly to aerosol radiation intensity, playing an important role in atmospheric energy balance.

Numerical simulation is an important means for researching the digestion and evolution of sea salt aerosol in the atmosphere, and in order to realize the simulation of sea salt, a sea salt discharge flux parameterization scheme needs to be constructed firstly, and then the sea salt discharge flux parameterization scheme is coupled into an atmospheric chemical transmission mode to realize the simulation of the discharge, advection, diffusion, dry-wet sedimentation and the like of sea salt particles. In recent decades, scholars at home and abroad develop various sea salt discharge flux calculation schemes and develop simulation analysis research, the schemes are based on empirical formulas, and the sea salt flux is mainly related to the wind speed on the sea surface. Researches show that the calculation schemes can better simulate the basic distribution characteristics of sea salt, but the fineness, the calculation efficiency, the simulation applicability of different space-time scales and the like of various calculation schemes have respective limitations, and no optimal scheme can be suitable for various application scenes. In addition, most of these sea salt flux calculation schemes are empirical formulas obtained by analyzing and summarizing or fitting observation data, and parameters of the formulas have large uncertainties and need to be continuously verified, improved and localized.

In recent years, with the deep progress of air pollution control, the concentration of particulate matters generated by artificial emission in various regions around the world is continuously reduced. Under the condition of low particulate matter concentration, the influence of natural source aerosol on the global and regional particulate matter concentration is more prominent. Developed cities or urban communities around the world are mostly distributed in coastal areas, and sea salt aerosol is the most important natural aerosol source. Therefore, the research on the growth and digestion evolution of the sea salt aerosol and the influence of the sea salt aerosol on the concentration, the composition, the regional climate feedback and the like of particulate matters in coastal regions is of great significance.

Disclosure of Invention

In order to solve the problems in the prior art, the embodiment of the invention provides a method and a device for determining sea salt discharge flux. The technical scheme is as follows:

according to an aspect of the present invention, there is provided a method of determining sea salt discharge flux, the method comprising:

receiving a sea salt discharge flux calculation request, wherein the sea salt discharge flux calculation request at least comprises calculation scheme configuration parameters;

determining a target calculation scheme according to the calculation scheme configuration parameters;

acquiring the required meteorological information and geographic information according to the target calculation scheme;

and determining sea salt discharge flux according to the target calculation scheme, the meteorological information and the geographic information.

Optionally, the sea salt discharge flux calculation request further includes a space configuration parameter;

the acquiring of the required weather information and the geographic information according to the target calculation scheme comprises:

determining a plurality of geographical grids of a target area according to the spatial configuration parameters;

acquiring meteorological information and geographic information of each geographic grid according to the target calculation scheme;

said determining sea salt discharge flux from said target calculation scenario, said meteorological information, said geographic information, comprising:

determining the underlying surface attribute of each geographic grid according to the geographic information of each geographic grid;

calculating through the target calculation scheme based on the meteorological information and underlying surface attributes of each geographic grid, and determining sea salt discharge flux of each geographic grid;

and determining the sea salt discharge flux of the target area according to the sea salt discharge flux of each geographical grid.

Optionally, the target calculation scheme is a first calculation scheme, and the first calculation scheme is a calculation scheme based on geographic attributes; the underlying surface attribute used by the first computing scheme is a geographic attribute comprising a land attribute, a sea ice attribute, and a sea attribute;

said determining sea salt discharge flux for each said geographic grid calculated by said target calculation scheme based on said meteorological information and underlying surface properties for said each geographic grid comprises:

setting a sea salt discharge flux of a first geographic grid to a first preset value when a geographic attribute of the first geographic grid is a land attribute or a sea ice attribute;

and when the geographic attribute of the second geographic grid is the marine attribute, determining the sea salt discharge flux of the second geographic grid according to the meteorological information of the second geographic grid.

Optionally, the determining the sea salt discharge flux of the second geographic grid according to the weather information of the second geographic grid includes:

acquiring the wind speed and the particle radius of sea salt at the sea surface of the second geographic grid at 10 meters;

determining sea salt discharge flux of the second geographic grid by:

wherein F is sea salt discharge flux, r is the particle radius of sea salt, and U₁₀Wind speed at 10 meters above sea surface, B₁Is an empirical operator related to the radius of the particle.

Optionally, the target calculation scheme is a second calculation scheme, where the second calculation scheme is a calculation scheme based on geographic attributes and particle radii; the underlying surface attribute used by the second calculation scheme is a geographic attribute comprising a land attribute, a sea ice attribute, and a sea attribute;

said determining sea salt discharge flux for each said geographic grid calculated by said target calculation scheme based on said meteorological information and underlying surface properties for said each geographic grid comprises:

when the geographic attribute of the third geographic grid is a land attribute or a sea ice attribute, setting the sea salt discharge flux of the third geographic grid to be a second preset value;

when the geographic attribute of the fourth geographic grid is the marine attribute, respectively determining a first sea salt discharge flux F with the particle radius smaller than a first radius threshold value according to the meteorological information of the fourth geographic grid₁And a second sea salt discharge flux F having a particle radius equal to or greater than a first radius threshold₂；

Determining sea salt discharge flux for the fourth geographic grid by:

wherein F is the sea salt discharge flux and r is the particle radius of the sea salt.

Optionally, determining a first sea salt discharge flux F with a particle radius less than a first radius threshold based on meteorological information of said fourth geographic grid₁The method comprises the following steps:

acquiring the wind speed and the particle radius of sea salt at the sea surface of the fourth geographical grid of 10 meters;

determining said first sea salt discharge flux F by the following formula₁：

Wherein, U₁₀Wind speed at 10 meters above sea surface, B₂Is an empirical operator related to the radius of the particle.

Optionally, determining a second sea salt discharge flux F with a particle radius greater than or equal to a first radius threshold value according to meteorological information of the fourth geographic grid₂The method comprises the following steps:

acquiring the wind speed and the particle radius of sea salt at the sea surface of the fourth geographical grid of 10 meters;

determining the second sea salt discharge flux F by the following formula₂：

Wherein, U₁₀The wind speed at 10 meters above sea surface.

Optionally, the method further includes:

determining the hygroscopic increased sea salt emission flux of the fourth geographical grid by:

wherein M is the sea salt discharge flux after moisture absorption and growth, r_dIs the dry particle radius of sea salt, f_rhIs a moisture absorption growth factor, ρ₁Is the sea salt dry particle density.

Optionally, the target calculation scheme is a third calculation scheme, where the third calculation scheme is a calculation scheme based on geographic attributes, particle radii, and sea area attributes; the underlay surface attribute used by the third calculation scheme is a geographical attribute comprising a land attribute, an ice attribute and a sea attribute, wherein the sea attribute is a sea area attribute comprising a sea area and a sea area;

said determining sea salt discharge flux for each said geographic grid calculated by said target calculation scheme based on said meteorological information and underlying surface properties for said each geographic grid comprises:

when the geographic attribute of a fifth geographic grid is a land attribute or a sea ice attribute, setting the sea salt discharge flux of the fifth geographic grid to a third preset value;

when the geographic attribute of the sixth geographic grid is the marine attribute and the sea area attribute is the open sea area, respectively determining a third sea salt discharge flux F with the particle radius smaller than a second radius threshold value according to the meteorological information of the sixth geographic grid_aAnd a fourth sea salt discharge flux F having a particle radius equal to or greater than a second radius threshold_b(ii) a According to the third sea salt discharge flux F_aAnd fourth sea salt discharge flux F_bDetermining sea salt discharge flux F of said sixth geogrid_{Far away}；

When the geographic attribute of the seventh geographic grid is the marine attribute and the sea area attribute is the offshore area, determining the sea salt discharge flux F of the seventh geographic grid according to the meteorological information of the seventh geographic grid_{Near to}。

Optionally, determining a third sea salt discharge flux F with a particle radius less than a second radius threshold based on meteorological information of said sixth geographic grid_aThe method comprises the following steps:

acquiring the wind speed, the relative humidity and the particle radius of sea salt at the sea surface of the sixth geographical grid of 10 meters;

determining said third sea salt discharge flux F by the following formula_a：

A₁＝-5.001E³+0.808E⁶r_RHC₈₀-1.98E⁷(r_RHC₈₀)²+2.188E⁸(r_RHC₈₀)³ -1.144E⁹(r_RHC₈₀)⁴+2.29E⁹(r_RHC₈₀)⁵

A₂＝3.854E³+1.168E⁴r_RHC₈₀-6.572E⁴(r_RHC₈₀)²+1.003E⁵(r_RHC₈₀)³ -6.407E⁴(r_RHC₈₀)⁴+1.493E⁴(r_RHC₈₀)⁵

A₃＝4.498E²+0.839E³r_RHC₈₀-5.394E²(r_RHC₈₀)²+1.218E²(r_RHC₈₀)³ -1.213E¹(r_RHC₈₀)⁴+4.514E^-1(r_RHC₈₀)⁵

ρ₂＝1000×(3.8033-16.248RH+46.085RH²-68.317RH³+50.932RH⁴-15.261RH⁵)

x＝3.1657-19.079RH+55.72RH²-83.998RH³+63.436RH⁴-19.248RH⁵

Wherein A is_iThe value range of i is {1, 2, 3} which is an empirical coefficient; r is_iIs the particle radius of sea salt, r₀Has a value range of (0,0.1), r₁Has a value range of [0.1,1), r₂Has a value range of [1,2.5), r₃The value range of (1) is [2.5,5 ]; RH is relative humidity, r_RHIs the aerosol particle radius at RH humidity; u shape₁₀The wind speed at 10 meters above sea level; c₈₀Is a correction factor for salinity; AB₁The conversion coefficient of the aerosol concentration and the mass flux is obtained; rho₂Is sea salt particle density; x is the mass fraction of solute; eⁿIs 10ⁿ。

Optionally, determining a fourth sea salt emission flux F with a particle radius greater than or equal to a radius threshold according to the meteorological information of the sixth geographic grid_bThe method comprises the following steps:

acquiring the wind speed, the relative humidity and the particle radius of sea salt at the sea surface of the sixth geographical grid of 10 meters;

determining the fourth sea salt discharge flux F by the following formula_b：

ρ₂＝1000×(3.8033-16.248RH+46.085RH²-68.317RH³+50.932RH⁴-15.261RH⁵)

x＝3.1657-19.079RH+55.72RH²-83.998RH³+63.436RH⁴-19.248RH⁵

Wherein A is₄Is an empirical coefficient; r is₃、r₄Is the particle radius of sea salt, r₃Has a value range of [, ]2.5,5)，r₄Has a value range of [5,10 ]](ii) a RH is relative humidity, r_RHIs the particle radius of the sea salt at RH humidity; u shape₁₀The wind speed at 10 meters above sea level; c₈₀Is a correction factor for salinity; AB₁The conversion coefficient of the aerosol concentration and the mass flux is obtained; rho₂Is sea salt particle density; x is the mass fraction of solute; b is₃Is an empirical operator related to the radius of the particle.

Optionally, said determining sea salt emission flux F of said seventh geographical grid based on meteorological information of said seventh geographical grid_{Near to}The method comprises the following steps:

acquiring the wind speed, the relative humidity, the particle radius of sea salt and the salinity of sea water at the sea surface of the seventh geographical grid of 10 meters; determining sea salt discharge flux F of the seventh geogrid by_{Near to}：

AB₂＝2×10^-15×4.188×r_RH ³×ρ₂x

ρ₂＝1000×(3.8033-16.248RH+46.085RH²-68.317RH³+50.932RH⁴-15.261RH⁵)

x＝3.1657-19.079RH+55.72RH²-83.998RH³+63.436RH⁴-19.248RH⁵

Wherein, F_cFlux F for sea salt discharge_{Near to}RH is the relative humidity, r_RHIs the particle radius, U, of sea salt at RH humidity₁₀The wind speed at 10 meters above sea surface, S is the salinity of seawater,C₀as a correction factor for salinity, AB₂Is the conversion coefficient of aerosol concentration to mass flux, rho₂Is the sea salt particle density and x is the solute mass fraction.

Optionally, for the sixth geographic grid and the seventh geographic grid, the method further comprises:

acquiring preset proportion of multiple chemical components;

and determining the sea salt discharge flux corresponding to each chemical component according to the ratio of each chemical component and the sea salt discharge flux of the geographical grid.

Optionally, the method further includes:

simulating a physical and chemical process of sea salt aerosol according to the sea salt discharge flux;

and outputting the concentration distribution of the sea salt aerosol according to the calculation result of the physical and chemical process of the sea salt aerosol and the output frequency setting.

According to another aspect of the present invention, there is provided a sea salt discharge flux determining apparatus, the apparatus comprising:

the sea salt discharge flux calculation module is used for receiving a sea salt discharge flux calculation request, and the sea salt discharge flux calculation request at least comprises calculation scheme configuration parameters;

the acquisition module is used for determining a target calculation scheme according to the calculation scheme configuration parameters; acquiring the required meteorological information and geographic information according to the target calculation scheme;

and the determining module is used for determining sea salt discharge flux according to the target calculation scheme, the meteorological information and the geographic information.

Optionally, the sea salt discharge flux calculation request further includes a space configuration parameter;

the acquisition module is configured to: determining a plurality of geographical grids of a target area according to the spatial configuration parameters; acquiring meteorological information and geographic information of each geographic grid according to the target calculation scheme;

the determination module is to: determining the underlying surface attribute of each geographic grid according to the geographic information of each geographic grid; and calculating through the target calculation scheme based on the meteorological information and the underlying surface attribute of each geographic grid, and determining the sea salt discharge flux of each geographic grid.

Optionally, the target calculation scheme is a first calculation scheme, and the first calculation scheme is a calculation scheme based on geographic attributes; the underlying surface attribute used by the first computing scheme is a geographic attribute comprising a land attribute, a sea ice attribute, and a sea attribute;

the determination module is to:

setting a sea salt discharge flux of a first geographic grid to a first preset value when a geographic attribute of the first geographic grid is a land attribute or a sea ice attribute;

and when the geographic attribute of the second geographic grid is the marine attribute, determining the sea salt discharge flux of the second geographic grid according to the meteorological information of the second geographic grid.

Optionally, the determining module is configured to:

acquiring the wind speed and the particle radius of sea salt at the sea surface of the second geographic grid at 10 meters;

determining sea salt discharge flux of the second geographic grid by:

wherein F is sea salt discharge flux, r is the particle radius of sea salt, and U₁₀Wind speed at 10 meters above sea surface, B₁Is an empirical operator related to the radius of the particle.

Optionally, the target calculation scheme is a second calculation scheme, where the second calculation scheme is a calculation scheme based on geographic attributes and particle radii; the underlying surface attribute used by the second calculation scheme is a geographic attribute comprising a land attribute, a sea ice attribute, and a sea attribute;

the determination module is to:

when the geographic attribute of the third geographic grid is a land attribute or a sea ice attribute, setting the sea salt discharge flux of the third geographic grid to be a second preset value;

when the geographic attribute of the fourth geographic grid is the marine attribute, respectively determining a first sea salt discharge flux F with the particle radius smaller than a first radius threshold value according to the meteorological information of the fourth geographic grid₁And a second sea salt discharge flux F having a particle radius equal to or greater than a first radius threshold₂；

Determining sea salt discharge flux for the fourth geographic grid by:

wherein F is the sea salt discharge flux and r is the particle radius of the sea salt.

Optionally, the determining module is configured to:

acquiring the wind speed and the particle radius of sea salt at the sea surface of the fourth geographical grid of 10 meters;

determining said first sea salt discharge flux F by the following formula₁：

Wherein, U₁₀Wind speed at 10 meters above sea surface, B₂Is an empirical operator related to the radius of the particle.

Optionally, the determining module is configured to:

acquiring the wind speed and the particle radius of sea salt at the sea surface of the fourth geographical grid of 10 meters;

determining the second sea salt discharge flux F by the following formula₂：

Wherein, U₁₀The wind speed at 10 meters above sea surface.

Optionally, the determining module is further configured to:

determining the hygroscopic increased sea salt emission flux of the fourth geographical grid by:

wherein M is the sea salt discharge flux after moisture absorption and growth, r_dIs the dry particle radius of sea salt, f_rhIs a moisture absorption growth factor, ρ₁Is the sea salt dry particle density.

Optionally, the target calculation scheme is a third calculation scheme, where the third calculation scheme is a calculation scheme based on geographic attributes, particle radii, and sea area attributes; the underlay surface attribute used by the third calculation scheme is a geographical attribute comprising a land attribute, an ice attribute and a sea attribute, wherein the sea attribute is a sea area attribute comprising a sea area and a sea area;

the determination module is to:

when the geographic attribute of a fifth geographic grid is a land attribute or a sea ice attribute, setting the sea salt discharge flux of the fifth geographic grid to a third preset value;

when the geographic attribute of the sixth geographic grid is the marine attribute and the sea area attribute is the open sea area, respectively determining a third sea salt discharge flux F with the particle radius smaller than a second radius threshold value according to the meteorological information of the sixth geographic grid_aAnd a fourth sea salt discharge flux F having a particle radius equal to or greater than a second radius threshold_b(ii) a According to the third sea salt discharge flux F_aAnd fourth sea salt discharge flux F_bDetermining sea salt discharge flux F of said sixth geogrid_{Far away}；

When the geographic attribute of the seventh geographic grid is the marine attribute and the sea area attribute is the offshore area, determining the sea salt discharge flux F of the seventh geographic grid according to the meteorological information of the seventh geographic grid_{Near to}。

Optionally, the determining module is configured to:

acquiring the wind speed, the relative humidity and the particle radius of sea salt at the sea surface of the sixth geographical grid of 10 meters;

determining said third sea salt discharge flux F by the following formula_a：

A₁＝-5.001E³+0.808E⁶r_RHC₈₀-1.98E⁷(r_RHC₈₀)²+2.188E⁸(r_RHC₈₀)³ -1.144E⁹(r_RHC₈₀)⁴+2.29E⁹(r_RHC₈₀)⁵

A₂＝3.854E³+1.168E⁴r_RHC₈₀-6.572E⁴(r_RHC₈₀)²+1.003E⁵(r_RHC₈₀)³ -6.407E⁴(r_RHC₈₀)⁴+1.493E⁴(r_RHC₈₀)⁵

A₃＝4.498E²+0.839E³r_RHC₈₀-5.394E²(r_RHC₈₀)²+1.218E²(r_RHC₈₀)³ -1.213E¹(r_RHC₈₀)⁴+4.514E^-1(r_RHC₈₀)⁵

ρ₂＝1000×(3.8033-16.248RH+46.085RH²-68.317RH³+50.932RH⁴-15.261RH⁵)

x＝3.1657-19.079RH+55.72RH²-83.998RH³+63.436RH⁴-19.248RH⁵

Wherein A is_iThe value range of i is {1, 2, 3} which is an empirical coefficient; r is_iIs the particle radius of sea salt, r₀Has a value range of (0,0.1), r₁Has a value range of [0.1,1), r₂Has a value range of [1,2.5), r₃The value range of (1) is [2.5,5 ]; RH is relative humidity, r_RHIs the aerosol particle radius at RH humidity; u shape₁₀The wind speed at 10 meters above sea level; c₈₀Is a correction factor for salinity; AB₁The conversion coefficient of the aerosol concentration and the mass flux is obtained; rho₂Is sea salt particle density; x is the mass fraction of solute; eⁿIs 10ⁿ。

Optionally, the determining module is configured to:

acquiring the wind speed, the relative humidity and the particle radius of sea salt at the sea surface of the sixth geographical grid of 10 meters;

determining the fourth sea salt discharge flux F by the following formula_b：

ρ₂＝1000×(3.8033-16.248RH+46.085RH²-68.317RH³+50.932RH⁴-15.261RH⁵)

x＝3.1657-19.079RH+55.72RH²-83.998RH³+63.436RH⁴-19.248RH⁵

Wherein A is₄Is an empirical coefficient; r is₃、r₄Is the particle radius of sea salt, r₃Has a value range of [2.5,5), r₄Has a value range of [5,10 ]](ii) a RH is relative humidity, r_RHIs the particle radius of the sea salt at RH humidity; u shape₁₀The wind speed at 10 meters above sea level; c₈₀Is a correction factor for salinity; AB₁The conversion coefficient of the aerosol concentration and the mass flux is obtained; rho₂Is sea salt particle density; x is the mass fraction of solute; b is₃Is an empirical operator related to the radius of the particle.

Optionally, the determining module is configured to:

acquiring the wind speed, the relative humidity, the particle radius of sea salt and the salinity of sea water at the sea surface of the seventh geographical grid of 10 meters; determining sea salt discharge flux F of the seventh geogrid by_{Near to}：

AB₂＝2×10^-15×4.188×r_RH ³×ρ₂x

ρ₂＝1000×(3.8033-16.248RH+46.085RH²-68.317RH³+50.932RH⁴-15.261RH⁵)

x＝3.1657-19.079RH+55.72RH²-83.998RH³+63.436RH⁴-19.248RH⁵

Wherein, F_cFlux F for sea salt discharge_{Near to}RH is the relative humidity, r_RHIs the particle radius, U, of sea salt at RH humidity₁₀The wind speed at 10m above sea surface, S is the salinity of seawater, C₀As a correction factor for salinity, AB₂Is the conversion coefficient of aerosol concentration to mass flux, rho₂Is the sea salt particle density and x is the solute mass fraction.

Optionally, for the sixth geographic grid and the seventh geographic grid, the determining module is further configured to:

acquiring preset proportion of multiple chemical components;

and determining the sea salt discharge flux corresponding to each chemical component according to the ratio of each chemical component and the sea salt discharge flux of the geographical grid.

Optionally, the determining module is further configured to:

simulating a physical and chemical process of sea salt aerosol according to the sea salt discharge flux;

and outputting the concentration distribution of the sea salt aerosol according to the calculation result of the physical and chemical process of the sea salt aerosol and the output frequency setting.

According to another aspect of the present invention, there is provided an electronic apparatus including:

a processor; and

a memory for storing a program, wherein the program is stored in the memory,

wherein the program comprises instructions which, when executed by the processor, cause the processor to carry out the above method of determining sea salt discharge flux.

According to another aspect of the present invention, there is provided a non-transitory computer readable storage medium having stored thereon computer instructions for causing a computer to execute the above-described method of determining sea salt emission flux.

In the embodiment of the invention, the sea salt discharge flux calculation request can carry the calculation scheme configuration parameters, and the server can select the target calculation scheme through the calculation scheme configuration parameters so as to determine the sea salt discharge flux. The method realizes integration and improvement of a plurality of sea salt calculation schemes with different complexities under the same mode frame, can be suitable for multi-scale simulation of the whole world, the region, the city and the like, and simultaneously gives consideration to the calculation efficiency and the simulation precision.

Drawings

Further details, features and advantages of the invention are disclosed in the following description of exemplary embodiments with reference to the accompanying drawings, in which:

fig. 1 shows a flow chart of a sea salt discharge flux determination method according to an exemplary embodiment of the present invention;

FIG. 2 shows a flow chart of a sea salt discharge flux determination method according to an exemplary embodiment of the present invention;

FIG. 3 shows a flow chart of a sea salt discharge flux determination method according to an exemplary embodiment of the present invention;

fig. 4 shows a schematic block diagram of a sea salt discharge flux determining apparatus according to an exemplary embodiment of the present invention;

FIG. 5 illustrates a block diagram of an exemplary electronic device that can be used to implement an embodiment of the invention.

Detailed Description

Embodiments of the present invention will be described in more detail below with reference to the accompanying drawings. While certain embodiments of the present invention are shown in the drawings, it should be understood that the present invention may be embodied in various forms and should not be construed as limited to the embodiments set forth herein, but rather are provided for a more thorough and complete understanding of the present invention. It should be understood that the drawings and the embodiments of the present invention are illustrative only and are not intended to limit the scope of the present invention.

It should be understood that the various steps recited in the method embodiments of the present invention may be performed in a different order and/or performed in parallel. Moreover, method embodiments may include additional steps and/or omit performing the illustrated steps. The scope of the invention is not limited in this respect.

The term "include" and variations thereof as used herein are open-ended, i.e., "including but not limited to". The term "based on" is "based, at least in part, on". The term "one embodiment" means "at least one embodiment"; the term "another embodiment" means "at least one additional embodiment"; the term "some embodiments" means "at least some embodiments". Relevant definitions for other terms will be given in the following description. It should be noted that the terms "first", "second", and the like in the present invention are only used for distinguishing different devices, modules or units, and are not used for limiting the order or interdependence relationship of the functions performed by the devices, modules or units.

It is noted that references to "a", "an", and "the" modifications in the present invention are intended to be illustrative rather than limiting, and that those skilled in the art will recognize that reference to "one or more" unless the context clearly dictates otherwise.

The names of messages or information exchanged between devices in the embodiments of the present invention are for illustrative purposes only, and are not intended to limit the scope of the messages or information.

The embodiment of the invention provides a method for determining sea salt discharge flux, which can be completed by a terminal, a server and/or other equipment with processing capacity. The method provided by the embodiment of the present invention may be implemented by any one of the above devices, or may be implemented by a plurality of devices, for example, the terminal may determine the target calculation scheme according to the sea salt discharge flux calculation request, and then the server may determine the sea salt discharge flux through the target calculation scheme, that is, the calculation process is placed at the server side, which is not limited in the present invention.

Taking a terminal as an example, the sea salt discharge flux determining method will be described below with reference to the flow chart of the sea salt discharge flux determining method shown in fig. 1.

Step 101, a server receives a sea salt discharge flux calculation request.

Wherein the sea salt discharge flux calculation request may include at least a calculation scheme configuration parameter. Optionally, the sea salt discharge flux calculation request may further include a spatial configuration parameter and a temporal configuration parameter. The calculation scheme configuration parameters may include an identification of the calculation scheme. The spatial configuration parameters may include geographic location information and resolution. The time configuration parameter may include a point in time or a period of time.

In a possible implementation, the sea salt discharge flux determination method provided by the embodiment can be set in an atmospheric chemical transmission mode, and a user can use the mode through a terminal. When a user needs to simulate the transmission of chemical substances in the atmosphere, the simulation parameters used can be configured in the atmospheric chemical transmission mode of the terminal. For example, the user may select the calculation method used and select the chemical substances to be simulated (e.g., sulfate, nitrate, ammonium salts, black carbon, organic carbon, sand dust, sea salt, etc.) and set the geographical area and period to be simulated. In this embodiment, a calculation method of sea salt aerosol is mainly described, and a calculation method of other chemical substances is not described in this embodiment.

When the user confirms the configuration parameters, the terminal can package the configuration parameters to obtain a corresponding sea salt discharge flux calculation request and send the sea salt discharge flux calculation request to the server. Further, the server may receive a corresponding sea salt discharge flux calculation request.

And step 102, the server determines a target calculation scheme according to the calculation scheme configuration parameters.

In a possible implementation, for the sea salt aerosol, a plurality of calculation schemes may be stored in the server in advance, and the calculation complexity and the applicable scenario of each calculation scheme are different. Different target calculation schemes have different calculation complexity, can be suitable for multi-scale simulation of the whole world, the region, the city and the like, simultaneously considers calculation efficiency and simulation precision, and simultaneously supports the refined simulation of the near sea area, the far sea area, different sea salt particle sizes and different sea salt chemical components.

The present embodiment mainly relates to 3 calculation schemes, and the features of each calculation scheme will be described below.

1. A first calculation scheme.

The first calculation scheme may refer to a calculation scheme based on geographic attributes.

The computational complexity is: low complexity.

Scheme description: and (3) simplifying an empirical formula to the utmost extent, only distinguishing non-marine attributes from marine attributes, and adopting the same calculation scheme for the particle sizes of different sea salts.

The advantages are that: the method is convenient to realize and high in calculation speed.

The disadvantages are as follows: the simulation precision is low, the near sea area and the far sea area are not distinguished, the generation mechanism of sea salt with different particle sizes is not distinguished, and the simulation of the chemical components of the sea salt particles is not supported.

Applicable scenarios are as follows: global scale long-term simulation; and an application scene with low requirement on sea salt simulation precision.

2. A second calculation scheme.

The second calculation scheme may refer to a calculation scheme based on geographic attributes and particle radii.

The computational complexity is: of medium complexity.

Scheme description: on the basis of a low-complexity scheme, the generation mechanism of sea salt with different particle sizes is considered, and the moisture absorption growth effect of the sea salt particles is also considered.

The advantages are that: the method is convenient to realize and high in calculation speed.

The disadvantages are as follows: the simulation precision is medium, offshore and open sea grids are not distinguished, and the simulation of the chemical components of the sea salt particles is not supported.

Applicable scenarios are as follows: global or regional scale simulation; and (3) an application scene with low requirement on sea salt simulation precision.

3. A third calculation scheme.

The third calculation scheme refers to a calculation scheme based on geographic attributes, particle radii, and sea area attributes.

The computational complexity is: high complexity.

Scheme description: distinguishing a far sea area and a near sea area, and adopting different calculation schemes; adopting different calculation schemes in the open sea area according to different particle sizes; different chemical compositions in sea salt are considered.

The advantages are that: the simulation precision is high, the far-sea area and the near-sea area are distinguished, different particle sizes of sea salt are distinguished, and the simulation of the fine chemical components of the sea salt can be supported.

The disadvantages are as follows: the realization is relatively complex and the calculation speed is slow.

Applicable scenarios are as follows: simulating regional or urban dimensions; and (3) an application scene with high requirements on sea salt simulation precision.

Each calculation scheme may correspond to a different identification, for example, the identification of the first calculation scheme may be "0", the identification of the second calculation scheme may be "1", and the identification of the third calculation scheme may be "2". This embodiment is not limited to this.

The server receives the sea salt discharge flux calculation request, and can analyze the sea salt discharge flux calculation request to acquire information in the sea salt discharge flux calculation request. After the server obtains the configuration parameters of the calculation scheme, the server can obtain the corresponding target calculation scheme according to the identifier of the calculation scheme.

Optionally, the server may further obtain geographic location information and resolution corresponding to the spatial configuration parameter, and a time point or a time period corresponding to the temporal configuration parameter in the sea salt discharge flux calculation request.

And 103, acquiring the required weather information and geographic information by the server according to the target calculation scheme.

The meteorological information may at least include atmospheric information and ocean information, the atmospheric information may at least include a wind speed U component, a wind speed V component, relative humidity, aerosol particle radius, temperature, air pressure, boundary layer height, precipitation, and the like, and the ocean information may at least include seawater salinity, seawater temperature, and the like. The source of the weather information may be global weather reanalysis data, global atmospheric circulation mode simulation and/or mesoscale weather mode simulation, and may also be obtained from other modes, which is not limited in this embodiment.

The geographic information may include at least geographic location information, land use type, and sea ice identification. The geographical location information may correspond to the geographical location information in the sea salt discharge flux calculation request. Land use types may include at least ocean and land. The sea ice flag may be used to indicate whether sea ice is present (e.g., sea ice flag "1" may indicate sea ice present and "0" may indicate no sea ice present).

In one possible implementation, the server may pre-process the acquired weather information to meet the input requirements of the atmospheric chemical transmission mode, and the pre-processed weather information may be stored in a database. After the server determines the target calculation scheme, the data required by the target calculation scheme, that is, the weather information and the geographic information, can be obtained from the database.

Alternatively, the geographic area may be divided into geographic grids, and the calculation may be performed in units of geographic grids. Correspondingly, as shown in the flow chart of the sea salt discharge flux determining method shown in fig. 2, the specific processing of the step 103 may be as follows:

and step 1031, the server determines a plurality of geographic grids of the target area according to the spatial configuration parameters.

In a possible implementation manner, the server may obtain the geographic location information in the spatial configuration parameter, and obtain data of the corresponding target area. Then, the server may obtain a resolution in the spatial configuration parameter, and divide the target area according to the resolution to obtain a plurality of geographic grids.

Optionally, the target region may further include a plurality of sub-regions, each of which may have a corresponding resolution. For example, the first sub-region may be a global region, and the resolution may be 1 × 1 degree; the second subregion can be a Chinese region and a peripheral region, and the resolution can be 0.33 × 0.33 degrees; the first sub-area and the second sub-area are in a nested relation, the second sub-area is contained in the first sub-area, and the resolution of the second sub-area is higher than that of the first sub-area. That is, the areas in china and around can be the areas of the key research, and the resolution with higher precision can be adopted; and the global region may be a region for auxiliary research, and a resolution with lower precision may be adopted. By adopting the mode, the computing resources can be saved, and the processing efficiency is improved.

Step 1032, the server obtains the weather information and the geographic information of each geographic grid according to the target calculation scheme.

In one possible implementation, the server may obtain the data needed to compute each geographic grid on a targeted basis for different computing scenarios.

And step 104, the server determines sea salt discharge flux according to the target calculation scheme, the meteorological information and the geographic information.

In one possible embodiment, the server may use the weather information and the geographic information as input, perform calculation through a target calculation scheme, and determine the sea salt discharge flux in the target area corresponding to the geographic information.

Alternatively, as shown in the flow chart of the sea salt discharge flux determination method shown in fig. 2, the specific processing of step 104 may be as follows:

step 1041, the server determines the underlying surface attribute of each geographic grid according to the geographic information of each geographic grid.

In one possible implementation, the server may identify the underlying surface attribute according to the land use type and the sea ice identification in the geographic information. The underlying surface attribute may be a geographical attribute including a land attribute, a sea ice attribute, and a sea attribute, the sea attribute is subdivided, and the sea attribute may further include at least a sea area attribute divided into a distant sea area and a near sea area.

When the land use type is land, the server may identify an underlying surface attribute of the geographic grid as a land attribute, that is, the geographic attribute of the geographic grid is a land attribute.

When the land use type is sea, the server may further determine a sea ice identification. When the sea ice identification indicates that sea ice is present, the server may identify an underlying surface attribute of the geographic grid as a sea ice attribute, that is, the geographic attribute of the geographic grid is a sea ice attribute. When the sea ice identification indicates that no sea ice exists, the server may identify an underlying surface attribute of the geographic grid as a sea attribute, that is, the geographic attribute of the geographic grid is a sea attribute.

And further dividing the sea attributes, and determining whether a geographical grid with land attributes exists in a plurality of adjacent geographical grids for any geographical grid with the sea attributes. If the geographic grid exists, identifying the underlying surface attribute of the geographic grid as an offshore area, namely, the geographic attribute of the geographic grid is an ocean attribute, and the sea area attribute is the offshore area; and if not, identifying the underlying surface attribute of the geographic grid as the open sea area, namely, the geographic attribute of the geographic grid is the sea attribute, and the sea area attribute is the open sea area.

For example, in practical applications, the identifier of the underlying surface attribute may be set to "distant sea region", "near sea region" or "other". "open sea region" may be used to indicate that the geographic attribute is an ocean attribute, and the sea domain attribute is an open sea region; "offshore region" may be used to indicate that the geographic attribute is an ocean attribute and the sea domain attribute is an offshore region; "other" may be used to indicate that the geographic attribute is a terrestrial attribute or a sea ice attribute. Of course, the underlying surface attribute may also adopt other forms of identifiers to represent the geographic grids of the above various attributes, and the specific identifier of the underlying surface attribute is not limited in this embodiment.

And 1042, calculating by the server through a target calculation scheme based on the meteorological information and the underlying surface attribute of each geographic grid, and determining the sea salt discharge flux of each geographic grid.

The first calculation scheme, i.e. the low complexity calculation scheme, is first described below.

The underlying surface attributes used by the first computing scheme are geographic attributes including land attributes, sea ice attributes, and sea attributes. For convenience of description, herein, in the first computing scheme, a geographic grid whose geographic attribute is a terrestrial attribute or an ice-sea attribute is referred to as a first geographic grid; the geographic grid whose geographic attribute is a marine attribute is referred to as a second geographic grid.

When the target calculation scheme is the first calculation scheme, the process of step 1042 may be as follows: when the geographic attribute of the first geographic grid is a land attribute or an ice attribute, the server sets the sea salt discharge flux of the first geographic grid to a first preset value; and when the geographic attribute of the second geographic grid is the marine attribute, the server determines the sea salt discharge flux of the second geographic grid according to the meteorological information of the second geographic grid.

In one possible implementation, when the server identifies the geographic grid as a first geographic grid in step 1041, the sea salt discharge flux of the geographic grid may be set to a first preset value. The first preset value may be set by a technician, for example may be set to 0, indicating that no sea salt is discharged on land or in sea ice areas. Of course, the first preset value may also be a constant, which is not limited in this embodiment.

When the server identifies the geographic grid as the second geographic grid in step 1041, the server may invoke a calculation function in the first calculation scheme to calculate the weather information of the geographic grid and determine the corresponding sea salt discharge flux.

Specifically, the server may obtain the wind speed at 10 meters of the sea surface of the second geographic grid, and the particle radius of the sea salt.

Determining sea salt discharge flux of the second geographic grid by:

wherein F is sea salt discharge flux, r is the particle radius of sea salt, and U₁₀Wind speed at 10 meters above sea surface, B₁Is an empirical operator related to the radius of the particle.

The calculation process of the first calculation scheme is integrated, so that the first calculation scheme only distinguishes an ocean area from a non-ocean area, and the same calculation scheme is adopted for different sea salt particle sizes.

A second calculation scheme, a medium complexity calculation scheme, will be described below.

The underlying surface attribute used by the second computing scheme may be a geographic attribute including a land attribute, a sea ice attribute, and a sea attribute. For convenience of description, herein, in the second calculation scheme, a geographic grid whose geographic attribute is a terrestrial attribute or a sea ice attribute is referred to as a third geographic grid; the geographic grid whose geographic attribute is a sea attribute is referred to as a fourth geographic grid.

When the target calculation scheme is the second calculation scheme, the process of step 1042 may be as follows:

the server sets the sea salt discharge flux of the third geographic grid to a second preset value when the geographic attribute of the third geographic grid is a land attribute or a sea ice attribute. The processing of the third geographic grid is the same as the first geographic grid, and is not described herein again.

When the geographic attribute of the fourth geographic grid is the marine attribute, the server respectively determines a first sea salt discharge flux F with the particle radius smaller than a first radius threshold value according to the meteorological information of the fourth geographic grid₁And a second sea salt discharge flux F having a particle radius equal to or greater than a first radius threshold₂。

Determining sea salt discharge flux for the fourth geogrid by:

wherein F is the sea salt discharge flux and r is the particle radius of the sea salt.

The first threshold radius may be used to classify fine particle sea salt aerosols and coarse particle sea salt aerosols. For example, the threshold radius may be 2.5 μm, and a particle radius smaller than 2.5 μm is a fine particle sea salt aerosol, and a particle radius greater than or equal to 2.5 μm is a coarse particle sea salt aerosol. The generation mechanism of the fine particle sea salt aerosol is a direct mechanism, and the direct mechanism means that when the wind speed exceeds 10m/s, the top of the spray is directly broken by strong turbulence to generate foam drops. The generation mechanism of the coarse particle sea salt aerosol is an indirect mechanism, and the indirect mechanism means that bubbles generated by wave breaking are broken on the water surface, and countless liquid drops are sprayed after the bubbles are broken.

Specifically, for the fine particle sea salt aerosol, the server may obtain the wind speed at 10 meters of the sea surface of the fourth geographic grid and the corresponding fine particle radius of the sea salt.

Determining a first sea salt discharge flux F of a fine particle sea salt aerosol by the following equation₁：

Wherein, U₁₀Wind speed at 10 meters above sea surface, B₂Which is an empirical operator related to the particle radius, r is the fine particle radius of the corresponding sea salt.

That is, the calculation scheme for fine particle sea salt aerosols may be consistent with a low complexity calculation scheme.

For the coarse particle sea salt aerosol, the server may obtain the wind speed at 10 meters of the sea surface of the fourth geographic grid, the corresponding coarse particle radius of the sea salt.

Determining coarse particle sea salt aerosols by the following formulaSecond sea salt discharge flux F₂：

Wherein, U₁₀The wind speed at 10 meters above the sea surface, r is the coarse particle radius of the corresponding sea salt.

Alternatively, in an actual ocean boundary layer, the relative humidity changes greatly, and the particle radius of the sea salt needs to be corrected as necessary, so as to consider the moisture absorption growth effect of the sea salt, namely, r_d×f_rh。

Determining the hygroscopic increased sea salt emission flux of the fourth geographical grid by:

wherein F is sea salt discharge flux of the fourth geogrid, M is sea salt discharge flux after moisture absorption and growth, and r is_dIs the dry particle radius of sea salt, f_rhIs a moisture absorption growth factor, ρ₁Is the sea salt dry particle density. Exemplary, f_rhThe value may be 1.8, ρ₁The value can be 2200kg/m^-3。

The calculation process of the second calculation scheme is integrated, so that the second calculation scheme distinguishes an ocean area from a non-ocean area, distinguishes fine-particle sea salt aerosol and coarse-particle sea salt aerosol, and different calculation schemes are adopted for different sea salt particle sizes. Also, the second calculation scheme may also take into account the hygroscopic growth effect of sea salt.

A third calculation scheme, a high complexity calculation scheme, will be described below.

The underlying surface attribute used by the third calculation scheme may be a geographical attribute including a land attribute, an ice attribute, and a sea attribute, and the sea attribute may be a sea attribute including a sea area and a sea area. For convenience of description, herein, in the third calculation scheme, a geographic grid whose geographic attribute is a terrestrial attribute or an ice-sea attribute is referred to as a fifth geographic grid; the geographic attribute is an ocean attribute, and the sea domain attribute is a geographic grid of a sea area, which is called a sixth geographic grid; the geographic attribute is the sea attribute, and the sea domain attribute is the geographic grid of the offshore area, which is called as a seventh geographic grid.

When the target computing solution is the third computing solution, the process of step 1042 may be as follows:

the server may set the sea salt discharge flux of the fifth geographic grid to a third preset value when the geographic attribute of the fifth geographic grid is a land attribute or a sea ice attribute. The processing of the fifth geographic grid is the same as the first geographic grid, and is not described herein again.

When the geographic attribute of the sixth geographic grid is the marine attribute and the sea area attribute is the open sea area, the server may respectively determine the third sea salt discharge flux F with the particle radius smaller than the second radius threshold according to the meteorological information of the sixth geographic grid_aAnd a fourth sea salt discharge flux F having a particle radius equal to or greater than a second radius threshold_b(ii) a According to the third sea salt discharge flux F_aAnd fourth sea salt discharge flux F_bDetermining sea salt discharge flux F for a sixth geogrid_{Far away}. In a manner similar to the second calculation scheme described above, the sea salt discharge flux in the open sea area can distinguish between fine particle sea salt aerosols and coarse particle sea salt aerosols, wherein the third sea salt discharge flux F_aNamely the sea salt discharge flux of the fine particle sea salt aerosol and the fourth sea salt discharge flux F_bI.e. the sea salt discharge flux of the coarse particle sea salt aerosol.

When the geographic attribute of the seventh geographic grid is the marine attribute and the sea area attribute is the offshore area, the server may determine the sea salt discharge flux F of the seventh geographic grid according to the meteorological information of the seventh geographic grid_{Near to}。

Specifically, for the open sea area, the server may obtain the wind speed, the relative humidity, and the particle radius of the sea salt at 10 meters of the sea surface of the sixth geographic grid.

Determining the third sea salt discharge flux F of the fine particle sea salt aerosol by the following equation_a：

A₁＝-5.001E³+0.808E⁶r_RHC₈₀-1.98E⁷(r_RHC₈₀)²+2.188E⁸(r_RHC₈₀)³ -1.144E⁹(r_RHC₈₀)⁴+2.29E⁹(r_RHC₈₀)⁵

A₂＝3.854E³+1.168E⁴r_RHC₈₀-6.572E⁴(r_RHC₈₀)²+1.003E⁵(r_RHC₈₀)³ -6.407E⁴(r_RHC₈₀)⁴+1.493E⁴(r_RHC₈₀)⁵

A₃＝4.498E²+0.839E³r_RHC₈₀-5.394E²(r_RHC₈₀)²+1.218E²(r_RHC₈₀)³ -1.213E¹(r_RHC₈₀)⁴+4.514E^-1(r_RHC₈₀)⁵

ρ₂＝1000×(3.8033-16.248RH+46.085RH²-68.317RH³+50.932RH⁴-15.261RH⁵)

x＝3.1657-19.079RH+55.72RH²-83.998RH³+63.436RH⁴-19.248RH⁵

Wherein A is_iThe value range of i is {1, 2, 3} which is an empirical coefficient; r is_iIs the particle radius (in μm) of sea salt, r₀Has a value range of (0,0.1), r₁Has a value range of [0.1,1), r₂Has a value range of [1,2.5), r₃The value range of (1) is [2.5,5 ]; RH is relative humidity, r_RHIs the aerosol particle radius at RH humidity; u shape₁₀The wind speed at 10 meters above sea level; c₈₀Is a correction factor for salinity; AB₁The conversion coefficient of the aerosol concentration and the mass flux is obtained; rho₂Is sea salt particle density; x is the mass fraction of solute; eⁿIs 10ⁿ。

Determining the fourth sea salt discharge flux F of the coarse particle sea salt aerosol by the following formula_b：

ρ₂＝1000×(3.8033-16.248RH+46.085RH²-68.317RH³+50.932RH⁴-15.261RH⁵)

x＝3.1657-19.079RH+55.72RH²-83.998RH³+63.436RH⁴-19.248RH⁵

Wherein A is₄Is an empirical coefficient; r is₃、r₄Is the particle radius (in μm) of sea salt, r₃Has a value range of [2.5,5), r₄Has a value range of [5,10 ]](ii) a RH is relative humidity, r_RHIs the particle radius of the sea salt at RH humidity; u shape₁₀The wind speed at 10 meters above sea level; c₈₀Is a correction factor for salinity; AB₁The conversion coefficient of the aerosol concentration and the mass flux is obtained; rho₂Is sea salt particle density; x is the mass fraction of solute; b is₃Is an empirical operator related to the radius of the particle.

The server may then discharge a third sea salt flux F_aAnd fourth sea salt discharge flux F_bAdding the sea salt discharge flux of the fine particle sea salt aerosol and the coarse particle sea salt aerosol to obtain the sea salt discharge flux of the open sea area, namely F_{Far away}＝F_a+F_b。

For the offshore region, the server may obtain wind speed, relative humidity, particle radius of sea salt, and sea water salinity at 10 meters of the sea surface of the seventh geographic grid.

Determining sea salt discharge flux F for the seventh geogrid by the following formula_{Near to}：

AB₂＝2×10^-15×4.188×r_RH ³×ρ₂x

ρ₂＝1000×(3.8033-16.248RH+46.085RH²-68.317RH³+50.932RH⁴-15.261RH⁵)

x＝3.1657-19.079RH+55.72RH²-83.998RH³+63.436RH⁴-19.248RH⁵

Wherein, F_cFlux F for sea salt discharge_{Near to}RH is the relative humidity, r_RHIs the particle radius, U, of sea salt at RH humidity₁₀The wind speed at 10m above sea surface, S is the salinity of seawater, C₀As a correction factor for salinity, AB₂Is the conversion coefficient of aerosol concentration to mass flux, rho₂Is the sea salt particle density and x is the solute mass fraction.

Optionally, in order to further refine different chemical components in the simulated sea salt, according to comprehensive research of documents and observation results in the near sea area of china, the component ratio of the sea salt particles is defined as follows:

TABLE 1 ratio of different chemical components in sea salt particles

For the sixth geographic grid and the seventh geographic grid, the processing of the server may further include: acquiring preset proportion of multiple chemical components; and determining the corresponding discharge flux of each chemical component in the sea salt according to the ratio of each chemical component and the sea salt discharge flux of the geographical grid.

In a possible embodiment, after the server determines the sea salt discharge flux of the geographical grids of the open sea region and the offshore region by the above method, the sea salt discharge flux may be multiplied by the ratio of each chemical component to obtain the discharge flux corresponding to each chemical component in the sea salt.

Optionally, in order to further simulate the physicochemical process of calculating the sea salt aerosol and obtain the concentration distribution of the sea salt aerosol in the atmosphere, referring to the flow charts of the sea salt emission flux determination methods shown in fig. 2 and fig. 3, after performing step 104, the processing of the server may further be as follows:

and 105, simulating the physical and chemical process of the sea salt aerosol by the server according to the sea salt discharge flux.

In a possible implementation manner, the server can calculate the physical and chemical processes of the sea salt aerosol in the atmosphere in a target area range by grid integration on the basis of obtaining the sea salt discharge flux with different particle sizes in the above steps, wherein the physical and chemical processes mainly comprise advection, diffusion, convection, dry sedimentation, wet sedimentation, gravity sedimentation and the like, and the concentrations of the sea salt aerosol and chemical components thereof with different particle sizes at different time and different spatial positions can be obtained.

In the above steps, according to the time configuration parameters in the sea salt discharge flux calculation request, for different geographic grids and/or time steps, the server may cyclically execute steps 1032, 1041, 1042, and 105, dynamically simulate the physicochemical process of the sea salt aerosol, that is, the processes of discharge, advection, diffusion, convection, dry-wet precipitation, and the like of the sea salt aerosol in the atmosphere, and finally obtain the three-dimensional space-time distribution of the sea salt aerosol concentration in different time and space ranges. In the atmospheric chemical model simulation process, the simulation can be realized by parallel calculation on the geographic grids.

And 106, outputting the concentration distribution of the sea salt aerosol by the server according to the calculation result of the physical and chemical process of the sea salt aerosol and the output frequency setting.

In one possible embodiment, the server may output three-dimensional distributions of concentrations of sea salt aerosols with different particle sizes and chemical components thereof according to a preset simulation result output frequency (for example, per hour) in the simulation period. The output results can be used to analyze the temporal-spatial distribution characteristics of the sea salt aerosol and its impact on the total aerosol concentration and composition in the atmosphere.

Of course, the server may output the concentration distribution of the sea salt aerosol in the form of a graph, or may output the concentration distribution of the sea salt aerosol in the form of a two-dimensional map, which is not limited in this embodiment.

The invention has good performance in the experimental stage. By adopting a high-complexity sea salt discharge flux calculation scheme, the simulation time period covers the whole year of 2010, and on a high-performance computing cluster, 128-core parallel simulation is adopted, so that 1-year simulation is completed in about 8 days, and the total calculation speed is moderate. On a global scale, the high-value sea salt concentration area is mainly located in a medium-high latitude sea area in north-south latitude and is close to the wind zone distribution of the south-north hemisphere. The annual average concentration of sea salt can reach 20 mu g/m on open sea³The above. In coastal region of China, the annual average concentration of sea salt is between 3 and 7 mu g/m³Generally, the sea salt concentration gradually increases from coastal areas to open sea areas. The simulated spatial distribution characteristics of sea salt concentration may reasonably reflect the global and regional distribution of sea salt as a whole.

In the embodiment of the invention, the sea salt discharge flux calculation request can carry the calculation scheme configuration parameters, and the server can select the target calculation scheme through the calculation scheme configuration parameters so as to determine the sea salt discharge flux. The method realizes integration and improvement of a plurality of sea salt calculation schemes with different complexities under the same mode frame, can be suitable for multi-scale simulation of the whole world, the region, the city and the like, and simultaneously gives consideration to the calculation efficiency and the simulation precision.

The embodiment of the invention provides a device for determining sea salt discharge flux, which is used for realizing the method for determining the sea salt discharge flux. A schematic block diagram of a sea salt discharge flux determining apparatus as shown in fig. 4, the apparatus comprising:

a receiving module 701, configured to receive a sea salt discharge flux calculation request, where the sea salt discharge flux calculation request at least includes a calculation scheme configuration parameter;

an obtaining module 702, configured to determine a target calculation scheme according to the calculation scheme configuration parameters; acquiring the required meteorological information and geographic information according to the target calculation scheme;

a determining module 703, configured to determine the sea salt discharge flux according to the target calculation scheme, the meteorological information, and the geographic information.

Optionally, the sea salt discharge flux calculation request further includes a space configuration parameter;

the obtaining module 702 is configured to: determining a plurality of geographical grids of a target area according to the spatial configuration parameters; acquiring meteorological information and geographic information of each geographic grid according to the target calculation scheme;

the determining module 703 is configured to: determining the underlying surface attribute of each geographic grid according to the geographic information of each geographic grid; and calculating through the target calculation scheme based on the meteorological information and the underlying surface attribute of each geographic grid, and determining the sea salt discharge flux of each geographic grid.

Optionally, the target calculation scheme is a first calculation scheme, and the first calculation scheme is a calculation scheme based on geographic attributes; the underlying surface attribute used by the first computing scheme is a geographic attribute comprising a land attribute, a sea ice attribute, and a sea attribute;

the determining module 703 is configured to:

setting a sea salt discharge flux of a first geographic grid to a first preset value when a geographic attribute of the first geographic grid is a land attribute or a sea ice attribute;

and when the geographic attribute of the second geographic grid is the marine attribute, determining the sea salt discharge flux of the second geographic grid according to the meteorological information of the second geographic grid.

Optionally, the determining module 703 is configured to:

acquiring the wind speed and the particle radius of sea salt at the sea surface of the second geographic grid at 10 meters;

determining sea salt discharge flux of the second geographic grid by:

wherein F is sea salt discharge flux, r is the particle radius of sea salt, and U₁₀Wind speed at 10 meters above sea surface, B₁Is an empirical operator related to the radius of the particle.

Optionally, the target calculation scheme is a second calculation scheme, where the second calculation scheme is a calculation scheme based on geographic attributes and particle radii; the underlying surface attribute used by the second calculation scheme is a geographic attribute comprising a land attribute, a sea ice attribute, and a sea attribute;

the determining module 703 is configured to:

when the geographic attribute of the third geographic grid is a land attribute or a sea ice attribute, setting the sea salt discharge flux of the third geographic grid to be a second preset value;

when the geographic attribute of the fourth geographic grid is the marine attribute, respectively determining a first sea salt discharge flux F with the particle radius smaller than a first radius threshold value according to the meteorological information of the fourth geographic grid₁And a second sea salt discharge flux F having a particle radius equal to or greater than a first radius threshold₂；

Determining sea salt discharge flux for the fourth geographic grid by:

wherein F is the sea salt discharge flux and r is the particle radius of the sea salt.

Optionally, the determining module 703 is configured to:

acquiring the wind speed and the particle radius of sea salt at the sea surface of the fourth geographical grid of 10 meters;

determining said first sea salt discharge flux F by the following formula₁：

Wherein, U₁₀Wind speed at 10 meters above sea surface, B₂Is an empirical operator related to the radius of the particle.

Optionally, the determining module 703 is configured to:

acquiring the wind speed and the particle radius of sea salt at the sea surface of the fourth geographical grid of 10 meters;

determining the second sea salt discharge flux F by the following formula₂：

Wherein, U₁₀The wind speed at 10 meters above sea surface.

Optionally, the determining module 703 is further configured to:

determining the hygroscopic increased sea salt emission flux of the fourth geographical grid by:

wherein M is the sea salt discharge flux after moisture absorption and growth, r_dIs the dry particle radius of sea salt, f_rhIs a moisture absorption growth factor, ρ₁Is the sea salt dry particle density.

Optionally, the target calculation scheme is a third calculation scheme, where the third calculation scheme is a calculation scheme based on geographic attributes, particle radii, and sea area attributes; the underlay surface attribute used by the third calculation scheme is a geographical attribute comprising a land attribute, an ice attribute and a sea attribute, wherein the sea attribute is a sea area attribute comprising a sea area and a sea area;

the determining module 703 is configured to:

when the geographic attribute of a fifth geographic grid is a land attribute or a sea ice attribute, setting the sea salt discharge flux of the fifth geographic grid to a third preset value;

when the geographic attribute of the sixth geographic grid is the marine attribute and the sea area attribute is the open sea area, respectively determining a third sea salt discharge flux F with the particle radius smaller than a second radius threshold value according to the meteorological information of the sixth geographic grid_aAnd a fourth sea salt discharge flux F having a particle radius equal to or greater than a second radius threshold_b(ii) a According to the third sea salt discharge flux F_aAnd fourth sea salt discharge flux F_bDetermining sea salt discharge flux F of said sixth geogrid_{Far away}；

When the geographic attribute of the seventh geographic grid is the marine attribute and the sea area attribute is the offshore area, determining the sea salt discharge flux F of the seventh geographic grid according to the meteorological information of the seventh geographic grid_{Near to}。

Optionally, the determining module 703 is configured to:

acquiring the wind speed, the relative humidity and the particle radius of sea salt at the sea surface of the sixth geographical grid of 10 meters;

determining said third sea salt discharge flux F by the following formula_a：

A₁＝-5.001E³+0.808E⁶r_RHC₈₀-1.98E⁷(r_RHC₈₀)²+2.188E⁸(r_RHC₈₀)³ -1.144E⁹(r_RHC₈₀)⁴+2.29E⁹(r_RHC₈₀)⁵

A₂＝3.854E³+1.168E⁴r_RHC₈₀-6.572E⁴(r_RHC₈₀)²+1.003E⁵(r_RHC₈₀)³ -6.407E⁴(r_RHC₈₀)⁴+1.493E⁴(r_RHC₈₀)⁵

A₃＝4.498E²+0.839E³r_RHC₈₀-5.394E²(r_RHC₈₀)²+1.218E²(r_RHC₈₀)³ -1.213E¹(r_RHC₈₀)⁴+4.514E^-1(r_RHC₈₀)⁵

ρ₂＝1000×(3.8033-16.248RH+46.085RH²-68.317RH³+50.932RH⁴-15.261RH⁵)

x＝3.1657-19.079RH+55.72RH²-83.998RH³+63.436RH⁴-19.248RH⁵

Wherein A is_iThe value range of i is {1, 2, 3} which is an empirical coefficient; r is_iIs the particle radius of sea salt, r₀Has a value range of (0,0.1), r₁Has a value range of [0.1,1), r₂Has a value range of [1,2.5), r₃The value range of (1) is [2.5,5 ]; RH is relative humidity, r_RHIs the aerosol particle radius at RH humidity; u shape₁₀The wind speed at 10 meters above sea level; c₈₀Is a correction factor for salinity; AB₁The conversion coefficient of the aerosol concentration and the mass flux is obtained; rho₂Is sea salt particle density; x is the mass fraction of solute; eⁿIs 10ⁿ。

Optionally, the determining module 703 is configured to:

acquiring the wind speed, the relative humidity and the particle radius of sea salt at the sea surface of the sixth geographical grid of 10 meters;

determining the fourth sea salt discharge flux F by the following formula_b：

ρ₂＝1000×(3.8033-16.248RH+46.085RH²-68.317RH³+50.932RH⁴-15.261RH⁵)

x＝3.1657-19.079RH+55.72RH²-83.998RH³+63.436RH⁴-19.248RH⁵

Wherein A is₄Is an empirical coefficient; r is₃、r₄Is the particle radius of sea salt, r₃Has a value range of [2.5,5), r₄Has a value range of [5,10 ]](ii) a RH is relative humidity, r_RHIs the particle radius of the sea salt at RH humidity; u shape₁₀The wind speed at 10 meters above sea level; c₈₀Is a correction factor for salinity; AB₁As aerosol concentrationConversion coefficient to mass flux; rho₂Is sea salt particle density; x is the mass fraction of solute; b is₃Is an empirical operator related to the radius of the particle.

Optionally, the determining module 703 is configured to:

acquiring the wind speed, the relative humidity, the particle radius of sea salt and the salinity of sea water at the sea surface of the seventh geographical grid of 10 meters; determining sea salt discharge flux F of the seventh geogrid by_{Near to}：

AB₂＝2×10^-15×4.188×r_RH ³×ρ₂x

ρ₂＝1000×(3.8033-16.248RH+46.085RH²-68.317RH³+50.932RH⁴-15.261RH⁵)

x＝3.1657-19.079RH+55.72RH²-83.998RH³+63.436RH⁴-19.248RH⁵

Wherein, F_cFlux F for sea salt discharge_{Near to}RH is the relative humidity, r_RHIs the particle radius, U, of sea salt at RH humidity₁₀The wind speed at 10m above sea surface, S is the salinity of seawater, C₀As a correction factor for salinity, AB₂Is the conversion coefficient of aerosol concentration to mass flux, rho₂Is the sea salt particle density and x is the solute mass fraction.

Optionally, for the sixth geographic grid and the seventh geographic grid, the determining module 703 is further configured to:

acquiring preset proportion of multiple chemical components;

and determining the sea salt discharge flux corresponding to each chemical component according to the ratio of each chemical component and the sea salt discharge flux of the geographical grid.

Optionally, the determining module 703 is further configured to:

simulating a physical and chemical process of sea salt aerosol according to the sea salt discharge flux;

and outputting the concentration distribution of the sea salt aerosol according to the calculation result of the physical and chemical process of the sea salt aerosol and the output frequency setting.

In the embodiment of the invention, the sea salt discharge flux calculation request can carry the calculation scheme configuration parameters, and the server can select the target calculation scheme through the calculation scheme configuration parameters so as to determine the sea salt discharge flux. The method realizes integration and improvement of a plurality of sea salt calculation schemes with different complexities under the same mode frame, can be suitable for multi-scale simulation of the whole world, the region, the city and the like, and simultaneously gives consideration to the calculation efficiency and the simulation precision.

An exemplary embodiment of the present invention also provides an electronic device including: at least one processor; and a memory communicatively coupled to the at least one processor. The memory stores a computer program executable by the at least one processor, the computer program, when executed by the at least one processor, is for causing the electronic device to perform a method according to an embodiment of the invention.

Exemplary embodiments of the present invention also provide a non-transitory computer-readable storage medium storing a computer program, wherein the computer program, when executed by a processor of a computer, is operable to cause the computer to perform a method according to an embodiment of the present invention.

Exemplary embodiments of the present invention also provide a computer program product comprising a computer program, wherein the computer program is operative, when executed by a processor of a computer, to cause the computer to perform a method according to an embodiment of the present invention.

Referring to fig. 5, a block diagram of a structure of an electronic device 800, which may be a server or a client of the present invention, which is an example of a hardware device that may be applied to aspects of the present invention, will now be described. Electronic device is intended to represent various forms of digital electronic computer devices, such as laptops, desktops, workstations, personal digital assistants, servers, blade servers, mainframes, and other suitable computers. The electronic device may also represent various forms of mobile devices, such as personal digital processing, cellular phones, smart phones, wearable devices, and other similar computing devices. The components shown herein, their connections and relationships, and their functions, are meant to be exemplary only, and are not meant to limit implementations of the inventions described and/or claimed herein.

As shown in fig. 5, the electronic device 800 includes a computing unit 801 that can perform various appropriate actions and processes according to a computer program stored in a Read Only Memory (ROM) 802 or a computer program loaded from a storage unit 808 into a Random Access Memory (RAM) 803. In the RAM 803, various programs and data required for the operation of the device 800 can also be stored. The calculation unit 801, the ROM 802, and the RAM 803 are connected to each other by a bus 804. An input/output (I/O) interface 805 is also connected to bus 804.

A number of components in the electronic device 800 are connected to the I/O interface 805, including: an input unit 806, an output unit 807, a storage unit 808, and a communication unit 809. The input unit 806 may be any type of device capable of inputting information to the electronic device 800, and the input unit 806 may receive input numeric or character information and generate key signal inputs related to user settings and/or function controls of the electronic device. Output unit 807 can be any type of device capable of presenting information and can include, but is not limited to, a display, speakers, a video/audio output terminal, a vibrator, and/or a printer. The storage unit 808 may include, but is not limited to, a magnetic disk, an optical disk. The communication unit 809 allows the electronic device 800 to exchange information/data with other devices via a computer network, such as the internet, and/or various telecommunications networks, and may include, but is not limited to, modems, network cards, infrared communication devices, wireless communication transceivers and/or chipsets, such as bluetooth devices, WiFi devices, WiMax devices, cellular communication devices, and/or the like.

Computing unit 801 may be a variety of general and/or special purpose processing components with processing and computing capabilities. Some examples of the computing unit 801 include, but are not limited to, a Central Processing Unit (CPU), a Graphics Processing Unit (GPU), various dedicated Artificial Intelligence (AI) computing chips, various computing units running machine learning model algorithms, a Digital Signal Processor (DSP), and any suitable processor, controller, microcontroller, and the like. The calculation unit 801 executes the respective methods and processes described above. For example, in some embodiments, the method of determining sea salt discharge flux may be implemented as a computer software program tangibly embodied in a machine-readable medium, such as storage unit 808. In some embodiments, part or all of the computer program can be loaded and/or installed onto the electronic device 800 via the ROM 802 and/or the communication unit 809. In some embodiments, the computing unit 801 may be configured to perform the determination method of sea salt discharge flux by any other suitable means (e.g., by means of firmware).

Program code for implementing the methods of the present invention may be written in any combination of one or more programming languages. These program codes may be provided to a processor or controller of a general purpose computer, special purpose computer, or other programmable data processing apparatus, such that the program codes, when executed by the processor or controller, cause the functions/operations specified in the flowchart and/or block diagram to be performed. The program code may execute entirely on the machine, partly on the machine, as a stand-alone software package partly on the machine and partly on a remote machine or entirely on the remote machine or server.

In the context of the present invention, a machine-readable medium may be a tangible medium that can contain, or store a program for use by or in connection with an instruction execution system, apparatus, or device. The machine-readable medium may be a machine-readable signal medium or a machine-readable storage medium. A machine-readable medium may include, but is not limited to, an electronic, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus, or device, or any suitable combination of the foregoing. More specific examples of a machine-readable storage medium would include an electrical connection based on one or more wires, a portable computer diskette, a hard disk, a Random Access Memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disc read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the foregoing.

As used herein, the terms "machine-readable medium" and "computer-readable medium" refer to any computer program product, apparatus, and/or device (e.g., magnetic discs, optical disks, memory, Programmable Logic Devices (PLDs)) used to provide machine instructions and/or data to a programmable processor, including a machine-readable medium that receives machine instructions as a machine-readable signal. The term "machine-readable signal" refers to any signal used to provide machine instructions and/or data to a programmable processor.

To provide for interaction with a user, the systems and techniques described here can be implemented on a computer having: a display device (e.g., a CRT (cathode ray tube) or LCD (liquid crystal display) monitor) for displaying information to a user; and a keyboard and a pointing device (e.g., a mouse or a trackball) by which a user can provide input to the computer. Other kinds of devices may also be used to provide for interaction with a user; for example, feedback provided to the user can be any form of sensory feedback (e.g., visual feedback, auditory feedback, or tactile feedback); and input from the user may be received in any form, including acoustic, speech, or tactile input.

The systems and techniques described here can be implemented in a computing system that includes a back-end component (e.g., as a data server), or that includes a middleware component (e.g., an application server), or that includes a front-end component (e.g., a user computer having a graphical user interface or a web browser through which a user can interact with an implementation of the systems and techniques described here), or any combination of such back-end, middleware, or front-end components. The components of the system can be interconnected by any form or medium of digital data communication (e.g., a communication network). Examples of communication networks include: local Area Networks (LANs), Wide Area Networks (WANs), and the Internet.

The computer system may include clients and servers. A client and server are generally remote from each other and typically interact through a communication network. The relationship of client and server arises by virtue of computer programs running on the respective computers and having a client-server relationship to each other.

Claims (17)

1. A method for determining sea salt discharge flux, the method comprising:

receiving a sea salt discharge flux calculation request, wherein the sea salt discharge flux calculation request at least comprises calculation scheme configuration parameters;

determining a target calculation scheme according to the calculation scheme configuration parameters;

acquiring the required meteorological information and geographic information according to the target calculation scheme;

and determining sea salt discharge flux according to the target calculation scheme, the meteorological information and the geographic information.

2. The method of determining sea salt discharge flux according to claim 1, wherein said sea salt discharge flux calculation request further comprises spatial configuration parameters;

the acquiring of the required weather information and the geographic information according to the target calculation scheme comprises:

determining a plurality of geographical grids of a target area according to the spatial configuration parameters;

acquiring meteorological information and geographic information of each geographic grid according to the target calculation scheme;

said determining sea salt discharge flux from said target calculation scenario, said meteorological information, said geographic information, comprising:

determining the underlying surface attribute of each geographic grid according to the geographic information of each geographic grid;

and calculating through the target calculation scheme based on the meteorological information and the underlying surface attribute of each geographic grid, and determining the sea salt discharge flux of each geographic grid.

3. The method for determining sea salt discharge flux according to claim 2, wherein the target calculation scheme is a first calculation scheme, and the first calculation scheme is a calculation scheme based on geographic attributes; the underlying surface attribute used by the first computing scheme is a geographic attribute comprising a land attribute, a sea ice attribute, and a sea attribute;

said determining sea salt discharge flux for each said geographic grid calculated by said target calculation scheme based on said meteorological information and underlying surface properties for said each geographic grid comprises:

setting a sea salt discharge flux of a first geographic grid to a first preset value when a geographic attribute of the first geographic grid is a land attribute or a sea ice attribute;

and when the geographic attribute of the second geographic grid is the marine attribute, determining the sea salt discharge flux of the second geographic grid according to the meteorological information of the second geographic grid.

4. The method for determining sea salt discharge flux according to claim 3, wherein said determining sea salt discharge flux of the second geographic grid according to the weather information of the second geographic grid comprises:

acquiring the wind speed and the particle radius of sea salt at the sea surface of the second geographic grid at 10 meters;

determining sea salt discharge flux of the second geographic grid by:

wherein F is sea salt discharge flux, r is the particle radius of sea salt, and U₁₀Wind speed at 10 meters above sea surface, B₁Is an empirical operator related to the radius of the particle.

5. The method for determining sea salt discharge flux according to claim 2, wherein the target calculation scheme is a second calculation scheme, and the second calculation scheme is a calculation scheme based on geographic attributes and particle radii; the underlying surface attribute used by the second calculation scheme is a geographic attribute comprising a land attribute, a sea ice attribute, and a sea attribute;

said determining sea salt discharge flux for each said geographic grid calculated by said target calculation scheme based on said meteorological information and underlying surface properties for said each geographic grid comprises:

when the geographic attribute of the third geographic grid is a land attribute or a sea ice attribute, setting the sea salt discharge flux of the third geographic grid to be a second preset value;

when the geographic attribute of the fourth geographic grid is the marine attribute, respectively determining a first sea salt discharge flux F with the particle radius smaller than a first radius threshold value according to the meteorological information of the fourth geographic grid₁And a second sea salt discharge flux F having a particle radius equal to or greater than a first radius threshold₂；

Determining sea salt discharge flux for the fourth geographic grid by:

wherein F is the sea salt discharge flux and r is the particle radius of the sea salt.

6. The method of claim 5, wherein determining the first sea salt discharge flux F with a particle radius less than a first radius threshold is based on meteorological information of the fourth geographic grid₁The method comprises the following steps:

acquiring the wind speed and the particle radius of sea salt at the sea surface of the fourth geographical grid of 10 meters;

determining said first sea salt discharge flux F by the following formula₁：

Wherein, U₁₀Wind speed at 10 meters above sea surface, B₂Is an empirical operator related to the radius of the particle.

7. Method for determining sea salt discharge flux according to claim 5, characterized in that a second sea salt discharge flux F is determined with a particle radius greater than or equal to a first radius threshold value, based on meteorological information of said fourth geographical grid₂The method comprises the following steps:

acquiring the wind speed and the particle radius of sea salt at the sea surface of the fourth geographical grid of 10 meters;

determining the second sea salt discharge flux F by the following formula₂：

Wherein, U₁₀The wind speed at 10 meters above sea surface.

8. The method of determining sea salt discharge flux of claim 5, further comprising:

determining the hygroscopic increased sea salt emission flux of the fourth geographical grid by:

wherein M is the sea salt discharge flux after moisture absorption and growth, r_dIs the dry particle radius of sea salt, f_rhIs a moisture absorption growth factor, ρ₁Is the sea salt dry particle density.

9. The method of determining sea salt discharge flux according to claim 2, wherein the target calculation scheme is a third calculation scheme, the third calculation scheme being a calculation scheme based on geographical attributes, particle radius and sea area attributes; the underlay surface attribute used by the third calculation scheme is a geographical attribute comprising a land attribute, an ice attribute and a sea attribute, wherein the sea attribute is a sea area attribute comprising a sea area and a sea area;

said determining sea salt discharge flux for each said geographic grid calculated by said target calculation scheme based on said meteorological information and underlying surface properties for said each geographic grid comprises:

when the geographic attribute of a fifth geographic grid is a land attribute or a sea ice attribute, setting the sea salt discharge flux of the fifth geographic grid to a third preset value;

when the geographic attribute of the sixth geographic grid is the marine attribute and the sea area attribute is the open sea area, respectively determining a third sea salt discharge flux F with the particle radius smaller than a second radius threshold value according to the meteorological information of the sixth geographic grid_aAnd a fourth sea salt discharge flux F having a particle radius equal to or greater than a second radius threshold_b(ii) a According to the third sea salt discharge flux F_aAnd the fourth sea salt discharge channelQuantity F_bDetermining sea salt discharge flux F of said sixth geogrid_{Far away}；

When the geographic attribute of the seventh geographic grid is the marine attribute and the sea area attribute is the offshore area, determining the sea salt discharge flux F of the seventh geographic grid according to the meteorological information of the seventh geographic grid_{Near to}。

10. The method of claim 9, wherein a third sea salt discharge flux F is determined having a particle radius less than a second radius threshold based on meteorological information of the sixth geographic grid_aThe method comprises the following steps:

acquiring the wind speed, the relative humidity and the particle radius of sea salt at the sea surface of the sixth geographical grid of 10 meters;

determining said third sea salt discharge flux F by the following formula_a：

A₁＝-5.001E³+0.808E⁶r_RHC₈₀-1.98E⁷(r_RHC₈₀)²+2.188E⁸(r_RHC₈₀)³-1.144E⁹(r_RHC₈₀)⁴+2.29E⁹(r_RHC₈₀)⁵

A₂＝3.854E³+1.168E⁴r_RHC₈₀-6.572E⁴(r_RHC₈₀)²+1.003E⁵(r_RHC₈₀)³-6.407E⁴(r_RHC₈₀)⁴+1.493E⁴(r_RHC₈₀)⁵

A₃＝4.498E²+0.839E³r_RHC₈₀-5.394E²(r_RHC₈₀)²+1.218E²(r_RHC₈₀)³-1.213E¹(r_RHC₈₀)⁴+4.514E^-1(r_RHC₈₀)⁵

ρ₂＝1000×(3.8033-16.248RH+46.085RH²-68.317RH³+50.932RH⁴-15.261RH⁵)

x＝3.1657-19.079RH+55.72RH²-83.998RH³+63.436RH⁴-19.248RH⁵

Wherein A is_iThe value range of i is {1, 2, 3} which is an empirical coefficient; r is_iIs the particle radius of sea salt, r₀Has a value range of (0,0.1), r₁Has a value range of [0.1,1), r₂Has a value range of [1,2.5), r₃The value range of (1) is [2.5,5 ]; RH is relative humidity, r_RHIs the aerosol particle radius at RH humidity; u shape₁₀The wind speed at 10 meters above sea level; c₈₀Is a correction factor for salinity; AB₁The conversion coefficient of the aerosol concentration and the mass flux is obtained; rho₂Is sea salt particle density; x is the mass fraction of solute; eⁿIs 10ⁿ。

11. The method of claim 9, wherein a fourth sea salt discharge flux F having a particle radius equal to or greater than a radius threshold is determined based on the weather information of the sixth geographic grid_bThe method comprises the following steps:

acquiring the wind speed, the relative humidity and the particle radius of sea salt at the sea surface of the sixth geographical grid of 10 meters;

determining the fourth sea salt discharge flux F by the following formula_b：

ρ₂＝1000×(3.8033-16.248RH+46.085RH²-68.317RH³+50.932RH⁴-15.261RH⁵)

x＝3.1657-19.079RH+55.72RH²-83.998RH³+63.436RH⁴-19.248RH⁵

Wherein A is₄Is an empirical coefficient; r is₃、r₄Is the particle radius of sea salt, r₃Has a value range of [2.5,5), r₄Has a value range of [5,10 ]](ii) a RH is relative humidity, r_RHIs the particle radius of the sea salt at RH humidity; u shape₁₀The wind speed at 10 meters above sea level; c₈₀For correcting salinityA positive factor; AB₁The conversion coefficient of the aerosol concentration and the mass flux is obtained; rho₂Is sea salt particle density; x is the mass fraction of solute; b is₃Is an empirical operator related to the radius of the particle.

12. The method of claim 9, wherein said determining sea salt discharge flux F of said seventh geographic grid is based on meteorological information of said seventh geographic grid_{Near to}The method comprises the following steps:

acquiring the wind speed, the relative humidity, the particle radius of sea salt and the salinity of sea water at the sea surface of the seventh geographical grid of 10 meters;

determining sea salt discharge flux F of the seventh geogrid by_{Near to}：

AB₂＝2×10^-15×4.188×r_RH ³×ρ₂x

ρ₂＝1000×(3.8033-16.248RH+46.085RH²-68.317RH³+50.932RH⁴-15.261RH⁵)

x＝3.1657-19.079RH+55.72RH²-83.998RH³+63.436RH⁴-19.248RH⁵

Wherein, F_cFlux F for sea salt discharge_{Near to}RH is the relative humidity, r_RHIs the particle radius, U, of sea salt at RH humidity₁₀The wind speed at 10m above sea surface, S is the salinity of seawater, C₀As a correction factor for salinity, AB₂Is the conversion coefficient of aerosol concentration to mass flux, rho₂Is the sea salt particle density and x is the solute mass fraction.

13. The method of determining sea salt discharge flux of claim 9, wherein for the sixth geographical grid and the seventh geographical grid, the method further comprises:

acquiring preset proportion of multiple chemical components;

and determining the corresponding discharge flux of each chemical component in the sea salt according to the ratio of each chemical component and the sea salt discharge flux of the geographical grid.

14. The method of determining sea salt discharge flux of claim 1, further comprising:

simulating a physical and chemical process of sea salt aerosol according to the sea salt discharge flux;

and outputting the concentration distribution of the sea salt aerosol according to the calculation result of the physical and chemical process of the sea salt aerosol and the output frequency setting.

15. An apparatus for determining sea salt discharge flux, the apparatus comprising:

the sea salt discharge flux calculation module is used for receiving a sea salt discharge flux calculation request, and the sea salt discharge flux calculation request at least comprises calculation scheme configuration parameters;

the acquisition module is used for determining a target calculation scheme according to the calculation scheme configuration parameters; acquiring the required meteorological information and geographic information according to the target calculation scheme;

and the determining module is used for determining sea salt discharge flux according to the target calculation scheme, the meteorological information and the geographic information.

16. An electronic device, comprising:

a processor; and

a memory for storing a program, wherein the program is stored in the memory,

wherein the program comprises instructions which, when executed by the processor, cause the processor to carry out the method according to any one of claims 1-14.

17. A non-transitory computer readable storage medium having stored thereon computer instructions for causing a computer to perform the method of any one of claims 1-14.

Images