A method and system for inverting the height distribution of gravity wave disturbance parameters at the top of the intermediate layer

Abstract

The present invention proposes a method and system for inverting the height distribution of gravity wave disturbance quantity parameters at the top of the mesosphere. The present invention constructs an event sequence and height sequence of multiple meteors, and constructs a corresponding data matrix and drift matrix according to the parameter set of each event; the top of the mesosphere is equidistantly divided into multiple height intervals, and each meteor event is placed therein; in each height interval, the data matrix and drift matrix corresponding to all meteor events in the interval are sequentially connected in series to construct a corresponding joint data matrix and joint drift matrix; in each height interval, according to the corresponding two joint matrices, the corresponding gravity wave disturbance quantity parameter matrix is ​​obtained by the least squares method; all gravity wave disturbance quantity parameter matrices are sequentially connected in series to obtain their height distribution. The height distribution of gravity wave disturbance quantity parameters at the top of the mesosphere can be used to quantitatively analyze the changes in gravity waves and solve the dynamic effects of gravity waves on the background atmosphere. The present invention can be widely used in space physics related industries.

Description

Technical Field

[0001] The present invention relates to the field of space detection technology, and in particular to a method and system for inverting the height distribution of gravity wave disturbance parameters at the top of an intermediate layer. Background Art

[0002] The top of the mesosphere is an important region where matter and energy are exchanged between the Earth and the outer space, and gravity waves are one of the most important atmospheric fluctuations active at the top of the mesosphere. Most gravity waves originating from the lower atmosphere propagate upward to the top of the mesosphere, break and deposit the energy they carry into the background, thereby driving the structural evolution of the local space environment. Gravity waves can also lead to adiabatic cooling of the polar mesosphere in summer and adiabatic heating in winter, resulting in a residual circulation from the summer pole to the winter pole at the top of the mesosphere. Therefore, in order to quantitatively evaluate the effects of gravity waves on the dynamic and thermal structures of the background atmosphere and understand the vertical coupling process of the atmosphere and its circulation mechanism at the top of the mesosphere, it is necessary to establish more observations of gravity wave activities. On the other hand, in order to achieve a more realistic simulation of the global climate at the subgrid scale, more observational constraints are also needed to improve the existing gravity wave parameterization scheme.

[0003] The parameters of gravity wave disturbance and their evolution in altitude are crucial for quantitatively characterizing the impact of gravity waves on the background atmosphere. Their acquisition depends on continuous observation of the wind field at the top of the mesosphere. However, it is challenging to detect the atmosphere at the top of the mesosphere, especially in terms of wind field detection, which leads to the current scarcity of research on gravity waves and their related effects. Currently, only a limited number of lidars that generally observe only at night can detect the wind field at the top of the mesosphere. Airglow imagers can usually only provide wavelength and frequency information of gravity waves at a specific altitude, making it difficult to evaluate the evolution of gravity wave energy and momentum with altitude and its impact on the background. The development of satellite technology in the past few decades has made global detection of the atmosphere possible, but there are relatively few satellites capable of detecting the atmosphere at the top of the mesosphere; in addition, satellite observations can usually only resolve gravity waves with longer horizontal wavelengths (usually corresponding to lower natural frequencies); however, high-frequency gravity waves with shorter horizontal wavelengths contribute more than 70% to momentum transport at the top of the mesosphere. Therefore, continuous and reliable observations of wind fields and high-frequency, small-scale gravity waves at the top of the mesosphere are crucial for studying the effects of gravity waves and related mechanisms of atmospheric vertical coupling.

[0004] Radio radar (such as meteor radar, very high frequency radar, etc.) is an effective means of continuous measurement of atmospheric wind field in the mesosphere. Among them, very high frequency radar is generally a device with a large power aperture product, and the cost of operation and maintenance is relatively high.

[0005] Meteor radar is not limited by day and night and is not affected by weather. It is an advanced device that is widely deployed for continuous measurement of atmospheric wind fields and fluctuations. At present, there are related technologies based on meteor radar to invert the parameters of gravity waves at the top of the mesosphere. However, the processing of a large amount of meteor trail drift data and wind field data in the inversion process consumes too much computing resources, and there is currently a lack of mature and time-saving fast solution methods for gravity wave disturbance parameters and their height distribution. In practical applications, meteor radars have to collect and process massive amounts of observation data at all times, and also undertake important observation tasks such as meteor selection, meteor parameter solution and atmospheric wind field inversion. The original inversion method for the height distribution of gravity wave disturbance parameters contains a large number of cumulative operations, resulting in excessive occupation of computing resources, making it impossible for the radar system to further undertake the analysis of atmospheric dynamic processes related to high-frequency small-scale fluctuations. In order to reduce time and computing costs and realize real-time dynamic analysis of the impact of gravity waves on background atmospheric dynamic processes, it is necessary to develop a reliable inversion method and system for the height distribution of gravity wave disturbance parameters at the top of the mesosphere that can meet the requirements of real-time dynamic analysis. Summary of the invention

