Prediction of Target Detection Probability Based on Air-to-Air Long-Range Scenarios in Anomalous Atmospheric Environments

: We investigate a target detection probability (TDP) using path loss of an airborne radar based on air-to-air scenarios in anomalous atmospheric and weather environments. In the process of calculating the TDP, it is necessary to obtain the overall path loss including the anomalous atmospheric environment, gas attenuation, rainfall attenuation, and beam scanning loss. The path loss including the quad-linear refractivity model and other radar input parameters is simulated using the Advanced Refractive Effects Prediction System (AREPS) software along the range and the altitude. For the gas and rainfall attenuations, ITU-R models are used to consider the weather environment. In addition, the radar beam scan loss and a radar cross section (RCS) of the target are considered to estimate the TDP of the airborne long-range radar. The TDP performance is examined by employing the threshold evaluations of the total path loss derived from the detectability factor and the free-space radar range equation. Finally, the TDPs are obtained by assuming various air-to-air scenarios for the airborne radar in anomalous atmospheric and weather environments.


Introduction
With the dramatic advances of radar system design technologies, the use of long-range airborne radar systems, such as synthetic aperture radars (SARs), airborne early warning (AEW) radars, and active electronically scanned array (AESA) radars, has been growing extensively [1][2][3][4][5][6]. Such airborne radar systems essentially require high-performance specifications for each component to increase the detection probability of long-range targets. However, the TDP can often be degraded by environmental or external factors including noises, clutter, atmospheric gas attenuation, multipath interference, atmospheric refraction, and rainfall attenuation. To be specific, an anomalous atmospheric refractive index including an unusual or abnormal distributions of temperature and relative humidity along the altitude can cause the refraction of the electromagnetic (EM) wave propagation different from that in normal conditions, i.e., sub-refraction, super-refraction, and ducting [7][8][9]. These wave refractions can mislead the accuracy of target position predictions, and in the worst case, long-range targets can be missed under anomalous atmospheric conditions. Therefore, it is important to model the detailed atmospheric refractive index in order for the precise predictions of the propagation path loss and propagation factor. Many studies have been conducted to model the refractive index under anomalous atmospheric conditions through the radar signal measurement [10][11][12], global positioning system tropospheric delay observation [13], and statistical analysis of the stored meteorological observatory data [14,15]. In addition, various studies have been conducted to analyze the wave propagation characteristics in consideration of low-altitude actual atmospheric data in ground-to-ground and ground-to-air scenarios. For example, the real atmospheric data of specific coastal areas in the UK [16], United States [17], Greece [18], China [19], and Korea [20] are used to calculate the long-range path loss using the propagation models, i.e., the parabolic equation model, ray optics model, waveguide mode model, and hybrid model [21][22][23][24][25]. These aforementioned studies have demonstrated the high accuracy of the path loss estimation in low-altitude situations. However, their path loss calculation processes need more in-depth consideration of the weather or atmospheric environment to observe the air-to-air airborne radar propagation characteristics at a high altitude over 5 km.
In this paper, we propose a TDP calculation process for the airborne long-range radar based on air-to-air scenarios in anomalous atmospheric environments. In the proposed process, it is necessary to obtain the overall path loss including the anomalous atmospheric environment, gas attenuation, rainfall attenuation, and beam scanning loss. To observe the airborne radar path loss in an anomalous atmospheric environment, the refractive index as a function of altitude is modeled using four linear lines to represent various wave refractions. Then, the path loss, including the refractivity model along the range and the altitude, is simulated using the AREPS software [26], which is a commercial software based on a hybrid model with a parabolic equation model and raytracing model. In this AREPS simulation, the radar antenna beam pattern and the digital terrain elevation data (DTED) are also employed as input parameters to accurately predict the airborne radar propagation. The 32 × 32 airborne radar array antenna with a triangular array configuration is used to calculate the radiation pattern, and the DTED of the southwest region in South Korea from the National Geographic Information Institute is employed. For gas and rainfall attenuations, ITU-R models are used to consider the weather environment in accordance with the water vapor pressure and precipitation. In addition, the radar beam scan loss and RCS of targets are considered to estimate the TDP of the airborne long-range radar. Finally, the TDPs are examined by assuming various air-to-air scenarios in anomalous atmospheric and weather environments. Figure 1 shows a conceptual figure of the air-to-air airborne radar wave propagation in accordance with an anomalous atmospheric environment. The gray and red aircrafts indicate an airborne radar and a target, which are located at heights of h r and h t , having a distance of R t (slant range) and a horizontal range of R h . The target is placed under various atmospheric conditions according to the wave refraction. To observe the wave refraction of the airborne radar, the refractive index needs be calculated using the measured atmospheric data such as temperature, air pressure, and relative humidity based on Equations (1)-(3), which are usually obtained by employing a rawinsonde or a GPS sonde [27][28][29].

Anomalous Atmospheric Refractivity and Measurement
where p, T, rh indicate air pressure in hPa, temperature in K, and relative humidity in %. Then, the modified refractive index M at a height of h needs to be calculated using Equation (4) [30]. An example for a calculation of the modified refractive index M based on the measured data from Heuksando Meteorological Observatory [31] is listed in Table 1.   The gradient of the modified refractive index along the height M can determine four common wave propagation refractions: normal refraction (79 < M  157), superrefraction (0 < M  79), sub-refraction (157 < M), and ducting (M < 0). In accordance with these four refractions, the airborne radar wave propagation can become bent in the direction opposite or toward Earth's surface, and it can be trapped in the ducting region. In particular, elevated ducts sometimes occur at high altitudes, which affect the wave propagation with a high path loss near the elevated ducting region. Thus, it is necessary to precisely model the refractivity of the elevated ducting in the air-to-air situation of the anomalous environment using a quad-linear model [32]. Figure 2a shows the conceptual figure of the quad-linear refractivity model. Four linear lines with slopes of m1, m2, m3, and m4 are used to model the refractivity, and the heights where the slope changes are h1, h2, and h3. Figure 2b represents the modified refractivity according to height using the quadlinear refractivity model in comparison with the actual refractivity data from a meteorological observatory. The measured refractivity data is interpolated using the typical linear interpolation method due to the sparse data along the height. The solid and dashed lines indicate the comparisons between the actual data from Heuksando Meteorological Observatory [31] on 25 July 2017 and the quad-linear model.  The gradient of the modified refractive index along the height ∇M can determine four common wave propagation refractions: normal refraction (79 < ∇M ≤ 157), superrefraction (0 < ∇M ≤ 79), sub-refraction (157 < ∇M), and ducting (∇M < 0). In accordance with these four refractions, the airborne radar wave propagation can become bent in the direction opposite or toward Earth's surface, and it can be trapped in the ducting region. In particular, elevated ducts sometimes occur at high altitudes, which affect the wave propagation with a high path loss near the elevated ducting region. Thus, it is necessary to precisely model the refractivity of the elevated ducting in the air-to-air situation of the anomalous environment using a quad-linear model [32]. Figure 2a shows the conceptual figure of the quad-linear refractivity model. Four linear lines with slopes of m 1 , m 2 , m 3 , and m 4 are used to model the refractivity, and the heights where the slope changes are h 1 , h 2 , and h 3 . Figure 2b represents the modified refractivity according to height using the quad-linear refractivity model in comparison with the actual refractivity data from a meteorological observatory. The measured refractivity data is interpolated using the typical linear interpolation method due to the sparse data along the height. The solid and dashed lines indicate the comparisons between the actual data from Heuksando Meteorological Observatory calculating the array factor of a 32 × 32 airborne radar array antenna using a triangular array configuration [31]. The array distances along the xand y-axes are 0.475λ and 0.538λ at 10 GHz. This radiation pattern has a half-power beamwidth of 2.85 • and a side lobe level of 13.5 dB. The DTED of the southwest region in South Korea are obtained from the Korean National Geographic Information Institute [33]. Then, the path loss computed using the AREPS software can be defined as L AREPS (m 1 , m 2 , m 3 , m 4 , h 1 , h 2 , h 3 , R h , h r , h t ) considering the anomalous environment. For example, with specific input parameters (m 1 = m 2 = m 3 = m 4 = 85, h 1 = h 2 = h 3 = 15 km, R h = 150 km, h r = 11 km, and h t = 5 km), we can obtain the path loss L AREPS of 163.5 dB in the elevated ducting of the anomalous atmospheric environment. . These quad-linear refractivity models have 2001 data points along the altitude, which are employed as input parameters for the AREPS software to obtain the path loss according to the range and altitude. Furthermore, the other input parameters of the antenna beam pattern and DTED are also applied to the AREPS simulation to precisely estimate the air-toair wave propagation path loss. The antenna beam pattern is derived by calculating the array factor of a 32 × 32 airborne radar array antenna using a triangular array configuration [31]. The array distances along the x-and y-axes are 0.475 and 0.538at 10 GHz. This radiation pattern has a half-power beamwidth of 2.85° and a side lobe level of 13.5 dB. The DTED of the southwest region in South Korea are obtained from the Korean National Geographic Information Institute [33]. Then, the path loss computed using the AREPS software can be defined as LAREPS (m1, m2, m3, m4, h1, h2, h3, Rh, hr, ht) considering the anomalous environment. For example, with specific input parameters (m1 = m2 = m3 = m4 = 85, h1 = h2 = h3 = 15 km, Rh = 150 km, hr = 11 km, and ht = 5 km), we can obtain the path loss LAREPS of 163.5 dB in the elevated ducting of the anomalous atmospheric environment. Figure 3 illustrates the attenuation of the wave propagation by the atmospheric gases and rainfall. Such weather conditions can have a serious impact on the path loss, which is directly related to the TDP performance of the airborne radar. Hence, it is necessary to have the detailed attenuation models to increase the accuracy of the TDP estimation in the air-to-air airborne radar wave propagation. The atmospheric gas attenuation can be modeled based on ITU-R P.676-12 as follows [34]:  Figure 3 illustrates the attenuation of the wave propagation by the atmospheric gases and rainfall. Such weather conditions can have a serious impact on the path loss, which is directly related to the TDP performance of the airborne radar. Hence, it is necessary to have the detailed attenuation models to increase the accuracy of the TDP estimation in the air-to-air airborne radar wave propagation. The atmospheric gas attenuation can be modeled based on ITU-R P.676-12 as follows [34]:

Weather Environment Models
where γ g is the specific gas attenuation in dB/km, and f is the frequency in GHz. N Ox and N Wat are the imaginary parts of complex refractivities in terms of oxygen and water vapor. S i,Ox and S i,Wat are the strengths of the ith oxygen and water vapor lines, and N D (f ) is the dry continuum. P and e are dry air and water vapor partial pressures in hPa, and T is the absolute temperature in K. F i is the oxygen or water vapor line shape. a 1 , a 2 , b 1 , and b 2 are the line strength coefficients. The detailed coefficient values are listed in tables in [34]. Thus, the gas attenuation according to the range R t can be defined as L g in Equation (7).
where γg is the specific gas attenuation in dB/km, and f is the frequency in GHz. N ″ Ox and N ″ Wat are the imaginary parts of complex refractivities in terms of oxygen and water vapor. Si,Ox and Si,Wat are the strengths of the ith oxygen and water vapor lines, and N ″ D(f) is the dry continuum. P and e are dry air and water vapor partial pressures in hPa, and T is the absolute temperature in K. Fi is the oxygen or water vapor line shape. a1, a2, b1, and b2 are the line strength coefficients. The detailed coefficient values are listed in tables in [34]. Thus, the gas attenuation according to the range Rt can be defined as Lg in equation (7). In addition, the rainfall attenuation can be modeled based on ITU-R P.530-17, P.837-7, and P.838-3 as follows [35][36][37]: where γr is the specific rainfall attenuation in dB/km, and  is the rainfall precipitation in mm/hr. α and κ are the rainfall attenuation parameters that depend on the frequency, polarization state, and angle of the signal path to the target. deff is the effective propagation distance considering the target range Rt multiplied by the scale factor of r. The rainfall attenuation in terms of the range can be defined as Lr in Equation (11). In general, the precipitation is measured using rain gauge, weather radar, and satellite, and the measured weather information can be employed to calculate the rainfall attenuation [38][39][40][41][42][43]. To verify these attenuation models, we have investigated the attenuation measurements [41] and then some of the measurements according to the frequency with different precipitations are compared to the ITU-R attenuation models, as shown in Figure 4a. Furthermore, Figure 4b depicts the attenuation results compared with those of the measurement in accordance with the precipitation [42,43]. These comparison results demonstrate that the ITU-R model agrees well with the measurement, and this model can be adopted to increase the accuracy of the TDP in real atmospheric and weather environments. To predict the TDP along the range, we calculate atmospheric gas and rainfall attenuation results according to the range from 0 to 190 km, as presented in Figure 4c. The solid and dashed lines indicate the atmospheric gas attenuations in dry and wet weathers, In addition, the rainfall attenuation can be modeled based on ITU-R P.530-17, P.837-7, and P.838-3 as follows [35][36][37]: where γ r is the specific rainfall attenuation in dB/km, and ρ is the rainfall precipitation in mm/h. α and κ are the rainfall attenuation parameters that depend on the frequency, polarization state, and angle of the signal path to the target. d eff is the effective propagation distance considering the target range R t multiplied by the scale factor of r. The rainfall attenuation in terms of the range can be defined as L r in Equation (11). In general, the precipitation is measured using rain gauge, weather radar, and satellite, and the measured weather information can be employed to calculate the rainfall attenuation [38][39][40][41][42][43]. To verify these attenuation models, we have investigated the attenuation measurements [41] and then some of the measurements according to the frequency with different precipitations are compared to the ITU-R attenuation models, as shown in Figure 4a. Furthermore, Figure 4b depicts the attenuation results compared with those of the measurement in accordance with the precipitation [42,43]. These comparison results demonstrate that the ITU-R model agrees well with the measurement, and this model can be adopted to increase the accuracy of the TDP in real atmospheric and weather environments. To predict the TDP along the range, we calculate atmospheric gas and rainfall attenuation results according to the range from 0 to 190 km, as presented in Figure 4c  In the atmospheric environment modeling, it is important to consider uncertainties of the measurement for the atmospheric refractivity [44] and rainfall attenuation, which can possibly affect the TDP results. For example, if we make an assumption that the measurement of the refractivity had the uniform uncertainty of 2.5 % along the height, then the deviations between the original measurement and the uncertainties of the TDP can be increased up to 4.7% in the long-range over 80 km.

TDP Calculation with Airborne Radar Parameters
In Section 2.1, the refractive index, gas attenuation, rainfall attenuation, and radar parameters were modeled to estimate the total path loss of the wave propagation in an anomalous atmospheric environment; however, more radar parameters and computations are needed to accurately predict the TDP of the airborne radar. We determine the TDP using the threshold evaluations derived from the free-space path loss according to the range. To examine the threshold for the total path loss, the detectability factor is required for a given probability of detection and false alarm as written in Equation (14). In addition, the fluctuation loss is modeled using Equation (15) to calculate the TDP performance for describing the fluctuating target. Finally, these equations are directly related to the radar range equation for the detectable range that is introduced by Blake as follows [45]: where P d is the detection probability ranging from 0 to 0.99, and P fa indicates the false alarm rate of the airborne radar. L f is the fluctuation loss in Swerling Case 1, and N p is the number of pulses integrated by the detector depending on the radar system. D is the detectability factor, and R fs is the free-space detectable range in meters. c 0 is the light speed of 3 × 10 8 m/s, and σ is the RCS of the target in m 2 . k is the Boltzmann constant, and τ is the pulse length of the radar in seconds. P t and G are transmitting power in watts and the array antenna gain, respectively. N f and L s are the noise figure and miscellaneous system loss of the radar. In the radar range Equation (16), it is required to set the system parameters of the airborne radar to specific values: P t = 1 MW, τ = 2 µs, P fa = 10 −8 , G = 10 3 , k = 1.38 × 10 −23 , f = 10 GHz, N f = 10 0.5 , and L s = 10 0.3 . Then, the total path loss of the threshold related to the TDP can be calculated by using the detectable range as follows: Figure 5 shows the calculated path loss threshold and detectability factor results according to the detection probability at a fixed RCS value of 8.5 m 2 . The path loss threshold varies from 151.4 to 159.2 dB when the P d increases from 70% to 100%. The maximum detectability factor is 17.9. The total path loss threshold is significantly affected by the detectable range R fs , and the detectability factor is dominantly decided by the fluctuation loss for the fluctuating type of the target. To further observe the TDP performance with the ranges and angles, we additionally compute the theoretical scan losses in terms of the scan angle φ scan for the beam scanning of the radar array antenna, as written in Equation (18) [46,47]. where L scan is the scan loss of the radar, and n is the power of the cosine function set to be −2.5 in this research. Figure 6a presents the scan loss according to the scan angle. Due to the beam scanning, the array antenna gain can be gradually attenuated until 6.8 dB at the scanning angle of 70 • . We also model a target aircraft using FEKO full EM software to obtain the RCSs at the azimuth angles. Figure 6b presents the RCS for the rear direction of the target aircraft. In the scanning angle, the maximum and averaged RCS values are 67.2 and 2.8 dBsm at the scan angles of 0 • and 70 • . With the environment models and the radar parameters, we can predict the total wave propagation loss for the airborne radar in the air-to-air situation. This total loss L tot is then calculated to examine the TDP using Equation (19). For example, in the normal atmospheric environment, if the total path loss L tot is 147.5 dB at R t of 50 km, then the TDP can be obtained over 60% regardless of the target RCS. On the other hand, the TDP for very small targets with the RCS of lower than 2 m 2 (typical military airborne) dramatically decrease for the long-range target over 100 km (L tot > 154.2 dB). L tot = L AREPS + L g + L r + L scan [dB] (22) loss for the fluctuating type of the target. To further observe the TDP performance with the ranges and angles, we additionally compute the theoretical scan losses in terms of the scan angle scan for the beam scanning of the radar array antenna, as written in Equation (18) [46,47]. where Lscan is the scan loss of the radar, and n is the power of the cosine function set to be −2.5 in this research. Figure 6a presents the scan loss according to the scan angle. Due to the beam scanning, the array antenna gain can be gradually attenuated until 6.8 dB at the scanning angle of 70°. We also model a target aircraft using FEKO full EM software to obtain the RCSs at the azimuth angles. Figure 6b presents the RCS for the rear direction of the target aircraft. In the scanning angle, the maximum and averaged RCS values are 67.2 and 2.8 dBsm at the scan angles of 0° and 70°. With the environment models and the radar parameters, we can predict the total wave propagation loss for the airborne radar in the air-to-air situation. This total loss Ltot is then calculated to examine the TDP using Equation (19). For example, in the normal atmospheric environment, if the total path loss Ltot is 147.5 dB at Rt of 50 km, then the TDP can be obtained over 60% regardless of the target RCS. On the other hand, the TDP for very small targets with the RCS of lower than 2 m 2 (typical military airborne) dramatically decrease for the long-range target over 100 km (Ltot > 154.2 dB).  loss for the fluctuating type of the target. To further observe the TDP performance with the ranges and angles, we additionally compute the theoretical scan losses in terms of the scan angle scan for the beam scanning of the radar array antenna, as written in Equation (18) [46,47]. where Lscan is the scan loss of the radar, and n is the power of the cosine function set to be −2.5 in this research. Figure 6a presents the scan loss according to the scan angle. Due to the beam scanning, the array antenna gain can be gradually attenuated until 6.8 dB at the scanning angle of 70°. We also model a target aircraft using FEKO full EM software to obtain the RCSs at the azimuth angles. Figure 6b presents the RCS for the rear direction of the target aircraft. In the scanning angle, the maximum and averaged RCS values are 67.2 and 2.8 dBsm at the scan angles of 0° and 70°. With the environment models and the radar parameters, we can predict the total wave propagation loss for the airborne radar in the air-to-air situation. This total loss Ltot is then calculated to examine the TDP using Equation (19). For example, in the normal atmospheric environment, if the total path loss Ltot is 147.5 dB at Rt of 50 km, then the TDP can be obtained over 60% regardless of the target RCS. On the other hand, the TDP for very small targets with the RCS of lower than 2 m 2 (typical military airborne) dramatically decrease for the long-range target over 100 km (Ltot > 154.2 dB).
[]   Figure 7 shows three air-to-air scenarios considered for the airborne radar TDP simulation in anomalous atmospheric environments. The detailed scenario explanations are as follows:

TDP Calculation with Airborne Radar Parameters
(1) TDPs when azimuth beam scanning within 90 • in an anomalous atmospheric environment of the refractivity, as shown in Figure 7a.
(2) TDPs when azimuth beam scanning within 90 • in a rainy weather environment, as shown in Figure 7b. (3) TDPs according to the distance (R t < 190 km) to the target in a heavy rain environment, as shown in Figure 7c. Figure 7 shows three air-to-air scenarios considered for the airborne radar TDP simulation in anomalous atmospheric environments. The detailed scenario explanations are as follows:

TDP Calculation with Airborne Radar Parameters
(1) TDPs when azimuth beam scanning within 90° in an anomalous atmospheric environment of the refractivity, as shown in Figure 7a. (2) TDPs when azimuth beam scanning within 90° in a rainy weather environment, as shown in Figure 7b. In all scenarios, the height hr of the airborne radar is consistent at 11 km, and the target height ht is set to be 5 km. Additionally, these scenarios are in the long-range radar propagation situations, which assumes that the slant range is almost same with the horizontal range (Rt  Rh). The detailed explanations for the scenarios are listed in Table 2. Table 2. TDP simulation scenarios.  In all scenarios, the height h r of the airborne radar is consistent at 11 km, and the target height h t is set to be 5 km. Additionally, these scenarios are in the long-range radar propagation situations, which assumes that the slant range is almost same with the horizontal range (R t ≈ R h ). The detailed explanations for the scenarios are listed in Table 2.   Figure 8a illustrates the simulation results of the TDP along the range and the azimuth angle for the first scenario, where the deep red and blue colors indicate the TDP of 100% and 0%. The TDP over 90% is obtained at the scan angle of 0 • due to the high RCS level of the target aircraft, and the target is detected over the probability of 80% in the scan angle between −20 • and 60 • and the range from 50.5 to 153.8 km. On the other hand, the range from 0 to 50.2 km within the whole scan angle has the low TDP because this is the area below the airborne radar beam that cannot be reached. In addition, the scan loss of the airborne radar increases the path loss, which reduces the TDP in the wide angle from 60 • to 90 • . Figure 8b depicts the contour plot of the TDP over the 90% region according to the different atmospheric conditions of the refractivity. The solid, dashed, dotted, and dash-dotted lines denote the normal refraction, super-refraction, sub-refraction, and elevated ducting cases. We defined a detectable area by calculating the area inside the contour region to intuitively examine the TDP. The resulting detectable area of the super-refraction is larger than those of the others, which are similar to each other. It is because the super-refraction makes the wave propagation bend toward Earth's surface, which results in low path loss levels over a large area. The summary of the detectable area results for TDPs over 80% and 90% is shown in Figure 8c.

