Fizyczne Podstawy Teledetekcji Wykład 6

1. Fizyczne Podstawy TeledetekcjiWykład 6 Krzysztof Markowicz Instytut Geofizyki Uniwersytet Warszawski kmark@igf.fuw.edu.pl

2. Albedo planetarne Ziemi wynosi około 30 % • Albedo samej powierzchni ziemi 7 % • Pozostały wkład do albeda pochodzi od atmosfery (chmury, aerozolowe)

3. Dlaczego kolor oceanu? Jego pomiary dostarczają informacji o koncentracji fitoplanktonu, chlorofilu oraz rozpuszczonych cząstek w wodzie. Jest to istotne z punktu widzenia • Produkcji biologicznej w oceanach (w oceanach zawarta jest około 1/3 części CO2) • Optyki oceanu • Interakcji wiatru i prądów morskich z biologią oceanu • Wpływu antropogenicznego na środowisko oceaniczne

4. Kolor oceanu • Jest zdefiniowany przez spektralna zależności promieniowania wychodzącego z wody (water leaving radiance) • Jest wynikiem absorpcji oraz rozpraszania przez chlorofil, pigment oraz rozpuszczoną materię oceaniczną. • Promieniowanie wychodzące z wody stanowi zaledwie kilka procent promieniowania dochodzącego do satelity. Dlatego ważna role odgrywa tu tzw. poprawka (korekcja) atmosferyczna i kluczowe staje się jej poprawne wyznaczenie.

5. Spektralna zależność promieniowania wychodzącego z wody dla wysokiej (linia kropkowana) oraz niskiej (dolna linia ciągła) koncentracji chlorofilu. Górna linia przedstawia promieniowanie wychodzące z wody, które dochodzi do górnej granicy atmosfery.

6. Widmo absorpcyjne chlorofilu Chlorofil absorbuje znaczna cześć promieniowania czerwonego oraz niebieskiego tym samym rozpraszając światło zielone.

7. Dygresja

8. Albedo powierzchni oceanu -odbicie Fresnela od płaskiej powierzchni wody i kat zenitalny promieniowania padającego, t kat zenitalny promieniowania załamanego. n – współczynnik załamania światła np. dla wody ok. 1.33 i i t

9. Współczynniki odbicia: Rp współ. odbicia (reflektancja) dla wektora elektrycznego równoległego do płaszczyzny padania (polaryzacja pionowa), Rs - współ. odbicia dla wektora elektrycznego prostopadłego do płaszczyzny padania (polaryzacja pozioma), Współczynnik odbicia dla promieniowania niespolaryzowanego

10. Dla wody współczynnik odbicia dla promieniowania padającego z kierunku zenitalnego wynosi; Rnormal=0.022 (dla obszaru widzialnego) Dla małych kątów elewacyjnych i 90o współczynnik odbicia R 1 Z wyjątkiem i =0 współczynnik odbicia pionowej składowej jest znacznie mniejsza niż składowej horyzontalnej. Dla i=0o

11. Kąt Brewstera Rp=0 • To kąt dla którego współczynnika odbicia składowej prostopadłej znika. w przepadku wody kąt ten wynosi ok. 53o

13. Odstępstwa współczynnika odbicia wody od wzorów Fresnela • Falowanie oceanu - powierzchnia wody nie jest plaska. Teoria Coxa-Munka z 1956 opisuje kształt powierzchni wody w zależności od kierunki i prędkości wiatru. • Załamywanie grzbietów fal i pojawianie się piany prowadzi do wzrostu albeda powierzchni oceanu. Przy prędkości wiatru 10 m/s piana pokrywa około 4% powierzchni, przy 15 m/s około 10-20% . Albedo samej piany zmienia się od 0.6-0.9. • Dostające się do wody bąble powietrza z jednej strony zmieniają własności optyczne wody z drugiej zaś w czasie ich pękania emitowana jest do atmosfery sól morska oraz ładunki elektryczne.

14. Elementy Teorii Coxa-Munka • Statystyczne podejście do problemu kształtu powierzchni oceanu. • Na podstawie fotografii sfalowanej powierzchni wody Cox i Munk stwierdzili, iż rozkład prawdopodobieństwa nachyleń powierzchni wody niewiele różni się od rozkładu gaussowskiego. • Jedna z najprostszych formuł określających rozkład prawdopodobieństwa nachylania powierzchni wody ma postać: gdzie w jest prędkością wiatru [m/s] Bardziej zaawansowane rozkłady uwzględniają zależność szerokości rozkładu od kierunku wiatru i mają postać rozbudowanych wzorów empirycznych

15. Odblask - Sun Glint Zjawisko związane odbiciem zwierciadlanym (Fresnelowskim) od płaskiej powierzchni wody. Zależy ono od wzajemnego położenia Słońca i satelity oraz stanu morza.

