An Empirical Assessment of Temporal Decorrelation Using the Uninhabited Aerial Vehicle Synthetic Aperture Radar over Forested Landscapes

We present an empirical assessment of the impact of temporal decorrelation on interferometric coherence measured over a forested landscape. A series of repeat-pass interferometric radar images with a zero spatial baseline were collected with UAVSAR (Uninhabited Aerial Vehicle Synthetic Aperture Radar), a fully polarimetric airborne L-band radar system. The dataset provided temporal separations of 45 minutes, 2, 7 and 9 days. Coincident airborne lidar and weather data were collected. We theoretically demonstrate that UAVSAR measurement accuracy enables accurate quantification of temporal decorrelation. Data analysis revealed precipitation events to be the main driver of temporal decorrelation over the acquisition period. The experiment also shows temporal decorrelation increases with canopy height, and this pattern was found consistent across forest types and polarization.


Introduction
Interferometric Synthetic Aperture Radar (InSAR) has been used to produce digital elevation models (DEM) of the Earth for several decades [1][2][3], and has matured to produce global elevation datasets [4].In radar interferometry, two radar acquisitions are combined to obtain coherence (i.e., similarity between the acquisitions) and phase difference between these acquisitions [5,6].When obtained from slightly different vantage points, this phase difference (or phase center) can be used to estimate the elevation of the target (e.g.surface topography).In most cases, DEMs derived from interferometry use two antennas that are mounted on an aircraft or spaceborne platform, acquiring data simultaneously.But early work has shown that a single sensor repeating nearly the same flight track can also be used to produce DEM's [3].It was soon realized that the interferometric coherence decreases with time between radar acquisitions because of intrinsic changes in surface reflectivity that increase phase noise [5] and reduce elevation measurement accuracy [3,[7][8][9].Such phenomena are a major unknown in the design of InSAR systems [5] and are still not well understood.
Recent advances that exploit the complementarity between radar interferometry and polarimetric diversity now enable estimates of forest canopy height [10][11][12][13].The height of the phase center within a forest canopy depends on the structural parameters of the canopy as well as the incident and received signal polarization.The so-called Polarimetric interferometric Synthetic Aperture Radar (PolinSAR) technique uses a simplified radar scattering model to invert canopy height or other structural parameters given a set of coherence and phase observations at various polarization configurations.Over short repeat-pass time (i.e., scale of minutes to a few hours), temporal decorrelation over forest usually may be neglected [14,11] but is otherwise significant at longer time scales [7,9].The motion of leaves and branches caused by wind as well as variations in moisture content (for example from precipitation events) reduce the observed interferometric correlation as do phenological changes of leaf emergence, growth, senescense and fall [7,8].In [14], it has been shown, using a physical model, that temporal decorrelation depends on the structural parameters of forests, such as canopy height, and changes with polarization.Because the observed coherence becomes dominated by temporal changes rather than the forest volume, it may inhibit estimation of canopy height.The impact of temporal decorrelation has been partially addressed in the PolinSAR model [12,13].In order to assess DEM accuracy and compensate for its impact in the PolinSAR model, it is necessary to quantify the contribution of temporal decorrelation on InSAR coherence.
In this paper we empirically measure the impact of vegetation structure and type, weather, and sensor configuration on temporal decorrelation of L-band interferometric coherence.Measurements of were acquired during a UAVSAR (Uninhabited Aerial vehicle Synthetic Aperture Radar) experiment specifically designed to assess temporal decorrelation.UAVSAR is an airborne L-band fully polarimetric SAR that enables accurate repeat-pass interferometry [15].The remainder of this paper is structured as follows.First, we describe the data collection strategy and parameters of the experiment.Next we develop a theoretical treatment of coherence and demonstrate, using this derivation, that UAVSAR performance enables accurate measurements of temporal decorrelation.We then analyze the dependence of InSAR coherence on forest type, canopy height, repeat-pass time interval, weather and polarization using UAVSAR data and conclude with a discussion of our major findings.

