Nothing Special   »   [go: up one dir, main page]

You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Sign in to use this feature.

Years

Between: -

Subjects

remove_circle_outline
remove_circle_outline
remove_circle_outline

Journals

Article Types

Countries / Regions

Search Results (18)

Search Parameters:
Keywords = dual polarization phased-array radar

Order results
Result details
Results per page
Select all
Export citation of selected articles as:
17 pages, 9836 KiB  
Article
An Algorithm to Retrieve Range Ocean Current Speed under Tropical Cyclone Conditions from Sentinel-1 Synthetic Aperture Radar Measurements Based on XGBoost
by Yuhang Zhou, Weizeng Shao, Ferdinando Nunziata, Weili Wang and Cheng Li
Remote Sens. 2024, 16(17), 3271; https://doi.org/10.3390/rs16173271 - 3 Sep 2024
Viewed by 590
Abstract
In this study, a novel algorithm to retrieve the current speed along the range direction under extreme sea states is developed from C-band synthetic aperture radar imagery. To this aim, a Sentinel-1 (S-1) dual-polarized synthetic aperture radar (SAR) dataset consisting of 2300 images [...] Read more.
In this study, a novel algorithm to retrieve the current speed along the range direction under extreme sea states is developed from C-band synthetic aperture radar imagery. To this aim, a Sentinel-1 (S-1) dual-polarized synthetic aperture radar (SAR) dataset consisting of 2300 images is collected during 200 tropical cyclones (TCs). The dataset is complemented with collocated wave simulations from the Wavewatch-III (WW3) model and reanalysis currents from the HYbrid Coordinate Ocean Model (HYCOM). The corresponding TC winds are officially released by IFRMER, while the Stokes drift following the wave propagation direction is estimated from the waves simulated by WW3. In this study, first the dependence of wind, Stokes drift, and range current on the Doppler centroid anomaly is investigated, and then the extreme gradient boosting (XGBoost) machine learning model is trained on 87% of the S-1 dataset for range current retrieval purposes. The rest of the dataset is used for testing the retrieval algorithm, showing a root mean square error (RMSE) and a correlation coefficient (r) of 0.11 m/s and 0.97, respectively, with the HYCOM outputs. A validation against measurements collected from two high-frequency (HF) phased-array radars is also performed, resulting in an RMSE and r of 0.12 m/s and 0.75, respectively. Those validation results are better than the 0.22 m/s RMSE and 0.28 r achieved by the empirical CDOP model. Hence, the experimental results confirm the soundness of the XGBoost, exhibiting a certain improvement over the empirical model. Full article
(This article belongs to the Special Issue SAR Monitoring of Marine and Coastal Environments)
Show Figures

Figure 1

Figure 1
<p>The frames of the 2300 S-1 SAR images overlaid on the TC tracks and maximum wind speeds, in which the red boxes indicate the training set and the black boxes indicate the test set.</p>
Full article ">Figure 2
<p>(<b>a</b>) S-1 VV-polarized NRCS image collected during the TC Infa (24 July 2021, 09:55 UTC); (<b>b</b>) S-1 VH-polarized NRCS image collected during the TC Infa (24 July 2021, 09:55 UTC); (<b>c</b>) corresponding CyclObs wind map.</p>
Full article ">Figure 3
<p>(<b>a</b>) SWH map simulated by WW3 at 10:00 UTC on 24 July 2021 relevant to the TC Infa, where the footprint of the HY-2B altimeter track is highlighted. (<b>b</b>) Validation of WW3-simulated SWH against the HY-2B products during the period of June–August 2022 in China seas.</p>
Full article ">Figure 4
<p>Current speed map simulated for 24 July 2021, 09:00 UTC, using HYCOM. The black box represents the footprint of the S-1 SAR image shown in <a href="#remotesensing-16-03271-f002" class="html-fig">Figure 2</a>a.</p>
Full article ">Figure 5
<p>Current speed maps measured by the HF phased-array radars over the TC Infa on (<b>a</b>) 24 July 2021, 09:55 UTC, and (<b>b</b>) 25 July 2021, 21:51 UTC. The black line and red dot are the typhoon track, and the red triangle is the current typhoon time position.</p>
Full article ">Figure 6
<p>DCA versus (<b>a</b>) wind speed, (<b>b</b>) Stokes drift, and (<b>c</b>) current speed, projected onto the range direction. Black lines stand for the linear regression results.</p>
Full article ">Figure 7
<p>The flow chart of the eXtreme gradient boosting (XGBoost).</p>
Full article ">Figure 8
<p>(<b>a</b>) The behavior of the XGBoost training process and (<b>b</b>) the SHAP value map.</p>
Full article ">Figure 9
<p>(<b>a</b>) Range current speed map by XGBoost method and (<b>b</b>) the empirical algorithm obtained from the S-1 SAR scene collected over TC Infa on 24 July 2021, 09:55 UTC; (<b>c</b>) the HYCOM range current speed map on 24 July 2021, 09:00.</p>
Full article ">Figure 10
<p>Latitude variation given the longitude—see the dashed red line in (<a href="#remotesensing-16-03271-f009" class="html-fig">Figure 9</a>a)—of the range current speed retrieved by XGBoost (red line), by CDOP model (green line), and simulated with HYCOM (black line).</p>
Full article ">Figure 11
<p>Validation of SAR-derived range current wind speeds by (<b>a</b>) XGBoost and (<b>b</b>) the CDOP model against HYCOM data.</p>
Full article ">Figure 12
<p>Validation of SAR-derived range current wind speeds by (<b>a</b>) XGBoost and (<b>b</b>) CDOP model against HF radar measurements.</p>
Full article ">Figure 13
<p>Taylor diagram with respect to current speed up to 3 m/s at intervals of 1 m/s, in which the red and blue symbols represent the result using XGBoost and the CDOP model.</p>
Full article ">Figure 14
<p>Variation in the bias (SAR-derived minus HYCOM data) of the range current speed with respect to (<b>a</b>) the HYCOM current speed for 0.375 m/s, (<b>b</b>) the wind speed for 10 m/s, (<b>c</b>) the Stokes drift for a 0.125 m/s bin, and (<b>d</b>) the DCA for a 12.5 Hz bin.</p>
Full article ">
18 pages, 11970 KiB  
Article
Contrasting the Effects of X-Band Phased Array Radar and S-Band Doppler Radar Data Assimilation on Rainstorm Forecasting in the Pearl River Delta
by Liangtao He, Jinzhong Min, Gangjie Yang and Yujie Cao
Remote Sens. 2024, 16(14), 2655; https://doi.org/10.3390/rs16142655 - 20 Jul 2024
Cited by 1 | Viewed by 793
Abstract
Contrasting the X-band phased array radar (XPAR) with the conventional S-Band dual-polarization mechanical scanning radar (SMSR), the XPAR offers superior temporal and spatial resolution, enabling a more refined depiction of the internal dynamics within convective systems. While both SMSR and XPAR data are [...] Read more.
Contrasting the X-band phased array radar (XPAR) with the conventional S-Band dual-polarization mechanical scanning radar (SMSR), the XPAR offers superior temporal and spatial resolution, enabling a more refined depiction of the internal dynamics within convective systems. While both SMSR and XPAR data are extensively used in monitoring and alerting for severe convective weather, their comparative application in numerical weather prediction through data assimilation remains a relatively unexplored area. This study harnesses the Weather Research and Forecasting Model (WRF) and its data assimilation system (WRFDA) to integrate radial velocity and reflectivity from the Guangzhou SMSR and nine XPARs across Guangdong Province. Utilizing a three-dimensional variational approach at a 1 km convective-scale grid, the assimilated data are applied to forecast a rainstorm event in the Pearl River Delta (PRD) on 6 June 2022. Through a comparative analysis of the results from assimilating SMSR and XPAR data, it was observed that the assimilation of SMSR data led to more extensive adjustments in the lower- and middle-level wind fields compared to XPAR data assimilation. This resulted in an enlarged convergence area at lower levels, prompting an overdevelopment of convective systems and an excessive concentration of internal hydrometeor particles, which in turn led to spurious precipitation forecasts. However, the sequential assimilation of both SMSR and XPAR data effectively reduced the excessive adjustments in the wind fields that were evident when only SMSR data were used. This approach diminished the generation of false echoes and enhanced the precision of quantitative precipitation forecasts. Additionally, the lower spectral width of XPAR data indicates its superior detection accuracy. Assimilating XPAR data alone yields more reasonable adjustments to the low- to middle-level wind fields, leading to the formation of small-to-medium-scale horizontal convergence lines in the lower levels of the analysis field. This enhancement significantly improves the model’s forecasts of composite reflectivity and radar echoes, aligning them more closely with actual observations. Consequently, the Threat Score (TS) and Equitable Threat Score (ETS) for heavy-rain forecasts (>10 mm/h) over the next 5 h are markedly enhanced. This study underscores the necessity of incorporating XPAR data assimilation in numerical weather prediction practices and lays the groundwork for the future joint assimilation of SMSR and XPAR data. Full article
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>) Location of Guangzhou SMSR (red pentagram) and nine XPARs in Guangdong Province (blue triangles), with the radar-coverage circles in 230 (solid red circle) and 60/40 (blue dashed circles) km ranges for SMSR and XPAR, respectively. The black dashed frame delineates the D03. (<b>b</b>) The position of automatic weather stations across Guangdong Province (green scatters).</p>
Full article ">Figure 2
<p>The potential geopotential height field (black contours, units: dagpm), temperature field (red contours, units: °C), wind field (wind barbs, units: m/s), and relative humidity field (shaded areas) from the ERA5 reanalysis data at 18:00 UTC on 6 June 2022 are depicted at different pressure levels: (<b>a</b>) 500 hPa; (<b>b</b>) 700 hPa; (<b>c</b>) 850 hPa; (<b>d</b>) 925 hPa. “D” in red represents the center of the cyclone. The brown line represents the wind shear line.</p>
Full article ">Figure 3
<p>Model grid configuration and topography (shaded). (<b>a</b>) Domain configuration; (<b>b</b>) D03 configuration.</p>
Full article ">Figure 4
<p>The flow chart for the DA experiments.</p>
Full article ">Figure 5
<p>Velocity spectrum width (SW) at 1.5° elevation angle for (<b>a</b>) SMSR and (<b>b</b>) XPARs. And (<b>c</b>) the spatial average velocity SW from the lowest to the highest of the first 9 elevation angles for both at 18:00 UTC on 6 June 2022. Unit: m/s.</p>
Full article ">Figure 6
<p>The average analysis increment in each model layer for the first assimilation cycle in the D03 region: (<b>a</b>) u (units: m/s), (<b>b</b>) v (units: m/s), (<b>c</b>) T (units: 10<sup>−3</sup> K), (<b>d</b>) q (units: g/kg).</p>
Full article ">Figure 7
<p>The first row shows the 850hPa horizontal wind field (vector) and wind speed (shaded, units: m/s) at 18:00 UTC on 6 June 2022, for (<b>a</b>) CTRL, (<b>b</b>) DA_S, (<b>c</b>) DA_X, and (<b>d</b>) DA_S_X. The second row depicts the incremental field of the horizontal wind field relative to CTRL, with wind speed greater than 5 m/s indicated by a red vector, for (<b>e</b>) DA_S, (<b>f</b>) DA_X, and (<b>g</b>) DA_S_X.</p>
Full article ">Figure 8
<p>Radar composite reflectivity-analysis field at 18:00 UTC in D03 on 6 June 2022, for (<b>a</b>) OBS-SMSR, (<b>b</b>) CTRL, (<b>c</b>) DA_S, (<b>d</b>) DA_X, and (<b>e</b>) DA_S_X; line AB is the profile in <a href="#remotesensing-16-02655-f009" class="html-fig">Figure 9</a> and <a href="#remotesensing-16-02655-f010" class="html-fig">Figure 10</a>.</p>
Full article ">Figure 9
<p>Vertical cross-sections of RF (shaded, units: dBZ) along the black solid line A (112.7°E, 21.6°N) and B (114.6°E, 23.2°N) in <a href="#remotesensing-16-02655-f008" class="html-fig">Figure 8</a> at 18:00 UTC on 6 June 2022. (<b>a</b>) OBS-SMSR, (<b>b</b>) CTRL, (<b>c</b>) DA_S, (<b>d</b>) DA_X, (<b>e</b>) DA_S_X. The height is from sea level.</p>
Full article ">Figure 10
<p>Vertical cross-sections of <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>q</mi> </mrow> <mrow> <mi>r</mi> </mrow> </msub> <mo> </mo> </mrow> </semantics></math>(first row), <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>q</mi> </mrow> <mrow> <mi>g</mi> </mrow> </msub> </mrow> </semantics></math> (second row), and <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>q</mi> </mrow> <mrow> <mi>s</mi> </mrow> </msub> </mrow> </semantics></math> (third row) along line AB for each experiment at 18:00 UTC on 6 June 2022. (<b>a</b>,<b>e</b>,<b>i</b>) CTRL, (<b>b</b>,<b>f</b>,<b>j</b>) DA_S, (<b>c</b>,<b>g</b>,<b>k</b>) DA_X, (<b>d</b>,<b>h</b>,<b>l</b>) DA_S_X. The black contours represent the distribution of hydrometeors in CTRL. The height is from sea level.</p>
Full article ">Figure 11
<p>Composite reflectivity at 19:00 UTC (first row), 20:00 UTC (second row), and 21:00 UTC (third row) on 6 June 2022 (units: dBZ). (<b>a</b>,<b>f</b>,<b>k</b>) OBS-SMSR, (<b>b</b>,<b>g</b>,<b>l</b>) CTRL, (<b>c</b>,<b>h</b>,<b>m</b>) DA_S, (<b>d</b>,<b>i</b>,<b>n</b>) DA_X, (<b>e</b>,<b>j</b>,<b>o</b>) DA_S_X.</p>
Full article ">Figure 12
<p>Vertical cross-sections of RF (shaded, units: dBZ) and wind (vector) along the black solid line A (113.0°E, 21.6°N) and B (114.5°E, 23.2°N) in <a href="#remotesensing-16-02655-f011" class="html-fig">Figure 11</a> at 19:00 UTC on 6 June 2022. (<b>a</b>) OBS-SMSR, (<b>b</b>) CTRL, (<b>c</b>) DA_S, (<b>d</b>) DA_X, (<b>e</b>) DA_S_X. The height is from sea level.</p>
Full article ">Figure 13
<p>Hourly precipitation from 18:00 UTC (first column) to 22:00 (last column) on 6 June 2022 (units: mm). (<b>a</b>–<b>d</b>) OBS, (<b>e</b>–<b>h</b>) CTRL, (<b>i</b>–<b>l</b>) DA_S, (<b>m</b>–<b>p</b>) DA_X, (<b>q</b>–<b>t</b>) DA_S_X.</p>
Full article ">Figure 14
<p>Threat Score (TS) (first row), False Alarm Rate (FAR) (second row) and Equitable Threat Score (ETS) (third row) for hourly precipitation from 18:00 to 22:00 on 6 June 2022 (<b>a</b>,<b>d,g</b>) for &gt;1 mm, (<b>b</b>,<b>e,h</b>) for &gt;5 mm, (<b>c</b>,<b>f,i</b>) for &gt;10 mm.</p>
Full article ">
24 pages, 15061 KiB  
Article
A Case Study on Two Differential Reflectivity Columns in a Convective Cell: Phased-Array Radar Observation and Cloud Model Simulation
by Gang Ren, Yue Sun, Hongping Sun, Yaning Dong, Yonglong Yang and Hui Xiao
Remote Sens. 2024, 16(3), 460; https://doi.org/10.3390/rs16030460 - 25 Jan 2024
Viewed by 1040
Abstract
A convective cell storm containing two differential reflectivity (ZDR) columns was observed with a dual-polarization phased-array radar (X-PAR) in Xixian County. Since a ZDR column is believed to correspond to a strong updraft and a single convective cell is considered [...] Read more.
A convective cell storm containing two differential reflectivity (ZDR) columns was observed with a dual-polarization phased-array radar (X-PAR) in Xixian County. Since a ZDR column is believed to correspond to a strong updraft and a single convective cell is considered to have a simple dynamic structure with one updraft core, how these two ZDR columns form and coexist is the focus of this study. The dynamic and microphysical structures around the two ZDR columns are studied under the mutual confirmation of the X-PAR observations and a cloud model simulation. The main ZDR column forms and maintains in an updraft whose bottom corresponds to a convergence of low-level and mid-level flow; it lasts from the early stages to the later stages. The secondary ZDR column emerges at the rear of the horizontal reflectivity (ZH) core relative to the moving direction of the cell; it forms in the middle stages and lasts for a shorter period, and its formation is under an air lifting forced by the divergent outflow of precipitation. Therefore, the secondary ZDR column is only a by-product in the middle stages of the convection rather than an indicator of a new or enhanced convection. Full article
Show Figures