[0006] In view of the technical defects of the prior art, the present invention proposes a method and system for inverting the height distribution of gravity wave disturbance parameters at the top of the intermediate layer, so as to provide a reliable and fast solution to the height distribution of gravity wave disturbance parameters in the top area of ​​the intermediate layer.

[0007] The technical solution of the method of the present invention is a method for inverting the height distribution of gravity wave disturbance parameters at the top of the intermediate layer, comprising the following steps:

[0008] Step 1: Acquire multiple meteor events to construct a meteor event sequence. According to the zenith angle, azimuth, radial drift velocity disturbance of the trail and the height of each meteor event, construct a data matrix, a drift matrix and a sequence of the height of each meteor event;

[0009] Step 2: Divide the mesosphere top area into multiple height intervals at equal intervals, and place each meteor event in the corresponding height interval according to the height of each meteor event;

[0010] Step 3: In each altitude interval, the data matrix and drift matrix obtained in step 1 corresponding to all meteor events in the same altitude interval are vertically connected in a one-to-one correspondence to construct a joint data matrix and a joint drift matrix in the altitude interval;

[0011] Step 4: In each altitude interval, according to the corresponding joint data matrix and joint drift matrix, the gravity wave disturbance parameter matrix in the corresponding altitude interval is obtained by the least square fitting method;

[0012] Step 5: Based on all gravity wave disturbance parameter matrices, connect them in series in order of height to obtain the distribution matrix of gravity wave disturbance parameters at different heights of the mesosphere top.

[0013] Step 6: The distribution matrix of gravity wave disturbance parameters at different heights at the top of the mesosphere is used to quantitatively analyze the changes in the momentum flux of gravity waves at different heights, and through further calculation of the divergence of the gravity wave momentum flux, the dynamic effect of gravity waves on the background atmosphere at the top of the mesosphere, especially the drag caused by gravity waves on the wind field, is solved.

[0014] Preferably, the meteor event sequence in step 1 is specifically:

[0015] Met=[Met1 Met2 ... Met x ... Met N ]

[0016] Among them, Met represents the meteor event sequence, Met x represents the xth meteor event in the meteor event sequence, N represents the total number of meteor events in the meteor event sequence, x∈[1,N];

[0017] Step 1 constructs a data matrix corresponding to each meteor event as follows:

[0018] The zenith angle of the xth meteor event is θ x , the azimuth of the x-th meteor event is φ x , the radial drift velocity disturbance of the wake of the x-th meteor event is Vr′ x ;

[0019] The data matrix A of the x-th meteor event x for:

[0020]

[0021] Among them, A x is the data matrix of the x-th meteor event, x∈[1,N];

[0022] Step 1 constructs the drift matrix corresponding to each meteor event as follows:

[0023]

[0024] Among them, V x is the drift matrix of the x-th meteor event, x∈[1,N];

[0025] Assume that the sequence of the height of the meteor event is Ht, specifically:

[0026] Ht=[Ht1 Ht2 ... Ht x ... Ht N ]

[0027] Among them, Ht x Indicates the xth meteor event Met x The height of the corresponding meteor event, N represents the total number of meteor events in the meteor event sequence, x∈[1,N];

[0028] Preferably, in step 2, the top area of ​​the intermediate layer is divided into a plurality of height intervals at equal intervals, specifically:

[0029] The height scale is marked at equal intervals in the middle layer top area. Let the height scale sequence be H, which is:

[0030] H=[H1 H2 ... H a ... H M ]

[0031] Among them, H a represents the ath height scale in the middle layer top area, M represents the total number of height intervals in the height scale sequence H, a∈[1,M-1];

[0032] Assume that the height interval sequence corresponding to the height scale sequence H is I, specifically:

[0033] I=[I1 I2 ... I a ... I M-1 ]

[0034] Among them, I a Indicates the ath height interval in the height interval sequence, located at the height scale H a and height scale H a+1 , M-1 represents the total number of height intervals in the height interval sequence I, a∈[1,M-1];

[0035] The xth meteor event Met x The corresponding height of the meteor event Ht x , x∈[1,N], N represents the total number of meteor events in the meteor event sequence, and the specific calculation is as follows:

[0036] if:

[0037] H a ≤Ht x <H a+1

[0038] Then it will be combined with Ht x Corresponding Met x Placed in height range I a within;

[0039] All N meteor events in the meteor event sequence Met are processed according to the above process until each meteor event is placed within the corresponding height interval;