Scenario I Scenario II Scenario III
range from 0 to 50.2 km within the whole scan angle has the low TDP because this is the area below the airborne radar beam that cannot be reached. In addition, the scan loss of the airborne radar increases the path loss, which reduces the TDP in the wide angle from 60° to 90°. Figure 8b depicts the contour plot of the TDP over the 90% region according to the different atmospheric conditions of the refractivity. The solid, dashed, dotted, and dash-dotted lines denote the normal refraction, super-refraction, sub-refraction, and elevated ducting cases. We defined a detectable area by calculating the area inside the contour region to intuitively examine the TDP. The resulting detectable area of the superrefraction is larger than those of the others, which are similar to each other. It is because the super-refraction makes the wave propagation bend toward Earth's surface, which results in low path loss levels over a large area. The summary of the detectable area results for TDPs over 80% and 90% is shown in Figure 8c. (c) detectable area results for TDPs over 80% and 90%. Figure 9 illustrates the contour plot of the TDP over 80% according to the different weather conditions in the sub-and normal refractions for the second scenario. The solid and dashed lines indicate the TDP results according to the dry and rainy weather in the normal atmospheric condition, while the dotted and dash-dotted lines denote the results in the amorous atmospheric condition (super-refraction). In the dry weather, the TDP for (c) detectable area results for TDPs over 80% and 90%. Figure 9 illustrates the contour plot of the TDP over 80% according to the different weather conditions in the sub-and normal refractions for the second scenario. The solid and dashed lines indicate the TDP results according to the dry and rainy weather in the normal atmospheric condition, while the dotted and dash-dotted lines denote the results in the amorous atmospheric condition (super-refraction). In the dry weather, the TDP for the normal and super-refractions is similar to the results of the first scenario, where the detectable areas are 4738 and 6298 km 2 , respectively. In contrast, in rainy weather, the TDP of the airborne radar is extremely deteriorated because of the rainfall attenuation. The resulting detectable areas for the normal and super-refraction atmospheric environments are 1670 and 1734 km 2 . Note that the rainfall loss considerably affects the TDP of the airborne radar. In the third scenario, we additionally simulate and calculate the TDP of the airborne radar using another wave propagation simulation software (PETOOL) based on one-way and two-way split-step parabolic equation [48] and the ITU-R P.528-4 model [49] to confirm the results of the proposed process. Figure 10 shows the comparison of the TDP results along the target range of R t using the different propagation models for the third scenario when the target and the airborne radar are gradually getting closer to each other. The airborne radar TDP becomes significantly low (almost zero TDP) when the heavy rainfall occurs. In the target range of R t below 100 km, the AREPS simulation result well agrees with those of the PETOOL and ITU-R model. These results demonstrate that the proposed process is feasible for observing the TDP of the airborne radar considering the anomalous atmospheric environments.
(c) Figure 8. Simulation results of the TDP along the range and the azimuth angle for the first scenario: (a) TDP in a super-refraction; (b) contour plots of TDP in an anomalous atmospheric environment; (c) detectable area results for TDPs over 80% and 90%. Figure 9 illustrates the contour plot of the TDP over 80% according to the different weather conditions in the sub-and normal refractions for the second scenario. The solid and dashed lines indicate the TDP results according to the dry and rainy weather in the normal atmospheric condition, while the dotted and dash-dotted lines denote the results in the amorous atmospheric condition (super-refraction). In the dry weather, the TDP for the normal and super-refractions is similar to the results of the first scenario, where the detectable areas are 4738 and 6298 km 2 , respectively. In contrast, in rainy weather, the TDP of the airborne radar is extremely deteriorated because of the rainfall attenuation. The resulting detectable areas for the normal and super-refraction atmospheric environments are 1670 and 1734 km 2 . Note that the rainfall loss considerably affects the TDP of the airborne radar. In the third scenario, we additionally simulate and calculate the TDP of the airborne radar using another wave propagation simulation software (PETOOL) based on one-way and two-way split-step parabolic equation [48] and the ITU-R P.528-4 model [49] to confirm the results of the proposed process. Figure 10 shows the comparison of the TDP results along the target range of Rt using the different propagation models for the third scenario when the target and the airborne radar are gradually getting closer to each other. The airborne radar TDP becomes significantly low (almost zero TDP) when the heavy rainfall occurs. In the target range of Rt below 100 km, the AREPS simulation result well agrees with those of the PETOOL and ITU-R model. These results demonstrate that the proposed process is feasible for observing the TDP of the airborne radar considering the anomalous atmospheric environments.