Figure 1

Figure 1
<p>The topography around the X-PAR and the nearest operational radar (C-band).</p>
Full article ">Figure 2
<p>Local weather conditions near the radar site from ERA5 reanalysis data on 8 July 2021: (<b>a</b>) wind and air temperature and (<b>b</b>) convective available potential energy (CAPE). The wind barb and contour lines in (<b>a</b>) represent horizontal wind and air temperature (with 10 °C intervals), respectively; the bold black line represents the 0 °C level (approximately 4.73 km height). The vertical red line in (<b>b</b>) represents the emerging moment (09:18 UTC) of the convective cell observed with the X-PAR. The local standard time is Beijing Time (UTC + 8), and the time used in this paper is UTC.</p>
Full article ">Figure 3
<p>Example of a comparison of X-PAR and C-band radar on the same PPI surface. Data near 09:51 are used from these two radars. (<b>a</b>) Z<sub>H</sub> in the original 4th-level PPI of the C-band radar at 3.4° elevation. (<b>b</b>) The same as (<b>a</b>) but zoomed at the target cloud, and the X and Y coordinates are converted to be relative to the X-PAR. (<b>c</b>) The scatter plot of data in (<b>b</b>,<b>d</b>), where both of them have valid data points. (<b>d</b>) Z<sub>H</sub> of X-PAR that was interpolated to the same PPI surface of (<b>b</b>). The dashed lines in (<b>a</b>) mark beam ranges of the azimuth of the C-band radar, which cover the target cloud.</p>
Full article ">Figure 4
<p>Examples of original RHIs observed with the X-PAR, which contain Z<sub>DR</sub> columns: (<b>a</b>) Z<sub>e</sub> and the direction mark of the RHIs; (<b>b</b>) Z<sub>H</sub> at 72° azimuth; and (<b>c</b>) Z<sub>H</sub> at 79.2° azimuth; (<b>d</b>–<b>f</b>) are the same as (<b>a</b>–<b>c</b>) but for Z<sub>DRC</sub> or Z<sub>DR</sub>. The dashed lines in (<b>a</b>,<b>d</b>) represent 72° (upper) and 79.2° (lower) azimuths.</p>
Full article ">Figure 5
<p>Evolution of Z<sub>e</sub>, Z<sub>DRC</sub>, and Z<sub>DRw</sub> observed with X-PAR. X and Y represent west–east and south–north distances relative to the radar site. Lines A–F are the locations of selected typical vertical profiles, which will be analyzed in the following. The front side and rear side relative to the moving direction of the cloud are marked in the upper left sub-figure. The time is UTC.</p>
Full article ">Figure 6
<p>Evolution of Z<sub>e</sub>, Z<sub>DRC</sub>, and Z<sub>DRW</sub> simulated with the cloud model. X and Y represent west–east and south–north distances of the model domain. Lines A–F are the locations of selected typical vertical profiles, which will be analyzed in the following.</p>
Full article ">Figure 7
<p>Variables in profiles AA″ observed with the X-PAR and simulated with the cloud model. (<b>a</b>) Observed Z<sub>H</sub>, (<b>b</b>) observed Z<sub>DR</sub>, (<b>c</b>) observed RVD, (<b>d</b>) simulated Z<sub>H</sub>, (<b>e</b>) simulated Z<sub>DR</sub>, and (<b>f</b>) simulated HWD. The dashed lines represent the height where the background air temperature is 0 °C. The outlines of the simulated echo are limited to 10 dBZ. Black vectors are the simulated wind field.</p>
Full article ">Figure 8
<p>Same as <a href="#remotesensing-16-00460-f007" class="html-fig">Figure 7</a> but for profiles BB″.</p>
Full article ">Figure 9
<p>Same as <a href="#remotesensing-16-00460-f007" class="html-fig">Figure 7</a> but for profiles CC″.</p>
Full article ">Figure 10
<p>Same as <a href="#remotesensing-16-00460-f007" class="html-fig">Figure 7</a> but for profiles DD″.</p>
Full article ">Figure 10 Cont.
<p>Same as <a href="#remotesensing-16-00460-f007" class="html-fig">Figure 7</a> but for profiles DD″.</p>
Full article ">Figure 11
<p>Same as <a href="#remotesensing-16-00460-f007" class="html-fig">Figure 7</a> but for profiles EE″.</p>
Full article ">Figure 12
<p>Same as <a href="#remotesensing-16-00460-f007" class="html-fig">Figure 7</a> but for profiles FF″.</p>
Full article ">Figure 12 Cont.
<p>Same as <a href="#remotesensing-16-00460-f007" class="html-fig">Figure 7</a> but for profiles FF″.</p>
Full article ">Figure 13
<p>Typical horizontal distribution in the early stages. The selected height is where the background temperature is near 0 °C. (<b>a</b>) Observed Z<sub>H</sub> and Z<sub>DR</sub>; (<b>b</b>) simulated Z<sub>H</sub> and Z<sub>DR</sub>; and (<b>c</b>) simulated vertical air velocity. The thinner black line in (<b>c</b>) is 0 dBZ for Z<sub>H</sub>, and the thicker black line is 0.2 dB for Z<sub>DR</sub>, which are the same in (<b>b</b>).</p>
Full article ">Figure 14
<p>Same as <a href="#remotesensing-16-00460-f013" class="html-fig">Figure 13</a> but for the middle stages. The black thin line in (<b>c</b>) is 0 dBZ for Z<sub>H</sub>, and the thicker black line is 0.2 dB for Z<sub>DR</sub>, which are the same in (<b>b</b>).</p>
Full article ">Figure 15
<p>Demonstration of the secondary Z<sub>DR</sub> column splitting from the main Z<sub>DR</sub> column: (<b>a</b>) observed Z<sub>DR</sub> and (<b>b</b>) simulated Z<sub>DR</sub>. The unit of the shading is dB. The selected height is where the background temperature is near 0 °C. The depicted outlines of the cloud correspond to Z<sub>H</sub> = 0 dBZ.</p>
Full article ">Figure 16
<p>Relative distribution of different hydrometeors in profiles (<b>a</b>) AA″, (<b>b</b>) BB″, (<b>c</b>) CC″, (<b>d</b>) DD″, (<b>e</b>) EE″, and (<b>f</b>) FF″. The hydrometeors include ice crystals (I), snow (S), graupel (G), hail (H), frozen drops (FD), cloud drops (C)l and raindrops (R). The bold black lines are the outlines of Z<sub>H</sub> = 10 dBZ. Colored contour lines are based on simulated water content and are automatically produced according to their own value range in a specific time and profile instead of a fixed value. The scales of these contour lines are seen in <a href="#remotesensing-16-00460-t003" class="html-table">Table 3</a>.</p>
Full article ">Figure 17
<p>Same as <a href="#remotesensing-16-00460-f016" class="html-fig">Figure 16</a> but for the simulated vertical flux of the water content of raindrops.</p>
Full article ">Figure 18
<p>Same as <a href="#remotesensing-16-00460-f016" class="html-fig">Figure 16</a> but for the hydrometeor classification retrieved from X-PAR data.</p>
Full article ">
14 pages, 4621 KiB  
Article
Attenuation Correction of the X-Band Dual-Polarization Phased Array Radar Based on Observed Raindrop Size Distribution Characteristics
by Jiabao Feng, Xiantong Liu, Feng Xia, Yu Zhang and Xiaona Rao
Atmosphere 2023, 14(6), 1022; https://doi.org/10.3390/atmos14061022 - 14 Jun 2023
Viewed by 1368
Abstract
X-band dual-polarization phased array radar (XPAR-D) possesses high resolution and plays a significant role in detecting meso- and micro-scale convective systems. However, the precipitation attenuation it endures necessitates an effective correction method. This study selected radar data from XPAR-D at the peak of [...] Read more.
X-band dual-polarization phased array radar (XPAR-D) possesses high resolution and plays a significant role in detecting meso- and micro-scale convective systems. However, the precipitation attenuation it endures necessitates an effective correction method. This study selected radar data from XPAR-D at the peak of Maofeng Mountain in Guangzhou during 16–17 May 2020 from three precipitation stages after quality control. Attenuation coefficients were calculated for different precipitation types through scattering simulations of raindrop size distribution (RSD) data. Next, an attenuation correction algorithm (MZH-KDP method) was proposed for the radar reflectivity factor (ZH) according to different raindrop types and compared to the ZH-KDP method currently in use. The results indicate that the attenuation amount of XPAR-D echoes depends on the attenuation path and echo intensity. When the attenuation path is shorter and the echo intensity is weaker, the amount of attenuation and correction is smaller. Difficulties arise when there are noticeable deviations, which are challenging to resolve using attenuation correction methods. Longer attenuation paths and stronger echoes highlight the advantages of the MZH-KDP method, while the ZH-KDP method tends to overcorrect the bias. The MZH-KDP method outperforms the ZH-KDP method for different precipitation types. The superior correction capability of the MZH-KDP method provides a significant advantage in improving the performance of XPAR-D for the detection of extreme weather. Full article
(This article belongs to the Special Issue Monsoon and Typhoon Precipitation in Asia: Observation and Prediction)
Show Figures

Figure 1

