ESA Earth Home Missions Data Products Resources Applications
    22-Jul-2014
EO Data Access
How to Apply
How to Access
FRINGE '96 Workshop: ERS SAR Interferometry, 30 September - 2 October 1996
Services
Site Map
Frequently asked questions
Glossary
Credits
Terms of use
Contact us
Search


 
 
 

Fringe 96

Forest Monitoring over Hilly Terrain Using ERS INSAR Data

André BEAUDOIN
Thierry CASTEL
LCT Cemagref-ENGREF
500, rue J. F. Breton, 34039 Montpellier cedex 05, France
andre.beaudoin@teledetection.fr
http://www.cemagref.fr/

Thierry RABAUTE

SCOT Conseil
1 rue Hermes, 31526 Ramonville, France
Tel : 61.39.46.00 / Fax : 61.39.46.10
tr@scot.cnes.fr

Abstract

ERS INSAR data proved recently to be of interest for forest applications. The interferometric coherence was found related to various land-uses and forest types, whereas the interferometric phase has been linked to the forest height, in the case of relatively flat terrain. This paper reports an investigation on the information content of the interferometric coherence over a hilly terrain supporting large pine plantations. The approach includes the use of a DEM-differential interferometric technique to remove the topographic effects on the interferogram. Then, coherence statistics for different acquisition intervals and baselines are analyzed as a function of land-use and forest types, forest parameters related to growth stage and wind conditions over the topography. Main results show that the coherence is linked to the land-use type and therefore, can be useful for discrimination purposes. Furthermore, the coherence is related indirectly to 1) the forest growth stage through variable ground contribution and tree stiffness under the wind influence and 2) the local wind speed over the topography. Assuming that these influences could be separated, these results suggest that the coherence could be used as an indicator of forest biomass or local wind speed.


Keywords: ERS repeat-pass interferometry, differential interferemetry, coherence, land-use and forest discrimination, topography, wind conditions.

Introduction

Recent studies using ERS SAR data showed that its configuration is somewhat limited for forest applications, due to ERS small incidence angle (23deg.) and low penetration depth in forest covers at C-band. In addition to low-level cover discrimination such as forest/non-forest mapping, possible applications are mainly related to the monitoring of natural or man-made forest perturbations such as clear-cutting (tropical deforestation), regrowth (low biomass levels) and environmental changes (thaw/freeze events). Fortunately, applications can be broadened significantly when interferometric phase and correlation derived from repeat-pass ERS INSAR data are considered, in addition to the usual backscatter information.

On one hand, the interferometric correlation, which is an indicator of the temporal stability of the target in terms of geometric and dielectric properties, proved to be a good discrimination in cultivated areas (Wegmuller and Werner 1995a) and forested landscapes (Wegmuller and Werner 1995b, Herland, 1995 ); multi-pairs correlation appeared to be particulary powerful (Wegmuller et al. 1995). On the other hand, it has been showed that the interforemetric phase, usually used to derive the terrain altitude, can also be linked to the height of the forest canopy (Askne et al. 1995, Hagberg et al. 1995, Ulander et al. 1995).

However, these are very recent results and the potential of such new data type is far from being fully explored. The interferometric, topographic and environmental conditions, for which INSAR data can provide useful thematic applications, have to be assessed furthermore. In this paper, we address the potential use of INSAR data over hilly terrain, towards the generalization of SAR-retrieved forest information in various conditions. The approach consists in using a differential interferometric technique developed at CNES (French Space Agency), to remove the topographic effects and keep the residual correlation and fringes. Then, differential interfrometric data are analyzed on a well-documented forest test-site featuring large plantations of coniferous trees found over hilly terrain.

Experimental set-up

The test-site

The test-site is situated in the central part of the Département de la Lozère, in Southeastern Massif Central (France). It is centered approximately at 44.5deg. N and 3.5deg. E. The area is characterized by large and gently rolling limestone plateaus culminating around 1200m, which are intersected by gorges with 300-500m depth and steep slopes up to 50deg. (see DEM in Fig.1a). Main land-uses are natural short grasslands, cultivated areas (small pot holes on the plateaus and valleys) and coniferous plantations. These latter are found on gorges slopes and partially on plateaus tops for erosion prevention and pulpwood production.

The planted species adapted to the difficult site conditions (soil thinness and summer drought) is almost exclusively Austrian pine (Pinus nigra). They are found in state-owned forest plantations made of even-aged and relatively homogeneous stands. The two main forest test-sites cover respectively 5 400 and 1200 ha with more than 500 stands (average area 5-15 ha), offering a large range of growth stages (0 to 140 year-old) as well as topographic situations.

Different ground data concerning these stands have been collected and entered in a GIS: 1) stand limits, 2) age classes, 3) detailed measurements of forest parameters and 4) a Digital Elevation Model (DEM) with a 50 m grid size in Lambert III projection (IGN (c)). Main stand limits and the DEM are included in Fig.1a. More than 300 homogeneous stands of different age classes and with area greater than 3ha were kept for the SAR data analysis. In addition, 103 stands were sampled for stem density, diameter at breast height (dbh) and stand height. Then, bole volume were estimated from simple allometric equations, with an estimation error around 10 to 20% (see Beaudoin et al. 1995 for more details).

INSAR acquisition and processing

Different ERS INSAR pairs, including one tandem acquisition, were acquired on this test-site. The various acquisition parameters that must be accounted for in the analysis are presented in Table 1. Only pairs with relatively small baselines (< 300m) were kept to prevent from excessive spatial decorrelation. Ignoring the tandem pair, it can be seen that the time acquisition interval is large, providing contrasted seasons and thus, variability in ground, forest and weather conditions. In particular, near-freezing conditions with a shallow snow layer on the ground were present for one winter scene included in pairs 1 & 3. However, low wind conditions (10 to 30 km/h) prevailed for all pairs.



Baseline
Acquisition Time
Wind speed
Air Temp.
ERS

Pair

(m)
Date
day

interval

Season
Date 1

(m/s)

Date 2

(m/s)

Date 1

(C deg.)

Date 2

(C deg.)

1
126
25/08/92

28/01/93

210
summer

winter

3.5
8
26.4
2.9
2
124
25/08/92

17/06/93

296
summer

summer

3.5
5
26.4
18.3
3
65
28/01/93

17/06/93

140
winter

summer

8
5
2.9
18.3
4
133
08/04/93

04/11/93

156
spring

fall

4
5
13.1
11.3
Tandem
98
19/08/95

20/08/95

1
summer

summer

4
3
28
25
Table 1 : Acquisition parameters of various ERS INSAR pairs over the forest test-sites in Lozère, South France.



Fig. 1a : DEM and pine plantations limits (red lines) over the Lozère test-sites.


Fig. 1b : Differential interferometric correlation (tandem mission) and pine plantations limits (red lines) over the Lozère test-sites. See Fig. 1a for visual comparison.


Several differential interferometric techniques, involving either SAR data pairs or triplets, have already been considered, mainly for the detection of ground displacements related to various types of hazards such as earthquakes, landslides or volcanic activity (Carnec et al. 1995, Massonnet et al. 1993, Massonnet et al. 1995). The concept of differential interferometry can also be applied to detect relative height differences within the interferogram : in that case, fringes patterns corresponding to these variations are superposed to those related to topography.

In our case, INSAR data were processed using a DEM differential technique (Massonnet 1994) : the site DEM is used in combination with the orbit geometry to simulate the topographic component and to remove it directly from the interferogram. Another approach would consist in unwrapping first the raw interferogram - including thus the terrain altitude and the forest height - and then substracting the terrain altitude provided by the DEM. The first method is more interesting in the way that there is no need for phase unwrapping if the average height of the forest cover is small compared to the altitude of ambiguity, that is to say the altitude variation corresponding to one interferometric fringe. However, the results may be corrupted if atmospheric perturbations occurred over the site area.

After this processing, interferometric coherence (or correlation) and residual phase are estimated from spatial averages within forest stands (more than 150 pixels) and are compared with parameters of interest. At this stage, we will focus only on the information content of the interferometric coherence, which must be analysed before investigating the phase, as the relevance of the latter directly depends on the value of the former.


2) EXPERIMENTAL RESULTS

2.1) General correlation behavior

Fig. 1b presents the geocoded correlation image including the forest stand limits, while Fig. 2 presents mean correlation values for the 3 main land-use types, plotted against time interval. Globally, brighter areas correspond to grasslands (Causses area in Fig. 1b ), which always exhibit a good correlation in Fig. 2 (>0.85), even for the pairs including the winter scene. Forested areas are darker and correspond very well to the stand limits. However, bright areas are also found in these areas (see under Mende in Fig. 1b). In Fig. 2, young stands have a higher correlation than mature ones; in addition, the correlation decreases for both stand ages with acquisition interval. Finally, cultivated areas also present in the scene can be found as small dark patches on the Causses (cultivated pot holes) or dark stripes in the main valley (Florac area). In most cases, agricultural areas are confused with mature forests. Beside this, the potential of interferometric correlation can be foreseen for general land-use discrimination, which was not possible using ERS backscatter intensity alone.

Fig. 2: Mean correlation as a function of acquisition interval, for 3 main land-use types.

2.2)Forest correlation behavior

Looking at Fig. 1b , large correlation variations are found within forested areas. Correlation was estimated on 103 stands presenting different biomass levels and topographic situations. Fig. 3a presents correlation estimates for the tandem acquisition, as a function of biomass level. A decrease with biomass is observed with r2=0.52 from linear regression, but the trend is too weak to estimate correctly the biomass. Lower r2 were observed for other pairs, but the same trend appeared with a more or less pronounced negative slope. Such decrease was partially observed also by Herland (1995) between some forest stands.

Fig. 3 : Correlation behavior with A) bole volume as an indicator of the growth
stage, and wind exposure angle for B) young pines and C) mature pines.

2.3)Preliminary interpretation

To interpret qualitatively above observations, simple theoretical considerations are addressed in the following. Considering that most of the spatial decorrelation due to terrain slope was remove during the DEM-differential interferometric processing, the total decorrelation rho can be decomposed in its different sources (Askne et al. 1995, Zebker and Villasenor, 1992):

rho total = rho thermal * rho spatial * rho temporal (1)

including the thermal noise decorrelation which is neglected here, the spatial decorrelation due to the baseline that has been removed after Zebker and Villasenor (1992) and finally the temporal decorrelation. This later is mainly an indicator of the temporal unstability of the target geometric and dielectric properties (Wegmuller and Werner 1995a).

For interpretation purposes, we first consider the backscatter from a discontinuous forest layer over a ground surface. At C-band, the backscatter mainly originates from needles and twigs volume scattering sigma forest, in addition to a direct ground scattering sigma ground which contributes more or less with soil moisture and roughness, along with transmissivity T and forest fraction cover F, which in turn are linked to the growth stage. This can be approximated by:

sigma = F (sigma forest + T2 sigma ground) + (1- F) sigma ground (2)

Now, we consider the temporal decorrelation of these 2 main scattering sources. In our case, for non-forested areas (F=0), the ground scatter should be geometrically a stable surface, considering natural short grasslands without mechanical cultivation. However, dielectric seasonal changes occur due to surface soil moisture changes (snow and rain), but should affect to a less extent the correlation (Wegmuller and Werner, 1995a). Therefore, this explains the high correlation observed for such surfaces in Fig. 2.

When adding a forest layer (here, F within [0.6-1.0] and crown layer depth from 3 to 8 meters), the ground contribution will decrease with growth stage. In addition, temporal decorrelation due to volume scattering at C-band will add and should be considered at 2 different time spans:

· On a monthly basis during the growing season, the crown geometry changes due to annual new shoots, which is probably a source of significant decorrelation.

· On a very short time basis (maybe down to few minutes), decorrelation occurs due to random displacements of the scatterers (needles and twigs) under the wind influence.

Zebker and Villasenor (1992) linked this decorrelation to the RMS displacement of the scatterers. For ERS configuration, only few RMS cm are enough to decorrelate the signal, which we qualitatively observed even for weak winds. However, the link between the displacement and the wind and forest parameters has been merely investigated experimentally (Askne et al. 1995). It should depend on:

· local wind speed over the canopy, in addition to the topographic exposure to wind which will be described in the following;

· tree stiffness, expected to be related to tree species and structure (Hagberg et al. 1993).

No model exists to relate the scatterers displacement to above parameters. Askne et al.( 1995) use a parameter k in a simple exponential model to describe the temporal decorrelation within the crown layer, k being related in an unknown way to above parameters. This will be investigated in a near future on our data set.

At this stage, we will focus on the topographic effect on the correlation, knowing the hourly wind speed |W| and the topography through the DEM. We consider that the topography attenuates the speed of the wind vector W, giving a local wind vector Wloc through:

Wloc = W*T(thetaw) (3)

where T is an unknown attenuation function assumed to depend on a wind exposure angle, that we defined as:

thetaw = cos-1[sin(alpha) * cos(beta-betaw)] (4)

with alpha and beta being the terrain slope and aspect angles, and betaw the wind direction. thetaw is simply the angle between the wind vector W and the normal vector to the terrain N, as sketched in Fig. 4.

Fig.4 : Modification of the wind vector W due to the topographic exposure to wind thetaw, and effect on the scatterers displacement for young and mature trees. For simplicity, aspect angle beta is here equal to betaw -180deg. on the left side and betaw on the right side.

First, from all above considerations, we can interpret Fig. 2 as follows :

· as usually observed, forests have lower correlation due to volume decorrelation;

· higher correlation for young pines compared to mature ones can be explained by higher ground contribution and/or higher tree stiffness, as observed in the field;

· lower correlation for larger time spans, for both pine ages, could be explained by an additional volume decorrelation due to new annual shoots (about 30 cm increase).

Second, we can interpret in more details the correlation decrease with bole volume in Fig. 3 :

· Correlation decrease with bole volume is mainly due to decreased ground contribution and increased displacement with stand age due to decreased stiffness; this can be understood as for mature trees, needles are mainly holded by branches of higher order (3 or 4) which are long and thin, compared to young trees where needles are holded by low order branches (1 or 2) which are short and relatively big, thus stiffer;

· The high dispersion around the trend in Fig. 3a can be explained by the topographic exposure to wind. If we consider in Fig. 3b&c two data subsets, that is young and mature pines respectively, and plot the correlation against the wind exposure angle (eq. 4), we see that:

· for mature trees in Fig. 3c, the correlation is highly correlated to thetaw, with a quasi-linear decrease (r2=0.78); this explains also the correlation variation visually observed just below the Mende town in Fig. 1b, which seems to be linked to the topography (2 small parallel valleys ( NE-SW) with old pine stands);

· for young pines in Fig. 3b (see also Fig. 1b on the Causses Méjean), no relationship with thetaw is observed; probably due to the correlation mainly driven by the ground contribution, in addition to the higher stability of young trees under the local wind influence. Therefore, dispersion in Fig. 3b could be explained in the variability of the transmissivity and the cover fraction, in addition to soil moisture and roughness that all modulate the ground contribution.


Conclusions

In this paper, ERS DEM-differential interferometric correlation was investigated over a hilly forested test-site. This correlation was found to be a good land-use discriminator for major land-use types, especially using shorter time intervals. On forest covers, the correlation is partially linked to the bole volume (growth stage of the forest in fact), due to decreased ground contribution and increased crown volume scattering, originating mainly from needles and twigs. With ERS characteristics (C-band 23deg.), such volume scattering rapidly decorrelates due to centimetric random displacements easily generated by low winds (10km/h). As is it rare that there is no wind at all, and as the wind speed can change rapidly, we can expect that very high correlation on forests can not be achieved using repeat-pass interferometry, towards a performant DEM generation over forested areas.

Although the wind conditions were quite similar for the different ERS pairs investigated, the correlation behavior with the wind speed was highlighted using the tandem scene, through the variation of the local wind speed due to the topographic exposure to wind. Therefore, in a sense, for similar mature pine forest, we can say in our case that the correlation provides us with an indicator of the local wind speed at the canopy level.

Globally, this study confirmed that the correlation is linked to many parameters, including those related to the interferometric configuration (baseline and incidence angle), the weather conditions (wind speed and direction), the topography and the ground and forest covers. Therefore, for generalization purposes, to which point this correlation can be linked robustly to characteristics of interest is a new and exciting challenge. Carying on such work is probably worthwhile, considering the potential of interferometric SAR data in addition to the usual backscatter information, especially at L-band.

To this aim, future works will include the use of semi-empirical models as proposed by Askne et al. (1995) for a better understanding of the data set. Also, other tandem couples will be processed. Finally, ERS differential interferometric pairs with sufficient correlation will be investigated for the phase information content, which should be linked to the forest height.

Acknowledgments

This work was performed within a joint LCT-SCOT project funded by the CNES (contract # CNES/94/0239), that we would like to warmly thank. Particularly, we thank Didier Massonnet for the access to the CNES interferometric processor. Thanks also go to ESA for providing the SAR data within the project-pilot PP2-F132, and P. Gigord from IGN who provided us with the DEM (BD CARTO IGN (c)) within a collaboration on the correction of SAR bacskcatter data acquired over hilly terrain. Final thanks are for Luce Castagnas, from LCT, for the html editing.

References

Askne, J., Dammert, P., Fransson, J., Israelsson, H. and Ulander, L.M.H, 1995:
Retrieval of forest parameters using intensity and repeat-pass interferometric SAR information, Proceeding of the International Symposium on Retrieval of bio- and geophysical parameters from SAR data for land applications, Toulouse, France, 10-13 october 1995, pp. 119-129.

Beaudoin, A., Castel, T., Deshayes, M., Stachs, N., Stussi, N. and Le Toan, T., 1995:
Biomass retrieval over hilly terrain from spaceborne SAR data, Proceeding of the International Symposium on Retrieval of bio- and geophysical parameters from SAR data for land applications, Toulouse, France, 10-13 october 1995, pp. 131-140.

Carnec C., Massonnet D. and King C., 1995 :
Measurement of land subsidence by differential SAR interferometry, Proceedings of the Fifth International Symposium on Land Subsidence, The Hague, The Netherlands, 16-20 October 1995.

Hagberg, J.O., Ulander, L.M.H. and Askne, J., 1995:
Repeat-pass SAR interferometry over forested terrain, IEEE Trans. on Geoscience and Remote Sensing, Vol. 33, No.2, pp. 331-340, 1995.

Herland, E.A., 1995:
SAR interferometry with ERS-1 in forested areas, Proc. of IGARSS'95 Symposium, Florence, Italy, 10-14 July 1995, pp. 202-204.

Massonnet D., Rossi M.,Carmona C.,Adragna F.,Peltzer G., Feigl K., Rabaute T.,1993:
The displacement field of the Landers earthquake mapped by radar interferometry, Nature, vol. 364, pp 138-142, 8 juillet 1993.

Massonnet D., 1994 :
Giving an operational status to SAR interferometry, Proceedings of the First
Workshop on ERS-1 Pilot Projects,
Toledo, Spain, 22-24 June 1994.

Massonnet D., Briole P. and Arnaud A., 1995 :
Deflation of Mount Etna monitored by spaceborne radar interferometry,
Nature, vol. 375, 5 juin 1995 pp, 138-142.

Ulander, L., Dammert, P.B.G. and Hagberg, J.O., 1995:
Measuring tree height with ERS-1 SAR interferometry, Proc. of IGARSS'95 Symposium, Florence, Italy, 0-14 July 1995, pp. 2189-2191.

Wegmuller, U. and Werner, C.L., 1995a:
Farmland monitoring with SAR interferometry, Proc. of IGARSS'95 Symposium, Florence, Italy, 10-14 July 1995, pp. 544-546.

Wegmuller, U. and Werner, C.L., 1995b:
SAR interferometric signatures of forest, IEEE Trans. on Geoscience and Remote Sensing, Vol. 33, No.5, Sept. 1995, pp. 950-959.

Wegmuller, U., Werner, C.L. and Nuesch, D., 1995:
Retrieval of vegetation parameters with SAR interferometry, Proceeding of the International Symposium on Retrieval of bio- and geophysical parameters from SAR data for land applications, Toulouse, France, 10-13 October 1995, pp. 11-18.

Zebker, H.A. and Villasenor, J., 1992:
Decorrelation interferometric radar echoes, IEEE Trans. on Geoscience and Remote Sensing, Vol. 30, No.5, Sept. 1992, pp. 950-959.



Keywords: ESA European Space Agency - Agence spatiale europeenne, observation de la terre, earth observation, satellite remote sensing, teledetection, geophysique, altimetrie, radar, chimique atmospherique, geophysics, altimetry, radar, atmospheric chemistry