[0040] Preferably, step 3 constructs a joint data matrix and a joint drift matrix located in the height interval, specifically:

[0041] Let correspond to the height interval I a The meteor event interval set is MetS a , specifically:

[0042] MetS a =[Met a，1 Met a，2 ... Met a，m ... Met a，n ]

[0043] Among them, Met a，m Indicates that it is in altitude interval I a The mth meteor event within the range, n represents the meteor event in the height interval I a The total number of meteor events within , m∈[1,n];

[0044] Let the interval set corresponding to the meteor event MetS a The data matrix set is AS a , specifically:

[0045] AS a =[A a，1 A a，2 ... A a，m ... A a，n ]

[0046] Among them, A a，m Indicates that it is in altitude interval I a The mth meteor event Met a，m The corresponding data matrix, n represents the height interval I a The total number of meteor events within , m∈[1,n];

[0047] Let the interval set corresponding to the meteor event MetS a The drift matrix set is MS a , specifically:

[0048] VS a =[V a，1 V a，2 ... V a，m ... V a，n ]

[0049] Among them, V a，m Indicates that it is in altitude interval I a The drift matrix corresponding to the mth meteor event within the range of a The total number of meteor events within , m∈[1,n];

[0050] Suppose that the meteor event interval set MetS a The meteor events in the AS are vertically connected in a one-to-one correspondence manner. a The obtained height interval I a The joint data matrix within is ComA a , specifically:

[0051]

[0052] Suppose that the meteor event interval set MetS a The meteor events in the vertical series VS a The obtained height interval I a The joint drift matrix within is ComV a , specifically:

[0053]

[0054] All M-1 height intervals in the height interval sequence I are processed according to the above process to obtain all joint data matrices and joint drift matrices corresponding to each height interval. The sequences of the two are represented by ComA and ComV, respectively, as follows:

[0055] ComA=[ComA1 ComA2 ... ComA a ... ComA M-1 ]

[0056] ComV=[ComV1 ComV2 ... ComV a ... ComV M-1 ]

[0057] Preferably, in step 4, the least squares fitting method is used to obtain the gravity wave disturbance parameter matrix in the corresponding height interval, specifically:

[0058] Let correspond to the height interval I a The gravity wave disturbance parameter matrix is ​​GW a , specifically:

[0059] GW a =(Com a TCom a ) -1 Com a T ComV a

[0060] All M-1 height intervals in the height interval sequence I are processed according to the above process to obtain the parameter matrix of all gravity wave disturbances corresponding to each height interval. The sequence is represented by GW, specifically:

[0061] GW=[GW1 GW2 ... GW a ... GW M-1 ]

[0062] Preferably, the step 5 is as follows:

[0063] Let the discrete height scale sequence corresponding to the height interval sequence I be H′, specifically:

[0064] H′=[H′1 H′2 ... H′ a ... H′ M-1 ]

[0065] Among them, H′ a Indicates the height interval I corresponding to the ath height interval I in the height interval sequence I a The discrete height scale is:

[0066]

[0067] Assume that the distribution matrix of gravity wave disturbance parameters at different heights of the mesosphere top is GW_MAT, then

[0068]

[0069] Among them, H′ a Indicates the height interval I corresponding to the ath height interval I in the height interval sequence I a Discrete height scale, GW a corresponds to the height interval I a The gravity wave perturbation parameter matrix, a∈[1,M-1];

[0070] The present invention also provides a system for inverting the height distribution of the parameters of the gravity wave disturbance at the top of the intermediate layer, comprising:

[0071] The data matrix and drift matrix calculation module of the meteor event is used to obtain multiple meteor events to construct a meteor event sequence, and construct the data matrix and drift matrix corresponding to each meteor event, as well as the sequence of the height of the meteor event according to the zenith angle, azimuth angle, radial drift velocity disturbance of the trail and the height of each meteor event;

[0072] A meteor event height interval division module is used to divide the mesosphere top area into multiple height intervals at equal intervals, and place each meteor event in a corresponding height interval according to the height of each meteor event;

[0073] The joint data matrix and joint drift construction calculation module of the meteor event is used to vertically connect the data matrices and drift matrices corresponding to all meteor events in the same height interval in a one-to-one correspondence manner in each height interval to construct the joint data matrix and joint drift matrix in the height interval;

[0074] A gravity wave disturbance parameter matrix calculation module is used to obtain the gravity wave disturbance parameter matrix in the corresponding height interval by the least square fitting method according to the corresponding joint data matrix and joint drift matrix;

[0075] A distribution matrix construction module is used to connect all gravity wave disturbance parameter matrices in series in order of height to obtain the distribution matrix of gravity wave disturbance parameter at different heights of the mesosphere top;

[0076] The distribution matrix application module is used for the distribution matrix of gravity wave disturbance parameters at different heights at the top of the mesosphere, and is used to quantitatively analyze the changes in the momentum flux of gravity waves at different heights. Through subsequent further calculation of the divergence of the gravity wave momentum flux, the dynamic effect of gravity waves on the background atmosphere at the top of the mesosphere, especially the drag caused by gravity waves on the wind field, is solved.

[0077] Compared with the prior art, the present invention has the following advantages and beneficial effects:

[0078] This method can be widely used in space exploration, space physics and other related industries.

[0079] This method can invert the distribution matrix of gravity wave disturbance parameters at different heights at the top of the mesosphere, and the calculation time is significantly shortened. The numerical value is within the internationally recognized reasonable range and is very close to the calculation results of traditional methods. Although there is an error between the two, the error percentage is small enough to be negligible. BRIEF DESCRIPTION OF THE DRAWINGS

[0080] Figure 1 : is a schematic flow chart of the method of the present invention;

[0081] Figure 2 : A comparison of the gravity wave momentum flux height distribution obtained by this method and the traditional method; DETAILED DESCRIPTION

[0082] The specific embodiments of the present invention will be described in detail below. In order to avoid too much unnecessary details, the well-known structures or functions will not be described in detail in the following examples. The approximate language used in the following examples can be used for quantitative expression, indicating that the quantity can be allowed to have a certain change without changing the basic function. Unless otherwise defined, the technical and scientific terms used in the following examples have the same meaning as those generally understood by those skilled in the art to which the present invention belongs.

[0083] In order to overcome the limitations of the prior art, the present invention proposes a method for rapidly inverting the height distribution of gravity wave disturbance parameters at the top of the mesosphere. This method constructs corresponding data matrices and drift matrices for all meteor events with good data quality and complete detection information identified by radar according to their respective parameter sets; places each meteor event in the corresponding height interval of the top of the mesosphere according to its height; connects the data matrices and drift matrices of all meteor events in the same height interval in series in a one-to-one correspondence to construct a joint data matrix and a joint drift matrix in the height interval; obtains the gravity wave disturbance parameter matrix in the corresponding height interval by the least squares fitting method based on the two joint matrices obtained in the previous step; connects all gravity wave disturbance parameter matrices in series in height order to obtain the height distribution matrix of gravity wave disturbance parameters at the top of the mesosphere. The present invention effectively saves computing resources while obtaining reliable results through inversion.

[0084] An embodiment of the present invention is an inversion of the height distribution of the parameters of the mesosphere top gravity wave disturbance quantity based on the meteor detection data of the meteor radar at the YC station at 40° north latitude on September 30, 2020.

[0085] Combine the following Figure 1 to Figure 2 The specific implementation modes of the present invention are introduced.

[0086] The specific embodiment of the method of the present invention is a method for inverting the height distribution of gravity wave disturbance parameters at the top of the intermediate layer, comprising the following steps:

[0087] Step 1: Acquire multiple meteor events to construct a meteor event sequence. According to the zenith angle, azimuth, radial drift velocity disturbance of the trail and the height of each meteor event, construct a data matrix, a drift matrix and a sequence of the height of each meteor event;

[0088] The meteor event sequence described in step 1 is specifically:

[0089] Met=[Met1 Met2 ... Met x ... Met N ]

[0090] Among them, Met represents the meteor event sequence, Metx represents the xth meteor event in the meteor event sequence, N represents the total number of meteor events in the meteor event sequence (here, N = 8038), x∈[1, N];

[0091] Step 1 constructs a data matrix corresponding to each meteor event as follows:

[0092] The zenith angle of the xth meteor event is θ x , the azimuth of the x-th meteor event is φ x , the radial drift velocity disturbance of the wake of the x-th meteor event is Vr′ x ;

[0093] The data matrix A of the x-th meteor event x for:

[0094]

[0095] Among them, A x is the data matrix of the x-th meteor event, x∈[1,N];

[0096] Step 1 constructs the drift matrix corresponding to each meteor event as follows:

[0097]

[0098] Among them, V x is the drift matrix of the x-th meteor event, x∈[1,N];

[0099] Assume that the sequence of the height of the meteor event is Ht, specifically:

[0100] Ht=[Ht1 Ht2 ... Ht x ... Ht N ]

[0101] Among them, Ht x Indicates the xth meteor event Met x The height of the corresponding meteor event, N represents the total number of meteor events in the meteor event sequence, x∈[1,N];

[0102] Step 2: Divide the mesosphere top area into multiple height intervals at equal intervals, and place each meteor event in the corresponding height interval according to the height of each meteor event;

[0103] Step 2 divides the middle layer top area into multiple height intervals at equal intervals, specifically:

[0104] The height scale is marked at equal intervals in the middle layer top area. Let the height scale sequence be H, which is:

[0105] H=[H1 H2 ... H a ... H M ]

[0106] Among them, H a represents the ath height scale in the middle layer top area, M represents the total number of height intervals in the height scale sequence H (here, M = 8), a∈[1, M-1];

[0107] Assume that the height interval sequence corresponding to the height scale sequence H is I, specifically:

[0108] I=[I1 I2 ... I a ... I M-1 ]

[0109] Among them, I a Indicates the ath height interval in the height interval sequence, located at the height scale H a and height scale H a+1 , M-1 represents the total number of height intervals in the height interval sequence I, a∈[1,M-1];

[0110] The xth meteor event Met x The corresponding height of the meteor event Ht x , x∈[1,N], N represents the total number of meteor events in the meteor event sequence, and the specific calculation is as follows:

[0111] if:

[0112] H a ≤Ht x <H a+1

[0113] Then it will be combined with Ht x Corresponding Met x Place in height range I a within;

[0114] All N meteor events in the meteor event sequence Met are processed according to the above process until each meteor event is placed within the corresponding height interval;

[0115] Step 3: In each altitude interval, the data matrix and drift matrix obtained in step 1 corresponding to all meteor events in the same altitude interval are vertically connected in a one-to-one correspondence to construct a joint data matrix and a joint drift matrix in the altitude interval;

[0116] Step 3 constructs the joint data matrix and the joint drift matrix in the height interval, specifically:

[0117] Let correspond to the height interval Ia The meteor event interval set is MetS a , specifically:

[0118] MetS a =[Met a，1 Met a，2 ... Met a，m ... Met a，n ]

[0119] Among them, Met a，m Indicates that it is in altitude interval I a The mth meteor event within the range, n represents the meteor event in the height interval I a The total number of meteor events within , m∈[1,n];

[0120] Let the interval set corresponding to the meteor event MetS a The data matrix set is AS a , specifically:

[0121] AS a =[A a，1 A a，2 ... A a，m ... A a，n ]

[0122] Among them, A a，m Indicates that it is in altitude interval I a The mth meteor event Met a，m The corresponding data matrix, n represents the height interval I a The total number of meteor events within , m∈[1,n];

[0123] Let the interval set corresponding to the meteor event MetS a The drift matrix set is VS a , specifically:

[0124] VS a =[V a，1 V a，2 ... V a，m ... V a，n ]

[0125] Among them, V a，m Indicates that it is in altitude interval I a The drift matrix corresponding to the mth meteor event within the range of a The total number of meteor events within , m∈[1,n];

[0126] Suppose that the meteor event interval set MetS aThe meteor events in the AS are vertically connected in a one-to-one correspondence manner. a The obtained height interval I a The joint data matrix within is ComA a , specifically:

[0127]

[0128] Suppose that the meteor event interval set MetS a The meteor events in the vertical series VS a The obtained height interval I a The joint drift matrix within is ComV a , specifically:

[0129]

[0130] All M-1 height intervals in the height interval sequence I are processed according to the above process to obtain all joint data matrices and joint drift matrices corresponding to each height interval. The sequences of the two are represented by ComA and ComV, respectively, as follows:

[0131] ComA=[ComA1 ComA2 ... ComA a ... ComA M-1 ]

[0132] ComV=[ComV1 ComV2 ... ComV a ... ComV M-1 ]

[0133] Step 4: In each altitude interval, according to the corresponding joint data matrix and joint drift matrix, the gravity wave disturbance parameter matrix in the corresponding altitude interval is obtained by the least square fitting method;

[0134] In step 4, the least squares fitting method is used to obtain the gravity wave disturbance parameter matrix in the corresponding height interval, which is specifically:

[0135] Let correspond to the height interval I a The gravity wave disturbance parameter matrix is ​​GW a , specifically:

[0136] GW a =(Com a T Com a ) -1 Com a T ComV a

[0137] All M-1 height intervals in the height interval sequence I are processed according to the above process to obtain the parameter matrix of all gravity wave disturbances corresponding to each height interval. The sequence is represented by GW, specifically:

[0138] GW=[GW1 GW2 ... GW a ... GW M-1 ]

[0139] Step 5: Based on all gravity wave disturbance parameter matrices, connect them in series in order of height to obtain the distribution matrix of gravity wave disturbance parameters at different heights of the mesosphere top.

[0140] The step 5 is specifically as follows:

[0141] Let the discrete height scale sequence corresponding to the height interval sequence I be H′, specifically:

[0142] H′=[H′1 H′2 ... H′ a ... H′ M-1 ]

[0143] Among them, H′ a Indicates the height interval I corresponding to the ath height interval I in the height interval sequence I a The discrete height scale is:

[0144]

[0145] Assume that the distribution matrix of gravity wave disturbance parameters at different heights of the mesosphere top is GW_MAT, then

[0146]

[0147] Among them, H′ a Indicates the height interval I corresponding to the ath height interval I in the height interval sequence I a Discrete height scale, GW a corresponds to the height interval I a The gravity wave perturbation parameter matrix, a∈[1,M-1];

[0148] The solution is:

[0149]

[0150] The calculation time is 1.16×10 -2 Second;

[0151] The GW_MAT′ obtained by the traditional method is:

[0152]

[0153] The calculation time is 2.15×10 -2 Second;

[0154] Step 6: The distribution matrix of gravity wave disturbance parameters at different heights at the top of the mesosphere can be used to quantitatively analyze the changes in the momentum flux of gravity waves at different heights, and through subsequent further calculation of the divergence of the gravity wave momentum flux, the dynamic effect of gravity waves on the background atmosphere at the top of the mesosphere, especially the drag caused by gravity waves on the wind field, can be solved.

[0155] The system embodiment of the present invention provides a system for inverting the height distribution of gravity wave disturbance parameters at the top of the intermediate layer, including:

[0156] The data matrix and drift matrix calculation module of the meteor event is used to obtain multiple meteor events to construct a meteor event sequence, and construct the data matrix and drift matrix corresponding to each meteor event, as well as the sequence of the height of the meteor event according to the zenith angle, azimuth angle, radial drift velocity disturbance of the trail and the height of each meteor event;

[0157] A meteor event height interval division module is used to divide the mesosphere top area into multiple height intervals at equal intervals, and place each meteor event in a corresponding height interval according to the height of each meteor event;

[0158] The joint data matrix and joint drift construction calculation module of the meteor event is used to vertically connect the data matrices and drift matrices corresponding to all meteor events in the same height interval in a one-to-one correspondence manner in each height interval to construct the joint data matrix and joint drift matrix in the height interval;

[0159] A gravity wave disturbance parameter matrix calculation module is used to obtain the gravity wave disturbance parameter matrix in the corresponding height interval by the least square fitting method according to the corresponding joint data matrix and joint drift matrix;

[0160] A distribution matrix construction module is used to connect all gravity wave disturbance parameter matrices in series in order of height to obtain the distribution matrix of gravity wave disturbance parameter at different heights of the mesosphere top;

[0161] The distribution matrix application module is used for the distribution matrix of gravity wave disturbance parameters at different heights above the mesosphere. It can be used to quantitatively analyze the changes in the momentum flux of gravity waves at different heights, and through the subsequent calculation of the divergence of the gravity wave momentum flux, it can solve the dynamic effects of gravity waves on the background atmosphere above the mesosphere, especially the drag caused by gravity waves on the wind field.

[0162] The data matrix of the meteor event, the drift matrix calculation module, the meteor event height interval division module, the joint data matrix of the meteor event, the joint drift construction calculation module, the gravity wave disturbance parameter matrix calculation module, the distribution matrix construction module, and the distribution matrix application module are all deployed on the server.

[0163] The comparison of gravity wave momentum flux height distribution obtained by the two methods is as follows: Figure 2 As shown;

[0164] In summary, the method in the embodiment can invert and obtain the distribution matrix of gravity wave disturbance parameters at different heights of the intermediate layer top that is very close to the calculation result of the traditional method, and the calculation time is significantly shortened by 46%.

[0165] The errors of the gravity wave disturbance parameter matrices obtained by the two methods are shown in Table 1.

[0166] The error percentages of the gravity wave perturbation parameter matrices obtained by the two methods are shown in Table 2.

[0167] Table 1: Errors in the gravity wave perturbation parameter matrix obtained by the two methods

[0168]

[0169] Table 2: Percent error of the gravity wave disturbance parameter matrix obtained by the two methods

[0170]

[0171] Combining Table 1 and Table 2, we can see that below 92 km, the error is less than 10 m. 2 / s 2 , and 92% of the error percentages are concentrated within 10%.

[0172] During specific implementation, the present invention can use computer software technology to realize automatic operation process.

[0173] The embodiments of the present invention are described in detail above, but the contents are only preferred embodiments of the present invention and are not intended to limit the present invention. Any modifications, equivalent substitutions and improvements made within the scope of the present invention shall be included in the protection scope of the present invention.

Claims

