Measuring thickness change of glaciers from satellite data: a case
Mean thickness change of Alaskan glaciers as a function of altitude (from. Arendt et al., 2002). 1. Scientifical background. 2. Methodology. 3. Results. 4. Outlooks ...
to estimate the contribution of glaciers to sea level rise. to improve our knowledge of the regional impact of climate change on glaciers and water ressources But: Only 250 glaciers on earth with mass balance measurements (2% of the total glacierazied area, Braithwaite, 2002)
" !
! " # $%&&%
'#
(
"
%#
)
*#+ "
# "
Field measurements “glaciological method” Time consuming and expensive, harsh environment
Aerial photography “geodetic method” (Funk et al., 1997) Precise but only cover limited area (30 km2)
Airborne laser altimetry (Arendt et al., 2002) Precise, rapid, but expensive and sampling along a few profiles
Aim of this study: to demonstrate that we can use satellite images to derive glacier mass balance. A first step: the measurement of thickness change
'#
(
"
%#
)
*#+ "
# "
Data : - SPOT images - 30 precise GCPs (D-GPS) - DEM computed using PCI geomatica - same GCPs for all images
Validation data : - DEM derived from aerial photographs (same process) - Transverse profile measured each year with GPS on glaciers ),
'#
(
%#
"
)
*#+ "
year #1
# "
year #2
Stereoscopic effect
DEM #1
DEM #2 DEM #2 - DEM #1 differential DEM
Thickness change of glaciers
(
"
%#
)
*#+ "
# "
! Using the fact that, out of the glaciers, no changes are expected ! 6
Planimetric adjustment based on the cross correlation of orthorectified images
Altimetric adjustment based on a polynomial fit to the bias with altitud between the 2 DEM
4 Bias between the two DEM
'#
2 0 -2 -4 -6 -8 -10 700
Before correction After correction 1200
1700
2200
2700
3200
altitude
%&&&
-" %&&*
( .
'#
(
%#
"
)
*#+ "
# "
" 2000
,-
.
/0
Aerial DEM
.
,$ 1 2
Satellite DEM
/,
2003
,0 1 2
GPS profile
.
/$- 1 2
,$ 1 2 GPS profile
Satellite DEM
Ref.
Aerial DEM
Cumulated ablation calculated from a degree day model at 3 different altitudes from the 1st August to the 15th September 2000
%&'(
%
!"# $ ) $$$ *&'+' /
$%&&%
(
Ablation (m we)
4 3,5
2550 m
3
3550 m
2,5
1550 m
2 1,5 1 0,5 0 1/8/00
10/8/00
19/8/00
Time
28/8/00
6/9/00
15/9/00
'#
(
"
%#
)
*#+ "
# "
#
( .
%&&*
%&&&
'0
'#
(
"
%#
*#+ "
)
$
# "
% &'
( ) $'
Accumulation area
!" " %&&&1 Comparison with aerial DEM
'0 $'&
!" " %&&*1
$%#
Out glacier : N = 4734
-= 1.5 m σ = 15.9 m
Out glacier : N = 1549446
-= -0.3 m σ = 6.2 m
On glacier : N = 3190
-= 6.2 m σ=6m
On glacier : N = 806630
-= 1.4 m σ = 2.5 m
'#
(
%#
"
*#+ "
)
# "
*
1994 - 2000
2000 - 2003
5
1600
1700
1800
1900
2000
2100
2200
2300
0 1500
-10
thickness change (m)
thickness change (m)
0 1500 -5
5
-15 -20 -25 -30 -35
GPS profile Satellite DEM
-40
1600
1700
1800
1900
2000
2100
-10 -15
GPS profile Satellite DEM Aerial Photo. DEM
-20
-45
-30 altitude (m.a.s.l.)
altitude (m.a.s.l.)
.
2300
-5
-25
-50
2200
2
2400
'#
(
%#
"
*#+ "
)
# "
" Larger coverage Good resolution
3
3 2
( .
%&&&
%&&*
'#
(
%#
"
)
*#+ "
# "
%' 0
thickness change (m/a)
-1 -2 -3
+
4 2
-4 -5
79-94 94-00
-6
00-03
-7 1500
1800
Altitude (m)
2100
2400
thinning at low altitudes Acceleration of the thinning
'#
(
%#
"
)
*#+ "
# "
" This method is suitable to measure thickness change occuring within a few years (signal over ~5 meters) Acceleration of the melting below 2400m Extend the study to the whole Mont Blanc range (effect of orientation, slope) Production of a new DEM from SPOT5 images during summer 2004 to estimate the annual mass balance With the advent of satellites with submetric resolution (0.7m for PLEIADE), a new way to semi-automatically estimate the annual mass balance of worldwide glaciers.
$
+
Project funded by: ACI-OT Glaciers ISIS program at CNES
Ed. Belin, 31401, Toulouse Cedex 9, France. bCNES, 18 av. ... Ed. Belin, 31400, Toulouse, France .... First, the surface features of the glaciers may change due.
followed in developing the Brussels UrbIS data dictionary, while the fourth section. 4. D. Pantazis .... charge in digital format or in paper (more expensive option).
A historical view of twentieth-century accounting developments, for instance, .... 21st of March, a son of Luis Medina, Rafael Medina was appointed as managing .... activity... For achieve these needs, accounting requires more resources to its.
[10,19,20]. The present day kinematics and strain ... (extension strain rates) [19], Engdahl et al. [18] (0â70 km .... a good wavelength coherence for the sources. A.
to contour maps (50s- ... ~1970s (USGS, YT maps) and 2003-2007 (ASTER, SPOT5 DEM). Thickening ... Google Earth images of glaciers illustrating that debris-.
edge slope (tes) adjustment variables; H0 is the mean altitude; sx and sy are the ..... computed by fitting a plane to the surface within a disc with a 20 km radius. Up to km ..... venice06/participants/181/paper_181_legrsy.pdf. Liu H, Jezek K, Li B
Despite progress in de¢ning the surface topog- raphy .... only degrees 0^60 when computing the free air ..... Center Key to Geophysical Records Documentation,.
Aug 18, 2010 - change pycnocline thickness dz through control by Richardson number Ri. Shear ..... In this case dz ¼ 4.0 m, du ¼ 0.4 m s21, Ëg, or du/dz, ¼ 0.1 ...
Received 13 July 2007; received in revised form 26 November 2007; accepted 1 December 2007 ... dependence of the microwave scattering on the snow grain size has been ... dynamic retrieval algorithm has been developed and validated ...... inclusion of
Simulation Experiment Project of the 5th EU Framework program, and in part by the INTAS ... environmental management. Nowadays ... times series of ice cover parameters necessary for practical ..... Science and co-ordination plan. Version 1 ... Sensin
with Department of Rehabilitation Sciences, Faculty of Physical. Education and ... Five complete crank cycles were analyzed for all parameters. MATLAB ...
also generate visible albedo (Diner 1996), and the. MODIS team is to produce all three broadband albedo products: visible, near-infrared, and total shortwave.
The ESRC Research Centre for Analysis of Social Exclusion (CASE) was established in ... stageâ indices â those that are applied directly to a sample from a.
authors introduced network morphogenesis models with diverse modes of preferen- tial link creation depending on various node properties (attractiveness [17], ...
structure. As a consequence, the structure studied with standard part a property of the shape of the protein itself. One of. NMA is always a distorted one. Finally ...
heavy-duty rail car with a specially built cover. Each nozzle has ...... self test. If available, check concentration on SM SYS SUMM 1. 0001/. /078. SM SYS SUMM ...
Abstract. The present-day distribution of glaciers in Africa is limited to three specific ..... then become progressively obscured by clouds that do not dissipate until.
the radius and if all airfoils would have a flat lower side, then ... edge of the propeller, before the felt tip color has dried. Doing the same at the trailing edge leaves you with two traces on the ... method - the error is around 1 percent of the.
Clinical data warehouse, Clinical information system, Hetero- geneous data ... ADT, healthcare record (EHR), and act management (CPOE) components from ...
June 12, 1998. NGK Insulators, Ltd. ... Specifications of measuring tools are shown in Table 1. Portable tools are ... Measuring and calculation procedure is shown as follows. ... g) Attention should be taken not to loose the quantity of the water (i
aInstitute of Electronic Commerce, National Chung-Hsing University, 250 Kuo Kuang Road, 403 Taichung, Taiwan, ROC. bDepartment of ..... Supply chain reengineering using core process .... Case-based reasoning for production scheduling.