16. Albedo oceanu a albedo powierzchni oceanu. • Ze względu na promieniowanie wychodzące w wody albedo oceanu jest inną wielkością niż albedo powierzchni oceanu. • Chlorofil i zawarte w wodzie substancje zwiększają albedo oceanu. • Ponadto zmieniają jego spektralną zależność – kolor oceanu

17. Pomiary koncentracji chlorofilu – pomiary naziemneprzy użyciu przyrządu SIMBAD. • Pomiar promieniowania bezpośredniego • Pomiar promieniowania odbitego od powierzchni oceanu dla kąta odbicia około 45o (blisko kąta Brewstera) • Pomiary w wąskich (10 nm) kanałach spektralnych: 443, 490, 560, 670, 870 nm.

18. Definiujemy wielkości zwaną Remote Sensing Reflectance (Reflektancja zdalna) • Lw(z=a,,,) – water leaving radiance (radiancja promieniowania wychodzącego z wody) z kierunku (,) na wysokości a=z po wyżej powierzchni wody • F - strumień promieniowana słonecznego • Remote Sensing Reflectance (Reflektancja zdalna) określa jaka cześć promieniowania padającego na powierzchnie wody powraca do atmosfery pod katem (,)

19. Znormalizowana radiancja promieniowania wychodzącego z wody zdefiniowana jest jako Zakładając ze powierzchnia jest typu Lamberta współ. odbicia związany jest z znormalizowaną reflektancją zdalnąrelacja: Wielkość silnie skorelowany z koncentracja chlorofilu w wodzie

20. Koncentracja chlorofilu można wyznaczyć z następującego wzoru empirycznego C wyrażone jest w [mg/m3] Aby wyznaczyć Rw()Nmusimy znać wpływ atmosfery na promieniowanie odbite od wody. W przypadku pomiarów naziemnych sytuacja jest znacznie prostsza w porównaniu z pomiarami satelitarnymi

21. LSKY Lmw i Lwater w Spektralne pomiary reflektancji zdalnej Mierzone wielkości: LSKY– radiancja nieba Lmw – radiancja od powierzchni wody Uwzględniając efekt odbicia promieniowania nieba mamy: RSKY – Fresnelowski współczynnik odbicia od wody

22. Opisana powyżej metoda jest poprawna dla płaskiego oceanu. W innym przepadku będzie prowadzić do istotnych błędów związanych z niewłaściwym oszacowaniem odbicia związanego z promieniowaniem nieba. • Redukcja promieniowania nieba pojawiająca się w wielkości Lw może być efektywnie uzyskana poprzez pomiar poziomej składowej promieniowania odbitego dla kąta Brewstera. • Ustawienie polaryzatora w płaszczyźnie poziomej sprawia, że efekt odbicia Fresnelowskiego jest minimalny. Wkład od promieniowania nieba związany jest tylko ze słaba składową pionowa pola elektrycznego związaną z pofalowaną powierzchnią oceanu.

23. Spektralna zmienność RSR

25. Poprawka atmosferyczna • Radiancja obserwowana przez detektor umieszczony na satelicie: Ipath() - radiancja promieniowania rozproszonego Ig() - radiancja związana z odblaskiem słonecznym (sun glint) Iwc() - radiancja związana z załamywaniem się fal (white caps), odbicie promieniowania całkowitego Iw() - water leaving radiance T() – całkowita transmisja atmosferyczna

26. Przyczynki do radiancji na górnej granicy atmosfery.

27. Analogiczne równanie możemy zapisać dla współczynnika odbicia R • Następujące wyrażenia musza być właściwie wyznaczone: • Rpath(), Tdir ()Rg (), T ()Rwc (), T ()Rw () • Tdir()Rg() usuwa się zakładając maskę na odblask słoneczny (nie rozpatruje się obszaru nim objętego) • Rwc() szacuje się z wzoru empirycznego: gdzie V jest prędkością wiatru w [m/s] na wysokości 10m 3) Rr i Ra współ. odbicia związane z rozpraszaniem Rayleigha oraz z aerozolami w przybliżeniu pojedynczego rozpraszania

28. Rra– współ. odbicia związany z wielokrotnym rozpraszaniem na aerozolach o molekułach powietrza. 4) W bliskiej podczerwieni promieniowanie wychodzące z wody jest zaniedbywanie małe a wiec: Rw(750 nm)=0, Rw(850 nm)=0 5) Znikające wielkości Lw w bliskiej podczerwieni pozwalają obliczyć poprawkę atmosferyczną. Różnica pomiędzy obliczona i mierzona na satelicie radiancja określa promieniowanie wychodzące z wody która dociera do górnej granicy atmosfery.

