COMPUTERIZED IONOSPHERIC TOMOGRAPHY METHOD BASED ON VERTICAL BOUNDARY TRUNCATION RAYS

Abstract

The disclosed technology includes a computerized ionospheric tomography method based on vertical boundary truncation rays, which relates to the technical field of computerized ionospheric tomography (CIT). The method includes: obtaining an initial ionospheric electron density (IED) of each voxel in a target region and an ionospheric total electron content (TEC) value along a propagation path from a global navigation satellite system (GNSS) satellite; extending the target region so that GNSS stations within a certain range beyond the target region are encompassed within the target region; for GNSS stations within a certain range in the target region, calculating a vertical boundary truncation TEC value; for the GNSS stations within the target region, calculating a vertical boundary truncation TEC value; and building a three-dimensional CIT model based on the vertical boundary truncation TEC values P.sub.rTEC and P.sub.sTEC.

Claims

1. A computerized ionospheric tomography method based on vertical boundary truncation rays, comprising: obtaining an initial ionospheric electron density (IED) of each voxel in a target region and an ionospheric total electron content (TEC) value along a propagation path from global navigation satellite system (GNSS) stations; extending the target region so that GNSS stations within a certain range beyond the target region are encompassed within the target region; calculating a vertical boundary truncation TEC value, denoted by P.sub.rTEC for GNSS stations within a predetermined range in the target region; calculating a vertical boundary truncation TEC value, denoted by P.sub.sTEC for GNSS stations within the target region; building a three-dimensional CIT model based on the vertical boundary truncation TEC values denotated by P.sub.rTEC and P.sub.sTEC; and resolving the CIT model for inversion of the IED of the target region.

2. The computerized ionospheric tomography method based on vertical boundary truncation rays according to claim 1, wherein obtaining the initial IED of each voxel in the target region and the ionospheric TEC value along the propagation path from GNSS stations comprises determining a longitude, a latitude, an altitude and time of the target region.

3. The computerized ionospheric tomography method based on vertical boundary truncation rays according to claim 1, wherein calculating the vertical boundary truncation TEC value denoted by P.sub.rTEC for GNSS stations within the predetermined range in the target region comprises: obtaining an initial IED value denotated by N.sub.e.sup.0(r,t) by a NeQuick model and a TEC value denotated by TEC.sub.0 with the NeQuick model as prior information; calculating a truncation factor λ.sub.r with the initial IED value denotated by N.sub.e.sup.0(r,t) and the TEC value denotated by TEC.sub.0 for GNSS stations within a predetermined range beyond the target region; and calculating the vertical boundary truncation TEC value denotated by P.sub.rTEC with the truncation factor λ.sub.r and an ionospheric TEC value denotated by TEC.sub.G along the propagation path from GNSS stations.

4. The computerized ionospheric tomography method based on vertical boundary truncation rays according to claim 1, wherein calculating the vertical boundary truncation TEC value denoted by P.sub.rTEC for GNSS stations within the predetermined range in the target region comprises: obtaining an initial IED value denotated by N.sub.e.sup.0(r,t) by a NeQuick model and a TEC value denotated by TEC.sub.0 with the NeQuick model as prior information; calculating a truncation factor λ.sub.r with the initial IED value denotated by N.sub.e.sup.0(r,t) and the TEC value denotated by TEC.sub.0 for GNSS stations within a predetermined range beyond the target range according to a first formula:
P.sub.rTEC.sub.0=∫.sub.lN.sub.e.sup.0(r,t)ds; and a second formula: $\begin{array}{cc}{\lambda }_r=\frac{P_r{\mathrm{TEC}}_0}{{\mathrm{TEC}}_0},& \left[\left[\left(2\right)\right]\right]\end{array}$ wherein TEC.sub.0 is the TEC value obtained with the NeQuick model as prior information; and calculating the vertical boundary truncation TEC value denotated by P.sub.rTEC with the truncation factor λ.sub.r and an ionospheric TEC value denotated by TEC.sub.G along the propagation path from GNSS stations according to a third formula:
P.sub.rTEC=λ.sub.r.Math.TEC.sub.G wherein TEC.sub.G is the ionospheric TEC value extracted from GNSS observation data.

