1 Krzysztof Markowicz [email protected]Metody teledetekcyjne w badaniach atmosfery. Wykład 15. Techniki fotometryczne Krzysztof Markowicz

2 Fotometry słoneczne CIMEL, PREDE (POM-01, POM-02) PGS-100 MicrotopsMFR-7

3 Pomiary promieniowania rozproszonegoPromieniowanie rozproszone docierające do powierzchni ziemi czy też do górnych granic atmosfery zależy silnie od własności optycznych aerozoli W przeciwieństwie do promieniowania bezpośredniego promieniowanie rozproszone rośnie ze wzrostem grubości optycznej aerozoli. Wykorzystanie promieniowania rozproszonego w metodach teledetekcyjnych jest jednak znacznie trudniejsze, gdyż wymaga rozwiązania pełnego równania transferu w atmosferze. Aby to zrobić musimy dokonać wielu założeń np. założyć profil pionowy parametrów stanu atmosfery czy fizyczno-optycznych własności aerozoli. Jeśli w pierwszym przypadku założenie standardowych profili termodynamicznych jest uzasadnione to w przypadku aerozolu już tak nie jest.

4 W przypadku pomiarów radiancji nieba nie jesteśmy wstanie określić dystrybucji aerozolu a jedynie szacować wielkości uśrednione w pionowej kolumnie powietrza. Rozkład radiancji dla małych grubości optycznych możemy określać przy użyciu przybliżenia pojedynczego rozpraszania. W tym przypadku wzór na promieniowanie rozproszone ma analityczna postać o =o gdzie wielkości optyczne takie jak funkcja fazowa P czy albedo pojedynczego rozpraszania są wartościami uśrednionymi w pionowej kolumnie atmosfery.

5 gdzie  oraz ó są katami zenitalnymi satelity oraz Słońca zaś  względnym kątem azymutalnym Słońca oraz satelity. Promieniowanie rozproszone docierające do powierzchni ziemi zależy wiec od grubości optycznej atmosfery, albeda pojedynczego rozpraszania oraz funkcji fazowej. Zauważmy, że grubość optyczna jest w przybliżeniu suma grubości optycznej związanej z rozpraszaniem Rayleigha oraz z rozpraszaniem i absorpcją na aerozolu aerozolu. Albedo pojedynczego rozpraszania jest zaś średnią wartością rozpraszania molekularnego oraz aerozolowego

6 Ponieważ ok. 90% całego aerozolu znajduje się warstwie granicznej (od powierzchni ziemi do ok. 2 km) zaś ok. 80% rozpraszania molekularnego odbywa się powyżej tej warstwy to radiancję promieniowania rozproszonego możemy zapisać przy pomocy przybliżonego wzoru gdzie IR , IA jest radiancja związaną odpowiednio z rozpraszaniem na molekułach powietrza oraz aerozolach. Radiancja związana z rozpraszaniem Rayleigha ma postać gdzie wykorzystano ponownie przybliżenie pojedynczego rozpraszania dla kata zenitalnego Słońca i rozpraszania Rayleigha

7 Zauważmy że przybliżenie pojedynczego rozpraszania ma dwa rozwiązania 1)  o , przypadek tzw. principal plane 2) = o , przypadek almucantar plane W przyrządach pomiarowych wykonuje się dwa różne skanowania obszaru nieba. Są one wykonywane w: płaszczyźnie horyzontalnej (stały kąt zenitalny) płaszczyźnie prostopadłej (stały kąt azymutalny) W przypadku płaszczyzny horyzontalnej radiancja docierająca do powierzchni ziemi ma postać . Po przekształceniach mamy

8 Dla małych grubości optycznych aerozolu funkcja fazowa jest proporcjonalna do różnicy pomiędzy radiancją nieboskłonu a radiancją rozpraszania molekularnego Lewa strona wcześniejszego równania jest iloczynem dwóch poszukiwanych wielkości, a wiec nie można wyznaczyć każdej z ich na podstawie tego równania Zauważmy jednak, że gdy mamy pomiary radiancji dla dwóch kątów azymutalnych jedną ze zmiennych możemy uprościć Wielkości po prawej stronie równania są znane gdyż grubość optyczna aerozolu jest obliczana na podstawie pomiarów promieniowania bezpośredniego, zaś IR jest obliczane. Ostatecznie więc