Figure 1
<p>Spatial distribution of the 2D Video Disdrometer and radars.</p>
Full article ">Figure 2
<p>Flow chart for the quality control of the differential phase (Φ<sub>DP</sub>) and obtaining the specific differential phase (K<sub>DP</sub>).</p>
Full article ">Figure 3
<p>Flow chart for the attenuation correction algorithm of the MZ<sub>H</sub>-K<sub>DP</sub> method.</p>
Full article ">Figure 4
<p>Frequency distributions of Z<sub>H</sub>′-K<sub>DP</sub>′ from scattering simulations for different raindrop types, where SR represents a small raindrop type, MR represents a moderate raindrop type, and LR represents a large raindrop type.</p>
Full article ">Figure 5
<p>Attenuation coefficients and their fitting for (<b>a</b>) small raindrop type, (<b>b</b>) moderate raindrop type, (<b>c</b>) large raindrop type, and (<b>d</b>) all raindrop types.</p>
Full article ">Figure 6
<p>The radar reflectivity factor (Z<sub>H</sub>) at the 1.5° elevation of the S-band dual-polarization radar (SDPR) at (<b>a</b>) 19:06 (UTC+8) on 16 May, (<b>b</b>) 21:12 (UTC+8) on 17 May, and (<b>c</b>) 21:54 (UTC+8) on 17 May 2020.</p>
Full article ">Figure 7
<p>Comparison of the Z<sub>H</sub> at the 2.7° elevation angle from (<b>a</b>,<b>e</b>,<b>i</b>) the SDPR data and from the XPAR-D data (<b>b</b>,<b>f</b>,<b>j</b>) before correction and after correction with (<b>c</b>,<b>g</b>,<b>k</b>) the Z<sub>H</sub>-K<sub>DP</sub> method and (<b>d</b>,<b>h</b>,<b>l</b>) the MZ<sub>H</sub>-K<sub>DP</sub> method at (<b>a</b>–<b>d</b>) 19:06 (UTC+8) on 16 May, (<b>c</b>–<b>h</b>) 21:12 (UTC+8) on 17 May, and (<b>i</b>–<b>l</b>) 21:54 (UTC+8) on 17 May 2020.</p>
Full article ">Figure 8
<p>Comparisons of azimuth averaged Z<sub>H</sub> at the 2.7° elevation angle in (<b>a</b>) the range of 110°–140° at 19:06 (UTC+8) on 16 May, (<b>b</b>) the range of 40°–65° average Z<sub>H</sub> at 21:12 (UTC+8) on 17 May, and (<b>c</b>) the range of 80°–110° at 21:54 (UTC+8) on 17 May. The red line represents the Z<sub>H</sub> from the SDPR, the green line indicates the original Z<sub>H</sub> of the XPAR-D, and the blue and black lines denote the XPAR-D Z<sub>H</sub> corrected by the Z<sub>H</sub>-K<sub>DP</sub> and MZ<sub>H</sub>-K<sub>DP</sub> methods, respectively.</p>
Full article ">Figure 9
<p>Comparisons of the Z<sub>H</sub> probability density distribution of the SDPR data with the XPAR-D data (<b>a</b>,<b>d</b>,<b>g</b>) before and after the corrections of (<b>b</b>,<b>e</b>,<b>h</b>) the Z<sub>H</sub>-K<sub>DP</sub> method and (<b>c</b>,<b>f</b>,<b>i</b>) the MZ<sub>H</sub>-K<sub>DP</sub> method in terms of (<b>a</b>–<b>c</b>) local heavy precipitation, (<b>d</b>–<b>f</b>) stratocumulus precipitation, and (<b>g</b>–<b>i</b>) stratiform precipitation from 16 to 17 May 2020. Only colored areas with frequencies greater than 5% are displayed. “R”, “RMSE”, “NAE”, and “NRE” denote the correlation coefficients, root mean square errors, normalized absolute errors, and normalized relative errors of the Z<sub>H</sub> greater than 20 dBZ between the SDPR and the XPAR-D. The white line range shows the area with a maximum frequency standardization greater than 70%.</p>
Full article ">
16 pages, 10028 KiB  
Article
1-Bit Hexagonal Meander-Shaped Wideband Electronically Reconfigurable Transmitarray for Satellite Communications
by Qasim Ali, Yu Xiao, Shozab Shafiq, Wenhao Tan, Waseem Shahzad, Syed Muzahir Abbas and Houjun Sun
Electronics 2023, 12(9), 1957; https://doi.org/10.3390/electronics12091957 - 22 Apr 2023
Cited by 1 | Viewed by 1619
Abstract
This paper proposes a hexagonal meander-shaped wideband electronically reconfigurable transmitarray (HMRTA) at Ku band for satellite communications and radar applications. The proposed transmitarray offers high gain, low profile, and wideband characteristics with beam-scanning and beam-forming features. The cascaded structure is a low-profile and [...] Read more.
This paper proposes a hexagonal meander-shaped wideband electronically reconfigurable transmitarray (HMRTA) at Ku band for satellite communications and radar applications. The proposed transmitarray offers high gain, low profile, and wideband characteristics with beam-scanning and beam-forming features. The cascaded structure is a low-profile and compact transmitarray. The transmitter (Tx) layer has an angular hexagonal patch with a meandered shape and resonating parasitic patches to enhance the bandwidth. The receiver (Rx) layer comprises a two-part hexagonal receiver patch and a dual ring impedance-matching receiver layer. The current reversal phenomena have executed the 180° phase shift by integrating two diodes in opposite directions. The measured results of a unit cell achieve a minimum insertion loss of 0.86 dB and 0.92 dB for state I and state II. The maximum insertion loss is 2.58 dB from 14.12 GHz to 18.02 GHz and is about 24.83% at 16.5 GHz. The full-wave simulations of a 20 × 20 space-fed reconfigurable transmitarray were performed. Good radiation patterns at all scanning angles of two principal planes are achieved, and the cross-polarization level remains less than −20 dB. The simulated 3 dB gain fluctuation bandwidth of the array is 15.85~18.35 GHz, and the wideband characteristics are verified. The simulation results show that the array can perform beam scanning ±60° in the elevation (y-z) plane and obtain the beam-scanning characteristics for ±60° in the Azimuth (x-z) plane. Full article
(This article belongs to the Special Issue Sparse Array Design, Processing and Application)
Show Figures

Figure 1

Figure 1
<p>Three-dimensional exploded view: (<b>a</b>) top view, (<b>b</b>) bias circuit, (<b>c</b>) ground, (<b>d</b>) receiver layer, (<b>e</b>) impedance-matching receiving layer.</p>
Full article ">Figure 2
<p>Operating Principle of the Proposed Unit Cell.</p>
Full article ">Figure 3
<p>Electric field distribution (E-field) in different states: (<b>a</b>) state I, (<b>b</b>) state II.</p>
Full article ">Figure 4
<p>(<b>a</b>) Theoretical scattering parameter of the proposed unit cell. (<b>b</b>) Phase shift.</p>
Full article ">Figure 5
<p>(<b>a</b>) E-plane oblique incident performance. (<b>a</b>) Transmission coefficient (S21). (<b>b</b>) Reflection coefficients S11, S22.</p>
Full article ">Figure 6
<p>(<b>a</b>) E-plane oblique incident performance of transmission coefficient (S21). (<b>b</b>) E-plane oblique incident performance of transmission coefficient (S21).</p>
Full article ">Figure 7
<p>Bias lines in the proposed unit cell.</p>
Full article ">Figure 8
<p>Scattering parameters with and without bias lines.</p>
Full article ">Figure 9
<p>Fabricated unit cell: (<b>a</b>) top view, (<b>b</b>) bottom view, (<b>c</b>) size reference unit cell.</p>
Full article ">Figure 9 Cont.
<p>Fabricated unit cell: (<b>a</b>) top view, (<b>b</b>) bottom view, (<b>c</b>) size reference unit cell.</p>
Full article ">Figure 10
<p>Theoretical and measured scattering parameters: (<b>a</b>) magnitude for S<sub>11</sub>, S<sub>22</sub>, and S<sub>21</sub>; (<b>b</b>) measured phase difference, corresponding to states I and II.</p>
Full article ">Figure 11
<p>Horn feed: (<b>a</b>) horn model (<b>b</b>) 17 GHz simulated radiation pattern.</p>
Full article ">Figure 12
<p>Three-dimensional and side view of proposed 20 × 20 transmitarray.</p>
Full article ">Figure 13
<p>The 20 × 20 frequency versus gain graph.</p>
Full article ">Figure 14
<p>The 20 × 20 hexagonal meander-shaped transmitarray: (<b>a</b>) 24 pin connectors for biasing, (<b>b</b>) full scheme of the bias line.</p>
Full article ">Figure 15
<p>Wideband hexagonal unit cell: (<b>a</b>) E-plane phase distribution and gain spectrum, (<b>b</b>) H-plane phase distribution and gain spectrum.</p>
Full article ">Figure 15 Cont.
<p>Wideband hexagonal unit cell: (<b>a</b>) E-plane phase distribution and gain spectrum, (<b>b</b>) H-plane phase distribution and gain spectrum.</p>
Full article ">Figure 16
<p>A 20 × 20 gain radiation pattern for 20 × 20 array at 17.5 GHz. (<b>a</b>) E-plane max gain (0–60 degrees). (<b>b</b>) E-plane co-pol and cross-pol patterns.</p>
Full article ">
19 pages, 12389 KiB  
Article
Design of Broadband Low-RCS Array Antennas Based on Characteristic Mode Cancellation
by Jialiang Han, Dan Jia, Biao Du, Guodong Han, Yongtao Jia and Zekang Zhao
Electronics 2023, 12(7), 1536; https://doi.org/10.3390/electronics12071536 - 24 Mar 2023
Cited by 3 | Viewed by 1910
Abstract
In this letter, a design method for low radar cross section (RCS) array antennas based on characteristic mode cancellation (CMC) is presented. Based on the characteristic mode theory (CMT), two novel microstrip elements are designed by introducing rectangular slots and cross slots, which [...] Read more.
In this letter, a design method for low radar cross section (RCS) array antennas based on characteristic mode cancellation (CMC) is presented. Based on the characteristic mode theory (CMT), two novel microstrip elements are designed by introducing rectangular slots and cross slots, which produce 180° scattering phase difference by adjusting the size of slots. The dominant characteristic modes of the two elements achieve broadband dual-linear polarization CMC and similar radiation performances. The 4 × 4 array antenna consisting of these two antenna elements is designed. The operating band is from 4.55 GHz to 5.49 GHz (relative bandwidth 18.7%). The gain loss of the proposed array is about 0.1 dB compared to the reference array. The monostatic RCS is reduced for dual−linear polarized waves, and the 6 dB radar cross section reduction (RCSR) bandwidths are 62.3% and 35.7%, respectively. The prototype is fabricated and measured. The measured results of radiation pattern and RCS are in good agreement with the simulated results. Full article
(This article belongs to the Topic Antennas)
Show Figures

Figure 1