5. The computerized ionospheric tomography method based on vertical boundary truncation rays according to claim 1, wherein calculating the vertical boundary truncation TEC value denoted by P.sub.sTEC for the GNSS stations within the target region comprises: obtaining an initial IED value denotated by N.sub.e.sup.0(r,t) by a NeQuick model and a TEC value denotated by TEC.sub.0 with the NeQuick model as prior information; calculating a truncation factor λ.sub.s with the initial IED value denotated by N.sub.e.sup.0(r,t) and the TEC value denotated by TEC.sub.0 for the GNSS stations within the target region; and calculating the vertical boundary truncation TEC value denotated by P.sub.sTEC with the truncation factor λ.sub.s and an ionospheric TEC value TEC.sub.G along the propagation path from GNSS stations.

6. The computerized ionospheric tomography method based on vertical boundary truncation rays according to claim 1, wherein calculating the vertical boundary truncation TEC value denoted by P.sub.sTEC for the GNSS stations within the target region comprises: obtaining an initial IED value denotated by N.sub.e.sup.0(r,t) by a NeQuick model and a TEC value denotated by TEC.sub.0 with the NeQuick model as prior information; calculating a truncation factor λ.sub.s with the initial IED value denotated by N.sub.e.sup.0(r,t) and the TEC value denotated by TEC.sub.0 for GNSS stations within an extended section of the target region according to a first formula:
P.sub.sTEC.sub.0=∫.sub.lN.sub.e.sup.0(r,t)ds; and a second formula: $\begin{array}{cc}{\lambda }_s=\frac{P_s{\mathrm{TEC}}_0}{{\mathrm{TEC}}_0}& \left[\left[\left(5\right)\right]\right]\end{array}$ wherein TEC.sub.0 is the TEC value obtained with the NeQuick model as prior information; and calculating the vertical boundary truncation TEC value denotated by P.sub.sTEC with the truncation factor λ.sub.s and an ionospheric TEC value TEC.sub.G along the propagation path from GNSS stations according to third formula:
P.sub.sTEC=λ.sub.s.Math.TEC.sub.G wherein TEC.sub.G is the ionospheric TEC value extracted from GNSS observation data.

7. The computerized ionospheric tomography method based on vertical boundary truncation rays according to claim 1, wherein the inversion of the IED of the target region in the resolving the CIT model for inversion of the IED of the target region is carried out according to a first formula:
y.sub.m×1=A.sub.m×n.Math.x.sub.n×1+e.sub.m×1 wherein y is an m-dimensional column vector constituted by observed ionospheric TEC values, A.sub.m×n is dimensional row vectors constituted by intercepts of rays in respective grids, x is an n-dimensional column vector constituted by unknown parameters, and e is an m-dimensional column vector constituted by observation noise and a discretization error.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0038] FIG. 1 is a flowchart of a computerized ionospheric tomography method based on vertical boundary truncation rays according to an embodiment of the present disclosure.

[0039] FIG. 2 is a schematic distribution diagram of GNSS rays based on vertical boundaries according to an embodiment of the present disclosure.

[0040] FIG. 3 is a schematic diagram of calculations of vertical boundary truncation TECs of a target region according to an embodiment of the present disclosure.

[0041] FIG. 4 is a distribution map of GNSS stations in an experimental region according to an embodiment of the present disclosure.

[0042] FIG. 5 shows comparison of average numbers of available rays per day between the TRM method and a traditional method according to an embodiment of the present disclosure.

[0043] FIG. 6 shows comparison of IED profiles between the TRM method and a traditional method, and an ionosonde station at 9:00 UT according to an embodiment of the present disclosure.

[0044] FIG. 7 shows comparison of IED profiles between the TRM method and a traditional method, and an ionosonde station at 19:00 UT according to an embodiment of the present disclosure.

DETAILED DESCRIPTION

[0045] To make the objective, technical solution and advantages of the present disclosure clearer and more comprehensible, the present disclosure will be further described below in detail by way of example. It will be understood that the specific embodiment described herein is merely used to explain, rather than limit, the present disclosure.

[0046] In the embodiment of the present disclosure, the ionospheric TEC value refers to ionospheric total electron content, and GNSS denotes global navigation satellite system.

[0047] Referring to FIG. 1, a computerized ionospheric tomography method based on vertical boundary truncation rays includes:

[0048] S1: obtain an initial IED of each voxel in a target region and an ionospheric TEC value along a propagation path from a GNSS satellite.

[0049] The target region may be selected based on data of a longitude ranging from 99° E to 115° E, a latitude ranging from 26° N to 34° N, and an altitude ranging from 100 km to 1000 km for inversion. Grid intervals in the directions of longitude and latitude are set to 1°, while an interval in the direction of altitude is set to 50 km, as shown in FIG. 4.

