Post-Flood UAV-Based Free Space Optics Recovery Communications with Spatial Mode Diversity

: The deployment of unmanned aerial vehicles (UAVs) for free space optical communications is an attractive solution for forwarding the vital health information of victims from a ﬂood-stricken area to neighboring ground base stations during rescue operations. A critical challenge to this is maintaining an acceptable signal quality between the ground base station and UAV-based free space optics relay. This is largely unattainable due to rapid UAV propeller and body movements, which result in ﬂuctuations in the beam alignment and frequent link failures. To address this issue, linearly polarized Laguerre–Gaussian modes were leveraged for spatial mode diversity to prevent link failures over a 400 m link. Spatial mode diversity successfully improved the bit error rate by 38% to 55%. This was due to a 10% to 19% increase in the predominant mode power from spatial mode diversity. The time-varying channel matrix indicated the presence of nonlinear deterministic chaos. This opens up new possibilities for research on state-space reconstruction of the channel matrix. capacity and wide coverage. This requires the design of satellite constellations comprising multiple satellites across several countries for reliable coverage in the ﬂood-stricken area, and to backhaul unﬂooded cells via satellite [12–14]. High-altitude platforms (HAPs) such as airships, aerostats, and balloons have been used to ferry antennas in the affected area, for relaying data to neighboring ground base stations to a large number of users [15]. HAPs are located in the stratosphere, enabling a wide coverage and accommodating a large number of users [16]. Compared to satellites, HAPs provide a higher area throughput and resource utilization [15]. HAPs can also distribute the recording of the orbital paths of satellites and monitor the probability of a collision between satellites [17]. Data transmission from a HAP involves the control of the HAP ﬂight trajectory and the transmit power of the HAP antenna [18–20]. diversity. Under heavy rain, the power in the dominant mode was approximately 56% to 63% without mode diversity reception and increased to 71% to 75% with mode diversity reception. Uplink and downlink data transmissions demonstrated similar power coupling coefﬁcient characteristics, whereas uplink and downlink data transmissions demonstrated similar power coupling coefﬁcient characteristics.


Introduction
The risk of floods is rising worldwide, due to an increase in the intensity and frequency of rainfalls as a consequence of global warming [1]. Floods cause fatalities, damage to buildings, deterioration of health conditions, severe economic losses, losses of livelihood, and disruption of global trade [2]. Globally, it is estimated that floods have directly affected 2.3 billion people and caused USD 662 billion in damages between 1995 and 2015 [3]. In Malaysia, destructive floods occur frequently in the three eastern states of the peninsular during the seasonal monsoon between October and March, affecting more than 4.8 million people annually [4]. In the recent December 2021 massive nationwide flood, as many as eight states were struck by a 1-in-a-100-year heavy rainfall spanning two weeks, which displaced thousands of residents and strained emergency services [5][6][7]. The aftermath of the floods also witnessed unscheduled water cuts and disruptions to the electricity supply [8,9]. In addition, telecommunication base stations were severely damaged, leading to the disruption of communication services [10,11]. This prevented the exchange of situational awareness and hampered the coordination of search and rescue operations.
To connect emergency responders to flood rescue centers and to victims, several recovery communications technologies are being explored, as shown in Figure 1. To compensate for flood-impaired ground radio base stations, satellites have been considered due its large To connect emergency responders to flood rescue centers and to victims, several recovery communications technologies are being explored, as shown in Figure 1. To compensate for flood-impaired ground radio base stations, satellites have been considered due its large capacity and wide coverage. This requires the design of satellite constellations comprising multiple satellites across several countries for reliable coverage in the flood-stricken area, and to backhaul unflooded cells via satellite [12][13][14]. High-altitude platforms (HAPs) such as airships, aerostats, and balloons have been used to ferry antennas in the affected area, for relaying data to neighboring ground base stations to a large number of users [15]. HAPs are located in the stratosphere, enabling a wide coverage and accommodating a large number of users [16]. Compared to satellites, HAPs provide a higher area throughput and resource utilization [15]. HAPs can also distribute the recording of the orbital paths of satellites and monitor the probability of a collision between satellites [17]. Data transmission from a HAP involves the control of the HAP flight trajectory and the transmit power of the HAP antenna [18][19][20].   Recently, the deployment of unmanned aerial vehicles (UAVs) during natural disasters has become increasingly prevalent for forwarding information from the affected area to undamaged ground base stations during rescue operations [21]. UAVs generally have limited onboard energy [22]. Typically, several UAVs are connected in topologies such as a mesh, star, or cluster. UAVs are linked to each other in order to reach the undamaged ground base station, in conjunction with routing and positioning algorithms [23][24][25][26].
The remainder of the paper organized as follows. Section 2 provides an overview of free space optics (FSO) for post-flood recovery communications and related work on spatial mode diversity using FSO. Section 3 elucidates the novelty and contributions of the paper. Section 4 describes the design of a UAV-based spatial mode diversity FSO system for post-flood recovery communications. Section 5 reports on the performance analysis of the proposed system.