Data Collection Strategy
The UAVSAR campaign was designed specifically to identify factors that impact InSAR coherence and quantify temporal decorrelation.UAVSAR data were acquired in the Réserve Faunique des Laurentides (Québec, Canada), located between Québec City and Saguenay (47.7°N, −71.3°W) in Eastern Canada, in August 2009.The site is characterized by a 1,000-m elevational gradient with a transition from deciduous temperate to boreal coniferous forest, and includes mature forests in conservation units, as well as managed stands in forest extraction zones.Thus, this site enables the study of temporal decorrelation across broad gradients in forest composition, structure and successional stage.
In general, repeat-pass InSAR coherence estimated over forested areas is dominated by volume decorrelation, temporal decorrelation and thermal noise decorrelation [4].For an arbitrary spatial and temporal baseline, the volume and temporal decorrelation can be of the same magnitude.To isolate the impact of temporal decorrelation from the volume decorrelation, we performed a set of UAVSAR flights with nominal zero spatial baselines.The UAVSAR system is designed to fly within 5 m of the planned flight track [16].The residual motion of the platform leads to local deviations of the baseline along the flight trajectory.In the next section, we show UAVSAR's zero-baseline repeat-pass InSAR acquisitions are sufficiently accurate to ensure that the InSAR coherence measurement is dominated by temporal rather than other system and volume effects.
To measure temporal decorrelation, UAVSAR data were collected several times repeating the same flight track (i.e., zero baseline) to obtain a set of temporal intervals from minutes up to days.Flights were performed twice on 5, 7 and 14 August 2009.This strategy provided temporal periods of 45 min, 2 days, 7 days and 9 days by combining all possible pairs and generating radar interferograms.