[0050] The initial IED of each voxel in the target region is obtained by an ionospheric empirical model NeQuick.

[0051] S2: extend the target region in step S1 so that GNSS stations within a certain range beyond the target region are encompassed within the target region.

[0052] S3: for GNSS stations within a certain range in the target region, calculate a vertical boundary truncation TEC value, denoted by P.sub.rTEC. The range denoted by P.sub.rTEC is as shown in FIG. 2 and FIG. 3.

[0053] Specific building process of P.sub.rTEC model may include:

[0054] S311: obtain an initial IED value N.sub.e.sup.0(r,t) by the NeQuick model and a TEC value TEC.sub.0 with the NeQuick model as prior information.

[0055] S321: for GNSS stations within a certain range beyond the target region, calculate a truncation factor λ.sub.r with the initial IED value N.sub.e.sup.0(r,t) and TEC.sub.0 according to formulas (1) and (2):

[00003]$\begin{array}{cc}{P}_r{\mathrm{TEC}}_0={\int }_l{N}_e^0\left(r,t\right)\mathrm{ds}& \left(1\right)\\ {\lambda }_r=\frac{P_r{\mathrm{TEC}}_0}{{\mathrm{TEC}}_0},& \left(2\right)\end{array}$

[0056] where TEC.sub.0 is the TEC value obtained with the NeQuick model as prior information; and

[0057] S331: calculate the value of P.sub.rTEC with the truncation factor λ.sub.r and the ionospheric TEC value TEC.sub.G along the propagation path from the GNSS satellite according to formula (3):

P.sub.rTEC=λ.sub.r.Math.TEC.sub.G  (3)

[0058] where TEC.sub.G is the ionospheric TEC value extracted from GNSS observation data.

[0059] S4: for the GNSS stations within the target region, calculate a vertical boundary truncation TEC value, denoted by P.sub.sTEC. The range denoted by P.sub.sTEC is as shown in FIG. 2 and FIG. 3.

[0060] Specific building process of P.sub.sTEC model may include:

[0061] S411: obtain an initial IED value N.sub.e.sup.0(r,t) by the NeQuick model and a TEC value TEC.sub.0 with the NeQuick model as prior information;

[0062] S421: for the GNSS stations within the target region, calculate a truncation factor λ.sub.s with the initial IED value N.sub.e.sup.0(r,t) and TEC.sub.0; and

[0063] S431: calculate the value of P.sub.sTEC with the truncation factor λ.sub.s and the ionospheric TEC value TEC.sub.G along the propagation path from the GNSS satellite.

[0064] S5: build a three-dimensional CIT model based on the vertical boundary truncation TEC values P.sub.rTEC and P.sub.sTEC.

[0065] S6: resolve the CIT model in step S5 for inversion of the IED of the target region.

[0066] The inversion of the IED of the target region in step S6 may be carried out according to formula (7):

y.sub.m×1=A.sub.m×n.Math.x.sub.n×1+e.sub.m×1  (7)

[0067] where y is an m-dimensional column vector constituted by observed ionospheric TEC values, while A is an m n-dimensional row vector constituted by intercepts of rays in respective grids, x is an n-dimensional column vector constituted by unknown parameters, and e is an m-dimensional column vector constituted by observation noise and a discretization error.

[0068] FIG. 5 shows comparison of average numbers of available rays per day between the method in the present disclosure and a traditional method, where the 152th day to the 181th day of 2015, 30 days in total, were selected for experimental analysis. As can be seen from the figure, the average number of available rays per day based on the TRM method is obviously greater than that based on the traditional method. This indicates that the TRM based on vertical boundaries in the present disclosure may greatly increase the number of available rays, so that information for modeling is increased, thus being conducive to improving the CIT inversion accuracy.

[0069] FIG. 6 and FIG. 7 show comparisons of IED profiles measured by an ionosonde at different moments with inversion results based on the TRM method and a traditional method.

[0070] Besides, FIG. 6 and FIG. 7 show the results of comparisons of electron density profiles at 9:00 UT and 19:00 UT on Jun. 19, 2015, respectively. From the results, the F2 peak IED of the inversion results based on the TRM method is closer to the data of the ionosonde. Experimental results indicate that the TRM based on vertical boundaries may effectively improve the CIT inversion accuracy.

[0071] The foregoing are merely descriptions of the preferred embodiment of the present disclosure, which are not intended to limit the present disclosure. Any modifications, equivalent replacements and improvements made within the spirit and principle of the present disclosure shall fall within the protection scope of the present disclosure.