Figure 1
<p>Illustration of the array antenna composed of two kinds of elements.</p>
Full article ">Figure 2
<p>Geometry of the reference antenna and antenna A. (<b>a</b>) Reference antenna; (<b>b</b>) antenna A.</p>
Full article ">Figure 3
<p>Radiation performance simulation results of the reference antenna and antenna A. (<b>a</b>) S<sub>11</sub>; (<b>b</b>) far−field patterns.</p>
Full article ">Figure 4
<p><math display="inline"><semantics> <mrow> <msub> <mrow> <mi>A</mi> </mrow> <mrow> <mi>n</mi> </mrow> </msub> </mrow> </semantics></math> of the radiation characteristic modes of the reference antenna and antenna A. (<b>a</b>) Reference antenna; (<b>b</b>) antenna A.</p>
Full article ">Figure 5
<p>Radiation characteristic current distribution on the radiation patch at 5.1 GHz. (<b>a</b>) Reference antenna; (<b>b</b>) antenna A.</p>
Full article ">Figure 6
<p><math display="inline"><semantics> <mrow> <msub> <mrow> <mi>A</mi> </mrow> <mrow> <mi>n</mi> </mrow> </msub> </mrow> </semantics></math> of scattering characteristic modes of the reference antenna and antenna A. (<b>a</b>) Reference antenna; (<b>b</b>) antenna A.</p>
Full article ">Figure 7
<p>The <math display="inline"><semantics> <mrow> <msub> <mrow> <mi mathvariant="sans-serif">ϕ</mi> </mrow> <mrow> <mi>n</mi> </mrow> </msub> </mrow> </semantics></math> and the difference of <math display="inline"><semantics> <mrow> <msub> <mrow> <mi mathvariant="sans-serif">ϕ</mi> </mrow> <mrow> <mi>n</mi> </mrow> </msub> </mrow> </semantics></math> of the dominant scattering characteristic modes of the reference antenna and antenna A. (<b>a</b>) The <math display="inline"><semantics> <mrow> <msub> <mrow> <mi mathvariant="sans-serif">ϕ</mi> </mrow> <mrow> <mi>n</mi> </mrow> </msub> </mrow> </semantics></math> of the reference antenna and antenna A; (<b>b</b>) difference of <math display="inline"><semantics> <mrow> <msub> <mrow> <mi mathvariant="sans-serif">ϕ</mi> </mrow> <mrow> <mi>n</mi> </mrow> </msub> </mrow> </semantics></math> of the dominant scattering characteristic modes of the reference antenna and antenna A.</p>
Full article ">Figure 8
<p>Scattering characteristic current distribution and far−field patterns comparison at 4.4 GHz. (<b>a</b>) Mode 1 of the reference antenna; (<b>b</b>) mode 1 of antenna A.</p>
Full article ">Figure 9
<p>Scattering characteristic currents distribution and far−field patterns comparison for at 5.0 GHz. (<b>a</b>) Mode 2 of the reference antenna; (<b>b</b>) mode 2 of antenna A.</p>
Full article ">Figure 10
<p>Geometry of antenna B.</p>
Full article ">Figure 11
<p>Radiation characteristic modes analysis results of antenna B. (<b>a</b>) <math display="inline"><semantics> <mrow> <msub> <mrow> <mi mathvariant="normal">A</mi> </mrow> <mrow> <mi>n</mi> </mrow> </msub> </mrow> </semantics></math>; (<b>b</b>) characteristic current distribution at 5.1 GHz.</p>
Full article ">Figure 12
<p>Radiation performances of antenna A and antenna B. (<b>a</b>) S<sub>11</sub>; (<b>b</b>) far−field patterns.</p>
Full article ">Figure 13
<p>Scattering characteristic mode analysis results of antenna B. (<b>a</b>) <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>A</mi> </mrow> <mrow> <mi>n</mi> </mrow> </msub> </mrow> </semantics></math> of antenna B; (<b>b</b>) <math display="inline"><semantics> <mrow> <msub> <mrow> <mi mathvariant="sans-serif">ϕ</mi> </mrow> <mrow> <mi>n</mi> </mrow> </msub> </mrow> </semantics></math> of antenna A and antenna B; (<b>c</b>) phase difference of antenna A and antenna B.</p>
Full article ">Figure 14
<p>Scattering characteristic currents distribution and far−field patterns comparison for at 7.0 GHz. (<b>a</b>) Mode 2 of antenna A; (<b>b</b>) mode 2 of antenna B.</p>
Full article ">Figure 15
<p><math display="inline"><semantics> <mrow> <msub> <mrow> <mi>A</mi> </mrow> <mrow> <mi>n</mi> </mrow> </msub> </mrow> </semantics></math> of scattering characteristic modes of the reference antenna and antenna A. (<b>a</b>) Antenna A; (<b>b</b>) antenna B.</p>
Full article ">Figure 16
<p><math display="inline"><semantics> <mrow> <msub> <mrow> <mi mathvariant="sans-serif">ϕ</mi> </mrow> <mrow> <mi>n</mi> </mrow> </msub> </mrow> </semantics></math> and difference of <math display="inline"><semantics> <mrow> <msub> <mrow> <mi mathvariant="sans-serif">ϕ</mi> </mrow> <mrow> <mi>n</mi> </mrow> </msub> </mrow> </semantics></math> of the dominant scattering characteristic modes of the reference antenna and antenna A. (<b>a</b>) <math display="inline"><semantics> <mrow> <msub> <mrow> <mi mathvariant="sans-serif">ϕ</mi> </mrow> <mrow> <mi>n</mi> </mrow> </msub> </mrow> </semantics></math> of antenna A and antenna B; (<b>b</b>) difference of <math display="inline"><semantics> <mrow> <msub> <mrow> <mi mathvariant="sans-serif">ϕ</mi> </mrow> <mrow> <mi>n</mi> </mrow> </msub> </mrow> </semantics></math> of the dominant scattering characteristic modes of the reference antenna and antenna A.</p>
Full article ">Figure 17
<p>Four forms of proposed arrays. (<b>a</b>) ABAB; (<b>b</b>) ABBA; (<b>c</b>) AABB; (<b>d</b>) checkerboard.</p>
Full article ">Figure 18
<p>RCS simulation results of different array forms.</p>
Full article ">Figure 19
<p>Geometry of reference array.</p>
Full article ">Figure 20
<p>Radiation far−field patterns of the reference and proposed array antenna.</p>
Full article ">Figure 21
<p>Monostatic RCS simulated results of the reference and proposed array antennas. (<b>a</b>) <span class="html-italic">x</span>−polarized; (<b>b</b>) <span class="html-italic">y</span>−polarized.</p>
Full article ">Figure 22
<p>Fabricated prototype of the proposed array antenna. (<b>a</b>) Structural diagram; (<b>b</b>) radiation far-field patterns test scene.</p>
Full article ">Figure 23
<p>Measured and simulated radiation performances of the proposed patch array antenna. (<b>a</b>) S<sub>11</sub>; (<b>b</b>) far−field patterns.</p>
Full article ">Figure 24
<p>The RCS test scenario of the array antenna prototype. (<b>a</b>) Array antenna prototype; (<b>b</b>) the RCS test environment.</p>
Full article ">Figure 25
<p>Measured and simulated monostatic RCS of the proposed array antenna for vertical incident plane waves. (<b>a</b>) <span class="html-italic">x</span>−polarized; (<b>b</b>) <span class="html-italic">y</span>−polarized.</p>
Full article ">Figure 26
<p>RCS of the array antenna in conditions of gap and displacement errors. (<b>a</b>) The space gap between the double PCB layers; (<b>b</b>) displacement between the double PCB layers.</p>
Full article ">
19 pages, 964 KiB  
Article
Study on the Backscatter Differential Phase Characteristics of X-Band Dual-Polarization Radar and its Processing Methods
by Fei Geng and Liping Liu
Remote Sens. 2023, 15(5), 1334; https://doi.org/10.3390/rs15051334 - 27 Feb 2023
Viewed by 1721
Abstract
The differential propagation phase (ΦDP) of X-band dual-polarization weather radar (including X-band dual-polarization phased-array weather radar, X-PAR) is important for estimating precipitation and classifying hydrometeors. However, the measured differential propagation phase contains the backscatter differential phase (δ), which [...] Read more.
The differential propagation phase (ΦDP) of X-band dual-polarization weather radar (including X-band dual-polarization phased-array weather radar, X-PAR) is important for estimating precipitation and classifying hydrometeors. However, the measured differential propagation phase contains the backscatter differential phase (δ), which poses difficulties for the application of the differential propagation phase from X-band radars. This paper presents the following: (1) the simulation and characteristics analysis of the backscatter differential phase based on disdrometer DSD (raindrop size distribution) measurement data; (2) an improved method of the specific differential propagation phase (KDP) estimation based on linear programming and backscatter differential phase elimination; (3) the effect of backscatter differential phase elimination on the specific differential propagation phase estimation of X-PAR. The results show the following: (1) For X-band weather radar, the raindrop equivalent diameters D > 2 mm may cause a backscatter differential phase between 0 and 20°; in particular, the backscatter differential phase varies sharply with raindrop size between 3.2 and 4.5 mm. (2) Using linear programming or smoothing filters to process the differential propagation phase could suppress the backscatter differential phase, but it is hard to completely eliminate the effect of the backscatter differential phase. (3) Backscatter differential phase correction may improve the calculation accuracy of the specific differential propagation phase, and the optimization was verified by the improved self-consistency of polarimetric variables, correlation between specific differential propagation phase estimations from S- and X-band radar and the accuracy of quantitative precipitation estimation. The X-PAR deployed in Shenzhen showed good observation performance and the potential to be used in radar mosaics with S-band weather radar. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>The relationship between <span class="html-italic">δ</span> and raindrops’ equivalent diameter for (<b>a</b>) S-, (<b>b</b>) C-, and (<b>c</b>) X-band at different temperatures.</p>
Full article ">Figure 2
<p>Scatterplots of the <span class="html-italic">δ</span> calculated from DSD measurement with <span class="html-italic">Z<sub>DR</sub></span> for the (<b>a</b>) C-band and (<b>b</b>) X- band, and of <span class="html-italic">δ</span> calculated from gamma assumptions of DSD with <span class="html-italic">Z<sub>DR</sub></span> for the (<b>c</b>) C-band and (<b>d</b>) X- band. The red line represents the polynomial-fitted <span class="html-italic">δ–Z<sub>DR</sub></span> relationship curve.</p>
Full article ">Figure 3
<p>Simulated (<b>a</b>) <span class="html-italic">Z<sub>H</sub></span> and Z<sub>DR</sub> curves with distance using DSD measurements from 12:12 to 15:58 on 10 June 2019, (<b>b</b>) <span class="html-italic">Φ<sub>DP</sub></span> and <span class="html-italic">K<sub>DP</sub></span> calculated from DSD with <span class="html-italic">δ</span> effects (red curve) and without <span class="html-italic">δ</span> effects (black curve), (<b>c</b>) Rain intensity calculated by <span class="html-italic">K<sub>DP</sub></span> without <span class="html-italic">δ</span> effects (R<sub>1</sub>, blue curve) and with <span class="html-italic">δ</span> effects (R<sub>2</sub>, red curve), and (<b>d</b>) Backscatter differential phase.</p>
Full article ">Figure 4
<p>Comparison of the actual attenuation and attenuation correction for (<b>a</b>) <span class="html-italic">Z<sub>H</sub></span> and (<b>b</b>) <span class="html-italic">Z<sub>DR</sub></span>. The black line shows <span class="html-italic">Z<sub>H</sub></span> or <span class="html-italic">Z<sub>DR</sub></span>, the red solid line shows the actual PIA for reflectivity and path integrated differential attenuation (PIDA) for differential reflectivity, and the blue line shows the attenuation correction (<span class="html-italic">PIA<sub>δ</sub></span> and <span class="html-italic">PIDA<sub>δ</sub></span>); the differences between the red and blue curves are the errors in the attenuation correction caused by <span class="html-italic">δ</span>.</p>
Full article ">Figure 5
<p>(<b>a</b>) <span class="html-italic">K<sub>DP</sub></span> (unaffected by <span class="html-italic">δ</span>) and <span class="html-italic">K<sub>DPδ</sub></span> calculated from <span class="html-italic">Φ<sub>DPR</sub></span> (affected by <span class="html-italic">δ</span>); (<b>b</b>) <span class="html-italic">K<sub>DP</sub></span><sub>1</sub> after the elimination of <span class="html-italic">δ</span> from <span class="html-italic">Φ<sub>DPR</sub></span> using the <span class="html-italic">Z<sub>DR</sub></span> without attenuation effect, and <span class="html-italic">K<sub>DP</sub></span><sub>2</sub> using the corrected <span class="html-italic">Z<sub>DR</sub>, K<sub>DP</sub></span><sub>1</sub>, and <span class="html-italic">K<sub>DP</sub></span><sub>2</sub> denote the corrected <span class="html-italic">K<sub>DP</sub></span> for the two tests, respectively; (<b>c</b>) rain intensity estimated from <span class="html-italic">K<sub>DP1</sub> (R<sub>1</sub>)</span> and <span class="html-italic">K<sub>DP2</sub> (R<sub>2</sub>)</span>; (<b>d</b>) the actual rainfall intensity.</p>
Full article ">Figure 6
<p>Colored density plot of the simulated <span class="html-italic">K<sub>DP</sub></span> versus (<b>a</b>) the <span class="html-italic">K<sub>DPδ</sub></span> affected by <span class="html-italic">δ</span> and (<b>b</b>) the corrected <span class="html-italic">K<sub>DP1</sub></span> based on the <span class="html-italic">δ</span>-elimination using <span class="html-italic">Z<sub>DR</sub></span>.</p>
Full article ">Figure 7
<p>Typical radar sounding of (<b>a</b>) <span class="html-italic">Φ<sub>DP</sub></span>, <span class="html-italic">Z<sub>DR</sub></span>, and <span class="html-italic">δ</span> obtained by <span class="html-italic">Z<sub>DR</sub></span>, and <span class="html-italic">Φ<sub>DP2</sub></span> obtained by <span class="html-italic">δ</span> correction; (<b>b</b>) LP-processed results obtained before (blue curve) and after (red curve) <span class="html-italic">δ</span> correction; (<b>c</b>) <span class="html-italic">Z<sub>H</sub></span> after attenuation correction (dBZ) and <span class="html-italic">K<sub>DP</sub></span> before and after <span class="html-italic">δ</span> correction using LS; and (<b>d</b>) <span class="html-italic">Z<sub>H</sub></span> after attenuation correction (dBZ) and <span class="html-italic">K<sub>DP</sub></span> before and after <span class="html-italic">δ</span> correction using LP.</p>
Full article ">Figure 8
<p>Radar PPI of (<b>a</b>) <span class="html-italic">Z<sub>H</sub></span> from S-POL, (<b>b</b>) <span class="html-italic">Z<sub>DR</sub></span> from S-POL, (<b>c</b>) <span class="html-italic">K<sub>DP</sub></span> from S-POL, (<b>d</b>) <span class="html-italic">Z<sub>H</sub></span> from X-PAR, (<b>e</b>) <span class="html-italic">Z<sub>DR</sub></span> from X-PAR, (<b>f</b>) <span class="html-italic">K<sub>DP</sub></span> from X-PAR calculated by Exp1 (LS without <span class="html-italic">δ</span>-elimination), (<b>g</b>) <span class="html-italic">K<sub>DP</sub></span> from X-PAR calculated by Exp2 (LS with <span class="html-italic">δ</span>-elimination), (<b>h</b>) <span class="html-italic">K<sub>DP</sub></span> from X-PAR calculated by Exp3 (LP without <span class="html-italic">δ</span>-elimination), (<b>i</b>) <span class="html-italic">K<sub>DP</sub></span> from X-PAR calculated by Exp4 (LP with <span class="html-italic">δ</span>-elimination).</p>
Full article ">Figure 9
<p>Based on the <span class="html-italic">δ</span>-elimination processing of <span class="html-italic">Φ<sub>DP</sub></span>, this figure shows the distribution of (<b>a</b>) <span class="html-italic">K<sub>DPS</sub></span> from S-POL vs. <span class="html-italic">K<sub>DPX</sub></span> calculated by LS, (<b>b</b>) <span class="html-italic">K<sub>DPS</sub></span> from S-POL vs. <span class="html-italic">K<sub>DPX</sub></span> calculated by LP, (<b>c</b>) <span class="html-italic">K<sub>DPM</sub></span> calculated using fitting relation of polarimetric variables vs. <span class="html-italic">K<sub>DPX</sub></span> calculated by LS, (<b>d</b>) <span class="html-italic">K<sub>DPM</sub></span> calculated using fitting relation of polarimetric variables vs. <span class="html-italic">K<sub>DPX</sub></span> calculated by LP.</p>
Full article ">Figure 10
<p>Colored density plot of 1 h rainfall by rain gauges vs. (<b>a</b>) QPE by EXP1 <span class="html-italic">K<sub>DP</sub></span> (LS without δ-elimination) and (<b>b</b>) QPE by EXP4 <span class="html-italic">K<sub>DP</sub></span> (LP with <span class="html-italic">δ</span>-elimination).</p>
Full article ">
22 pages, 5580 KiB  
Article
Study on Attenuation Correction for the Reflectivity of X-Band Dual-Polarization Phased-Array Weather Radar Based on a Network with S-Band Weather Radar
by Fei Geng and Liping Liu
Remote Sens. 2023, 15(5), 1333; https://doi.org/10.3390/rs15051333 - 27 Feb 2023
Cited by 1 | Viewed by 1778
Abstract
X-band dual-polarization phased-array weather radars (X-PARs) have been used in South China extensively. Eliminating the attenuation and system bias of X-band radar data is the key to utilizing the advantage of X-PAR networks. In this paper, the disdrometer raindrop-size distribution (DSD) measurements are [...] Read more.
X-band dual-polarization phased-array weather radars (X-PARs) have been used in South China extensively. Eliminating the attenuation and system bias of X-band radar data is the key to utilizing the advantage of X-PAR networks. In this paper, the disdrometer raindrop-size distribution (DSD) measurements are used to calculate the radar polarimetric variables and analyze the characteristics of precipitation attenuation. Furthermore, based on the network of S-band dual-polarization Doppler weather radar (S-POL) and X-PARs, an attenuation-correction method for X-PAR reflectivity is proposed with S-POL constraints in view of the radar-mosaic requirements of a multi-radar network. Linear programming is used to calculate the attenuation-correction parameters of different rainfall areas, which realizes the attenuation correction for X-PAR. The results show that the attenuation-correction parameters simulated based on the disdrometer DSD vary with different precipitation classification; the attenuation-corrected reflectivity of X-PARs is consistent with S-POL and can realize a more precise observation of the evolution of the convective system. Compared with previous attenuation-correction methods with constant correction parameters, the improved method can reduce the deviation between X-PAR reflectivity and that of S-POL in heavy rainfall areas and areas of strong attenuation. The method proposed in this paper is stable and effective. After effective quality control, it is found that the X-PAR network deployed in South China observes data accurately and is consistent with S-POL; thus, it is expected to achieve high temporal–spatial resolution within a radar mosaic. Full article
Show Figures