Free Space Optics for Post-Flood Recovery Communications
Free space optics is an attractive solution during floods when conventional base stations spanning several cities have been damaged, providing connectivity to the floodstricken area and ensuring minimal disruption in unaffected areas. Free space optics provides a high data bandwidth and rapid deployment, without spectrum licensing costs [27,28]. In addition, the light signals can penetrate easily through water droplets [29].
Although satellites and HAPs cover a larger area, the deployment time is longer [30]. Motivated by the rapid deployment and agility of UAVs for recovery communications [22,30], a UAV-based FSO is proposed as a wireless access technology for rapid recovery communications during floods.
Several UAV-based FSO strategies have previously been developed for recovery communications between base stations for recovery communications, focusing on downlink scheduling [31], and the placement of FSO transceivers [32,33], divergence angle [34], or hybrid radio frequency (RF)/FSO links [35,36]. An unveiled challenge to UAV-based FSO systems is to maintain a high SNR between the base station and UAV-based FSO relay. This is largely unattainable due to minute but rapid UAV propeller and body movements, which result in swift fluctuations in the beam alignment and water scattering from flood waters. Consequently, this would lead to frequent link failures. To address this issue, spatial mode diversity is proposed for improving the resilience of the system from UAV movements.
Harnessing spatial modes as independent information carriers using space division multiplexing (SDM) techniques has recently gained traction for tackling the impending data capacity crunch [37]. Spatial modes provide an additional degree of freedom in wireless communications. Characteristically, in a SDM system, independent data-carrying beams are structured on distinct spatial modes that can be multiplexed at the transmitter aperture for co-transmission and demultiplexed at the receiver aperture, with minimal interference [37].
To increase the resilience of an FSO system, multiplexing in the amplitude, frequency, polarization, and time domains have been employed. Recently, a new degree of freedom based on spatial modes is being explored. The simultaneous transmission of several spatial modes can be realized without the separation between apertures, thus reducing the device footprint. Individual spatial modes encounter different refractive index perturbations from atmospheric turbulence, despite propagating in the same path [38,39]. Spatial modes may be leveraged for transmission of independent data streams to improve link reliability and prevent network interruption.
Several approaches have been demonstrated for spatial mode diversity in FSO systems. In [40], spatial mode diversity was employed in conjunction with multiple transmitterreceiver aperture pairs to improve FSO link reliability under atmospheric turbulence. Three aperture pairs were used and each aperture pair utilized a Gaussian beam and an OAM beam in the uplink and downlink direction simultaneously for carrying the same data stream. In [41], spatial mode diversity from a three-mode photonic lantern coupling FSO receiver was designed in conjunction with a digital maximal ratio combining to enhance the worst signal-to-noise ratio (SNR) by more than 10 dB and to mitigate the interruption probability under atmospheric turbulence. In [42], the three-mode diversity reception in free space optical links from low-Earth orbital satellites to a ground terminal was demonstrated by mode coupling to a few-mode fiber for tracking power variations under different elevation angles. In [43], to reduce the outage probability, aperture diversity was realized through several transmitter apertures transmitting the same data stream on fundamental Gaussian beams whilst mode diversity was realized through multiple receiver apertures that decompose the incoming beam to several orbital angular momentum (OAM) modes, in conjunction with digital signal processing. In [44], OAM modes were used for spatial mode diversity through a 2 km FSO link to improve the link resilience and system capacity. The overall channel capacity was maximized by selecting optimal OAM mode numbers at each value of the SNR and turbulence strength. Spatial mode diversity has also been demonstrated using photonic crystal fibers [45][46][47][48]. Another approach for spatial mode diversity was demonstrated using Hermite-Gaussian (HG) and Laguerre-Gaussian (LG) modes of identical size [39]. This was shown to improve the bit error rate by up to 54% without an increase in the total transmit power or radius of the receive aperture. Unmanned aerial vehicle (UAV)-based FSO strategies have been explored for improving coverage between base stations for emergency communications during natural disasters, focusing on downlink scheduling [31] and the placement of FSO transceivers [32,33] and hybrid RF/FSO links [35,36].

Contributions
The disastrous aftermath of floods has reinforced the need for complementary access networks for critical recovery communications from the flood-stricken area to unaffected areas [21]. In the event where many conventional base stations spanning several cities are damaged, unlicensed, high-bandwidth wireless access technology is required for complementary post-flood recovery communications, to expedite communications to the floodstricken area. Connecting first responders and healthcare workers to flood victims and hospitals is an integral part of recovery communications. Several UAV-based FSO strategies have previously been developed for recovery communications between base stations for recovery communications, focusing on downlink scheduling [31], the placement of FSO transceivers [32,33], divergence angle [34], or hybrid RF/FSO links [35,36]. An unveiled challenge to UAV-based FSO systems is to maintain a high SNR between the base station and UAV-based FSO relay. This is largely unattainable due to minute but rapid UAV propeller and body movements, which result in swift fluctuations in the beam alignment and water scattering from flood waters. Consequently, this would lead to frequent link failures. To address this issue, spatial mode diversity is proposed using linearly polarized modes for improving the resilience of the system from UAV movements.

Methods
The proposed UAV-based spatial mode diversity FSO system for post-flood recovery communications is illustrated in Figure 1. Four vertical cavity surface-emitting lasers (VCSELs) generating continuous-wave optical signals on the fundamental modes at wavelengths of 850 nm, 880 nm, 910 nm, and 940 nm were structured into eight distinct linearly polarized Laguerre-Gaussian beams, LP lm, where l = 0, 1, 2, 3, 4 and m = 1, 3. The LG beams were intensity-modulated at 20 Gbps by pre-processed binary data collected from individual sensors on their respective channels. The LP modes and wavelengths used for the respective channels are shown in Table 1. The ground transmitter was enclosed in a case for protection from rain. LP lm modes were generated from the first diffraction order in the Fourier plane using the binarized electric field displayed on individual liquid crystal gratings, as follows. It is known that the Fourier transform of a linear translation is a complex phase shift [49]. Thus, to produce a translation in the Fourier plane, the transverse modal field is first multiplied by a complex phase shift: where τ x and τ y are linear tilt constants in the horizontal and vertical directions, respectively. The complex tilted field of the LP mode, a + jb, is binarized such that, in regions represented by an expanded laser beam is transmitted through a liquid crystal display, and in other regions, no light passes through. The binarized beam may be expressed as a Fourier series expansion [45][46][47][48]50]: where m o is the constant term and n is the n-th diffraction order. In the Fourier plane, the diffraction orders are translated linearly and separated: where x 2 and y 2 are spatial coordinates in the Fourier plane, * is the complex conjugate, and M n (x 2 , y 2 ) is the n-th diffraction order in the Fourier plane. Each binarized field was then Fourier transformed by a 400 cm focal-length achromatic convex lens, and the linearly polarized Laguerre-Gaussian mode, LP lm, was extracted from the first diffraction order in the Fourier plane, M1.
The transverse modal field distribution of a linearly polarized LP lm mode is expressed as [51]: where R is the normalized radius, φ is the azimuthal angle of the transmitted transverse modal field, V is the normalized frequency, L l m−1 is the generalized Laguerre polynomial whereby l is the azimuthal mode number, and m is the radial mode number.
For channel diversity, each sensor transmitted data on two LP modes so that they experience distinct refractive index fluctuations and modal power coupling. For any data stream, the first mode was used as the primary channel, and a second mode was used as the backup channel during adverse weather conditions, such as rain. Different combinations of the azimuthal mode number and radial mode numbers were used for various primary and backup channels, as shown in Table 1. The backup channels operated on different LP modes from the primary channels, thus experiencing distinct refractive index fluctuations and mode coupling.
The optical signals from all channels were then amplified, multiplexed, aligned, and transmitted to the UAV, which was flown at a height of 7 m above ground, in a straight line toward the ground receiver, located 40 m away from the transmitter. The horizontal distance travelled and height of the flying drone is shown in Figure 1. Following this, the optical beams are reflected from the UAV to the ground receiver. System parameters are provided in Table 2. The beam waist is the location along the propagation direction where the beam radius is minimum. The Rayleigh length is the distance from the beam waist, in the direction of beam propagation, where the beam radius is increased by a factor of the square root of 2. The beam divergence is an angular measure of the increase in beam diameter or radius with distance from the grating. The photodetector responsivity is a measure of optical-to-electrical conversion efficiency of a photodetector. The transverse electric field of the set of eight LP modes used for the transmission may be modeled as: where X i is the transverse electric field of the transmitted i-th mode and T is the transpose of the matrix. The estimated received electric field after propagating through the FSO channel is where Y i is estimated electric field of the i-th received mode and T is the transpose of the matrix. Y is related to the channel matrix, H by [52] where the channel matrix H contains the power coupling coefficients, A charge-coupled device was used to collect the intensity distribution of the received optical signals. The power coupling coefficients were computed for all channels using a modal decomposition method [53].
Multiplicative slow fading and inter-symbol interference from adjacent modes i = j were assumed. N is additive noise from the surroundings.
The received electric field of the i-th mode may be estimated as the summation of the individual components of the transmitted electric field, X j [54]: The power coupling coefficient h i,j is the overlap integral between the received transverse electric field, Y i , and transmitted transverse electric field, X j [55]: where A is the plane of the transverse electric field. When the set of modes propagate through the atmosphere, it encounters spatially and temporally varying refractive indices, due to random pressure temperature variations. Spatial modes individually experience different refractive index perturbations, even when propagating in the same path [39]. This causes unique random wavefront aberrations and power spreading for each mode into adjacent modes, evident by the degradation of power coupled into the original mode h i, j . Under atmospheric turbulence, these power fluctuations are slowly time-varying and vary significantly slower than the signal [56,57]. The Kolmogorov model for turbulent flow is the basis for many contemporary theories and models for emulating atmospheric turbulence for spatial modes [52,[58][59][60]. However, these models of turbulence typically only provide statistical averages for the random variations of the atmosphere. Thus, they may be considered insufficient to represent temporal intensity fluctuations and modal crosstalk for a UAV-based system. Hence, in our work, in order to demonstrate the performance of mode diversity for maintaining the signal quality for composite channels, the time-varying optical signal-to-interference-and-noise ratio (SINR) was evaluated under various weather conditions over a time interval of 30 min and sampling interval of 1 s. The optical SINR for the i-th channel is given by: where P i is the received power from the i-th desired mode (j = i), either in the primary channel or the backup channel. The measured interference power I i arises from crosstalk from undesirable spatial modes (j = i). The measured additive noise power, N, constituting thermal noise, ambient light from sunlight, and surrounding objects is minimal compared to the interfering modal power. For optimal signal detection, the ground receiver requires knowledge of the current channel matrix, H. The main challenge lies in the alignment between the ground transmitter to the UAV and the alignment between the UAV to the ground receiver, in order to prevent tip-tilt errors, despite a line-of-sight link between the ground transmitter and ground receiver. To address this issue at the ground transmitter, the conjugates of a known sequence of complex LP mode fields was transmitted toward the UAV and then reflected back to the ground transmitter to determine the optimum angle for pointing the optical beam from the ground transmitter toward the UAV. Adjustments to the pointing angle were made on the ground transmitter using a rotating shaft, translation stages, and goniometers until the exact LP mode was acquired on the reflected optical signal. Similarly, at the ground receiver, beam tracking of the signals was performed to avoid tip-tilt errors using a rotation shaft, translation stages, and goniometers. Using the conjugate of a known LP mode, the direction of arrival of the received optical beam was determined for maximizing spatial acquisition. After the optimal angles were determined, the pointing angle and angle of arrivals were fixed, as the UAV travelled in a straight line from the initial point. An alternative method for beam tracking is to measure the displacement offsets of the arriving at the detector, which translate into angles of arrival variations [61].
The SINR threshold SINRγ was set at the minimum level for turbulence-free transmission under heavy rain, such that: where ∑ 8 i=1 P i refers to the sum of the transmitted power in all desired modes and r i is the average receiver noise power per mode. The calculated SINRγ value was 25 dB.
As shown in Table 1, two different LP modes were used, for the primary channel and backup channels, under the same weather condition. The time-varying primary and backup channels were used for obtaining the composite channel, with the aid of a microelectromechanical systems (MEMS) optical switch for temporal switching between the two channels. The MEMS switch is controlled dynamically using an algorithm designed LabVIEW software by National Instruments (Austin, TX, USA), for automated switching between the two channels based on the comparison of the current SINR to the threshold SNR value, SINRγ. While the received SINR of the primary channel is equal to greater than SINRγ, the received data from the primary channel are forwarded to the composite channel by the MEMS switch. On the other hand, when the received SINR of the primary channel is less than SINRγ, the received data from the backup channel is forwarded to the composite channel by the MEMS switch.
A key parameter for indicating the strength of the atmospheric turbulence is the refractive index structure parameter, C 2 n , which was evaluated by transmitting the concerned mode and then calculating the scintillation index at the receiver using the intensity profile measurements from a charge-coupled-device camera. The scintillation index is given by the Rytov variance [62]: where I is the on-axis beam intensity and the angle brackets <> denote the ensemble average. The approach for measuring the scintillation index is similar to that in [61]. The refractive index structure parameter C 2 n is then evaluated by substituting the scintillation index into: where wave number k = 2π/λ and L is the link distance. Using Equation (13) to Equation (14), the C 2 n values were computed for various weather conditions based on the data transmission experiments using the designed FSO transceiver. The turbulence strength and turbulence fluctuation are indicated by C 2 n and the Rytov variance σ 2 , respectively. The values of C 2 n in the atmosphere are typically in the range from 10 −17 m −2/3 for weaker turbulence strength to 10 −13 m −2/3 for stronger turbulence strength. The value of σ 2 is typically related to the refractive index inhomogeneities C 2 n , optical wavelength λ, and the propagating distance L. A higher value of σ 2 generally represents stronger turbulence fluctuation [62]. Weak turbulence is associated with σ 2 < 1, and moderate fluctuation conditions are characterized by σ 2 ≈ 1. Strong fluctuations are associated with σ 2 > 1 [62].

Results and Discussion
To evaluate the effectiveness of spatial mode diversity using linearly polarized modes for maintaining the signal quality of a UAV-based free space optical system, the timevarying SINR of the primary, backup, and composite channels were measured under various weather conditions: (i) light rain, (ii) medium rain, (iii) heavy rain, (iv) clear day. The rain level was classified as 'light', 'medium', and 'heavy' for precipitation rates of accumulation of about 0 to 2 mm per h, 2 to 10 mm per h, and 10 to 50 mm per h, respectively [63]. For each weather condition, the minimum and maximum C 2 n values were computed to determine the range of turbulence strengths.
Although the signals were received at 20 Gbps, the acquisition and evaluation of the channel matrix and interfering modes at high data rates were computationally intensive and required time. Thus, the SINR was sampled at 1 Gbps. In addition, the hovering time of the UAV was 45 min, so a 30 min interval was selected to give the UAV sufficient time to return to the ground station.
For each case, the time-varying SINR is plotted for the primary channel, backup channel, and the composite channel. The composite channel comprises various temporal segments of data from the primary channel and the backup channel such that when the SINR is equal to or more than the threshold value of 25 dB, the received data from the primary channel are forwarded to the composite channel. On the other hand, when the SINR on the primary channel drops below 25 dB, the received data from backup channel are forwarded to the composite channel instead. Figure 2 shows typical plots of the SINR versus time under light rain. The average SINR was 28.5 dB, with a standard deviation of 5.6 dB. Channel 1 was the primary channel (first graph) and Channel 2 was the backup channel (second graph), utilized when the signal in Channel 1 declined below 25 dB. The third graph in Figure 2 shows that the utilization of Channel 2 as a backup channel successfully prevented the SINR from dropping below 25 dB, especially between 1157 s and 1172 s. various weather conditions: (i) light rain, (ii) medium rain, (iii) heavy rain, (iv) clear day. The rain level was classified as 'light', 'medium', and 'heavy' for precipitation rates of accumulation of about 0 to 2 mm per h, 2 to 10 mm per h, and 10 to 50 mm per h, respectively [63]. For each weather condition, the minimum and maximum 2 values were computed to determine the range of turbulence strengths.
Although the signals were received at 20 Gbps, the acquisition and evaluation of the channel matrix and interfering modes at high data rates were computationally intensive and required time. Thus, the SINR was sampled at 1 Gbps. In addition, the hovering time of the UAV was 45 min, so a 30 min interval was selected to give the UAV sufficient time to return to the ground station.
For each case, the time-varying SINR is plotted for the primary channel, backup channel, and the composite channel. The composite channel comprises various temporal segments of data from the primary channel and the backup channel such that when the SINR is equal to or more than the threshold value of 25 dB, the received data from the primary channel are forwarded to the composite channel. On the other hand, when the SINR on the primary channel drops below 25 dB, the received data from backup channel are forwarded to the composite channel instead. Figure 2 shows typical plots of the SINR versus time under light rain. The average SINR was 28.5 dB, with a standard deviation of 5.6 dB. Channel 1 was the primary channel (first graph) and Channel 2 was the backup channel (second graph), utilized when the signal in Channel 1 declined below 25 dB. The third graph in Figure 2 shows that the utilization of Channel 2 as a backup channel successfully prevented the SINR from dropping below 25 dB, especially between 1157 s and 1172 s.    (first graph) and Channel 4 was the backup channel (second graph), utilized when the SINR in Channel 3 decreased below 25 dB. The third graph in Figure 3 illustrates that employing Channel 4 as a backup channel successfully prevented the composite signal from depreciating below an SINR of 25 dB, preventing link failures to the cloud, especially between 670 s and 701 s, and 1352 s and 1396 s.  Figure 4 shows typical plots of the SINR versus time under heavy rain. The average SINR was 25.9 dB, with a standard deviation of 6.5 dB, showing more pronounced fluctuations in the SINR as compared to those under medium rain. Channel 7 was the primary channel (first graph) whilst Channel 8 was the backup channel (second graph). The number of regions with compromised SINR in heavy rain was higher than that under medium rain. The SINR for the composite channel, the third plot in Figure 4, shows that leveraging Channel 7 as a backup channel successfully alleviated a large portion of the composite signal from deteriorating below 25 dB, including during a predominant dip between 667 s and 710 s.  Figure 4 shows typical plots of the SINR versus time under heavy rain. The average SINR was 25.9 dB, with a standard deviation of 6.5 dB, showing more pronounced fluctuations in the SINR as compared to those under medium rain. Channel 7 was the primary channel (first graph) whilst Channel 8 was the backup channel (second graph). The number of regions with compromised SINR in heavy rain was higher than that under medium rain. The SINR for the composite channel, the third plot in Figure 4, shows that leveraging Channel 7 as a backup channel successfully alleviated a large portion of the composite signal from deteriorating below 25 dB, including during a predominant dip between 667 s and 710 s.  On a clear day, the average SINR was 30.4 dB. Thus, the backup channel was not required, as the SINR remained above 25 dB, as shown in Figure 5.  On a clear day, the average SINR was 30.4 dB. Thus, the backup channel was not required, as the SINR remained above 25 dB, as shown in Figure 5. On a clear day, the average SINR was 30.4 dB. Thus, the backup channel w required, as the SINR remained above 25 dB, as shown in Figure 5. Uplink and downlink transmissions demonstrated similar SINR characte Overall, a significant improvement in the SINR and outage performance was ob due to spatial mode diversity.
At the receiver, the received optical beam for each channel is demultiplexed t rate spatial modes. It is well-established that refractive index fluctuations in cloud pressure and temperature fluctuations cause rain [64]. This leads to aberration of th wavefront from random power and phase perturbations [65]. This results in the initially launched in a particular mode to spread to other modes [66]. The amo spreading is influenced by the strength of the atmospheric turbulence.
On a clear day, values were found to be between 7 × 10 −16 m −2/3 and 4 × 10 − peaking in the afternoon. On a clear night, the values were found to be lower, b 8 × 10 −17 m −2/3 and 1 × 10 −16 m −2/3 . For light rain, the values were in the range of 4 m −2/3 to 7 × 10 −16 m −2/3 . For moderate rain, the values were slightly higher, in th Uplink and downlink transmissions demonstrated similar SINR characteristics. Overall, a significant improvement in the SINR and outage performance was observed due to spatial mode diversity.
At the receiver, the received optical beam for each channel is demultiplexed to separate spatial modes. It is well-established that refractive index fluctuations in clouds from pressure and temperature fluctuations cause rain [64]. This leads to aberration of the beam wavefront from random power and phase perturbations [65]. This results in the signal initially launched in a particular mode to spread to other modes [66]. The amount of spreading is influenced by the strength of the atmospheric turbulence.
The power coupling coefficients were computed for all channels using a noninterferometric modal decomposition method [53], under various weather conditions, at the minimum and maximum C 2 n values. With knowledge of the power coupling coefficients, the modal crosstalk matrices were constructed, without and with spatial mode diversity. The power coupling coefficients were normalized to the total received optical power for the set of modes present.
Samples of inter-modal crosstalk matrices before mode diversity are shown in Figures 6  and 7, at the minimum and maximum C 2 n values, respectively. Samples of inter-modal crosstalk matrices after mode diversity are shown in Figures 8 and 9. For clear weather conditions in the day and at night, the power in the dominant mode was more than 75%. Under rain, power from the transmitted mode leaked into adjacent modes. The heavier the rain, the more substantial the spreading of power into adjacent modes due to increased atmospheric turbulence.
Under light rain, the power in the dominant mode was in the range of 68% to 75% before channel diversity reception and 79% to 83% after mode diversity and power were successfully maintained in the dominant modes through the switch to the backup channels during SINR dips. Under medium rain, the power in the dominant mode was in the range of 64% to 68% prior to mode diversity and improved to 75% to 78% in the presence of mode diversity. Under heavy rain, the power in the dominant mode was approximately 56% to 63% without mode diversity reception and increased to 71% to 75% with mode diversity reception. Uplink and downlink data transmissions demonstrated similar power coupling coefficient characteristics, whereas uplink and downlink data transmissions demonstrated similar power coupling coefficient characteristics. modal crosstalk matrices were constructed, without and with spatial mode diversit power coupling coefficients were normalized to the total received optical power f set of modes present.
Samples of inter-modal crosstalk matrices before mode diversity are shown ures 6 and 7, at the minimum and maximum values, respectively. Samples of modal crosstalk matrices after mode diversity are shown in Figures 8 and 9. Fo weather conditions in the day and at night, the power in the dominant mode was than 75%. Under rain, power from the transmitted mode leaked into adjacent mode heavier the rain, the more substantial the spreading of power into adjacent modes increased atmospheric turbulence.   power coupling coefficients were normalized to the total received optical power f set of modes present.
Samples of inter-modal crosstalk matrices before mode diversity are shown ures 6 and 7, at the minimum and maximum values, respectively. Samples of modal crosstalk matrices after mode diversity are shown in Figures 8 and 9. Fo weather conditions in the day and at night, the power in the dominant mode was than 75%. Under rain, power from the transmitted mode leaked into adjacent mode heavier the rain, the more substantial the spreading of power into adjacent modes increased atmospheric turbulence.   panel (a-e) shows the power coupling coefficients at low-refractive-index structure values, and tom panel shows (f-j) those at high-refractive-index structure values. Color legend of norma power coupling coefficients for all plots is on the bottom right corner.  Under light rain, the power in the dominant mode was in the range of 68% to before channel diversity reception and 79% to 83% after mode diversity and power w successfully maintained in the dominant modes through the switch to the backup c nels during SINR dips. Under medium rain, the power in the dominant mode was in range of 64% to 68% prior to mode diversity and improved to 75% to 78% in the pres of mode diversity. Under heavy rain, the power in the dominant mode was approxima 56% to 63% without mode diversity reception and increased to 71% to 75% with m panel (a-e) shows the power coupling coefficients at low-refractive-index structure values, and tom panel shows (f-j) those at high-refractive-index structure values. Color legend of norma power coupling coefficients for all plots is on the bottom right corner.  Under light rain, the power in the dominant mode was in the range of 68% to before channel diversity reception and 79% to 83% after mode diversity and power w successfully maintained in the dominant modes through the switch to the backup c nels during SINR dips. Under medium rain, the power in the dominant mode was in range of 64% to 68% prior to mode diversity and improved to 75% to 78% in the pres of mode diversity. Under heavy rain, the power in the dominant mode was approxima 56% to 63% without mode diversity reception and increased to 71% to 75% with m Overall, 25% to 44% of power was leaked to other modes during rainfall, prior to mode diversity. The power leakage was more significant for adjacent modes to the dominant modes and gradually decreased for modes further away from the dominant modes. With increased turbulence strength, the crosstalk was more prevalent and the power distribution of the modes was more dispersed. After mode diversity reception, an increase of 10% to 19% in the dominant mode power was observed. A larger amount of power was successfully maintained in the dominant modes through the switch to the backup channels during SINR dips.
The inter-modal crosstalk matrices for the primary channel, secondary channel, and composite channel after mode diversity is shown Figure 10. The power in the dominant mode of each channel was improved under all weather conditions. Without the mode diversity scheme, there were interruptions in the data transmission from sensors. By incorporating mode diversity, uninterrupted data transmission from sensors was observed, even under heavy rain due to the backup channels on different LG modes experiencing different channel degradations from the primary channels.
Overall, 25% to 44% of power was leaked to other modes during rainfall, prior to mode diversity. The power leakage was more significant for adjacent modes to the dominant modes and gradually decreased for modes further away from the dominant modes. With increased turbulence strength, the crosstalk was more prevalent and the power distribution of the modes was more dispersed. After mode diversity reception, an increase of 10% to 19% in the dominant mode power was observed. A larger amount of power was successfully maintained in the dominant modes through the switch to the backup channels during SINR dips.
The inter-modal crosstalk matrices for the primary channel, secondary channel, and composite channel after mode diversity is shown Figure 10. The power in the dominant mode of each channel was improved under all weather conditions. Without the mode diversity scheme, there were interruptions in the data transmission from sensors. By incorporating mode diversity, uninterrupted data transmission from sensors was observed, even under heavy rain due to the backup channels on different LG modes experiencing different channel degradations from the primary channels. Figure 10. Inter-modal crosstalk matrices for primary channel and secondary channel before mode diversity, and composite channel after mode diversity. The transmitted mode number and received mode number, i, are given in Table 1. The color legend shows the normalized power coupling coefficient.
The average diversity gain from spatial mode diversity was 0.97 dB. The average bit error rate (BER) versus SINR for all channels under various weather conditions was measured and compared, without and with spatial mode diversity, as shown in Figure 11. It Figure 10. Inter-modal crosstalk matrices for primary channel and secondary channel before mode diversity, and composite channel after mode diversity. The transmitted mode number and received mode number, i, are given in Table 1. The color legend shows the normalized power coupling coefficient.
The average diversity gain from spatial mode diversity was 0.97 dB. The average bit error rate (BER) versus SINR for all channels under various weather conditions was measured and compared, without and with spatial mode diversity, as shown in Figure 11. It was observed that spatial mode diversity successfully enhanced the BER by 38% to 55% from the original BER. After spatial mode diversity, a satisfactory BER below the 7% forward error correction (FEC) limit of 3.8 × 10 −3 was attained, at SINRs higher than 20 dB, under all levels of rain. An FEC of 7% was used to provide a 0.3 dB coding gain to better mitigate atmospheric turbulence. was observed that spatial mode diversity successfully enhanced the BER by 38% to 55% from the original BER. After spatial mode diversity, a satisfactory BER below the 7% forward error correction (FEC) limit of 3.8 × 10 −3 was attained, at SINRs higher than 20 dB, under all levels of rain. An FEC of 7% was used to provide a 0.3 dB coding gain to better mitigate atmospheric turbulence. Through UAV-FSO spatial diversity, despite mobility restrictions due to floods, vital health statistics from injured flood victims could be forwarded to a patient cloud-based health monitoring system accurately and timely, for hospitals to make critical preparations for relevant medical treatment upon the victims' arrival.
To verify the experimental measurements, the average values of were computed using [67]: where T is the average air temperature, W is the average relative humidity, S is the average wind speed from the local weather station, and w is a scaling coefficient. The average values of using [67] were 3 × 10 −15 on a clear day, 1 × 10 −16 on a clear night, 6 × 10 −16 under light rain, 8 × 10 −16 under moderate rain, and 2 × 10 −15 under heavy rain. This concurred with earlier experimental values computed using the scintillation index within 1 × 10 −15 . The refractive index structure parameter was in the range of 1 × 10 −15 to 1 × 10 −13 . Based on the turbulence strength classification in [62], this indicated weak to moderate turbulence.
A probability distribution function of the values from the experiment could be computed to determine if this adhered to relevant channel probability distribution functions for FSO spatial mode diversity and FSO aerial communications such as Malaga, log normal, Gamma-Gamma, Nagakami, and others [68][69][70][71][72]. Log normal and Gamma-Gamma channel models are considered as special cases of the Malaga channel model, for weak turbulence and moderate-strong turbulence, respectively [72].
Nevertheless, comparison of the time-varying SINR plots from the experiment with existing channel models do not provide accurate time-varying error analysis, as existing channel models are probabilistic [73]. Existing channel models only provide the statistical averages for the random variations and do not describe the evolution of the channel char- Through UAV-FSO spatial diversity, despite mobility restrictions due to floods, vital health statistics from injured flood victims could be forwarded to a patient cloud-based health monitoring system accurately and timely, for hospitals to make critical preparations for relevant medical treatment upon the victims' arrival.
To verify the experimental measurements, the average values of C 2 n were computed using [67]: where T is the average air temperature, W is the average relative humidity, S is the average wind speed from the local weather station, and w is a scaling coefficient. The average values of C 2 n using [67] were 3 × 10 −15 on a clear day, 1 × 10 −16 on a clear night, 6 × 10 −16 under light rain, 8 × 10 −16 under moderate rain, and 2 × 10 −15 under heavy rain. This concurred with earlier experimental values computed using the scintillation index within 1 × 10 −15 . The refractive index structure parameter C 2 n was in the range of 1 × 10 −15 to 1 × 10 −13 . Based on the turbulence strength classification in [62], this indicated weak to moderat turbulence.
A probability distribution function of the C 2 n values from the experiment could be computed to determine if this adhered to relevant channel probability distribution functions for FSO spatial mode diversity and FSO aerial communications such as Malaga, log normal, Gamma-Gamma, Nagakami, and others [68][69][70][71][72]. Log normal and Gamma-Gamma channel models are considered as special cases of the Malaga channel model, for weak turbulence and moderate-strong turbulence, respectively [72].
Nevertheless, comparison of the time-varying SINR plots from the experiment with existing channel models do not provide accurate time-varying error analysis, as existing channel models are probabilistic [73]. Existing channel models only provide the statistical averages for the random variations and do not describe the evolution of the channel characteristics with time. In simulations modelling the FSO channel, the time-varying atmospheric turbulence is typically emulated by continuously upgrading the random phase patterns to model random irradiance fluctuations, based on a known channel probability distribution function. This is only valid when the atmospheric turbulence is stationery during the laser pulse width [74]. Expanding on the previous work, we performed the false-nearest-neighbors algorithm on the measured time-varying channel matrix, which revealed that the fraction of false nearest neighbors converged to 0. This indicated the presence of nonlinear deterministic chaos. This opens up new possibilities for research on state-space reconstruction and the dynamic attractor for modeling the nonlinear FSO-UAV channel matrix.

Conclusions
A 4 × 2-channel UAV-based spatial mode diversity FSO system was developed for postflood recovery communications, achieving 2.4 Mbps for a distance of 400 m, under various weather conditions. Spatial mode diversity was realized by using distinct linearly polarized spatial modes in the primary and backup channels. Backup channels were used when the SINR in the primary channels declined below 25 dB. SINR plots from the composite channel showed that dips in the SINR from atmospheric fluctuations were successfully mitigated through spatial mode diversity. The power in the dominant modes increased by 10% to 19% through temporal switching to backup channels during SINR dips in the primary channels. Hence, this provided uninterrupted transmission, even under heavy rain. The spatial mode diversity scheme improved the BER by 38% to 55% compared to the original BER. In the future, the number of modes used for spatial mode diversity may be extended in conjunction with the maximal-ratio combination for better signal quality. The UAVbased spatial mode diversity FSO system is valuable for rapid recovery communications after natural disasters, especially when conventional ground base stations are damaged. Further analysis of nonlinear time-varying deterministic channel characteristics would be valuable for emulating atmospheric turbulence under various weather conditions, toward the implementation of more efficient turbulence compensation approaches.