CN114386860A - Method for determining regional atmospheric precipitation pattern index - Google Patents

Method for determining regional atmospheric precipitation pattern index Download PDF

Info

Publication number
CN114386860A
CN114386860A CN202210044462.8A CN202210044462A CN114386860A CN 114386860 A CN114386860 A CN 114386860A CN 202210044462 A CN202210044462 A CN 202210044462A CN 114386860 A CN114386860 A CN 114386860A
Authority
CN
China
Prior art keywords
precipitation
water vapor
regional
index
rainfall
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202210044462.8A
Other languages
Chinese (zh)
Other versions
CN114386860B (en
Inventor
权全
邓嘉祥
董宇翔
李平治
杨思敏
许美娇
秦毅
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xian University of Technology
Original Assignee
Xian University of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Xian University of Technology filed Critical Xian University of Technology
Priority to CN202210044462.8A priority Critical patent/CN114386860B/en
Publication of CN114386860A publication Critical patent/CN114386860A/en
Application granted granted Critical
Publication of CN114386860B publication Critical patent/CN114386860B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063—Operations research, analysis or management
    • G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
    • G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/10—Services
    • G06Q50/26—Government or public services
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Engineering & Computer Science (AREA)
  • Strategic Management (AREA)
  • Development Economics (AREA)
  • Educational Administration (AREA)
  • Economics (AREA)
  • Tourism & Hospitality (AREA)
  • General Physics & Mathematics (AREA)
  • Marketing (AREA)
  • Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Theoretical Computer Science (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Game Theory and Decision Science (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention discloses a method for determining regional atmospheric precipitation pattern index R, belongs to the technical field of atmospheric precipitation index analysis, and can solve the problems that the existing precipitation pattern index cannot reflect the change of the spatial and temporal distribution of precipitation and cannot reflect the amount and supply condition of regional aerial water vapor resources. The determination method comprises the following steps: acquiring precipitation data of a time period to be measured in a research area, and obtaining an area precipitation index M according to the precipitation data; acquiring the total water vapor conveying flux Q of the whole layer of atmosphere of the unit gas column in the research area, and obtaining an area water vapor conveying index T according to the total water vapor conveying flux Q; and calculating to obtain a regional atmospheric precipitation pattern index R according to the regional precipitation index M and the regional water vapor transmission index T.

Description

Method for determining regional atmospheric precipitation pattern index
Technical Field
The invention relates to a method for determining regional atmospheric precipitation pattern indexes, and belongs to the technical field of atmospheric precipitation index analysis.
Background
The change in regional precipitation can be viewed as two aspects: on the one hand the change in the total amount and on the other hand the supply of water vapour. The change of the total amount is relatively easy to quantify, and the existing research on regional precipitation is mainly on the precipitation amount, and the research on the concentration degree and the concentration period of the precipitation is less. Sufficient water vapor supply is needed for regional precipitation, water vapor resources in the air serve as potential water resources, and the water vapor content in the atmosphere and the transmission of the water vapor content are important components of regional precipitation and are closely related to regional atmospheric precipitation patterns.
At present, indexes reflecting regional precipitation patterns comprise regional annual average precipitation amount and annual average precipitation days, the indexes reveal the evolution rule of precipitation under a changing environment to a certain extent, but precipitation is one of important links of water circulation, the change of precipitation comprises the change of precipitation amount and space-time distribution, the possibility of flood and drought damage is increased when the space-time distribution of precipitation is excessively concentrated, the safe operation and the design implementation of watershed hydraulic engineering are threatened, and water supply, irrigation, power generation and the like in a region can be influenced by different degrees. The existing precipitation pattern indexes cannot reflect the change of the space-time distribution of precipitation and cannot reflect the amount and supply condition of regional aerial water vapor resources.
Disclosure of Invention
The invention provides a method for determining regional atmospheric precipitation pattern indexes, which can solve the problems that the existing precipitation pattern indexes cannot reflect the change of the spatial-temporal distribution of precipitation and cannot reflect the amount and supply conditions of regional aerial water vapor resources.
The invention provides a method for determining regional atmospheric precipitation pattern indexes, which comprises the following steps:
acquiring precipitation data of a time period to be measured in a research area, and obtaining an area precipitation index M according to the precipitation data;
acquiring the total water vapor conveying flux Q of the whole layer of atmosphere of the unit gas column in the research area, and obtaining an area water vapor conveying index T according to the total water vapor conveying flux Q;
and calculating to obtain a regional atmospheric precipitation pattern index R according to the regional precipitation index M and the regional water vapor transmission index T.
Optionally, the precipitation data includes daily precipitation of the study area in the time period to be measured;
obtaining a regional precipitation index M according to the precipitation data, specifically comprising:
dividing the daily rainfall into a plurality of levels from small to large, and obtaining the cumulative rainfall daily percentage and the cumulative rainfall percentage of each level;
drawing a relation curve of the cumulative rainfall days and the cumulative rainfall according to the cumulative rainfall day percentage and the cumulative rainfall percentage of each level in a plane rectangular coordinate system;
calculating to obtain a regional precipitation concentration index CI according to the relation curve and a quadrant bisector of a first quadrant of the plane rectangular coordinate system;
and calculating the regional precipitation index M according to the regional precipitation concentration index CI.
Optionally, the dividing the daily rainfall into a plurality of levels from small to large, and obtaining the cumulative rainfall percentage and the cumulative rainfall percentage of each level specifically includes:
dividing the daily rainfall into a plurality of levels from small to large to obtain precipitation days corresponding to the daily rainfall of each level, and recording as precipitation days;
acquiring the product of the rainfall days and the daily rainfall amount of the corresponding level, and recording the product as the rainfall amount of each level;
and obtaining the cumulative percentage of the rainfall days and the cumulative percentage of the rainfall capacity of each level according to the rainfall days and the rainfall capacity of each level.
Optionally, the regional precipitation concentration index CI is an area of a region surrounded by the relationship curve and a quadrant bisector of the first quadrant of the planar rectangular coordinate system.
Optionally, the precipitation data further includes a total precipitation day number of the time period to be measured in the research area;
the calculating the regional precipitation index M according to the regional precipitation concentration index CI specifically comprises:
obtaining the rainfall intensity I of the research area according to the daily rainfall and the total rainfall days;
calculating the standard deviation sigma of the rainfall intensity IIAnd the standard deviation sigma of the regional precipitation concentration index CICI；
Calculating the regional precipitation indicator M by a first formula, wherein the first formula is as follows:
Figure BDA0003471568700000031
calculating the value of the parameter n by a second formula, the second formula being:
Figure BDA0003471568700000032
optionally, the acquiring a total water vapor transport flux Q of the research area specifically includes:
calculating the total water vapor conveying flux Q according to a third formula, wherein the third formula is as follows:
Figure BDA0003471568700000033
wherein the content of the first and second substances,
Figure BDA0003471568700000034
Qλ＝Wm×u，
Figure BDA0003471568700000035
for radial transport of water vapor flux, QλFor transporting the water-vapor flux in the weft direction, WmThe water vapor content in the air column per unit area in the atmosphere, v is the latitudinal wind speed of the atmosphere in each layer, and u is the velocity of the air in each layerAtmospheric radial wind speed;
calculating the water vapor content W in the air column per unit area according to a fourth formulamThe fourth formula is:
Figure BDA0003471568700000036
wherein q is specific humidity, g is acceleration of gravity, and psIs the ground pressure and p is the pressure at the top of the atmosphere.
Optionally, the obtaining of the regional water vapor transmission index T according to the total water vapor transmission flux Q specifically includes:
calculating the water vapor transmission flux information entropy N according to the total water vapor transmission flux Qi；
According to the total water vapor transmission flux Q and the water vapor transmission flux information entropy NiAnd calculating the regional water vapor transmission index T.
Optionally, the entropy N of the water vapor transport flux information is calculated according to the total water vapor transport flux QiThe method specifically comprises the following steps:
calculating the water vapor transmission flux information entropy N according to a fifth formulaiThe fifth formula is:
Figure BDA0003471568700000041
wherein s is the time length of the time period to be measured, qi，jContributing rate to single precipitation water vapor transmission;
calculating the contribution rate q of single precipitation water vapor transmission according to a sixth formulai，jThe sixth formula is:
Figure BDA0003471568700000042
wherein Q isi，jFor the total water vapor transport flux, Q, of month j of year iiThe total water vapor transport flux for the ith year.
Optionally, the entropy N is obtained according to the total water vapor transport flux Q and the water vapor transport flux informationiCalculating the regional water vapor transmission index T, specifically comprising:
calculating the regional water vapor transport index T according to a seventh formula, wherein the seventh formula is as follows:
Figure BDA0003471568700000043
calculating the value of the parameter m by an eighth formula, the eighth formula being:
Figure BDA0003471568700000044
wherein σQFor the standard deviation of the total water vapor transport flux Q, σNEntropy N of flux information for said vapor deliveryiStandard deviation of (2).
Optionally, the calculating according to the regional precipitation index M and the regional water vapor transmission index T to obtain a regional atmospheric precipitation pattern index R specifically includes:
calculating the regional atmospheric precipitation pattern index R according to a ninth formula, wherein the ninth formula is as follows:
Figure BDA0003471568700000045
calculating the value of the parameter s by a tenth formula:
Figure BDA0003471568700000046
wherein σTFor the standard deviation of the regional water vapor transport index T, σMAnd the standard deviation of the regional atmospheric precipitation pattern index R. .
The invention can produce the beneficial effects that:
the inventionCalculating a regional atmospheric precipitation pattern index R through the regional precipitation index M and the regional water vapor transmission index T; when the regional rainfall index M is calculated, the rainfall intensity I is calculated and is used for describing parameters of rainfall total amount change conditions, and a regional rainfall concentration index CI is also calculated and is used for reflecting regional rainfall total amount change and rainfall concentration degree change; when the regional water vapor transmission index T is calculated, the total water vapor transmission flux Q and the water vapor transmission flux information entropy N are calculatediTherefore, the regional water vapor transmission index T can well reflect the regional water vapor transmission amount condition and the water vapor transmission stability.
The regional atmospheric precipitation pattern index R combines the existing index reflecting the regional precipitation pattern with the precipitation concentration index reflecting the regional precipitation space-time distribution change rule, the atmospheric water vapor total delivery quantity reflecting the regional water vapor delivery condition and the information entropy reflecting the regional water vapor delivery condition, so that the regional atmospheric precipitation pattern index can reflect the regional precipitation total quantity change, the regional precipitation space-time distribution and the regional water vapor delivery condition.
Drawings
FIG. 1 is a flow chart of a method for determining a regional atmospheric precipitation pattern indicator according to an embodiment of the present invention;
fig. 2 is a graph plotting cumulative raining days versus cumulative raining amounts.
Detailed Description
The present invention will be described in detail with reference to examples, but the present invention is not limited to these examples.
The embodiment of the invention provides a method for determining regional atmospheric precipitation pattern indexes, which comprises the following steps of:
and S1, acquiring precipitation data of the research area in the time period to be measured, and obtaining an area precipitation index M according to the precipitation data.
And S2, acquiring the total water vapor transmission flux Q of the whole layer of atmosphere of the unit gas column in the research area, and obtaining an area water vapor transmission index T according to the total water vapor transmission flux Q.
And S3, calculating according to the regional precipitation index M and the regional water vapor transmission index T to obtain a regional atmospheric precipitation pattern index R.
In this embodiment, the study area is the south minor gully basin of west peak of the province of Gansu, and the time periods to be measured are 2016 years and 2017 years. The precipitation data comprises daily precipitation and total precipitation days of 2016 and 2017 sections of the west-peak minor sulcus basin in Gansu province.
Obtaining regional precipitation index M according to precipitation data in S1, specifically including:
and S11, eliminating data with the daily rainfall being less than 0.1mm, dividing the residual daily rainfall into a plurality of levels from small to large, obtaining the precipitation days corresponding to the daily rainfall of each level, and recording as the precipitation days.
And S12, acquiring the product of the rainfall days and the daily rainfall amount of the corresponding level, and recording the product as the rainfall amount of each level.
For example: a certain level of daily precipitation in 2016 is 2mm, and a certain level of daily precipitation in 2016 is 6 days of precipitation with 2mm of daily precipitation, then the level of precipitation is 12 mm.
And S13, obtaining the cumulative percentage of rainfall days and the cumulative percentage of rainfall capacity of each level according to the rainfall days and the rainfall capacity of each level.
Specifically, the rainfall days of each level and the levels before the level are accumulated to obtain the accumulated rainfall days of each level; and accumulating the rainfall of each level and the level before the level to obtain the accumulated rainfall.
Dividing the rainfall day of each level by the total rainfall day of the time period to be measured to obtain the rainfall day percentage of the level, and accumulating the rainfall day percentages of the level and the previous levels to obtain the cumulative rainfall day percentage of the level; and dividing the rainfall of each level by the total rainfall of the time period to be measured to obtain the rainfall percentage of the level, and accumulating the rainfall percentages of the level and the previous levels to obtain the accumulated rainfall percentage of the level.
In this embodiment, the divided levels and their corresponding raining days, cumulative raining days, rainfall, cumulative rainfall percentage, and cumulative rainfall percentage are shown in table 1.
TABLE 1 calculation table of relation curve between cumulative raining day and cumulative raining quantity of time period to be measured in research area
Rank of Rainy day Accumulated rainy day Amount of rainfall Cumulative rainfall Cumulative percentage of rainy day Cumulative percentage of rainfall
1 11 11 11 11 11.3 0.2
2 5 16 10 21 16.5 0.3
3 4 20 12 33 20.6 0.5
4 1 21 4 37 21.6 0.5
5 4 25 20 57 25.8 0.8
6 1 26 6 63 26.8 0.9
7 3 29 21 84 29.9 1.2
9 2 31 18 102 32.0 1.5
10 1 32 10 112 33.0 1.7
11 2 34 22 134 35.1 2.0
14 3 37 42 176 38.1 2.6
15 2 39 30 206 40.2 3.0
17 1 40 17 223 41.2 3.3
18 1 41 18 241 42.3 3.6
20 2 43 40 281 44.3 4.2
22 2 45 44 325 46.4 4.8
25 1 46 25 350 47.4 5.2
26 1 47 26 376 48.5 5.6
27 2 49 54 430 50.5 6.4
28 5 54 140 570 55.7 8.4
30 1 55 30 600 56.7 8.9
And S14, drawing a relation curve of the accumulated rainfall days and the accumulated rainfall according to the accumulated rainfall day percentage and the accumulated rainfall percentage of each level in a plane rectangular coordinate system.
Specifically, the cumulative rainfall percentage value is used as an X coordinate, the cumulative rainfall percentage value is used as a Y coordinate, points corresponding to each level are drawn in a planar rectangular coordinate system, and the points corresponding to each level are sequentially connected by using a smooth curve, so that a relation curve between the cumulative rainfall and the cumulative rainfall is obtained. The relation curve is shown in fig. 2, the curve in fig. 2 is a relation curve of the accumulated rainfall days and the accumulated rainfall, and the straight line is a quadrant bisector of the first quadrant of the plane rectangular coordinate system.
And S15, calculating according to the relation curve and a quadrant bisector of the first quadrant of the plane rectangular coordinate system to obtain the regional precipitation concentration index CI.
As shown in fig. 2, the regional precipitation concentration index CI is the area of a region surrounded by the relationship curve and a quadrant bisector of the first quadrant of the rectangular plane coordinate system. A larger area indicates that the precipitation in the area of interest is more dispersed and more susceptible to extreme precipitation or extreme drought events.
Specifically, the cumulative percentage of rainfall is denoted as Y, and the cumulative percentage of rainfall is denoted as X, and as can be seen from table 1 and fig. 2, X and Y satisfy the relationship of the exponential distribution curve, and the relationship curve can be expressed as: y ═ axexp (bx) dX.
Where a and b are constants determined by the least square method.
According to
Figure BDA0003471568700000071
The value of CI is determined.
And S16, calculating the regional precipitation index M according to the regional precipitation concentration index CI.
Calculating a regional precipitation index M according to the regional precipitation concentration index CI, and specifically comprising the following steps:
and S161, obtaining the rainfall intensity I of the research area according to the daily rainfall and the total rainfall days.
Specifically, the sum of daily precipitation in a target time period of a research area is obtained and recorded as total precipitation, and rainfall intensity I is the total precipitation divided by the total precipitation day number.
S162, calculating standard deviation sigma of rainfall intensity IIStandard deviation sigma of sum regional precipitation concentration index CICI；
Calculating the regional precipitation index M through a first formula, wherein the first formula is as follows:
Figure BDA0003471568700000081
calculating the value of the parameter n through a second formula, namely making the rainfall intensity I and the rainfall concentration index CI have equal influence on the regional rainfall index M, and calculating the value of the parameter n, wherein the second formula is as follows:
Figure BDA0003471568700000082
the calculation of the regional precipitation index M comprises the calculation of rainfall intensity I and a precipitation concentration index CI, the rainfall intensity I is used for describing the change condition of the total rainfall, and the precipitation concentration index CI is used for describing the precipitation concentration degree, so that the regional precipitation index M can well reflect the change of the total rainfall in the region and the change of the precipitation concentration degree.
Acquiring a total water vapor transport flux Q of the research area in S2, specifically including:
s21, calculating the total water vapor conveying flux Q according to a third formula, wherein the third formula is as follows:
Figure BDA0003471568700000083
wherein the content of the first and second substances,
Figure BDA0003471568700000084
Qλ＝Wm×u，
Figure BDA0003471568700000085
for radial transport of water vapor flux, QλFor transporting the water-vapor flux in the weft direction, WmThe water vapor content in the air column per unit area in the atmosphere, v is the latitudinal wind speed of each layer of atmosphere, and u is the radial wind speed of each layer of atmosphere.
Calculating the water vapor content W in the air column per unit area in the atmosphere according to a fourth formulamThe fourth formula is:
Figure BDA0003471568700000086
wherein q is specific humidity, g is acceleration of gravity, and psIs the ground pressure and p is the pressure at the top of the atmosphere.
In this example, ps1000hPa is taken, and 300hPa is taken as p.
The specific humidity data for 2016 and 2017 in this example are shown in table 2:
TABLE 2 specific humidity data sheet
Figure BDA0003471568700000091
The radial wind speed data for 2016 and 2017 in this example are shown in Table 3:
TABLE 3 radial wind speed data sheet
Figure BDA0003471568700000101
The latitudinal wind speed data in 2016 and 2017 in this example are shown in Table 4:
TABLE 4 latitudinal wind speed data sheet
Figure BDA0003471568700000102
Obtaining a regional water vapor transmission index T according to the total water vapor transmission flux Q in S2, specifically including:
s22, calculating water vapor transmission flux information entropy N according to total water vapor transmission flux Qi。
S23, carrying out entropy N according to total water vapor transmission flux Q and water vapor transmission flux informationiAnd calculating the regional water vapor transmission index T.
Wherein the content of the first and second substances,s22, calculating water vapor transmission flux information entropy N according to total water vapor transmission flux QiThe method specifically comprises the following steps:
calculating the entropy N of the water vapor transport flux information according to a fifth formulaiThe fifth formula is:
Figure BDA0003471568700000111
wherein s is the time length of the time period to be measured, i.e. the number of months of the time period to be measured, qi，jContributing rate to single precipitation water vapor transmission;
calculating the contribution rate q of single precipitation water vapor transmission according to a sixth formulai，jThe sixth formula is:
Figure BDA0003471568700000112
wherein Q isi，jFor the total water vapor transport flux, Q, of month j of year iiThe total water vapor transport flux for the ith year.
Wherein, the water vapor transport flux information entropy NiThe values of (A) are as follows: n is not less than 0i≤1。
S23, according to the total water vapor transmission flux Q and the water vapor transmission flux information entropy NiCalculating a regional water vapor transmission index T, specifically comprising:
calculating the regional water vapor transport index T according to a seventh formula, wherein the seventh formula is as follows:
Figure BDA0003471568700000113
calculating the value of the parameter m by an eighth formula, namely ordering the total atmospheric water vapor transmission flux Q and the water vapor transmission flux information entropy N of the whole layer of the unit gas columniThe influence on the regional water vapor transmission index T is equal, the value of the parameter m is calculated, and the eighth formula is as follows:
Figure BDA0003471568700000114
wherein σQFor the standard deviation, σ, of the total water vapour transport flux QNEntropy N of flux information for vapor transportiStandard deviation of (2).
S3, calculating according to the regional precipitation index M and the regional water vapor transmission index T to obtain a regional atmospheric precipitation pattern index R, and concretely comprising:
calculating the regional atmospheric precipitation pattern index R according to a ninth formula, wherein the ninth formula is as follows:
Figure BDA0003471568700000115
calculating the value of the parameter s through a tenth formula, namely making the influence of the regional precipitation index M and the regional water vapor transmission index T on the regional atmospheric precipitation pattern index R equal, and calculating the value of the parameter s, wherein the tenth formula is as follows:
Figure BDA0003471568700000121
wherein σTStandard deviation, sigma, of regional water vapor transport index TMThe standard deviation of the regional atmospheric precipitation pattern index R.
In this embodiment, the calculation results of the regional precipitation index M, the regional water vapor transport index T, and the regional atmospheric precipitation pattern index R in 2016 and 2017 are shown in table 5:
TABLE 5 results of calculation
Year of year Regional precipitation indicator Regional water vapor transport indicator Regional precipitation climate pattern index
2017 2843.44 52.25 559.0073138
2016 1900.79 89.26 665.3755716
Standard deviation of 14.23505451 10.8163807 0.896930851
According to the method, a regional atmospheric precipitation pattern index R is obtained through calculation of a regional precipitation index M and a regional water vapor transmission index T; when the regional rainfall index M is calculated, the rainfall intensity I is calculated and is used for describing parameters of rainfall total amount change conditions, and a regional rainfall concentration index CI is also calculated and is used for reflecting regional rainfall total amount change and rainfall concentration degree change; when the regional water vapor transmission index T is calculated, the total water vapor transmission flux Q and the water vapor transmission flux information entropy N are calculatediTherefore, the regional water vapor transmission index T can well reflect the regional water vapor transmission amount condition and the water vapor transmission stability.
The regional atmospheric precipitation pattern index R combines the existing index reflecting the regional precipitation pattern with the precipitation concentration index reflecting the regional precipitation space-time distribution change rule, the atmospheric water vapor total delivery quantity reflecting the regional water vapor delivery condition and the information entropy reflecting the regional water vapor delivery condition, so that the regional atmospheric precipitation pattern index can reflect the regional precipitation total quantity change, the regional precipitation space-time distribution and the regional water vapor delivery condition.
Although the present application has been described with reference to a preferred embodiment, it should be understood that various changes, substitutions and alterations can be made herein without departing from the spirit and scope of the application.

Claims (10)

1. A method for determining a regional atmospheric precipitation pattern indicator, the method comprising:
acquiring precipitation data of a time period to be measured in a research area, and obtaining an area precipitation index M according to the precipitation data;
acquiring the total water vapor conveying flux Q of the whole layer of atmosphere of the unit gas column in the research area, and obtaining an area water vapor conveying index T according to the total water vapor conveying flux Q;
and calculating to obtain a regional atmospheric precipitation pattern index R according to the regional precipitation index M and the regional water vapor transmission index T.
2. The method of claim 1, wherein the precipitation data comprises daily precipitation for a period of time to be measured for the area of interest;
obtaining a regional precipitation index M according to the precipitation data, specifically comprising:
dividing the daily rainfall into a plurality of levels from small to large, and obtaining the cumulative rainfall daily percentage and the cumulative rainfall percentage of each level;
drawing a relation curve of the cumulative rainfall days and the cumulative rainfall according to the cumulative rainfall day percentage and the cumulative rainfall percentage of each level in a plane rectangular coordinate system;
calculating to obtain a regional precipitation concentration index CI according to the relation curve and a quadrant bisector of a first quadrant of the plane rectangular coordinate system;
and calculating the regional precipitation index M according to the regional precipitation concentration index CI.
3. The method according to claim 2, wherein the step of dividing the daily rainfall into a plurality of levels from small to large and obtaining an accumulated rainfall percentage and an accumulated rainfall percentage of each level specifically comprises:
dividing the daily rainfall into a plurality of levels from small to large to obtain precipitation days corresponding to the daily rainfall of each level, and recording as precipitation days;
acquiring the product of the rainfall days and the daily rainfall amount of the corresponding level, and recording the product as the rainfall amount of each level;
and obtaining the cumulative percentage of the rainfall days and the cumulative percentage of the rainfall capacity of each level according to the rainfall days and the rainfall capacity of each level.
4. The method of claim 3, wherein the regional precipitation concentration index CI is an area of a region surrounded by the relationship curve and a quadrant bisector of the first quadrant of the rectangular plane coordinate system.
5. A method of determining according to any one of claims 2 to 4, wherein the precipitation data further includes a total number of precipitation days for a period of time to be measured for the area of investigation;
the calculating the regional precipitation index M according to the regional precipitation concentration index CI specifically comprises:
obtaining the rainfall intensity I of the research area according to the daily rainfall and the total rainfall days;
calculating the standard deviation sigma of the rainfall intensity IIAnd the standard deviation sigma of the regional precipitation concentration index CICI；
Calculating the regional precipitation indicator M by a first formula, wherein the first formula is as follows:
Figure FDA0003471568690000021
calculating the value of the parameter n by a second formula, the second formula being:
Figure FDA0003471568690000022
6. the determination method according to claim 1, wherein the acquiring of the total water vapor transport flux Q of the investigation region specifically comprises:
calculating the total water vapor conveying flux Q according to a third formula, wherein the third formula is as follows:
Figure FDA0003471568690000023
wherein the content of the first and second substances,
Figure FDA0003471568690000024
Qλ＝Wm×u，
Figure FDA0003471568690000025
for radial transport of water vapor flux, QλFor transporting the water-vapor flux in the weft direction, WmThe water vapor content in the air column per unit area in the atmosphere, v is the latitudinal wind speed of each layer of atmosphere, and u is the radial wind speed of each layer of atmosphere;
calculating the water vapor content W in the air column per unit area according to a fourth formulamThe fourth formula is:
Figure FDA0003471568690000031
wherein q is specific humidity, g is acceleration of gravity, and psIs the ground pressure and p is the pressure at the top of the atmosphere.
7. The determination method according to claim 6, wherein the obtaining a regional water vapor transport indicator T according to the total water vapor transport flux Q specifically comprises:
calculating the water vapor transmission flux information entropy N according to the total water vapor transmission flux Qi；
According to the total water vapor transmission flux Q and the water vapor transmission flux information entropy NiAnd calculating the regional water vapor transmission index T.
8. The method of determining as claimed in claim 7, wherein said calculating a water vapor transport flux information entropy N based on said total water vapor transport flux QiThe method specifically comprises the following steps:
calculating the water vapor transmission flux information entropy N according to a fifth formulaiThe fifth formula is:
Figure FDA0003471568690000032
wherein s is the time length of the time period to be measured, qi，jContributing rate to single precipitation water vapor transmission;
calculating the contribution rate q of single precipitation water vapor transmission according to a sixth formulai，jThe sixth formula is:
Figure FDA0003471568690000033
wherein Q isi，jFor the total water vapor transport flux, Q, of month j of year iiThe total water vapor transport flux for the ith year.
9. The method of determining as claimed in claim 7 wherein said entropy of said total vapor transport flux Q and said entropy of vapor transport flux information N is based on said total vapor transport flux Q and said entropy of said total vapor transport flux NiCalculating the regional water vapor transmission index T, specifically comprising:
calculating the regional water vapor transport index T according to a seventh formula, wherein the seventh formula is as follows:
Figure FDA0003471568690000041
calculating the value of the parameter m by an eighth formula, the eighth formula being:
Figure FDA0003471568690000042
wherein σQFor the standard deviation of the total water vapor transport flux Q, σNEntropy N of flux information for said vapor deliveryiStandard deviation of (2).
10. The method for determining according to claim 1, wherein the calculating a regional atmospheric precipitation pattern indicator R from the regional precipitation indicator M and the regional water vapor transport indicator T comprises:
calculating the regional atmospheric precipitation pattern index R according to a ninth formula, wherein the ninth formula is as follows:
Figure FDA0003471568690000043
calculating the value of the parameter s by a tenth formula:
Figure FDA0003471568690000044
wherein σTFor the standard deviation of the regional water vapor transport index T, σMAnd the standard deviation of the regional atmospheric precipitation pattern index R.
CN202210044462.8A 2022-01-14 2022-01-14 Method for determining regional atmospheric precipitation pattern index Active CN114386860B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210044462.8A CN114386860B (en) 2022-01-14 2022-01-14 Method for determining regional atmospheric precipitation pattern index

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210044462.8A CN114386860B (en) 2022-01-14 2022-01-14 Method for determining regional atmospheric precipitation pattern index

Publications (2)

Publication Number Publication Date
CN114386860A true CN114386860A (en) 2022-04-22
CN114386860B CN114386860B (en) 2022-11-29

Family

ID=81202777

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210044462.8A Active CN114386860B (en) 2022-01-14 2022-01-14 Method for determining regional atmospheric precipitation pattern index

Country Status (1)

Country Link
CN (1) CN114386860B (en)

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH10288674A (en) * 1997-04-16 1998-10-27 Nippon Telegr & Teleph Corp <Ntt> Prediction method and device for rainfall pattern change
US20080178659A1 (en) * 2007-01-30 2008-07-31 Spinelli Charles B Methods and systems for measuring atmospheric water content
US20090107214A1 (en) * 2007-10-31 2009-04-30 Honeywell International Inc., Terahertz sensor to measure humidity and water vapor
JP2013224884A (en) * 2012-04-23 2013-10-31 Japan Radio Co Ltd Water vapor observation device and meteorological radar
CN109992746A (en) * 2019-03-11 2019-07-09 中国气象科学研究院 Atmosphere water substance total amount, total precipitable water and its corresponding precipitation efficiency calculation method
US20200109540A1 (en) * 2018-10-09 2020-04-09 Joseph Aoun Methods, systems to facilitate atmospheric water generation, and regulation of an environment of atmospheric water generation
CN112418684A (en) * 2020-11-26 2021-02-26 清华大学 Method, device, equipment and medium for evaluating space-time distribution rule of air water resource
CN112508352A (en) * 2020-11-18 2021-03-16 河北工程大学 Method for quantitatively distinguishing contributions of different factors in water circulation evolution process
CN113627690A (en) * 2021-09-03 2021-11-09 中国人民解放军国防科技大学 Method for predicting seasonal precipitation in southern China
CN113806948A (en) * 2021-09-23 2021-12-17 西安理工大学 Method and system for determining air water circulation influence factors

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH10288674A (en) * 1997-04-16 1998-10-27 Nippon Telegr & Teleph Corp <Ntt> Prediction method and device for rainfall pattern change
US20080178659A1 (en) * 2007-01-30 2008-07-31 Spinelli Charles B Methods and systems for measuring atmospheric water content
US20090107214A1 (en) * 2007-10-31 2009-04-30 Honeywell International Inc., Terahertz sensor to measure humidity and water vapor
JP2013224884A (en) * 2012-04-23 2013-10-31 Japan Radio Co Ltd Water vapor observation device and meteorological radar
US20200109540A1 (en) * 2018-10-09 2020-04-09 Joseph Aoun Methods, systems to facilitate atmospheric water generation, and regulation of an environment of atmospheric water generation
CN109992746A (en) * 2019-03-11 2019-07-09 中国气象科学研究院 Atmosphere water substance total amount, total precipitable water and its corresponding precipitation efficiency calculation method
CN112508352A (en) * 2020-11-18 2021-03-16 河北工程大学 Method for quantitatively distinguishing contributions of different factors in water circulation evolution process
CN112418684A (en) * 2020-11-26 2021-02-26 清华大学 Method, device, equipment and medium for evaluating space-time distribution rule of air water resource
CN113627690A (en) * 2021-09-03 2021-11-09 中国人民解放军国防科技大学 Method for predicting seasonal precipitation in southern China
CN113806948A (en) * 2021-09-23 2021-12-17 西安理工大学 Method and system for determining air water circulation influence factors

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
D.I.COOPER: "Spatial and temporal properties of water vapor and", 《AGRICULTURAL AND FOREST METEOROLOGY》 *
SIMIN YANG ETAL: "Characteristics of Agricultural Droughts and Spatial Stratified", 《ATMOSPHERE》 *
巩宁刚等: "1979-2016年祁连山地区大气水汽含量时空特征及其与降水的关系", 《干旱区地理》 *
王莹等: "近50年江苏省月降水分配格局的时空变化特征", 《人民长江》 *

Also Published As

Publication number Publication date
CN114386860B (en) 2022-11-29

Similar Documents

Publication Publication Date Title
Wu et al. Quantitative assessment and spatial characteristics analysis of agricultural drought vulnerability in China
Zhang et al. Changes in multiple ecosystem services between 2000 and 2013 and their driving factors in the Grazing Withdrawal Program, China
Mallakpour et al. Investigating the relationship between the frequency of flooding over the central United States and large-scale climate
Iizumi et al. Leveraging drought risk reduction for sustainable food, soil and climate via soil organic carbon sequestration
Santhi et al. Spatial calibration and temporal validation of flow for regional scale hydrologic modeling 1
Rocha et al. Assessing the impacts of sustainable agricultural practices for water quality improvements in the Vouga catchment (Portugal) using the SWAT model
Noory et al. Distributed agro-hydrological modeling with SWAP to improve water and salt management of the Voshmgir Irrigation and Drainage Network in Northern Iran
Aguayo et al. The glass half-empty: climate change drives lower freshwater input in the coastal system of the Chilean Northern Patagonia
Winchell et al. Using SWAT for sub-field identification of phosphorus critical source areas in a saturation excess runoff region
Xiao et al. Optimizing hotspot areas for ecological planning and management based on biodiversity and ecosystem services
Senti et al. Soil erosion, sediment yield and conservation practices assessment on Lake Haramaya Catchment
Lizárraga-Mendiola et al. Hydrological design of two low-impact development techniques in a semi-arid climate zone of central Mexico
Bai et al. Spatial-temporal variations of ecological vulnerability in the Tarim River Basin, Northwest China
Park et al. Evaluation of MODIS NDVI and LST for indicating soil moisture of forest areas based on SWAT modeling
Xiao et al. Assessing changes in water flow regulation in Chongqing region, China
Fischer et al. Hydrologic effects of climate change in a sub-basin of the Western Bug River, Western Ukraine
Jeong et al. Assessing the effects of indirect wastewater reuse on paddy irrigation in the Osan River watershed in Korea using the SWAT model
CN113011993A (en) Method for measuring and calculating water-entering load of agricultural pollution source based on standard data
Sehgal et al. Integrating climate forecasts with the soil and water assessment tool (SWAT) for high-Resolution hydrologic simulations and forecasts in the Southeastern US
Mamat et al. Ecological effect of the riparian ecosystem in the lower reaches of the Tarim River in northwest China
Seong et al. Hydroclimatic variability and change in the Chesapeake Bay Watershed
Chen et al. Spatial variation pattern analysis of hydrologic processes and water quality in Three Gorges Reservoir Area
Alitane et al. Towards a decision-making approach of sustainable water resources management based on hydrological modeling: a case study in central Morocco
Shetty et al. Comparing the hydrological performance of an irrigated native vegetation green roof with a conventional Sedum spp. green roof in New York City
CN114386860B (en) Method for determining regional atmospheric precipitation pattern index

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant