Natural Variability of the Land Surface - Phenology

1
Natural Variability of the Land Surface -
Phenology

• Lecture 12
• CLIM 714
• Paul Dirmeyer

2
Climate Drivers for Natural Land Cover

• Climatic Constraints
• Temperature
• annual and seasonal means, extremes, timing of
  first frost, ice free days, growing degree days
• Moisture
• annual and seasonal means, extreme events (floods
  and droughts), precipitation, actual and
  potential evapotranspiration

3
Regional Patterns

• Climate Zones
• Köppen (1928) - map coauthored by student Geiger
  (P T)
• Holdridges Life Zones (P T)
• USDA Hardiness Zones (minimum T)
• Potential Natural Vegetation
• climax vegetation
• A.W. Kuchler (1964)
• EPA and other groups have updated

4
Köppen-Geiger Climate Classification System
5
Köppen Classification Scheme
6
Holdridges Life Zones
7
Predicting current vegetation
Henderson-Sellers (1990) used a GCM and
Holdridges Life Zones to estimate the
distribution of vegetation that would be
consistent with GCM climate.
8
USDA Plant Hardiness
Gardeners in the U.S. are familiar with the USDA
classification scheme, which is based solely on
minimum temperature extremes (assuming
precipitation is irrelevant because gardeners can
water their plants during dry weather).
9
Geographic Controls on Land Cover

• Geomorphology
• Landforms are a product of the interaction of
  geology and climate
• Topography can modify climate
• Soils develop as a result of weathering of
  geologic substrates

10
Geologic map modified from King and Biekman (1974)
http//tapestry.usgs.gov/
11
Terrain of the Conterminous United States shaded
relief map
Dry adiabatic lapse rate Ga g / cp 9.7 K /
1000 m Observed lapse rate G 6.5 K / 1000
m This correction can be used to estimate
local temperatures near observation stations that
are at a different altitude.
12
A Statistical Solution A Neural Data-Driven
Assessment of Global Vegetation Classes
Global distribution of feature types obtained
after simulation with a SOM (Self-Organizing Map)
and the topological arrangement of the categories
on the network. As training data the monthly
means of temperature, precipitation, insolation
and the water storage capacity of soils are used.
(Kropp 1999)
Neural nets are black boxes. They dont tell
why or how.
13
Based on observed 1971-2000 climatology, along
with soils and other information (e.g. native
species habitats).
14
Monitoring Vegetation Variations

• Satellite are appropriate to a range of dynamic
  monitoring tasks
• monitoring vegetation dynamics over course of a
  year
• link to (crop) growth models to provide yield
  estimates
• distinguish cover types (classification)

15
Issues

• temporal sampling
• reconcile requirements of monitoring task with
  sensor characteristics and external influences
• repeat cycle of sensor
• spatial resolution of sensor
• lifespan of mission / historical data
• cloud cover aerosol effects on optical /
  thermal data

16
Issues

• discriminating surface changes from external
  influences on remote sensing data
• Viewing and illumination conditions can change
  over time
• Viewing
• wide field of view sensors
• pointable sensors
• Illumination
• variations in Sun position (visible and nearIR
  channels)
• variations in atmospheric conditions

17
Issues

• cloud cover

Composite image from U. Wisconsin SSEC
18
Issues

• Sensor calibration
• degradation over time
• variations between instruments
• Co-registration of data
• effects of mis-registration (practical)
• Quantity of data
• volume of data can be very large
• preprocessing requirements can be very large
• move towards formation of databases of
  remote-sensing derived 'products' (e.g., EOS
  Earth Observing System)

19
Mean Annual Cycles - Snow
20
Mean Annual Cycles Soil Wetness
Magnitude of annual cycle (top) and standard
deviation of April soil wetness (bottom).
21
Albedo movie
22
PAR movie
23
LAI movie
24
Green-up and climate
Is the green-up causing a cool spell, or is this
an artifact of weather variability?
25
Dates of start and length of growing season
26
Climate regimes reflected in growing seasons
27
First green-up
DMA delayed moving average Detects the time of
the first increase in NDVI from base
level. Year-to-year variations are caused by
climate variations (mainly temperature anomalies,
but many factors play a role, including
precipitation/soil moisture, and solar radiation).
1995
1996
Schwartz et al. (2002 Int. J. Climatol.)
28
Peak greening
SMN Seasonal Midpoint NDVI Detects the time of
the crossing of NDVI over the (maxmin)/2 value.
1995
Satellite NDVI measurements
Remove cloud-contaminated points
Fit a spline curve
Determine midpoint
1996
DMA looks for this first positive trend.
Schwartz et al. (2002 Int. J. Climatol.)
29
Hysteresis means you cant just map weather to
vegetation physical and biological processes
must be modeled
30
North American Monsoon
http//www.ofps.ucar.edu/name/
31
Predictors of North American Monsoon
32
Soil moisture and diurnal cycle
33
Dealing with issues

• Vegetation Indices (VIs)
• measured reflectance / radiance sensitive to
  variations in vegetation amount
• BUT also sensitive to external factors
• want contiguous data (clouds)
• Typically take VI compositing approach
• Assume highest measured VI is actual VI.
• Interpolate across missing data.

34
Use of VIs

• no one ideal VI - NDVI used historically
• empirical relationships will vary spatially and
  temporally
• direct
• attempt to find (empirical) relationship to
  biophysical parameter (e.g. LAI)
• indirect
• look at timing of vegetation events (phenology)
• VI can still be sensitive to external factors
  (Especially bi-directional reflectance
  distribution function (BRDF) effects)

35
BRDF explained
BRDF gives the reflectance of a target as a
function of illumination geometry and viewing
geometry. The BRDF depends on wavelength and is
determined by the structural and optical
properties of the surface, such as
shadow-casting, multiple scattering, mutual
shadowing, transmission, reflection, absorption
and emission by surface elements, facet
orientation distribution and facet
density. However, it should not be overlooked
that the BRDF simply describes what we all
observe every day that objects look differently
when viewed from different angles, and when
illuminated from different directions.
36
Real-world examples of BRDF variations
This is a black spruce forest in the BOREAS
experimental region in Canada. Left
backscattering (sun behind observer), note the
bright region (hotspot) where all shadows are
hidden. Right forward-scattering (sun opposite
observer), note the shadowed centers of trees and
transmission of light through the edges of the
canopies.
A barren field with rough surface Left
backscattering (sun behind observer), note the
bright region (hotspot) where all shadows are
hidden. Right forward-scattering (sun opposite
observer), note the specular reflection.
37
VI Issues

• IDEAL
• Attempt to make VI sensitive to vegetation amount
  but not to external factors
• atmospheric variations
• topographic effects
• BRDF effects (view and illumination)
• soil background effects
• SAVI, ARVI etc.
• PRACTICE
• VIs maintain some sensitivity to external factors
  
• Be wary of variations in satellite calibration
  etc. for time series

38
VI Issues
39
VI Issues
40
Examples/Techniques

• land cover change detection
• Vegetation Indices eg
• change in VI - infer change in vegetation state
• NDVI variation in Mozambique (UN World Food
  Programme)

41
(No Transcript)
42
Variance of Boreal Summer LAI over 8 years
(1987-1994)
43
Mean annual NPP (1981-2000) estimated with a DVM
at 8km resolution
NPP (gC m-2 yr-1)
44
Inter-annual variation in global total NPP(Gt C/
yr)
Year
45
Interannual trends in mean land NDVI and in total
land NPP
Annual growth rate K, significance 95, 99.
1991-92 omitted because of the Mount Pinatubo
eruption.
Trends in mean land NDVI
Trends in total land NPP
46
GLO-PEM estimate of changes in annual NPP
1982-2000
47
Temporal changes in NPP for the globe and each
hemisphere 1981-2000.
Global
MEI is an indicator of the intensity of El Niño
(ve) and La Niña (-ve). The 10-day anomaly is
the deseasonalized change in NPP. NPP anomalies
in g C m-2 per 10 days.
Interannual variability and trend in global
terrestrial net primaryproductivity satellite
analysis 1980-2000 M. Cao, S.D. Prince, J. Small,
S.J. Goetz Department of Geography, University of
Maryland
Northern hemisphere
Southern hemisphere
48
The responses of NPP to ENSO
Difference between NPP for the stated year and
mean value for whole period 1981-2000.
Anomalies. Yr. x-(mean of 1981-2000) gCm-2 per
10 d
La Niña - 1998-99
El Niño - 1986-87
Normal - 1995-96
El Niño 1993-94
El Niño - 1982-83
El Niño 1997-98
49
AFRICAWhole continent
50
NPP vs ENSO in Africa
Interannual Variation in NPP and ENSO Cycles for
the Whole African Continent
51
El Niño
1997
1987
1983
Normal
Negative anomaly
1996
1990
No change
La Niña
Positive anomaly
1989
1999
52
Coefficient of Variation and Mean Annual NPP 1982
- 1999
Low High
Low High
Coefficient of variation
Mean Annual NPP
53
Seasonal and interannual variability
One of the most important outstanding aspects of
land-atmosphere interactions is the feedbacks
associated with seasonal and interannual
vegetation cycles (Hutjes, 1998).
54
Areas exhibiting model sensitivity to
phenologyNumber of months per yearLAI-Phen
statistically different from LAI-Mean
In GLACE, HadAM ranked 11th out of 12 models for
sensitivity of precipitation to land surface.
HadAM3MOSES2
55
Models of vegetation