9 Jednak równanie to zawiera dalej dwie niewiadomeJednak równanie to zawiera dalej dwie niewiadome. Jeśli jednak założyć postać funkcji fazowej to problem upraszcza się. Dla przykładu niech ma ona postać funkcji Henyey-Greenstaina Wówczas otrzymujemy Rozwiązując to równanie kwadratowe ze względu na g obliczamy parametr asymetrii a następnie wyznaczamy albedo pojedynczego rozpraszania. Ostatecznie otrzymujemy wszystkie trzy parametry optyczne dla danej długości fali.

10 Definiujemy wielkości:Z reguły mamy do dyspozycji pomiary dla większej liczby kątów niż tylko dwa co umożliwia dokładniejsze wyznaczenie parametru asymetrii g czy też funkcji fazowej. Definiujemy wielkości: Wykreślamy funkcje Fij względem różnicy kąta rozpraszania a następnie dofitowujemy funkcję: Otrzymuje się w ten sposób znacznie lepsze oszacowanie na parametr asymetrii. Niestety funkcja fazowa H-G nie jest dobrym przybliżeniem dla dużych aerozoli ze względu na silną anizotropię rozpraszania do przodu.

11 Inna metoda opiera się na bezpośrednim wykorzystaniu funkcji FijInna metoda opiera się na bezpośrednim wykorzystaniu funkcji Fij. Kształt ten funkcji odpowiada funkcji fazowej. Po znormalizowaniu jej do 4 otrzymujemy funkcję fazową. Znacznym ograniczeniem tej metody jest pomiarów blisko Słońca, a co za tym idzie dla małych kątów rozpraszania. Zmienność funkcji fazowej od 0 do 5o jest duża i nie może być ona bezpośrednio mierzona. Podeście modelowe Na podstawie pomiarów prom. bezpośredniego oblicza się AOT oraz wykładnika Angstroma a następnie zgaduje się model optycznych własności aerozolu. Używając go obliczamy iloczyn P() i porównujemy z pomiarami. Minimalizujemy różnicę poprzez zmianę własności optycznych aerozolu (zmianę wyboru modelu).

12 Wychodząc z równania transferu w postaci:Dotychczasowe rozważania odnosiły się jedynie do przybliżenia pojedynczego rozpraszania, którego często we względu na wysoką wartość AOT nie można stosować. Uwzględnienie rozproszeń wyższego rzędu jest możliwe jednak zasadniczo komplikuje wcześniejsze wzory. Wychodząc z równania transferu w postaci: Radiancję można zapisać w postaci szeregu Szereg ten jest szybko zbieżny tylko dla <<1 co raczej nie ma miejsce w atmosferze. Zaletą tego podejścia jest analityczna postać rozwiązań, które niestety szybko komplikują się wraz ze stopniem kolejnych rozproszeń

13 Aby uniknąć problemu zbieżności korzysta się z modeli transferu promieniowania w atmosferze.Jest to jednak kosztowne obliczeniowo dlatego aby nie używać za każdym razem modelu konstruuje się specjalne tablice wynikowe (lookup table). Zawierają one informacje o radiancji nieboskłonu dla rożnej geometrii oraz różnego typu aerozolu. Wyznaczanie własności optycznych aerozolu sprowadza się do minimalizacji różnicy pomiędzy radiancją obserwowaną oraz obliczaną poprzez wybór odpowiedniego modelu aerozolu.

14 AERONET- sieć obserwacji aerozolowychPoziomy danych w AERONECIE Level 1.0 wstępnie przetworzone dane Level 1.5 dane po odrzuceniu chmur Level 2.0 ostateczna wersja uwzględniająca poprawki kalibracyjne i manualne odrzucenie chmur.

15 CIMEL – fotometr słoneczny

16

17

18

19

20