Figure 1

Figure 1
<p>Location and Digital Elevation Model (DEM) information of the S-POL and two X-PARs. The center of the red and black circles shows the sites of the S-POL and X-PARs, respectively. The distance of the black concentric circle is 10 km, and the distance of the red concentric circle is 50 km.</p>
Full article ">Figure 2
<p>Calculation of <span class="html-italic">A<sub>H</sub></span> for X-PAR attenuation correction based on networked S-band radar and precipitation classification.</p>
Full article ">Figure 3
<p>Eight-gate interpolation diagram, where <span class="html-italic">P</span> is the X-PAR observation gate, (<span class="html-italic">a</span>, <span class="html-italic">e</span>, <span class="html-italic">r</span>) is the S-POL radar coordinate of <span class="html-italic">P</span>, and <math display="inline"><semantics> <mrow> <msub> <mi>f</mi> <mn>1</mn> </msub> </mrow> </semantics></math>–<math display="inline"><semantics> <mrow> <msub> <mi>f</mi> <mn>8</mn> </msub> </mrow> </semantics></math> are the eight gates adjacent to <span class="html-italic">P</span> in S-POL observation.</p>
Full article ">Figure 4
<p>Exponential-fitting results of <span class="html-italic">Z<sub>H</sub></span> and <span class="html-italic">A<sub>H</sub></span> (<b>a</b>); linear relationship (blue line) and exponential relationship (red line) between <span class="html-italic">K<sub>DP</sub></span> and <span class="html-italic">A<sub>H</sub></span> (<b>b</b>); probability distribution of the <span class="html-italic">γ</span> value calculated according to the linear relationship for <span class="html-italic">Z<sub>H</sub></span> &lt; 45 dBZ (<b>c</b>) and for <span class="html-italic">Z<sub>H</sub></span> &gt; 45 dBZ (<b>d</b>).</p>
Full article ">Figure 5
<p>Simulated polarimetric variables and <span class="html-italic">A<sub>H</sub></span> based on DSD measurement in typical rainfall process.</p>
Full article ">Figure 6
<p>Simulated reflectivity distribution at X-band (<span class="html-italic">Z<sub>HX</sub></span>) and S-band (<span class="html-italic">Z<sub>HS</sub></span>) using DSD measurement data.</p>
Full article ">Figure 7
<p>Precipitation classification (<b>a</b>,<b>d</b>), <math display="inline"><semantics> <mrow> <msub> <mi>φ</mi> <mrow> <mi>D</mi> <mi>P</mi> </mrow> </msub> </mrow> </semantics></math> (<b>b</b>,<b>e</b>), and <span class="html-italic">PIA</span> (<b>c</b>,<b>f</b>) of X-PAR1 at 4.5° elevation-angle PPI from 01:44 to 01:50, where (<b>a</b>–<b>c</b>) are the calculation results at 01:44, (<b>d</b>–<b>f</b>) are the calculation results at 01:50, (<b>c</b>) is the <span class="html-italic">PIA</span><sub>0</sub> calculated by Equation (16), and (<b>f</b>) is the <span class="html-italic">PIA<sub>t</sub></span> calculated by attenuation-correction method proposed in this paper.</p>
Full article ">Figure 8
<p>Reflectivity measured by X-PAR and S-POL from 01:44 to 01:50 (UTC) on 5 June 2020. (<b>a</b>–<b>e</b>) are <span class="html-italic">Z<sub>H</sub></span> of X-PAR without attenuation correction at 01:44:32, 01:46: 04, 01:47:38, 01:49:10, and 01:50:42; (<b>f</b>–<b>j</b>) are <span class="html-italic">Z<sub>H</sub></span> of X-PAR after attenuation correction; and (<b>k</b>,<b>l</b>) are <span class="html-italic">Z<sub>SX</sub></span> (<span class="html-italic">Z<sub>H</sub></span> from S-POL and convert to X-PAR) at 01:42 and 01:48, respectively; A, B and C are three rapidly evolving cells in the precipitation system.</p>
Full article ">Figure 9
<p>Frequency plots of reflectivity from S-POL (after wavelength conversion and system-bias correction) and X-PAR before (<b>a</b>) and after (<b>b</b>) attenuation correction and frequency plots of the gates with <span class="html-italic">φ<sub>DP</sub></span> &gt; 40° before (<b>c</b>) and after (<b>d</b>) attenuation correction at 01:50. <span class="html-italic">n/N</span> is the frequency of the X-PAR and S-POL at the corresponding reflectivity.</p>
Full article ">Figure 10
<p>Frequency plots of Z<span class="html-italic"><sub>H</sub></span>, and <span class="html-italic">A<sub>H</sub></span> before (<b>a</b>) and after (<b>b</b>) attenuation correction at 01:50 and frequency plots of gates with <span class="html-italic">φ<sub>DP</sub></span> &gt; 40° before (<b>c</b>) and after (<b>d</b>) attenuation correction, where the gray curve is the fitting results of DSD.</p>
Full article ">
17 pages, 14993 KiB  
Article
Study on the Quantitative Precipitation Estimation of X-Band Dual-Polarization Phased Array Radar from Specific Differential Phase
by Guo Zhao, Hao Huang, Ye Yu, Kun Zhao, Zhengwei Yang, Gang Chen and Yu Zhang
Remote Sens. 2023, 15(2), 359; https://doi.org/10.3390/rs15020359 - 6 Jan 2023
Cited by 4 | Viewed by 2439
Abstract
In this study, the quantitative precipitation estimation (QPE) capability of three X-band dual-polarization phased array radars (PAR) in Guangzhou, South China, was demonstrated, with an S-band operational dual-polarization radar as the benchmark. Rainfall rate (R) estimators based on the specific differential [...] Read more.
In this study, the quantitative precipitation estimation (QPE) capability of three X-band dual-polarization phased array radars (PAR) in Guangzhou, South China, was demonstrated, with an S-band operational dual-polarization radar as the benchmark. Rainfall rate (R) estimators based on the specific differential phase (KDP) for summer precipitation for both X-band and S-band radars were derived from the raindrop size distributions (DSDs) observed by a 2-dimensional video disdrometer (2DVD) in South China. Rainfall estimates from the radars were evaluated with gauge observations in three events, including pre-summer rainfall, typhoon precipitation, and local severe convective precipitation. Observational results showed that radar echoes from the X-band PARs suffered much more severely from attenuation than those from the S-band radar. Compared to S-band observations, the X-band echoes can disappear when the signal-to-noise ratio drops to a certain level due to severe attenuation, resulting in different estimated rainfall areas for X- and S-band radars. The attenuation corrected by KDP had good consistency with S-band observations, but the accuracy of attenuation correction was affected by DSD uncertainty and may vary in different types of precipitation. The QPE results demonstrated that the R(KDP) estimator produced better rainfall accumulations from the X-band PAR observations compared to the S-band observations. For both the X-band and S-band radars, the estimates of hourly accumulated rainfall became more accurate in heavier rainfall, due to the decreases of both the DSD uncertainty and the impact of measurement errors. In the heavy precipitation area, the estimation accuracy of the X-band radar was high, and the overestimation of the S-band radar was obvious. Through the analysis of the ZH-ZDR distribution in the three weather events, it was found that the X-band PAR with the capability of high spatiotemporal observations can capture minute-level changes in the microphysical characteristics, which help improve the estimation accuracy of ground rainfall. Full article
Show Figures

Figure 1

Figure 1
<p>Location of X-band dual-polarization phased array radar and S-band dual-polarization radar at Guangzhou on the digital terrain elevation map.</p>
Full article ">Figure 2
<p><span class="html-italic">R</span>(<span class="html-italic">K</span><sub>DP</sub>) relationship for X-band (<b>a</b>) and S-band (<b>b</b>) radars, derived by the least squares method of the observations.</p>
Full article ">Figure 3
<p>A comparison between the <span class="html-italic">R</span>(<span class="html-italic">K</span><sub>DP</sub>) and the intrinsic rainfall rate <span class="html-italic">R</span><sub>cal</sub> directly observed by the 2DVD for X-band (<b>a</b>) and S-band (<b>b</b>).</p>
Full article ">Figure 4
<p>The <span class="html-italic">K</span><sub>DP</sub>-<span class="html-italic">A<sub>H</sub></span> (<b>a</b>) and <span class="html-italic">K</span><sub>DP</sub>-<span class="html-italic">A</span><sub>DP</sub> (<b>b</b>) relationships for the X-band PAR fitted by DSD data.</p>
Full article ">Figure 5
<p>Measurements of <span class="html-italic">Z<sub>H</sub></span> from the Panyu PAR (<b>a</b>) without attenuation correction and (<b>b</b>) after attenuation correction at the 6.3 deg elevation at 19:30 on 2 May 2018; (<b>c</b>) Measurements of <span class="html-italic">Z<sub>H</sub></span> from the Guangzhou S-band radar at the 6 deg elevation at 19:30 on 2 May 2018.</p>
Full article ">Figure 6
<p>The <span class="html-italic">Z<sub>H</sub></span> of the X-band PAR before (blue) and after (red) attenuation correction at an azimuth of 5 deg and an elevation of 6.3 deg corresponding to <a href="#remotesensing-15-00359-f005" class="html-fig">Figure 5</a>. The black line is from the S-band radar, also at an azimuth of 5 deg and an elevation of 6 deg.</p>
Full article ">Figure 7
<p>The NDFs of estimated hourly accumulated rainfall from (<b>a</b>) X-band <span class="html-italic">R</span>(<span class="html-italic">K</span><sub>DP</sub>) and (<b>b</b>) S-band <span class="html-italic">R</span>(<span class="html-italic">K</span><sub>DP</sub>) versus the gauge observations for all precipitation cases.</p>
Full article ">Figure 8
<p><span class="html-italic">Z<sub>H</sub></span>, <span class="html-italic">Z</span><sub>DR</sub>, and <span class="html-italic">K</span><sub>DP</sub> observed by the X-band PAR radar at the 2.7-deg-elevation in the (<b>a</b>–<b>c</b>) pre-summer rainband at 10:30 on June 7, the (<b>d</b>–<b>f</b>) typhoon rainband at 13:30 on August 19, and the (<b>g</b>–<b>i</b>) severe convective rainband at 17:30 on September 4.</p>
Full article ">Figure 9
<p>Comparison of hourly accumulated rainfall obtained using <span class="html-italic">R</span>(<span class="html-italic">K</span><sub>DP</sub>) with gauge observations in the pre-summer rainfall between 10:00–11:00 on June 7, the typhoon precipitation between 13:00–14:00 on August 19, and the severe convective precipitation rainband between 17:00–18:00 on September 4. (<b>a</b>,<b>c</b>,<b>e</b>) for X-band band radar and (<b>b</b>,<b>d</b>,<b>f</b>) for S-band band radar. The color of the bubble charts shows bias ratios (<math display="inline"><semantics> <mrow> <msub> <mi>R</mi> <mrow> <msub> <mi>K</mi> <mrow> <mi>DP</mi> </mrow> </msub> </mrow> </msub> </mrow> </semantics></math>/<span class="html-italic">R</span><sub>gauge</sub>) between the hourly accumulated rainfall estimates and independent gauge observations.</p>
Full article ">Figure 10
<p>The normalized density of <span class="html-italic">Z<sub>H</sub></span>-<span class="html-italic">Z</span><sub>DR.</sub> in the pre-summer rainband observed by the Huadu X-band PAR at (<b>a</b>) 10:48, (<b>b</b>) 10:50, (<b>c</b>) 10:52, and (<b>d</b>) the S-band radar at 10:54 on 7 June 2020.</p>
Full article ">Figure 11
<p>The normalized density of <span class="html-italic">Z<sub>H</sub></span>-<span class="html-italic">Z</span><sub>DR.</sub> in the Typhoon Higos rainband observed by the Nansha X-band PAR at (<b>a</b>) 13:00, (<b>b</b>) 13:02, (<b>c</b>) 13:04, and (<b>d</b>) the S-band radar at 13:00 on 19 August 2020.</p>
Full article ">Figure 12
<p>The normalized density of <span class="html-italic">Z<sub>H</sub></span>-<span class="html-italic">Z</span><sub>DR</sub> in the local heavy precipitation rainband observed by the Panyu X-band PAR at (<b>a</b>) 17:30, (<b>b</b>) 17:31, (<b>c</b>) 17:33, and (<b>d</b>) the S-band radar at 17:30 on 4 September 2020.</p>
Full article ">
14 pages, 4499 KiB  
Article
Circularly Polarized Antenna Array with Decoupled Quad Vortex Beams
by Shuo Xu, He-Xiu Xu, Yanzhao Wang, Jian Xu, Chaohui Wang, Zhichao Pang and Huiling Luo
Nanomaterials 2022, 12(17), 3083; https://doi.org/10.3390/nano12173083 - 5 Sep 2022
Cited by 3 | Viewed by 1955
Abstract
Achieving multiple vortex beams with different modes in a planar microstrip array is pivotal, yet still extremely challenging. Here, a hybrid method combining both Pancharatnam−Berry (PB) phase that is induced by the rotation phase and excitation phase of a feeding line has been [...] Read more.
Achieving multiple vortex beams with different modes in a planar microstrip array is pivotal, yet still extremely challenging. Here, a hybrid method combining both Pancharatnam−Berry (PB) phase that is induced by the rotation phase and excitation phase of a feeding line has been proposed for decoupling two orthogonal circularly polarized vortex beams. Theoretical analysis is derived for array design to generate quad vortex beams with different directions and an arbitrary number of topological charges. On this basis, two 8 × 8 planar arrays were theoretically designed in an X band, which are with topological charges of l1 = −1, l2 = 1, l3 = −1, and l4 = 1 in Case I and topological charges of l1 = −1, l2 = 1, l3 = −1, and l4 = 1 in Case II. To further verify the above theory, the planar array in Case I is fabricated and analyzed experimentally. Dual-LP beams are realized by using rectangular patch elements with two orthogonally distributed feeding networks on different layers based on two types of feeding: proximity coupling and aperture coupling. Both the numerical simulation and experimental measurement results are in good agreement and showcase the corresponding quad-vortex-beam characteristics within 8~12 GHz. The array achieves a measured S11 < −10 dB and S22 < −10 dB bandwidth of more than 33.4% and 29.2%, respectively. In addition, the isolation between two ports is better than −28 dB. Our strategy provides a promising way to achieve large capacity and high integration, which is of great benefit to wireless and radar communication systems. Full article
(This article belongs to the Special Issue Metasurfaces for Photonic Devices: Theory and Applications)
Show Figures

Figure 1

Figure 1
<p>Schematic of planar arrays that are composed of dual LP radiating sub-elements.</p>
Full article ">Figure 2
<p>Geometry of the proposed dual-LP element. (<b>a</b>) Schematic diagram. (<b>b</b>) Layout and parametric illustration. The optimized dimensions are: <span class="html-italic">a</span> = 8.15 mm, <span class="html-italic">b</span> = 6 mm, <span class="html-italic">w</span><sub>1</sub> = 0.56 mm, <span class="html-italic">w</span><sub>2</sub> = 0.67 mm, <span class="html-italic">l</span><sub>1</sub> = 5.5 mm, <span class="html-italic">l</span><sub>2</sub> = 3.9 mm, <span class="html-italic">l</span><sub>3</sub> = 1.8 mm, <span class="html-italic">l</span><sub>4</sub> = 5.7 mm, <span class="html-italic">h</span><sub>1</sub> = 0.8 mm, <span class="html-italic">h</span><sub>2</sub> = 0.8 mm, and <span class="html-italic">h</span><sub>3</sub> = 0.5 mm.</p>
Full article ">Figure 3
<p>Numerically simulated S-parameters of the dual-LP patch element.</p>
Full article ">Figure 4
<p>Normalized co-polarization and the cross-polarization patterns for LP element when the (<b>a</b>) sub-element A and (<b>b</b>) sub-element B working in <span class="html-italic">ϕ</span> = 72°, 150°, 252°, 330°.</p>
Full article ">Figure 5
<p>Rotation angles and excitation phases of sub-elements A and B for (<b>a</b>) quad CP beams without OAM and (<b>b</b>,<b>c</b>) quad-vortex-beam with (<b>b</b>) <span class="html-italic">l<sub>(A</sub></span><sub>,</sub> <span class="html-italic"><sub>L)</sub></span>= <span class="html-italic">l<sub>(B</sub></span><sub>,</sub> <span class="html-italic"><sub>L)</sub></span>= −1, <span class="html-italic">l<sub>(A</sub></span><sub>,</sub> <span class="html-italic"><sub>R)</sub></span>= <span class="html-italic">l<sub>(B</sub></span><sub>,</sub> <span class="html-italic"><sub>R)</sub></span>= 1 in Case I and (<b>c</b>) <span class="html-italic">l<sub>(A</sub></span><sub>,</sub> <span class="html-italic"><sub>L)</sub></span>= −1, <span class="html-italic">l<sub>(B</sub></span><sub>,</sub> <span class="html-italic"><sub>L)</sub></span>= 0, <span class="html-italic">l<sub>(A</sub></span><sub>,</sub> <span class="html-italic"><sub>R)</sub></span>= 2, <span class="html-italic">l<sub>(B</sub></span><sub>,</sub> <span class="html-italic"><sub>R)</sub></span>= 1. in case II (<b>d</b>,<b>e</b>) Simulated far-field patterns (<b>d</b>) in Case I and (<b>e</b>) in case II.</p>
Full article ">Figure 6
<p>Layouts of the (<b>a</b>) upper and (<b>b</b>) bottom feeding network for the quad-vortex-beam array in Case I.</p>
Full article ">Figure 7
<p>3-D view of the planar microstrip quad-vortex patch array in Case I.</p>
Full article ">Figure 8
<p>Simulated 3-D radiation pattern at 10 GHz for the quad-vortex-beam in Case I with different OAM modes (<span class="html-italic">l<sub>(A</sub></span><sub>,</sub> <span class="html-italic"><sub>L)</sub></span>= <span class="html-italic">l<sub>(B</sub></span><sub>,</sub> <span class="html-italic"><sub>L)</sub></span>= −1, <span class="html-italic">l<sub>(A</sub></span><sub>,</sub> <span class="html-italic"><sub>R)</sub></span>= <span class="html-italic">l<sub>(B</sub></span><sub>,</sub> <span class="html-italic"><sub>R)</sub></span>= 1) at predicted deflection angles of (<span class="html-italic">θ</span>, <span class="html-italic">ϕ</span>) = (30°, 72°), (30°, 150°), (30°, 252°), and (30°, 330°).</p>
Full article ">Figure 9
<p>Photographs of the final fabricated quad-vortex-beam planar microstrip array. (<b>a</b>) Assembled model, (<b>b</b>) radiating elements, (<b>c</b>) upper feeding network, (<b>d</b>) ground plane with H slot, and (<b>e</b>) bottom feeding network.</p>
Full article ">Figure 10
<p>(<b>a</b>) Simulated and measured S-parameters and (<b>b</b>) simulated AR of the microstrip quad-vortex-beam planar array.</p>
Full article ">Figure 11
<p>Experimental setup for (<b>a</b>) far-field and (<b>b</b>) near-field measurement.</p>
Full article ">Figure 12
<p>Numerical far-field patterns of the LHCP and RHCP beams in the <span class="html-italic">ϕ</span><sub>(<span class="html-italic">A</span>,</sub> <span class="html-italic"><sub>L</sub></span><sub>)</sub> = 72°; <span class="html-italic">ϕ</span><sub>(<span class="html-italic">B</span>,</sub> <span class="html-italic"><sub>L</sub></span><sub>)</sub> = 150°; <span class="html-italic">ϕ</span><sub>(<span class="html-italic">B</span>,</sub> <span class="html-italic"><sub>R</sub></span><sub>)</sub> = 252°; and <span class="html-italic">ϕ</span><sub>(<span class="html-italic">A</span>,</sub> <span class="html-italic"><sub>R</sub></span><sub>)</sub> = 330° planes at 10 GHz.</p>
Full article ">Figure 13
<p>The measured magnitude and phase distributions at 10 GHz for four beams of different OAM modes at different directions. (<b>a</b>) <span class="html-italic">l</span><sub>(<span class="html-italic">A</span>,</sub> <span class="html-italic"><sub>L</sub></span><sub>)</sub> = 1 at <span class="html-italic">θ</span><sub>(<span class="html-italic">A</span>,</sub> <span class="html-italic"><sub>L</sub></span><sub>)</sub> = 30°, <span class="html-italic">ϕ</span><sub>(<span class="html-italic">A</span>,</sub> <span class="html-italic"><sub>L</sub></span><sub>)</sub> = 72°. (<b>b</b>) <span class="html-italic">l</span><sub>(<span class="html-italic">B</span>,</sub> <span class="html-italic"><sub>L</sub></span><sub>)</sub> = 1 at <span class="html-italic">θ</span><sub>(<span class="html-italic">B</span>,</sub> <span class="html-italic"><sub>L</sub></span><sub>)</sub> = 30°, <span class="html-italic">ϕ</span><sub>(<span class="html-italic">B</span>,</sub> <span class="html-italic"><sub>L</sub></span><sub>)</sub> = 150°. (<b>c</b>) <span class="html-italic">l</span><sub>(<span class="html-italic">B</span>,</sub> <span class="html-italic"><sub>R</sub></span><sub>)</sub> = −1 at <span class="html-italic">θ</span><sub>(<span class="html-italic">B</span>,</sub> <span class="html-italic"><sub>R</sub></span><sub>)</sub> = 30°, <span class="html-italic">ϕ</span><sub>(<span class="html-italic">B</span>,</sub> <span class="html-italic"><sub>R</sub></span><sub>)</sub> = 252°. (<b>d</b>) <span class="html-italic">l</span><sub>(<span class="html-italic">A</span>,</sub> <span class="html-italic"><sub>R</sub></span><sub>)</sub> = −1 at <span class="html-italic">θ</span><sub>(<span class="html-italic">A</span>,</sub> <span class="html-italic"><sub>R</sub></span><sub>)</sub> = 30°, <span class="html-italic">ϕ</span><sub>(<span class="html-italic">A</span>,</sub> <span class="html-italic"><sub>R</sub></span><sub>)</sub> = 330°.</p>
Full article ">
13 pages, 3974 KiB  
Article
Complementary Multi-Band Dual Polarization Conversion Metasurface and Its RCS Reduction Application
by Fengan Li and Baiqiang You
Electronics 2022, 11(10), 1645; https://doi.org/10.3390/electronics11101645 - 21 May 2022
Cited by 8 | Viewed by 2497
Abstract
In this paper, we present a metasurface composed of complementary units that can realize orthogonal linear and linear-to-circular polarization conversion in multi-band. Linear polarization conversion has seven high-conversion frequency bands: 9.1–9.7 GHz, 15.6–17.6 GHz, 19.4–19.7 GHz, 21.2–23.1 GHz, 23.5–23.8 GHz, 26.2 GHz, and [...] Read more.
In this paper, we present a metasurface composed of complementary units that can realize orthogonal linear and linear-to-circular polarization conversion in multi-band. Linear polarization conversion has seven high-conversion frequency bands: 9.1–9.7 GHz, 15.6–17.6 GHz, 19.4–19.7 GHz, 21.2–23.1 GHz, 23.5–23.8 GHz, 26.2 GHz, and 27.9 GHz. Linear-to-circular polarization conversion also has seven frequency bands with axial ratios (ARs) less than 3 dB: 8.9–9.0 GHz, 9.9–14.7 GHz, 19.1–19.3 GHz, 23.2–23.35 GHz, 23.4 GHz, 24.1–25.4 GHz, and 27.2–27.8 GHz, with the generation of multiple bands extended by the combination of complementary units. Then, we utilize the combined polarization conversion unit’s mirror placement to form a 4 × 4 array to realize the phase difference cancellation of the reflective field, giving the metasurface the radar cross section (RCS) reduction function and the dual-band 10-dB monostatic RCS reduction bandwidth: 8.9–9.7 GHz and 15.5–26.1 GHz. The measured and simulated results were essentially identical. Because the design uses the complementary units to form an array to expand the polarization conversion frequency bands, it provides a novel idea for future designs and can be applied to multiple microwave frequency bands. Full article
(This article belongs to the Topic Antennas)
Show Figures

Figure 1

Figure 1
<p>Polarization conversion unit design. (<b>a</b>) Complementary unit. (<b>b</b>) Structure of unit 1. (<b>c</b>) Structure of unit 2.</p>
Full article ">Figure 2
<p>Simulation results of unit. (<b>a</b>) Reflection coefficient of unit 1. (<b>b</b>) Reflection coefficient of unit 2. (<b>c</b>) PCR for units 1 and 2. (<b>d</b>) AR for units 1 and 2.</p>
Full article ">Figure 3
<p>Simulation results of combined unit. (<b>a</b>) Reflection coefficient. (<b>b</b>) PCR. (<b>c</b>) Phase difference. (<b>d</b>) AR.</p>
Full article ">Figure 4
<p>Polarization conversion at oblique incidence: (<b>a</b>) PCR and (<b>b</b>) AR.</p>
Full article ">Figure 5
<p>Schematic diagram of PCM.</p>
Full article ">Figure 6
<p>RCS reduction metasurface-based PCM.</p>
Full article ">Figure 7
<p>(<b>a</b>) Monostatic RCS. (<b>b</b>) RCS reduction.</p>
Full article ">Figure 8
<p>3D scattering field of metal plane and RCS plane at 9.2 GHz and 16.5 GHz.</p>
Full article ">Figure 9
<p>Bistatic RCS reduction: (<b>a</b>) f = 9.2 GHz and (<b>b</b>) f = 16.5 GHz.</p>
Full article ">Figure 10
<p>Physical fabrication: (<b>a</b>) PCM and (<b>b</b>) RCS plane.</p>
Full article ">Figure 11
<p>Schematic diagram of measured device and environment.</p>
Full article ">Figure 12
<p>Comparison of measured and simulated results. (<b>a</b>) Measured and simulated reflection coefficients. (<b>b</b>) Measured and simulated RCS reduction.</p>
Full article ">
28 pages, 9175 KiB  
Article
Improved k-NN Mapping of Forest Attributes in Northern Canada Using Spaceborne L-Band SAR, Multispectral and LiDAR Data
by André Beaudoin, Ronald J. Hall, Guillermo Castilla, Michelle Filiatrault, Philippe Villemaire, Rob Skakun and Luc Guindon
Remote Sens. 2022, 14(5), 1181; https://doi.org/10.3390/rs14051181 - 27 Feb 2022
Cited by 13 | Viewed by 3160
Abstract
Satellite forest inventories are the only feasible way to map Canada’s vast, remote forest regions, such as those in the Northwest Territories (NWT). A method used to create such inventories is the k-nearest neighbour (k-NN) algorithm, which spatially extends information [...] Read more.
Satellite forest inventories are the only feasible way to map Canada’s vast, remote forest regions, such as those in the Northwest Territories (NWT). A method used to create such inventories is the k-nearest neighbour (k-NN) algorithm, which spatially extends information from forest inventory (FI) plots to the entire forest land base using wall-to-wall features typically derived from Landsat data. However, the benefits of integrating L-band synthetic aperture radar (SAR) data, strongly correlated to forest biomass, have not been assessed for Canadian northern boreal forests. Here we describe an optimized multivariate k-NN implementation of a 151,700 km2 area in southern NWT that included ca. 2007 Landsat and dual-polarized Phased Array type L-band SAR (PALSAR) data on board the Advanced Land Observing Satellite (ALOS). Five forest attributes were mapped at 30 m cells: stand height, crown closure, stand/total volume and aboveground biomass (AGB). We assessed accuracy gains compared to Landsat-based maps. To circumvent the scarcity of FI plots, we used 3600 footprints from the Geoscience Laser Altimeter System (GLAS) as surrogate FI plots, where forest attributes were estimated using Light Detection and Ranging (LiDAR) metrics as predictors. After optimization, k-NN predicted forest attribute values for each pixel as the average of the 4 nearest (k = 4) surrogate FI plots within the Euclidian space of 9 best features (selected among 6 PALSAR, 10 Landsat, and 6 environmental features). Accuracy comparisons were based on 31 National Forest Inventory ground plots and over 1 million airborne LiDAR plots. Maps that included PALSAR HV backscatter resulted in forest attribute predictions with higher goodness of fit (adj. R2), lower percent mean error (ME%), and percent root mean square error (RMSE%), and lower underestimation for larger attribute values. Predictions were most accurate for conifer stand height (RMSE% = 32.1%, adj. R2 = 0.58) and AGB (RMSE% = 47.8%, adj. R2 = 0.74), which is much more abundant in the area than mixedwood or broadleaf. Our study demonstrates that optimizing k-NN parameters and feature space, including PALSAR, Landsat, and environmental variables, is a viable approach for inventory mapping of the northern boreal forest regions of Canada. Full article
(This article belongs to the Section Forest Remote Sensing)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>MVI phase 1 study area (red outline) within a broader area across 2 provinces and 2 territories (separated by thin black outline) that has as a backdrop a ca. 2007 landcover map that includes forest cover types (C: conifer; M: mixedwood; B: broadleaf) with 3 density classes (sparse, open, dense) along with the Geoscience Laser Altimeter System (GLAS) reference dataset of surrogate forest inventory (FI) plots and 2 validation sample sets. The top right zoomed-in inset shows a single GLAS FI plot surrounded by BT−ALS plots in a 500 m by 500 m area corresponding to an intersection between the BT−ALS transect and an ICESat track. Map is in Albers equal area conic projection.</p>
Full article ">Figure 2
<p>Percent occurrence of forest cover types (WT: wetland treed; C: conifer; M: mixedwood; B: broadleaf) with cover density classes (sparse, open, dense) across all forested pixels of the study area and the initial and final GLAS samples of surrogate forest inventory FI plots, respectively, according to the ca. 2007 landcover map (<a href="#remotesensing-14-01181-f001" class="html-fig">Figure 1</a>).</p>
Full article ">Figure 3
<p><span class="html-italic">k</span>-NN optimization and mapping workflow to generate the Satellite Vegetation Inventory (SVI) raster maps of five forest attributes and SVI map comparative accuracy assessment using Landsat-based map version (SVI_L) and published (PUB) maps. Numbers in brackets refer to related sections in the article.</p>
Full article ">Figure 4
<p>Multivariate global root mean square difference (GRMSD) metric (left Y axis) and the number of times a particular feature was selected in univariate feature selection across five attributes (right Y axis) for (<b>a</b>) initial selection features among the 20 candidate features and (<b>b</b>) the adjusted final selection of nine features (marked as * in panel (<b>a</b>)).</p>
Full article ">Figure 5
<p>Percent statistic values relative to their optimal values (100%) across a range of <span class="html-italic">k</span> values (rel<sub>stat_opt</sub>) for pseudo-R<sup>2</sup> (T<sup>2</sup>, Equation (2)), root mean square difference (RMSD, Equation (3)), and mean difference using the lower and upper 5% of distribution (MD<sub>5</sub>, MD<sub>95,</sub> Equations (5) and (6)), supporting the selection of the optimal <span class="html-italic">k</span> value of 4.</p>
Full article ">Figure 6
<p>SVI raster maps from <span class="html-italic">k</span>-NN predictions of (<b>a</b>) stand height and (<b>b</b>) AGB for the Phase 1 area. White pixels are non-forested lands, whereas light blue pixels are water bodies. Low and high attribute values are the 5% and 95% percentile, respectively. SVI maps for all five forest attributes are found in <a href="#app1-remotesensing-14-01181" class="html-app">Figure S2</a>.</p>
Full article ">Figure 6 Cont.
<p>SVI raster maps from <span class="html-italic">k</span>-NN predictions of (<b>a</b>) stand height and (<b>b</b>) AGB for the Phase 1 area. White pixels are non-forested lands, whereas light blue pixels are water bodies. Low and high attribute values are the 5% and 95% percentile, respectively. SVI maps for all five forest attributes are found in <a href="#app1-remotesensing-14-01181" class="html-app">Figure S2</a>.</p>
Full article ">Figure 7
<p>Comparison of scatterplots of observations versus predictions of stand height (left column) and AGB (right column) from (<b>a</b>) SVI maps, (<b>b</b>) Landsat-based SVI_L maps and (<b>c</b>) previously published Landsat-based maps using all NFI plots (blue dots) and BT−ALS LiDAR plots (density scatterplot). Dashed blue and black lines are regression lines along with equations and adj. R<sup>2</sup> values, respectively, based on NFI plots and BT−ALS LiDAR plots (see <a href="#app1-remotesensing-14-01181" class="html-app">Table S3</a> for linear regression statistics of all attributes). Distinctive symbols for the NFI plots (blue dots) distinguish the three forest cover types.</p>
Full article ">Figure 8
<p>Goodness of fit (adj. R<sup>2</sup>) for (<b>a</b>) stand height and (<b>b</b>) AGB relative to the validation datasets comprising NFI plots and BT−ALS LiDAR plots using all samples (ALL) then partitioned by forest cover type (C: conifer, M: mixedwood, B: broadleaf) for SVI maps compared to Landsat-based SVI_L maps and published (PUB) maps.</p>
Full article ">Figure 9
<p>Percent mean error (ME%) for (<b>a</b>) stand height and (<b>b</b>) AGB relative to the validation datasets comprising NFI plots and BT−ALS LiDAR plots using all samples (ALL) then partitioned by forest cover type (C: conifer, M: mixedwood, B: broadleaf) for SVI maps compared to Landsat-based SVI_L maps and published (PUB) maps.</p>
Full article ">Figure 10
<p>Percent root mean square error (RMSE%) for (<b>a</b>) stand height and (<b>b</b>) AGB relative to the validation datasets comprising NFI plots and BT−ALS LiDAR plots using all samples (ALL) then partitioned by forest cover type (C: conifer, M: mixedwood, B: broadleaf) for SVI maps compared to Landsat-based SVI_L maps and published (PUB) maps.</p>
Full article ">Figure 11
<p>Plot of mean prediction error (predicted minus observed) for (<b>a</b>) stand height and (<b>b</b>) AGB using all NFI plots (horizontal lines) and NFI plots grouped by quartiles Q1 to Q4 (mean error ± one standard deviation, dots) for SVI maps, Landsat-based SVI_L maps and published (PUB) maps. Dotted lines are added to highlight trends across four quartiles.</p>
Full article ">
20 pages, 8230 KiB  
Article
Sequential 90° Rotation of Dual-Polarized Antenna Elements in Linear Phased Arrays with Improved Cross-Polarization Level for Airborne Synthetic Aperture Radar Applications
by Diego Lorente, Markus Limbach, Bernd Gabler, Héctor Esteban and Vicente E. Boria
Remote Sens. 2021, 13(8), 1430; https://doi.org/10.3390/rs13081430 - 8 Apr 2021
Cited by 4 | Viewed by 3694
Abstract
In this work, a novel rotation approach for the antenna elements of a linear phased array is presented. The proposed method improves by up to 14 dB the cross-polarization level within the main beam by performing a sequential 90° rotation of the identical [...] Read more.
In this work, a novel rotation approach for the antenna elements of a linear phased array is presented. The proposed method improves by up to 14 dB the cross-polarization level within the main beam by performing a sequential 90° rotation of the identical array elements, and achieving measured cross-polarization suppressions of 40 dB. This configuration is validated by means of simulation and measurements of a manufactured linear array of five dual-polarized cavity-box aperture coupled stacked patch antennas operating in L-Band, and considering both uniform amplitude and phase distribution and beamforming with amplitude tapering. The analysis is further extended by applying and comparing the proposed design with the 180° rotation and non-rotation topologies. This technique is expected to be used for the next generation L-Band Airborne Synthetic Aperture Radar Sensor of the German Aerospace Center (DLR). Full article
(This article belongs to the Special Issue New Technologies for Earth Remote Sensing)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>F-SAR Antenna carrier.</p>
Full article ">Figure 2
<p>Antenna multilayer structure.</p>
Full article ">Figure 3
<p>Antenna feeding layout and constructed prototype.</p>
Full article ">Figure 4
<p>Simulated and measured S-parameters.</p>
Full article ">Figure 5
<p>Simulated and measured radiation pattern of the single element @ <math display="inline"><semantics> <msub> <mi>f</mi> <mn>0</mn> </msub> </semantics></math> = 1.325 GHz.</p>
Full article ">Figure 6
<p>Picture of the antenna array in the DLR’s Compact Test Range.</p>
Full article ">Figure 7
<p>Manufactured normal array.</p>
Full article ">Figure 8
<p>Comparison simulation and measurement @ <math display="inline"><semantics> <msub> <mi>f</mi> <mn>0</mn> </msub> </semantics></math> = 1.325 GHz: co-polar patterns. Normal array. Uniform amplitude and phase distribution.</p>
Full article ">Figure 9
<p>Comparison simulation and measurement @ <math display="inline"><semantics> <msub> <mi>f</mi> <mn>0</mn> </msub> </semantics></math> = 1.325 GHz: co-polar &amp; cross-polar patterns (<math display="inline"><semantics> <mrow> <mi>φ</mi> <mo>=</mo> <msup> <mn>90</mn> <mo>∘</mo> </msup> </mrow> </semantics></math>). Normal array. Uniform amplitude and phase distribution.</p>
Full article ">Figure 10
<p>Comparison simulation and measurement @ <math display="inline"><semantics> <msub> <mi>f</mi> <mn>0</mn> </msub> </semantics></math> = 1.325 GHz: normalized co-polar &amp; cross-polar patterns (<math display="inline"><semantics> <mrow> <mi>φ</mi> <mo>=</mo> <msup> <mn>90</mn> <mo>∘</mo> </msup> </mrow> </semantics></math>). Normal array. Beamforming.</p>
Full article ">Figure 11
<p>Manufactured array performing a 180° rotation of antenna elements 2 and 4.</p>
Full article ">Figure 12
<p>Comparison normal array and array with 180° rotation @ <math display="inline"><semantics> <msub> <mi>f</mi> <mn>0</mn> </msub> </semantics></math> = 1.325 GHz: normalized co-polar and cross-polar patterns (<math display="inline"><semantics> <mrow> <mi>φ</mi> <mo>=</mo> <mspace width="3.33333pt"/> <msup> <mn>90</mn> <mo>∘</mo> </msup> </mrow> </semantics></math>). Uniform amplitude and phase distribution.</p>
Full article ">Figure 13
<p>Comparison normal array and array with 180° rotation @ <math display="inline"><semantics> <msub> <mi>f</mi> <mn>0</mn> </msub> </semantics></math> = 1.325 GHz: normalized co-polar &amp; cross-polar patterns (<math display="inline"><semantics> <mrow> <mi>φ</mi> <mo>=</mo> <mspace width="3.33333pt"/> <msup> <mn>90</mn> <mo>∘</mo> </msup> </mrow> </semantics></math>). Beamforming.</p>
Full article ">Figure 14
<p>Sequential 90° rotation of the antenna elements.</p>
Full article ">Figure 15
<p>Comparison normal array and array with sequential 90° rotation @ <math display="inline"><semantics> <msub> <mi>f</mi> <mn>0</mn> </msub> </semantics></math> = 1.325 GHz: normalized co-polar &amp; cross-polar patterns (<math display="inline"><semantics> <mrow> <mi>φ</mi> <mo>=</mo> <mspace width="3.33333pt"/> <msup> <mn>90</mn> <mo>∘</mo> </msup> </mrow> </semantics></math>). Uniform amplitude and phase distribution.</p>
Full article ">Figure 16
<p>Measured cross-polar suppression @ <math display="inline"><semantics> <msub> <mi>f</mi> <mn>0</mn> </msub> </semantics></math> = 1.325 GHz. Uniform amplitude and phase distribution. Horizontal polarization.</p>
Full article ">Figure 17
<p>Comparison normal array and array with sequential 90° rotation @ <math display="inline"><semantics> <msub> <mi>f</mi> <mn>0</mn> </msub> </semantics></math> = 1.325 GHz: normalized co-polar &amp; cross-polar patterns (<math display="inline"><semantics> <mrow> <mi>φ</mi> <mo>=</mo> <mspace width="3.33333pt"/> <msup> <mn>90</mn> <mo>∘</mo> </msup> </mrow> </semantics></math>). Beamforming.</p>
Full article ">Figure 18
<p>Measured cross-polar suppression @ <math display="inline"><semantics> <msub> <mi>f</mi> <mn>0</mn> </msub> </semantics></math> = 1.325 GHz. Beamforming. Vertical polarization.</p>
Full article ">
13 pages, 1274 KiB  
Article
A Low Cross-Polarization Configuration Method for Phased Array Radar Antenna
by Yongzhen Li, Zhanling Wang, Chen Pang and Xuesong Wang
Electronics 2020, 9(3), 396; https://doi.org/10.3390/electronics9030396 - 27 Feb 2020
Cited by 5 | Viewed by 4146
Abstract
The cross-polarization isolation (CPI) is a key parameter to assess the dual-polarization antenna because the cross-polarization closely affects the antenna application. A polarization state configuration (PSC) approach is proposed to configure the polarization state of the polarimetric phased array radar antenna. Unlike the [...] Read more.
The cross-polarization isolation (CPI) is a key parameter to assess the dual-polarization antenna because the cross-polarization closely affects the antenna application. A polarization state configuration (PSC) approach is proposed to configure the polarization state of the polarimetric phased array radar antenna. Unlike the traditional fixed polarization states such as the linear polarization (LP) and the circular polarization (CP), the PSC method modulates the polarization state of the radiated wave continuously. In addition, the optimal excitation magnitude and phase of the dual-polarization element is calculated, thereby maximizing the CPI. Most of the configured polarization state is the elliptical polarization (EP), and a lower cross-polarization level and higher CPI could be obtained. This method could expand the acceptable angle range when compared with the LP and CP waves. Numerical simulations and comparisons are conducted to manifest the validity of the proposed method. Full article
(This article belongs to the Special Issue Recent Advances in Array Antenna and Array Signal Processing)
Show Figures