Conclusions
In this paper, we proposed the TDP calculation process using the airborne long-range radar path loss based on the air-to-air scenarios in anomalous atmospheric and weather environments. The refractive index along an altitude was modeled using a quad-linear model. Then, the path loss including the refractivity model and other radar input parameters was simulated using the AREPS software along the range and the altitude. ITU-R models were used to consider the weather environment for the gas and rainfall attenua-

Conclusions
In this paper, we proposed the TDP calculation process using the airborne longrange radar path loss based on the air-to-air scenarios in anomalous atmospheric and weather environments. The refractive index along an altitude was modeled using a quadlinear model. Then, the path loss including the refractivity model and other radar input parameters was simulated using the AREPS software along the range and the altitude. ITU-R models were used to consider the weather environment for the gas and rainfall attenuations according to the range. In addition, the radar beam scan loss and RCS of the target were considered to estimate the TDP of the airborne long-range radar. Various air-to-air scenarios were assumed for the airborne radar TDP simulation in the anomalous atmospheric and weather environments. In the first scenario, the target was detected over the probability of 80% in the scan angle between −20 • and 60 • and the range from 50.5 to 153.8 km. In the second scenario, the TDP of the airborne radar was extremely deteriorated in the rainfall weather because of the rainfall attenuation. The resulting detectable areas for the normal and super-refraction atmospheric environments were 1670 and 1734 km 2 . In the third scenario, the airborne radar TDP became significantly low due to the rainfall.