3.Accuracy of Temporal Coherence Measurement
In this section, we estimate the expected accuracy of the temporal decorrelation measurement δ τ using the UAVSAR system.The observed InSAR coherence can be expressed as the product of noise γ snr , geometric γ g , volume γ v and temporal coherences γ t such that [7]: Given estimates for the other coherence terms, measurements of temporal coherence can be obtained from: (2) In order to obtain the error stemming from the various decorrelation sources, we differentiate Equation ( 2) to obtain: . ( We obtain the individual relative error terms from differentiation of published definitions for γ snr , γ v and γ g.The SNR correlation γ snr with SNR being the signal-to-noise ratio in linear units as given by [5,7]: .
Differentiating with respect to SNR, we obtain: .
For volumetric coherence, we set an upperbound using a uniform scattering model with [17]: where h and k z are the canopy height and the vertical wavenumber respectively.Assuming a flat Earth, the latter is given by [5,6]: where is the perpendicular baseline, is the signal wavelength, π = 2 for repeat-pass observations, the look angle and h p the platform altitude.Differentiating, we obtain: and therefore Finally, the geometric coherence for a flat surface is given by [5,7] where is the range resolution.Correcting for geometric decorrelation using filtering [17,18] results in where is the terrain slope in the range direction.Differentiating equation (11) with respect to leads to: From Equations ( 5), ( 9) and ( 12), we obtain the relative error in our measurement of temporal decorrelation δ τ (δ τ =1− ): Therefore, for the UAVSAR system with a SNR better than 25 dB, maximum baseline of 2.5 m, wavelength of 0.24 m, range resolution of 1.66 m and a platform height of 12.We used weather data from Environment Canada's National Climate Data and Information Archive covering the UAVSAR acquisition period.Figure 2 shows the recorded precipitation starting from 1 August 2009 and wind speed at the time of the radar data acquisitions.
We used the forest map from Natural Resources Canada (canadaforests.nrcan.gs.ca) and a Landsat scene to identify regions of coniferous boreal and mixed deciduous forests.Forest pixels were selected with RH75 greater than 2 m and UAVSAR HV backscatter above −20 dB.While the first criteria indicated the presence of a canopy the latter threshold removed smaller patches (pixels) of bare ground and water.

Analysis of Temporal Decorrelation
We analyzed the zero baseline InSAR coherence (Equation ( 1)) as a function of time, weather, canopy height, forest type as well as polarization.The InSAR HH coherence images are shown in Figure 3 for the 45 min, 2-, 7-and 9-day time intervals.Overall, the highest coherences were observed for the 45 min and 9-day intervals.

Weather
Variations in precipitation and wind induce changes in the complex radar backscattering properties of the canopy and the underlying ground.These changes result in a modification of the InSAR coherence, which could significantly impact the measurement accuracy of canopy height retrieval using repeat pass InSAR [7][8][9][12][13][14].Figure 2 shows the precipitation and wind speed recorded during the UAVSAR campaign.A day prior to the first acquisition on 5 August, 4 mm of rain was recorded within the Reserve, causing the canopy and soil to be moist at the time of the acquisition.On 7 August, both the canopy and ground were saturated with water as a large rain event occurred during the flight.The last acquisition was the driest with no rain during the previous days.
Figure 4 shows the frequency distribution of HH temporal coherence as a function of temporal baseline for all forest pixels within the UAVSAR image.The largest observed InSAR coherence occurs for same day InSAR pairs while the lowest coherence is measured from the 7-day temporal baseline.Not surprisingly, any interferogram using the August 7th acquisition produces low coherence.The 9-day temporal baseline interferogram produced with the August 5th and August 14th acquisitions has higher coherence than the 2-day August 5th and 7th interferogram.The large precipitation event that occurred during the acquisition of August 7th is the main factor responsible for the decrease in measured coherence.Frequency (% per 0.01 bin) 45'  45' (7-7) 45' (14-14) 2 days 7 days 9 days

Canopy Height
The canopy height estimates from LVIS (RH75) were used to investigate potential relationships between canopy structure and temporal decorrelation δ τ .The canopy structure per se is not expected to affect coherence.However, the impact of certain environmental phenomena, such as wind and interception of precipitation will vary as a function of canopy height [14].We aggregated UAVSAR coherence pixels to 20 m to compare with the LVIS canopy height estimates to investigate this effect.Figure 5(a) shows the mean coherence as a function of canopy height for all time intervals for both temperate and boreal forests.It shows that changes in target characteristics due to precipitation (e.g., moisture content, surface water, resulting displacement of leaves and branches) are the most significant drivers of temporal decorrelation.In addition, temporal decorrelation δ τ increases with canopy height, i.e., decrease in as in Figure 5(a).This trend is observed across all time intervals and forest types and is consistent with our model for the impact of wind.There is, however, a significant variability in the observed coherence within a canopy height class with standard deviations of about 0.08.

Polarization
The dependence of coherence on polarization is shown in Figure 5

Conclusions
The impact of temporal decorrelation in repeat-pass interferometry has been repeatedly observed during the historical development of repeat-pass radar interferometry (e.g., [3,7,9]).However, there has not been sufficient data to identify and understand the drivers of this decorrelation.The lack of a quantitative comprehension of these effects has limited the application of repeat-pass InSAR for canopy height derivation.To address this issue, we analyzed the dependence of InSAR temporal coherence on forest type, canopy height, repeat-pass time interval, weather and polarization using UAVSAR data.We formulated the expected measurement accuracy for the UAVSAR system and found it was adequate (i.e., ) for measurement of temporal decorrelation.Our experiment has shown that for our study area the main drivers of decorrelation were related to precipitation and canopy height.The main decrease in coherence was observed when interferograms were produced from radar acquisitions collected with different weather conditions, especially recent precipitation events.The range of InSAR coherence measurements for any given InSAR pair was large and shows a relationship to canopy height.InSAR coherence decreased by about 10% to 20% over the range of canopy heights (~20 m), even at short repeat time intervals regardless of polarization diversity and forest type.The standard deviation of coherence measurement over a given copy height class was around 0.08.
In interferometry applications related to measurement of ground topography, temporal decorrelation increases phase noise and therefore ground measurement error [5].Although this effect cannot be compensated for, estimates of temporal decorrelation can be used to estimate elevation measurement accuracy through analytic expressions [6].PolinSAR has shown great potential for measurement of canopy heights.However, repeat-pass methods for canopies, in contrast to bare earth, must contend with the pliant and dynamic nature of trees; leaves and branches change on short time scales as a result of their interaction with the weather.By quantifying the magnitudes of these decorrelations and exploring their causes, our work may potentially be used to improve PolinSAR canopy height estimates [12,13].

Figure 4 .
Figure 4. Observed zero-baseline InSAR coherence for all forests taller than 2 m.The legend gives the temporal baseline.The 45 min temporal baselines are shown for acquisitions on the 5th, 7th and 14th of August.
(b) for the 45 min and 2 day periods.Standard deviation (~0.08) is only shown on Figure 5(b) but also apply to Figure 5(a).Temporal decorrelation trends and levels for a given time interval are similar at all polarizations.

Figure 5 .
Figure 5. Variation of termporal coherence with canopy height.(a) HH coherence for mixed termperate and boreal coniferous forests at various time intervals and (b) for all forest types at various polarizations and two time periods.The vertical bars in (b) give the standard deviations for discrete height intervals and are not shown in (a) for clarity.