Figure 1

Figure 1
<p>Spherical coordinate system for electric fields radiating from dual-polarization element.</p>
Full article ">Figure 2
<p>Copolar and cross-polar electric fields for horizontally and vertically polarized dipoles. (<b>a</b>) co-polarization pattern of a horizontally polarized dipole; (<b>b</b>) cross-polarization pattern of a vertically polarized dipole; (<b>c</b>) cross-polarization pattern of a horizontally polarized dipole; (<b>d</b>) co-polarization pattern of a vertically polarized dipole.</p>
Full article ">Figure 3
<p>Copolar and cross-polar electric fields for LHCP and RHCP dipoles. (<b>a</b>) co-polarization pattern of LHCP dipoles; (<b>b</b>) cross-polarization pattern of RHCP dipoles; (<b>c</b>) cross-polarization pattern of LHCP dipoles; (<b>d</b>) co-polarization pattern of RHCP dipoles.</p>
Full article ">Figure 4
<p>Copolar and cross-polar electric fields for LHEP and RHEP dipoles. (<b>a</b>) co-polarization pattern of LHEP dipoles; (<b>b</b>) cross-polarization pattern of RHEP dipoles; (<b>c</b>) cross-polarization pattern of LHEP dipoles; (<b>d</b>) co-polarization pattern of RHEP dipoles.</p>
Full article ">Figure 5
<p>The available angle range of the polarized wave when the CPI is greater than 40 dB. (<b>a</b>) LP wave; (<b>b</b>) CP wave; (<b>c</b>) EP wave.</p>
Full article ">Figure 6
<p>Polarization states distribution on the Poincaré sphere in different directions. (<b>a</b>) along with the azimuth angle when <math display="inline"><semantics> <mrow> <mi>θ</mi> <mo>=</mo> <mn>60</mn> <mo>°</mo> </mrow> </semantics></math>; (<b>b</b>) along with the elevation angle <math display="inline"><semantics> <mrow> <mi>ϕ</mi> <mo>=</mo> <mo>−</mo> <mn>30</mn> <mo>°</mo> </mrow> </semantics></math>.</p>
Full article ">Figure 7
<p>The bias of the linear depolarization ratio versus beam scanning angle. (<b>a</b>) LP wave; (<b>b</b>) CP wave; (<b>c</b>) EP wave.</p>
Full article ">Figure 8
<p>The bias of the differential reflectivity versus beam scanning angle. (<b>a</b>) LP wave; (<b>b</b>) CP wave; (<b>c</b>) EP wave.</p>
Full article ">
26 pages, 14735 KiB  
Article
Monitoring and Detecting Archaeological Features with Multi-Frequency Polarimetric Analysis
by Jolanda Patruno, Magdalena Fitrzyk and Jose Manuel Delgado Blasco
Remote Sens. 2020, 12(1), 1; https://doi.org/10.3390/rs12010001 - 18 Dec 2019
Cited by 11 | Viewed by 4567
Abstract
In remote sensing for archaeology, an unequivocal method capable of automatic detection of archaeological features still does not exists. Applications of Synthetic Aperture Radar (SAR) remote sensing for archaeology mainly focus on high spatial resolution SAR sensors, which allow the recognition of structures [...] Read more.
In remote sensing for archaeology, an unequivocal method capable of automatic detection of archaeological features still does not exists. Applications of Synthetic Aperture Radar (SAR) remote sensing for archaeology mainly focus on high spatial resolution SAR sensors, which allow the recognition of structures of small dimension and give information of the surface topography of sites. In this study we investigated the potential of combined dual and fully polarized SAR data and performed polarimetric multi-frequency and multi-incidence angle analysis of C-band Sentinel-1, L-band Advanced Land Observing Satellite Phased Array type L-band Synthetic Aperture Radar (ALOS PALSAR) and of C-band Radar Satellite-2 (RADARSAT-2) datasets for the detection of surface and subsurface archaeological structures over the United Nations Educational, Scientific and Cultural Organization (UNESCO) site of Gebel Barkal (Sudan). While PALSAR offers a good historical reference, Sentinel-1 time series provide recent and systematic monitoring opportunities. RADARSAT-2 polarimetric data have been specifically acquired in 2012/2013, and have been scheduled to achieve a multi-temporal observation of the archaeological area under study. This work demonstrated how to exploit a complex but significant dataset composed of SAR full polarimetric and dual polarimetric acquisitions, with the purpose of identifying the most suitable earth observation technique for the preservation and identification of archaeological features. The scientific potential of the illustrated analysis fits perfectly with the current delicate needs of cultural heritage; such analysis demonstrates how multi-temporal and multi-data cultural heritage monitoring can be applied not only for documentation purposes, but can be addressed especially to those areas exposed to threats of different nature that require a constant and prompt intervention plans. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Gebel Barkal archaeological site © Google Earth.</p>
Full article ">Figure 2
<p>Gebel Barkal and the sites of the Napatan region © UNESCO.</p>
Full article ">Figure 3
<p>The Holy Mountain and the temple of Amun (B1500). © Max Farrar.</p>
Full article ">Figure 4
<p>Pyramids at Gebel Barkal (Courtesy of L. Perotti).</p>
Full article ">Figure 5
<p>Archaeological area of Gebel Barkal, Sudan: wide view (<b>a</b>) and arcaheological area (<b>b</b>).</p>
Full article ">Figure 6
<p>Erosion pebble conglomerate © Max Farrar.</p>
Full article ">Figure 7
<p>Left image illustrates the Pauli decomposition RGB image: |S<sub>hh</sub> + S<sub>vv</sub>| |S<sub>hh</sub> − S<sub>vv</sub>| |S<sub>hv</sub>|. Right image shows the Google Earth image available for 2006.</p>
Full article ">Figure 8
<p>Model-based Freeman decomposition. <b>f<sub>SB</sub></b>, f<sub>DB</sub> and <b>f</b><sub>vs</sub> represent Freeman-Durden components for single bounce, double bounce and volume scattering respectively, represented in blue, red and green colors. Plane is contained in normal direction to the line-of-sight of the satellite.</p>
Full article ">Figure 9
<p>Model-based Yamaguchi 4 component decomposition.</p>
Full article ">Figure 10
<p>Yamaguchi Y4O (top), Yamaguchi Y4R (middle) and Yamaguchi G4U1 (bottom) decomposition RGB images.</p>
Full article ">Figure 11
<p>Diagram listing pre-processing steps of ground range detected (GRD) scenes in Google Earth Engine (GEE).</p>
Full article ">Figure 12
<p>ALOS PALSAR Pauli RGB (left: 2006; rigth: 2009) and KOMPSAT-2 overlay.</p>
Full article ">Figure 13
<p>ALOS PALSAR Yamaguchi RGB images ((<b>a</b>): 2006; (<b>b</b>): 2009). Top, middle and bottom figures correspond to Y4O, Y4R and GU41 Yamaguchi decomposition respectively.</p>
Full article ">Figure 14
<p>G4U1 decomposition single channels 2006 (<b>a</b>) and 2009 (<b>b</b>) acquisitions. Double bounce (top); single bounce (middle); volume scattering (bottom).</p>
Full article ">Figure 15
<p>Pauli RGB decomposition ((<b>a</b>): 04/2012, (<b>b</b>): 11/2012, (<b>c</b>): 01/2013, (<b>d</b>): 07/2013) overlaid with Google Earth image (11/2012).</p>
Full article ">Figure 16
<p>Yamaguchi G4U1 decomposition RGB image ((<b>a</b>): 2012/04/28, (<b>b</b>): 2012/11/06, (<b>c</b>): 2013/01/17, (<b>d</b>): 2013/07/07).</p>
Full article ">Figure 17
<p>Incident wave on the pyramids.</p>
Full article ">Figure 18
<p>Yamaguchi G4U1 decomposition: double bounce (top), single bounce (middle) and volume scattering (bottom) for each acquisition date: 04/2012 (<b>a</b>), 11/2012 (<b>b</b>), 01/2013 (<b>c</b>), 07/2013 (<b>d</b>).</p>
Full article ">Figure 18 Cont.
<p>Yamaguchi G4U1 decomposition: double bounce (top), single bounce (middle) and volume scattering (bottom) for each acquisition date: 04/2012 (<b>a</b>), 11/2012 (<b>b</b>), 01/2013 (<b>c</b>), 07/2013 (<b>d</b>).</p>
Full article ">Figure 19
<p>Regions of interest (ROIs) in the archaeological area of Jebel Barkal.</p>
Full article ">Figure 20
<p>Radar backscatter values (<math display="inline"> <semantics> <mrow> <msub> <mi mathvariant="bold">σ</mi> <mn mathvariant="bold">0</mn> </msub> </mrow> </semantics> </math> in VV) evolution for ROIs.</p>
Full article ">Figure 21
<p>Radar backscatter values (<math display="inline"> <semantics> <mrow> <msub> <mi>σ</mi> <mn>0</mn> </msub> <mrow> <mtext> </mtext> <mi>in</mi> <mtext> </mtext> <mi>VH</mi> </mrow> </mrow> </semantics> </math>) evolution for ROIs.</p>
Full article ">Figure 22
<p>Layout performed in the GIS project with the representation of excavations areas and the anomaly.</p>
Full article ">Figure 23
<p>Gebel Barkal investigated area. Courtesy of E. Ciampini (Dec, 2013).</p>
Full article ">Figure 24
<p>Yamaguchi decomposition. High amplitude values recorded for the anomaly (white arrow).</p>
Full article ">
Back to TopTop