1. A method for inverting the height distribution of gravity wave disturbance parameters at the top of the intermediate layer, characterized in that: The following steps are involved: Step 1: Acquire multiple meteor events to construct a meteor event sequence. According to the zenith angle, azimuth, radial drift velocity disturbance of the trail and the height of each meteor event, construct a data matrix, a drift matrix and a sequence of the height of each meteor event; Step 2: Divide the mesosphere top area into multiple height intervals at equal intervals, and place each meteor event in the corresponding height interval according to the height of each meteor event; Step 3: In each altitude interval, the data matrix and drift matrix obtained in step 1 corresponding to all meteor events in the same altitude interval are vertically connected in a one-to-one correspondence to construct a joint data matrix and a joint drift matrix in the altitude interval; Step 4: In each altitude interval, according to the corresponding joint data matrix and joint drift matrix, the gravity wave disturbance parameter matrix in the corresponding altitude interval is obtained by the least square fitting method; Step 5: According to all gravity wave disturbance parameter matrices, connect them in series in order of height to obtain the distribution matrix of gravity wave disturbance parameters at different heights of the mesosphere top; Step 6: The distribution matrix of gravity wave disturbance parameters at different heights of the mesosphere is used to quantitatively analyze the changes of gravity wave momentum flux at different heights, and further calculate the divergence of gravity wave momentum flux to solve the dynamic effect of gravity waves on the background atmosphere of the mesosphere, including the drag caused by gravity waves on the wind field; The meteor event sequence described in step 1 is specifically: With=[With1 With2 ... With x ... Of N ] Among them, Met represents the meteor event sequence, Met x represents the xth meteor event in the meteor event sequence, N represents the total number of meteor events in the meteor event sequence, x∈[1,N]; Step 1 constructs a data matrix corresponding to each meteor event as follows: The zenith angle of the xth meteor event is θ x , the azimuth of the x-th meteor event is φ x , the radial drift velocity disturbance of the wake of the xth meteor event is Vr' x ; The data matrix A of the x-th meteor event x for: Among them, A x is the data matrix of the x-th meteor event, x∈[1,N]; Step 1 constructs the drift matrix corresponding to each meteor event as follows: Among them, V x is the drift matrix of the x-th meteor event, x∈[1,N]; Assume that the sequence of the height of the meteor event is Ht, specifically: Ht=[Ht1 Ht2 ... Ht x ... Ht N ] Among them, Ht x Indicates the xth meteor event Met x The height of the corresponding meteor event, N represents the total number of meteor events in the meteor event sequence, x∈[1,N].

2. The method for inverting the height distribution of the gravity wave disturbance parameters at the top of the intermediate layer according to claim 1 is characterized by: Step 2 divides the middle layer top area into multiple height intervals at equal intervals, specifically: The height scale is marked at equal intervals on the top area of ​​the middle layer. Let the height scale sequence be H, which is: H=[H1 H2 ... H a ... H M ] Among them, H a represents the ath height scale in the middle layer top area, M represents the total number of height intervals in the height scale sequence H, a∈[1,M-1]; Assume that the height interval sequence corresponding to the height scale sequence H is I, specifically: I=[I1 I2 ... I a ... I M-1 ] Among them, I a Indicates the ath height interval in the height interval sequence, located at the height scale H a and height scale H a+1 , M-1 represents the total number of height intervals in the height interval sequence I, a∈[1,M-1]; The xth meteor event Met x The corresponding height of the meteor event Ht x ,x∈[1,N], N represents the total number of meteor events in the meteor event sequence, and the specific calculation is as follows: if: H a ≤Ht x <H a+1 Then it will be combined with Ht x Corresponding Met x Placed in height range I a within; All N meteor events in the meteor event sequence Met are processed according to the above process until each meteor event is placed within the corresponding height interval.