29. Przybliżenie pojedynczego rozpraszania Rpath()=Rr()+Ra() Rr– oblicza się z teorii Rayleigha dla danego ciśnienia na powierzchni ziemi z uwzględnieniem funkcji fazowej Ra(750nm) oraz Ra(750nm) są wyznaczane na podstawie pomiarów satelitarnych. Definiujemy iloraz: Zakładamy ponadto, iż atmosfera jest płasko równoległa i składa się z dwóch warstw; górnej gdzie mamy tylko molekuły i dolnej zawierającej cały aerozol.

30. Ostatecznie pozostaje wyznaczenie wielkości (i,865) dla =443 oraz 550 nm na podstawie pomiarów (750,865) • Po założeniu jednego z typów aerozolu (morski, troposferyczny, miejski, wiejski, pustynny) wyznacza się poszukiwana wielkość (i,865). • Wybór modelu aerozolu zależny jest od typu przyrządu i dostępnych kanałów spektralnych. To zagadnienie będzie analizowane w czasie kolejnych wykładów.

31. Kluczowe czynniki dla poprawki atmosferycznej • Rozproszenie wielokrotne • Wpływ aerozoli absorbujących (potrzeba znajomości profilu aerozolu) • Istnienie aerozoli stratosferycznych • Obecność chmur cirrus • Polaryzacja promieniowania • Współczynnik odbicia oceanu związany z jego falowaniem

32. Detektory satelitarne • CZCS (Coastal Zone Color Scanner) na satelicie NIMBUS 7, 1978-1986 • OCTS (Ocean Color and Temperature Scanner) na ADEOS 1996-1997 • SeaWIFS (Sea-viewing Field of View Scanner) na satelicie Orbitview-2 od 1997 roku. • MODIS (Moderate Resolution Imaging Spectroradiometer) na satelicie Terra od 1999 oraz Aqua 2002

33. CZCS

34. SeaWIFS

35. Poprawka atmosferyczna dla przyrządów MODIS oraz SeaWIFS jest wyznaczana na podstawie pomiarów własności optyczne aerozoli np. grubość optyczna aerozolu.

38. Cyrkulacja powierzchniowa oraz koncentracja chlorofilu oparta na danych z MODIS-a oraz SeaWIFS-a

39. Chlorofil a cyrkulacja oceaniczna • Produkcja fitiplanktionu związana jest z dwoma czynnikami: światłem słonecznym oraz substancjami odżywczymi (nutrients) • Za substancje odżywcze odpowiedzialne jest zjawisko wypływania chłodnych wód oceanicznych na powierzchnie (upwelling) Fitoplankton rozwija się w obszarach upwelling (stref wypływania chłodnych wód głębinowych). Zaliczany do nich: związki azotu oraz fosforan • Produktywność fitoplanktonu związana jest z ilością węgla organicznego zawartego w wodzie wyprodukowanego w czasie fotosyntezy. • Komórki fitoplanktonu zawierają chlorofil więc obszary o wysokiej produktywności fitoplanktonu odznaczają się wysoką koncentracja chlorofilu.

40. Upwellingzwiązane są z cyrkulacją Ekamana w warstwie granicznej oceanu. Jest procesem wznoszenia się chłodnych wód docierających do powierzchni oceanu. Powstaje w : • rejonie zachodnich wybrzeży kontynentów • strefach dywergencji prądów morskich oraz • rejonie równikowym • Upwellingjest bardzo dobrze widoczny na mapach temperatury powierzchni oceanu (SST) jako obszary o niższej temperaturze znajdujące się z reguły blisko linii brzegowej

41. Cyrkulacja Ekamana w wodach przypowierzchniowych Średni transport masy pod kątem 90 stopni w prawo (półkula północna) w lewo (półkula południowa) do kierunku wiatru

42. Upwelling dywergencyjny zachodnich wybrzeży równikowy

43. Upwelling równikowy powstający na prądach równikowych

45. Typowa cyrkulacja Walkera Cyrkulacja Walkera w czasie Al. Nino

46. Koncentracja chlorofilu, SeaWIFS Sep 2004

47. Dwa typy wód morskich • Klasyfikacja na podstawie stężenie chlorofilu w wodzie, który warunkuje produkcję materii organicznej i tym samym wpływa na stężenie innych składników wody odziaływujących ze światłem. • 1) wody pierwszego rodzaju to głównie otwarte wody oceaniczne, w których składniki optycznie aktywne (oprócz soli morskiej i samej wody) są autogeniczne czyli powstają lokalnie w wyniku funkcjonowania ekosystemów morskich. Wody charakteryzują się niskimi stężeniami chlorofilu (0.01-1 mg/m3). • 2) wody drugiego rodzaju To najczęściej wody przybrzeżne w których znaczną część substancji aktywnych optycznie jest pochodzenia alogenicznego tj. z dopływów ze źródeł zewnętrznych (transportowane z rzek i zanieczyszczonej atmosfery). Stężenie chlorofilu jest wysokie (1-30 mg/m3)