3. The method for inverting the height distribution of the parameters of the gravity wave disturbance at the top of the intermediate layer according to claim 2 is characterized by: Step 3 constructs the joint data matrix and the joint drift matrix in the height interval, specifically: Let correspond to the height interval I a The meteor event interval set is MetS a , specifically: WithS a =[With a,1 Of a,2 ... Of a,m ... Of a,n ] Among them, Met a , m Indicates that it is in altitude interval I a The mth meteor event within the range, n represents the meteor event in the height interval I a The total number of meteor events within , m∈[1,n]; Let the interval set corresponding to the meteor event MetS a The data matrix set is AS a , specifically: AS a =[A a,1 A a,2 ... A a,m ... A a,n ] Among them, A a,m Indicates that it is in altitude interval I a The mth meteor event Met a,m The corresponding data matrix, n represents the height interval I a The total number of meteor events within , m∈[1,n]; Let the interval set corresponding to the meteor event MetS a The drift matrix set is VS a , specifically: VS a =[V a,1 In a，2 ... In a,m ... In a,n ] Among them, V a,m Indicates that it is in altitude interval I a The drift matrix corresponding to the mth meteor event within the range of a The total number of meteor events within , m∈[1,n]; Suppose that the meteor event interval set MetS a The meteor events in the AS are connected vertically in a one-to-one correspondence manner. a The obtained height interval I a The joint data matrix within is ComA a , specifically: Suppose that the meteor event interval set MetS a The meteor events in the vertical series VS a The obtained height interval I a The joint drift matrix within is ComV a , specifically: All M-1 height intervals in the height interval sequence I are processed according to the above process to obtain all joint data matrices and joint drift matrices corresponding to each height interval. The sequences of the two are represented by ComA and ComV, respectively, as follows: ComA=[ComA1 ComA2 ... ComA a ... ComA M-1 ] ComV=[ComV1 ComV2 ... ComV a ... ComV M-1 ]。 4. The method for inverting the height distribution of the gravity wave disturbance parameters at the top of the intermediate layer according to claim 3 is characterized by: In step 4, the least squares fitting method is used to obtain the gravity wave disturbance parameter matrix in the corresponding height interval, which is specifically: Let correspond to the height interval I a The gravity wave disturbance parameter matrix is ​​GW a , specifically: GW a =(With a T With a ) -1 With a T ComV a All M-1 height intervals in the height interval sequence I are processed according to the above process to obtain the parameter matrix of all gravity wave disturbances corresponding to each height interval. The sequence is represented by GW, specifically: GW=「GW1 GW2 ... GW a ... GW M-1 ]。 5. The method for inverting the height distribution of the gravity wave disturbance parameters at the top of the intermediate layer according to claim 4 is characterized by: The step 5 is specifically as follows: Suppose the discrete height scale sequence corresponding to the height interval sequence I is H', specifically: H′=[H′1 H′2 ... H′ C ... H′ M-1 ] Among them, H' a Indicates the height interval I corresponding to the ath height interval I in the height interval sequence I a The discrete height scale is: Assume that the distribution matrix of gravity wave disturbance parameters at different heights of the mesosphere top is GW_MAT, then Among them, H' a Indicates the height interval I corresponding to the ath height interval I in the height interval sequence I a Discrete height scale, GW a corresponds to the height interval I a The gravity wave perturbation parameter matrix, a∈[1,M-1].

6. A system for inverting the height distribution of gravity wave disturbance parameters at the top of the intermediate layer, characterized in that: include: The data matrix and drift matrix calculation module of the meteor event is used to obtain multiple meteor events to construct a meteor event sequence, and construct the data matrix and drift matrix corresponding to each meteor event, as well as the sequence of the height of the meteor event according to the zenith angle, azimuth angle, radial drift velocity disturbance of the trail and the height of each meteor event; A meteor event height interval division module is used to divide the mesosphere top area into multiple height intervals at equal intervals, and place each meteor event in a corresponding height interval according to the height of each meteor event; The joint data matrix and joint drift construction calculation module of the meteor event is used to vertically connect the data matrices and drift matrices corresponding to all meteor events in the same height interval in a one-to-one correspondence manner in each height interval to construct the joint data matrix and joint drift matrix in the height interval; A gravity wave disturbance parameter matrix calculation module is used to obtain the gravity wave disturbance parameter matrix in the corresponding height interval by the least square fitting method according to the corresponding joint data matrix and joint drift matrix; A distribution matrix construction module is used to connect all gravity wave disturbance parameter matrices in series in order of height to obtain the distribution matrix of gravity wave disturbance parameter at different heights of the mesosphere top; The distribution matrix application module is used for the distribution matrix of gravity wave disturbance parameters at different heights of the mesosphere top, which is used to quantitatively analyze the changes of gravity wave momentum flux at different heights, and through the subsequent further calculation of the divergence of gravity wave momentum flux, the dynamic effect of gravity waves on the background atmosphere at the mesosphere top, including the drag caused by gravity waves on the wind field, is solved; The meteor event sequence is specifically: With=[With1 With2 ... With x ... Of N ] Among them, Met represents the meteor event sequence, Met x represents the xth meteor event in the meteor event sequence, N represents the total number of meteor events in the meteor event sequence, x∈[1,N]; The data matrix corresponding to each meteor event is constructed as follows: The zenith angle of the xth meteor event is θ x , the azimuth of the x-th meteor event is φ x , the radial drift velocity disturbance of the wake of the xth meteor event is Vr' x ; The data matrix A of the x-th meteor event x for: Among them, A x is the data matrix of the x-th meteor event, x∈[1,N]; Step 1 constructs the drift matrix corresponding to each meteor event as follows: Among them, V x is the drift matrix of the x-th meteor event, x∈[1,N]; Let the sequence of meteor event heights be H t , specifically: Ht=[Ht1 Ht2 ... Ht x ... Ht N ] Among them, Ht x Indicates the xth meteor event Met x The height of the corresponding meteor event, N represents the total number of meteor events in the meteor event sequence, x∈[1,N].

Images