Metodologia przeglądu: protokół wyszukiwania, kryteria kwalifikowalności i przesiewanie
Systematyczny protokół oparty na wytycznych PRISMA regulował proces identyfikacji, przesiewania i wyboru13. Przeszukiwania objęły bazy Scopus, Web of Science Core Collection, ScienceDirect, Wiley Online Library, MDPI i Google Scholar (pierwsze 200 wyników uporządkowanych według trafności), repozytoria HAL, IRD/Horizon, WASCAL i FAO/AQUASTAT oraz regionalne archiwa uniwersyteckie, w języku francuskim i angielskim; listy referencyjne włączonych rekordów zostały poddane analizie wstecznej. Zastosowano dwa łańcuchy operatorów Boole'a, jeden tematyczny i jeden zorientowany na metodologię; oba zostały przedstawione w całości, wraz z adaptacjami do poszczególnych platform oraz liczbą wyników dla każdego źródła, w Tabeli 1. Zakres czasowy objął lata od 1979 do 2026, przy czym wstępne wyszukiwanie przeprowadzono 12 stycznia 2026 roku, a końcową aktualizację 17 kwietnia 2026 roku.
Zapisy kwalifikowano do analizy tylko wtedy, gdy spełniały wszystkie cztery kryteria: dotyczyły zlewni w Afryce Zachodniej lub metody wyraźnie możliwej do przeniesienia na warunki Afryki Zachodniej, gdzie brakuje danych; zawierały jawne, a nie przypadkowe ujęcie zmian LULC lub urbanizacji; opisem metody była procedura możliwa do odtworzenia, obejmująca źródła danych, model i strategię kalibracji; oraz posiadały status publikacji recenzowanej lub możliwej do prześledzenia literatury szarej. Odrzucano zapisy o charakterze czysto koncepcyjnym, ograniczone wyłącznie do jakości wody lub ekologii, pochodzące spoza regionu, pozbawione przenoszalnej treści metodologicznej lub niedostępne w pełnym tekście. Tytuły i streszczenia były niezależnie przesiewane przez dwóch autorów, a ewentualne rozbieżności rozstrzygał trzeci autor. Ekstrakcja danych obejmowała zlewnię i kraj, powierzchnię, okres, rodzinę metodologiczną, model lub algorytm, dane wejściowe, zmienne kalibracyjne i walidacyjne, raportowane metryki, zmiany hydrologiczne oraz ujęcie niepewności. Liczba rekordów na każdym etapie została przedstawiona na Rysunku 1, a wyekstrahowane dowody w Tabeli 2. Badania pogrupowano w osiem kategorii analitycznych, przypisując każdą do rodziny, która wniosła główny wkład; rodziny wtórne odnotowano osobno, aby nie pominąć żadnego wkładu.
Ilościowe charakterystyki bazy dowodowej
Korpus obejmuje 51 rekordów, z czego 35 odnosi się do Afryki Zachodniej, a 16 stanowi odniesienia metodologiczne lub porównawcze z innych regionów. Wśród 35 rekordów rozkład geograficzny jest nierównomierny: Benin 8 rekordów (23%), Burkina Faso i Niger po 5 (po 14% każdy), Senegal 3 (9%), Ghana i Wybrzeże Kości Słoniowej po 2 (po 6% każdy), transgraniczne dorzecza Mono i Benue po 1 (po 3% każdy) oraz 8 badań regionalnych (23%) bez przypisania do konkretnego państwa. Mali, Gwinea, Sierra Leone, Liberia, Togo i Gambia nie dostarczyły żadnego dedykowanego badania zlewni. Publikacje są aktualne: 2 rekordy pochodzą sprzed 2010 roku, 9 powstało między 2010 a 2019 rokiem, a 24 (69%) datowane są od 2020 roku wzwyż (Tabela 3).
Rozkład według rodzin metodologicznych, przedstawiony w Tabeli 4, jest równie nierównomierny: statystyka i uczenie maszynowe 8 z 35 (23%), modelowanie koncepcyjne i półrozproszone 7 (20%), modelowanie rozproszone i zintegrowane 5 (14%), obserwacje in situ 3 (9%), teledetekcja jako podejście podstawowe 3 (9%), hybrydowe łańcuchy scenariuszy 3 (9%), syntezy regionalne 3 (9%) oraz metody zorientowane na wody podziemne 2 (6%). Uwzględniając również wkład pomocniczy, teledetekcja pojawia się w 20 rekordach (57%), a modelowanie hydrologiczne w 15 (43%) – ta para stanowi metodologiczny fundament tej dziedziny.
Dwie cechy warunkują, z jaką pewnością można formułować wnioski regionalne. Eksperymenty kontrolowane są niemal nieobecne jako podejście podstawowe, przy czym tylko 1 zapis opiera się na pomiarach w skali działki15, w związku z czym parametry hydrodynamiczne w modelach rozproszonych są w dużej mierze przenoszone z innych badań, a nie mierzone lokalnie. Korpus nie jest również w pełni niezależny: 8 zapisów (23%) pochodzi z jednej sieci skoncentrowanej wokół Uniwersytetu w Bonn i programu WASCAL16,17,18,19,20,21,22,23, zatem zbieżność między nimi częściowo odzwierciedla wspólne dane, łańcuchy modeli i parametryzację, a nie niezależną replikację.
W odniesieniu do głównych twierdzeń, ustalenie, że rozszerzanie gruntów ornych i degradacja okrywy roślinnej zwiększają spływ powierzchniowy oraz plon wody, jednocześnie redukując ewapotranspirację lub wkład wód podziemnych, jest potwierdzone przez 6 badań18,20,24,25,26,27 przeprowadzonych w 4 krajach i 4 systemach zlewni, ale w ramach zaledwie 3 wzajemnie niezależnych grup. Wzrost spływu w Sahelu związany z degradacją stanu powierzchni jest potwierdzony przez 3 niezależne grupy3,4,15. Twierdzenie, że urbanizacja w szczególności napędza zmiany w stresie wodnym, opiera się na 3 zapisach6,28,29 pochodzących od 3 grup, wszystkie skoncentrowane w Niamey, i jest zatem raportowane tutaj jako prawdopodobne, lecz słabo zreplikowane.
In situ obserwacje i instrumentowane zlewnie
Bezpośrednia obserwacja pozostaje najsolidniejszą podstawą wiązania transformacji krajobrazu ze stresem wodnym. Długie serie pomiarowe z obszaru Sahelu wykazały rolę stanów powierzchni, degradacji gleb oraz powiększania się obszarów odsłoniętej gleby w szybkim spływie obserwowanym w kilku basenach3,4. Pionierskie pomiary w zlewniach miejskich Niamey w latach 60. i pod koniec lat 70. XX wieku udokumentowały współczynniki spływu osiągające kilkadziesiąt procent oraz czasy reakcji wynoszące około godziny, co wiązano z nieprzepuszczalnością powierzchni, sieciami dróg i organizacją drenażu28. W basenach wiejskich i agropastoralnych instrumentacja służy głównie do kalibracji i walidacji modeli. Obszary Dano, Couffo, Ouémé, White Volta i Bonsa dostarczają danych dotyczących opadów, temperatury, wilgotności gleby, przepływu strumieni, a sporadycznie także wód podziemnych, co pozwala na określenie bilansów wodnych, zasilania i spływu16,24,25,26,30, a wspólna analiza opadów i przepływów w północnym Beninie pozwala oddzielić zmienność międzyroczną od trendów długoterminowych23. Tam, gdzie sieć stacji pomiarowych jest rzadka, wielokryterialna walidacja narzędzia Soil and Water Assessment Tool (SWAT) z wykorzystaniem zdalnie wykrywanej ewapotranspiracji, wilgotności gleby oraz zasobów wodnych pochodnych z danych GRACE konsoliduje szacunki bilansu wodnego i redukuje efekt ekwifinalności21.
In situ procesy przechwytywania danych, a nie wskaźniki zastępcze, pozwalają na określenie odpowiedzi zdarzeń, progów, sezonowości oraz opóźnień między opadami, spływem a zasilaniem. Ich ograniczeniami są niska gęstość sieci, heterogeniczność protokołów, przerwane serie, koszty konserwacji oraz trudna ekstrapolacja na większe systemy lub systemy transgraniczne. Tylko 3 z 35 rekordów z Afryki Zachodniej (9%) opierają się głównie na pomiarach in situ, co stanowi wąską podstawę dla całej grupy danych, w której wszystkie pozostałe rekordy zależą od nich w celu walidacji.
Kontrolowane eksperymenty polowe i laboratoryjne
Eksperymenty kontrolowane są mniej rozpowszechnione niż modelowanie, lecz są niezbędne do zrozumienia mechanizmów. Obejmują one testy infiltracji, pomiary gęstości objętościowej i stabilności strukturalnej, nasyconą przewodność hydrauliczną, działki spływowe oraz obserwacje porównawcze na działkach uprawnych, ugorach, sawannach, gołej glebie i powierzchniach zurbanizowanych; w regionie Sahelu wyjaśniają one wzrost spływu poprzez tworzenie się skorupy glebowej, utratę pokrywy zielnej oraz koncentrację przepływu3,4. Pomiary w układzie zagnieżdżonym od działki do zlewni w niewielkiej półpustynnej zlewni w Burkina Faso wykazują systematyczny spadek współczynników spływu wraz ze wzrostem skali pomiaru, co pozwala określić zależność od skali ograniczającą przeniesienie wyników z działek do modeli zlewni 15. W hydrologii miejskiej badania w Afryce Zachodniej pozostają opisowe, a nie eksperymentalne, opierając się na instrumentacji epizodycznej, opartej na poszczególnych zdarzeniach, zamiast na ciągłym monitoringu, co osłabia możliwość przypisania konkretnych skutków urbanizacji9,28. Doświadczenia z ciągle monitorowanych zlewni w fazie rozwoju w innych regionach pokazują, jak bardzo poprawia się precyzja przypisywania przyczyn, gdy monitoring przed i po zabudowie jest zaprojektowany od samego początku31. Takie eksperymenty izolują poszczególne procesy i umożliwiają realistyczną parametryzację, jednak pozostają zależne od lokalizacji i są najużyteczniejsze w połączeniu z ciągłym monitoringiem i modelowaniem.
Satelitarna teledetekcja i wieloczasowa analiza pokrycia terenu
Teledetekcja jest najczęstszym punktem wyjścia w analizie zmian LULC, wykorzystywaną w 20 z 35 zapisów dla Afryki Zachodniej (57%). Dominuje system Landsat ze względu na głębię czasową i bezpłatny dostęp; rozdzielczość Sentinel-2 poprawia dyskryminację klas, a dopełnieniem są radar Sentinel-1, MODIS, ALOS, pochodne cyfrowe modele wysokościowe (DEMs) oraz drony8,22,32. Łańcuchy przetwarzania łączą korekcję atmosferyczną, kompozytowanie sezonowe, indeksy spektralne, klasyfikację nadzorowaną, ocenę dokładności oraz analizę przejść zmian33, przy czym lasy losowe (random forest), maszyny wektorów nośnych (support vector machines) i głębokie sieci zastępują klasyfikatory największej wiarygodności. W dorzeczu Senegalu wykorzystanie danych Landsat wraz z lasami losowymi i wielowarstwowym perceptronem w modelu Markowa pozwoliło na prognozowanie przyszłego użytkowania gruntów w zlewniach Bafing i Falémé32; w dorzeczu Ouémé narzędzie Land Change Modeler pozwoliło ilościowo określić przekształcanie terenów w grunty orne i osiedla22. Platformy chmurowe przyspieszają przetwarzanie serii czasowych8,12.
Druga generacja produktów operacyjnych może być bezpośrednio wykorzystywana jako ograniczenie hydrologiczne, jednak w skali regionalnej wciąż pozostaje ona niewystarczająco eksploatowana. Portal FAO WaPOR dostarcza dane o rzeczywistej ewapotranspiracji dla Afryki co dziesięć dni w rozdzielczościach 250, 100 i 30 m; walidacja kontynentalna w odniesieniu do czternastu stacji kowariancji wirów wykazała ogólny współczynnik R2 na poziomie 0,61 oraz błąd średnio-kwadratowy wynoszący 1,04 mm d-1 34. Produkty dotyczące wilgotności gleby, takie jak SMAP i rejestry ESA Climate Change Initiative, a także produkty opadowe, takie jak CHIRPS, TAMSAT i IMERG, dostarczają zmiennych komplementarnych niezbędnych do wielozmiennej walidacji35,36. Tylko 1 z 35 rejestrów wykorzystuje satelitarny produkt ewapotranspiracji jako ograniczenie kalibracyjne21, co stanowi najwyraźniejszy zidentyfikowany tutaj obszar niewykorzystanego potencjału.
Znane ograniczenia obejmują zachmurzenie w porze deszczowej w strefach gwinejskich i subhumidnych, konfuzję spektralną pomiędzy ugorami, zdegradowanymi sawannami a uprawami oraz niewielką powierzchnię wtórnych obszarów miejskich. Raportowanie jest nierównomierne: ogólna dokładność na poziomie 80–90% i współczynnik kappa Cohena wynoszący około 0,75–0,90 są typowe dla regionalnych badań ujętych w niniejszym zestawieniu, jednak korpus zawiera również analizy wykrywania zmian bez raportowanej oceny dokładności, a dokładności dla poszczególnych klas dla obszarów zabudowanych i gołej gleby, które mają kluczowe znaczenie dla stresu wodnego w miastach, rzadko są wyszczególniane. Niepewności te przenoszą się na szacunki spływu, ewapotranspiracji i zasilania, dlatego walidacja siatkowych produktów opadowych jest krytycznym krokiem w skali regionalnej35,36,37.
Metody zorientowane na wody podziemne: grawimetria, znaczniki i mapowanie zasilania
Wody podziemne stanowią główne źródło zaopatrzenia dla większości obszarów wiejskich i podmiejskich w Afryce Zachodniej, a mimo to są najmniej reprezentowaną składową: tylko 2 z 35 zapisów uznają je za główny obiekt badań, a kalibracja oparta wyłącznie na przepływach strumieni pozostaje normą. Trzy grupy metod pozwalają wypełnić tę lukę. Pierwszą z nich jest grawimetria satelitarna. Misja Gravity Recovery and Climate Experiment (GRACE) oraz jej misja kontynuacyjna dostarczają miesięcznych anomalii całkowitego ziemskiego zapasu wody z rozdzielczością około 300 km, z których anomalie wód podziemnych są izolowane poprzez usunięcie modelowanych składników wilgotności gleby i wód powierzchniowych. Zastosowanie tej metody do warstw wodonośnych w południowym Nigrze pozwoliło na oszacowanie zmian zapasów i zasilania w miejscach, gdzie zapisy piezometryczne są rzadkie38; z kolei w zastosowaniu odwrotnym ten sam sygnał pozwala ograniczyć wewnętrzne składniki bilansu wodnego, których nie można zidentyfikować na podstawie przepływów strumieni21. Niedopasowanie obszaru zasięgu pomiaru do zlewni ogranicza tę metodę do dużych dorzeczy, chyba że zastosuje się skalowanie w dół (downscaling).
Wskaźniki środowiskowe tworzą drugą grupę: stabilne izotopy, bilans masy chlorków oraz, w razie dostępności, tryt lub węgiel promieniotwórczy pozwalają odróżnić zasilanie rozproszone poprzez profil glebowy od zasilania skoncentrowanego poprzez stawy, doliny śródlądowe i obniżenia urbanistyczne; rozróżnienie to ma kluczowe znaczenie dla paradoksu Sahelu oraz oceny, czy urbanizacja zwiększa, czy zmniejsza całkowite zasilanie; zintegrowane projekty izotopowe, teledetekcyjne i modelowe w półpustynnych warstwach wodonośnych strefy przybrzeżnej można bezpośrednio przenieść na ten region39. Po trzecie, mapowanie zasilania łączące teledetekcję z systemami informacji geograficznej jest opcją wymagającą najmniejszej ilości danych, jednak dostarcza ono informacji o potencjale względnym, a nie o wolumetrycznym zasilaniu i wymaga walidacji w oparciu o wydajność otworów wiertniczych12. Zintegrowane modele powierzchniowo-podziemne symulują zasilanie w sposób jawny, przy znacznie wyższych kosztach27,29.
Koncepcyjne i półrozproszone modelowanie hydrologiczne
Modelowanie hydrologiczne jest najpowszechniejszą drogą prowadzącą od dynamiki użytkowania terenu do mierzalnej odpowiedzi hydrologicznej, stanowiąc 15 z 35 rekordów (43%) po połączeniu wszystkich rodzin modeli. Modele koncepcyjne i półrozproszone, takie jak UHP, SWAT oraz model Agricultural Catchments Research Unit (ACRU), równoważą opis procesów, wykonalność kalibracji i wymagania dotyczące danych, a także pozwalają na testowanie scenariuszy użytkowania terenu i klimatu16,18,20,24,25,26,30,40. Efektywność jest raportowana zgodnie z danymi źródłowymi i interpretowana według standardowych wytycznych, według których miesięczny współczynnik efektywności Nasha-Sutcliffe'a (NSE) powyżej 0,50 jest satysfakcjonujący, powyżej 0,65 dobry, a powyżej 0,75 bardzo dobry, przy błędzie systematycznym (percent bias) w granicach ±25% dla przepływu w ciekach14. Cornelissen i in. wykazali na przykładzie tropikalnego basenu w Beninie, że sam wybór modelu jest głównym źródłem niepewności, a dobra wydajność historyczna nie gwarantuje rzetelnej symulacji prospektywnej16.
Studia przypadków określają ilościowo analizowane zmiany. Model SWAT wykazał, że zmniejszenie naturalnej okrywy w dorzeczu Białej Volty obniżyło poziom wód powierzchniowych i przepływ bazowy, zwiększając jednocześnie ewapotranspirację24. W zlewni rzeki Bonsa model ACRU powiązał utratę lasów na poziomie około 39% ze wzrostem przepływu w strumieniu do 37%25. W dorzeczu Couffo w latach 2000–2011 obszary upraw rolnych zwiększyły się o 34%, podczas gdy powierzchnia sawanny, systemów agroleśnych i lasów spadła o 24–60%, poziom wód powierzchniowych wzrósł o 11 mm, a udział wód podziemnych i ewapotranspiracja spadły odpowiednio o 3,2 i 10,6 mm26. Połączenie modelu użytkowania gruntów z Modelem Symulacji Bilansu Wodnego (WaSiM) w dolinach śródlądowych Burkina Faso wykazało rozpiętość modelu klimatycznego od -44% do +95% w całkowitym odpływie, co było wartością większą niż sam sygnał zmian użytkowania gruntów20. Atutem tych badań jest możliwość przeprowadzania eksperymentów numerycznych, co jest niemożliwe w skali całego zlewni. Ich słabością jest kwestia strukturalna: reprezentacja jednostek hydrologicznych, regionalizacja parametrów, jakość numerycznego modelu terenu (DEM) i danych glebowych oraz ekwifinalność kalibracji17,30. Niedobór obserwacji wód podziemnych i wilgotności gleby często wymusza kalibrację opartą wyłącznie na przepływie strumienia, a spośród opracowanych w Tabeli 2 studiów modelowania tylko mniejszość została skalibrowana w oparciu o więcej niż jedną zmienną.
Rozproszone, oparte na fizyce i zintegrowane modelowanie powierzchni i podpowierzchni
Gdy cel wykracza poza sam przepływ strumienia i obejmuje sprzężone zachowanie infiltracji, wilgotności gleby, zwierciadła wody gruntowej, spływu oraz strumieni bocznych, odpowiednie stają się modele zintegrowane i oparte na fizyce20,27,29. W przypadku Dano model WaSiM został wspólnie skalibrowany i zwalidowany pod kątem odpływu, wilgotności gleby i poziomu wód gruntowych, osiągając wartości R2, NSE oraz efektywności Klinga-Gupty (KGE) w przedziale od 0,6 do 0,9 dla wszystkich trzech zmiennych, co stanowi jedną z niewielu dostępnych dla tego regionu autentycznych walidacji wielozmiennych18; model SHETRAN w tym samym zlewni osiągnął wartości od 0,66 do 0,79 wyłącznie dla odpływu19. Wysoko rozdzielczy model HydroGeoSphere zastosowany w półpustynnej zlewni miejskiej w Niamey traktował przepływ powierzchniowy i podpowierzchniowy jako w pełni sprzężone i wykazał, że samo zwiększenie rozdzielczości modelu klimatycznego nie poprawia jego wydajności29.
W delcie rzeki Ouémé model ParFlow-CLM wspólnie wyznaczył spływ, ewapotranspirację, głębokość zwierciadła wody oraz wilgotność gleby, rozszerzając analizę o zmienne nieosiągalne za pomocą prostszych modeli27. Zyskiem jest spójność fizyczna; kosztem są dane wejściowe, moc obliczeniowa, kalibracja i interpretacja. Modele te reprezentują 5 z 35 rekordów (14%), wykazując znaczny potencjał w instrumentowanych basenach pilotażowych oraz obszarach miejskich pod silną presją hydrologiczną41.
Podejścia statystyczne i uczenia maszynowego
Podejścia statystyczne zajmują pośrednią pozycję między obserwacją a modelowaniem opartym na procesach, łącząc metryki krajobrazowe lub klasy pokrycia terenu ze wskaźnikami hydrologicznymi, wykrywając punkty przełomowe i szacując trendy; w dorzeczach Afryki Zachodniej zazwyczaj uzupełniają one, a nie zastępują modele hydrologiczne42. Uczenie maszynowe działa na trzech poziomach. Pierwszym jest mapowanie pokrycia terenu, gdzie lasy losowe (random forest) i pokrewne algorytmy poprawiają klasyfikację i wykrywanie zmian8,22,32. Drugim jest przygotowanie danych: imputacja brakujących serii meteorologicznych za pomocą drzew decyzyjnych, lasów losowych i ekstremalnego wzmacniania gradientu (extreme gradient boosting)43, a także statystyczne downscaling projekcji opadów44. Trzeci poziom, czyli sama prognoza hydrologiczna, jest w regionie najmniej rozwinięty. Sieci długiej krótkotrwałej pamięci (LSTM) trenowane wspólnie na dużych próbkach zlewni przewyższają modele koncepcyjne skalibrowane na zlewni docelowej w regionalizacji metodą leave-one-out, redefiniując prognozowanie w nieopomiarowanych zlewniach (PUB) jako agregację danych, a nie transfer parametrów45; globalny model LSTM typu enkoder-dekoder zrównał się lub przewyższył uznany system prognozowania w nieopomiarowanych zlewniach, przy czym największe zyski odnotowano w regionach o rzadkich danych, w tym w Sahelu46. Żadna z analizowanych prac nie stosuje takich modeli w zlewni Afryki Zachodniej, co sprawia, że regionalne eksperymenty LSTM i PUB na zagregowanych rekordach przepływów rzek Niger, Volta i Senegal stanowią konkretny priorytet, warunkowany dostępnością otwartego regionalnego zbioru danych o cechach zlewni i seriach przepływów, który obecnie nie istnieje. Uczenie maszynowe dobrze radzi sobie z nieliniowością, złożonymi interakcjami i fuzją wielu źródeł danych oraz jest skuteczne w precyzyjnej klasyfikacji obszarów miejskich. Jego interpretowalność jest niższa niż w przypadku modeli opartych na procesach, wydajność zależy od jakości danych treningowych, a ekstrapolacja w warunkach głębokich zmian użytkowania terenu lub zmian klimatycznych pozostaje niepewna12.
Podejścia hybrydowe i sprzężone scenariusze użytkowania terenu i klimatu
Najbardziej uderzającą tendencją ostatniej dekady jest hybrydyzacja: łączenie mapowania wieloczasowego, projekcji przyszłych zmian, skalibrowanych modeli hydrologicznych oraz regionalnych scenariuszy klimatycznych w taki sposób, aby można było oddzielić efekty klimatyczne, efekty użytkowania gruntów oraz ich wzajemne oddziaływania10,18,20,47. Yira i inni wykorzystali takie sprzężenie w Dano do oceny wydajności wodnej i przepływu strumieni w ramach pięciu scenariuszy użytkowania gruntów18. W dorzeczu Lobo model Land Change Modeler, modele CORDEX oraz CEQUEAU pozwoliły na ilościowe określenie wpływu na dopływy do zbiorników, co zilustrowało operacyjną użyteczność tych łańcuchów w zarządzaniu zbiornikami i zaopatrzeniem miejskim47. Model Water Evaluation and Planning (WEAP) został zastosowany w dorzeczu Mono do analizy wpływu klimatu, LULC oraz rozwoju na zasoby wodne i energetykę wodną48.
Hybrydowe systemy modelowania najlepiej nadają się do wspierania procesów podejmowania decyzji, ponieważ odpowiadają na konkretne pytania: co stanie się z odpływem powierzchniowym, jeśli urbanizacja będzie postępować, lub co stanie się z zasilaniem wód podziemnych, jeśli sawanny zostaną przekształcone w grunty orne. Ich wadą jest kumulacja niepewności wynikająca z klasyfikacji, scenariuszy przejść, modeli klimatycznych oraz struktury hydrologicznej, co wymaga przeprowadzenia analiz wrażliwości i zastosowania ansambli wielomodelowych17,40. W przypadkach, gdy raportowane są zakresy rozbieżności, komponent modelu klimatycznego zazwyczaj dominuje nad komponentem użytkowania terenu20, zatem średnia z ansamblu przedstawiona bez wskazania rozbieżności jest myląca z punktu widzenia planowania.
Analiza porównawcza podejść
Porównanie metod musi uwzględniać różnorodność pytań naukowych: niektóre techniki służą do pomiaru, inne do mapowania, wyjaśniania lub symulacji. Nie istnieje jedna najlepsza metoda, a jedynie zestawy narzędzi w mniejszym lub większym stopniu dostosowane do problemu, skali, danych oraz badanej formy stresu wodnego (Rysunek 2). Aby uniknąć etykiet jakościowych, porównanie w Tabeli 5 zostało w miarę możliwości wyrażone za pomocą raportowanych wskaźników: NSE, KGE, R2 oraz błędu procentowego (percent bias) dla symulacji, a także ogólnej dokładności i kapy Cohena dla klasyfikacji, interpretowanych w odniesieniu do ustalonych progów14. Bezpośrednia obserwacja jest najwiarygodniejsza w zakresie tego, co mierzy, jednak jej wiarygodność terytorialna zależy od reprezentatywności sieci, a współczynniki spływu w skali działek nie mogą być przeniesione na skalę zlewni15. Modele dostarczają wglądu w rozkład przestrzenny, przy czym kalibracja i walidacja w Afryce Zachodniej dla NSE, KGE i R2 mieszczą się głównie w zakresie 0,6–0,9, co według standardowych kryteriów oznacza wynik od dobrego do bardzo dobrego14,18,19; jednak w większości przypadków osiągnięto to dla przepływu strumienia, a walidacja wielozmienna wykazuje, że porównywalne wyniki przepływu mogą współistnieć z wyraźnie różnymi wewnętrznymi bilansami wodnymi16,21. Teledetekcja skutecznie wykrywa przejścia, osiągając dokładność na poziomie 80–90% i wartość kapy bliską 0,75–0,90 przy połączeniu serii wieloźródłowych, walidacji gruntowej i odpornych algorytmów, natomiast jej skuteczność spada, gdy klasy są bliskie spektralnie33. Wrażliwość modeli na sformułowanie ewapotranspiracji potencjalnej, wykazana dla 21 metod w dorzeczu rzeki Senegal, nakazuje zachować ostrożność przy porównywaniu scenariuszy49.
Najlepszy kompromis między kosztem a zakresem pokrycia obejmuje darmowe obrazy, otwarte zestawy danych klimatycznych, oszczędny model oraz minimalny zestaw obserwacji przepływu strumieni, co wyjaśnia regionalny sukces modelu SWAT z klasyfikacjami Landsat i Sentinel8,24,26,30. Gęste sieci i zintegrowane modele pozostają kosztowne, choć najlepiej nadają się do reprezentowania mechanizmów kaskadowych, takich jak rozproszona urbanizacja, zmniejszone zasilanie oraz przekształcone tereny podmokłe27,29. W przypadku dorzeczy transgranicznych, takich jak Niger, Senegal i Volta, podejścia lokalne są niewystarczające, a najodpowiedniejsze projekty łączą bazę satelitarną, rygorystyczną kalibrację, ukierunkowane informacje terenowe oraz systematyczne raportowanie niepewności11,39.
Operacyjny schemat decyzyjny dla wyboru metody
Powyższe porównanie staje się jawną procedurą wyboru na Rysunku 3, poprzez cztery sekwencyjne filtry, z których mogą korzystać zarówno agencje zarządzające zlewnią, jak i konsorcja badawcze. Pierwszym z nich jest pytanie. Dokumentowanie zmian wymaga jedynie analizy trajektorii teledetekcyjnych z formalną oceną dokładności. Przypisywanie zmian użytkowaniu i pokryciu terenu (LULC) zamiast klimatowi wymaga w minimalnym zakresie projektu kontrfaktycznego: uruchomienia modelu z ustalonym użytkowaniem terenu i zmiennym klimatem oraz odwrotnie. Testowanie przyszłych opcji wymaga łańcucha scenariuszy. Ilościowe określenie odpowiedzi wód podziemnych wymaga albo jawnej reprezentacji podpowierzchniowej, albo podejścia grawimetrycznego lub opartego na znacznikiach, ponieważ modele kalibrowane według przepływu w potokach nie ograniczają zasilania.
Drugim filtrem jest skala. Poniżej około 100 km2, a zwłaszcza w zlewniach miejskich, dominują procesy w skali pojedynczych zdarzeń, w związku z czym odpowiednie są zintegrowane modele powierzchniowo-podpowierzchniowe lub modele miejskie oparte na zdarzeniach, pod warunkiem, że drenaż i frakcja nieprzepuszczalna są zmapowane dokładniej niż zazwyczaj dostępne klasy pokrycia terenu. Pomiędzy około 100 a 10 000 km2 najlepszy stosunek nakładów danych do uzyskanych wyników oferują modele półrozproszone. Powyżej 10 000 km2 oraz w przypadku dorzeczy transgranicznych realistycznymi opcjami stają się podejścia regionalne lub zbiorcze, księgowanie wody w oparciu o dane satelitarne oraz uczenie maszynowe na dużych próbach45,46.
Trzecim filtrem jest dostępność danych, stosowana w sposób rzetelny, a nie życzeniowy. Przynajmniej trzy lata współistniejących zapisów odpływu i opadów umożliwiają konwencjonalną kalibrację; krótsze lub przerywane zapisy wymagają zastosowania satelitarnych funkcji celu dla ewapotranspiracji, wilgotności gleby i całkowitego zapasu wody21,34; w przypadku braku zapisów odpływu wybór sprowadza się do regionalizacji z zlewni donacyjnych lub zbiorczej predykcji z wykorzystaniem uczenia maszynowego, przy czym niepewność jest propagowana, a nie ukrywana. Czwartym filtrem jest horyzont decyzyjny i dostępne zasoby, gdzie każda ścieżka jest powiązana z orientacyjnym kosztem w Tabeli 5. Zastosowanie tych ram do analizowanego zbioru danych ujawniło dwa powracające niedopasowania: stosowanie modeli w pełni rozproszonych w sytuacjach, gdy dane wspierają jedynie model półrozproszony, oraz przypisywanie urbanizacji do klasy osadniczej zbyt ogólnej dla przeprowadzonej inferencji.
Wyzwania i możliwości
Pierwszym wyzwaniem jest jakość danych wejściowych. Przerwane szeregi hydrometeorologiczne, niska gęstość sieci stacji oraz rzadkie obserwacje zwierciadła wód gruntowych, wilgotności gleby, poborów i odpływów wymuszają zależność od produktów pochodnych i interpolacji, co satelitarne produkty opadów, imputacja oraz downscaling dostosowany do Afryki Zachodniej łagodzą tylko częściowo35,36,37,43,44. Drugim wyzwaniem jest atrybucja przyczynowa: trajektorie stresu wodnego jednocześnie odzwierciedlają zmienność klimatu, intensyfikację ekstremalnych opadów, zmiany LULC, urbanizację, pobory wody do irygacji, zamulanie oraz infrastrukturę, a jedynie niewielka część zebranych tutaj badań modelowych wykorzystuje projekt faktorialny, który jest niezbędny do ich rozdzielenia.
Trzecim wyzwaniem jest skala, a w tym przypadku bilans dowodów musi zostać sformułowany ostrożnie. Urbanizacja Afryki Zachodniej ma charakter rozproszony i nieciągły, łącząc formalne obszary zabudowy, nieformalne osiedla, nieutwardzone drogi i hybrydowy drenaż, i często przywołuje się ją jako czynnik stresu wodnego; jednak tylko 3 opracowania (9%) czynią z niej główny przedmiot badań, wszystkie w Niamey6,28,29, w przeciwieństwie do 15, w których dominującym sygnałem jest ekspansja rolnictwa. Zatem kategoryczne stwierdzenia dotyczące wpływu obszarów miejskich opierają się na znacznie słabszej podstawie niż stwierdzenia dotyczące ekspansji gruntów ornych: mechanizmy te są dobrze rozpoznane w ogólnej hydrologii miejskiej9,31, ale ich skala w zlewniach Afryki Zachodniej została określona ilościowo tylko w kilku przypadkach. Jest to zatem luka do wypełnienia, a nie wniosek, który można z całą pewnością sformułować, co wymaga bardziej szczegółowego mapowania obszarów miejskich, z rozróżnieniem powierzchni nieprzepuszczalnych, gęstości zabudowy, stanu dróg, obszarów depresyjnych oraz stref infiltracji8. Czwartym wyzwaniem jest walidacja wielozmienna: stres wodny wpływa na przepływ w strumieniach, wilgotność gleby, ewapotranspirację, zasilanie oraz głębokość zwierciadła wody gruntowej, zatem należy wykorzystać kilka zmiennych, aby ograniczyć zjawisko ekwifinalności50.
Wyróżniają się trzy możliwości: rosnący dostęp do bezpłatnych obrazów, danych topograficznych, produktów klimatycznych oraz operacyjnych produktów dotyczących ewapotranspiracji i wilgotności gleby, co poprawia powtarzalność34; fuzja danych łącząca obserwacje terenowe, produkty satelitarne, modele i uczenie maszynowe; oraz regionalizacja, która otwiera drogę do obserwatoriów hydro-terytorialnych39.
Zalecane kierunki działań obejmują wzmocnienie zlewni obserwacyjnych, z priorytetem dla krajów nieobecnych w tym zbiorze; uogólnienie analiz niepewności, porównywania wielu modeli oraz systematycznego raportowania metryk; integrację abstrakcji, infrastruktury i zastosowań; opracowanie modeli hydrologii miejskiej odzwierciedlających formalną i nieformalną infrastrukturę Afryki Zachodniej; jawne udostępnienie regionalnego zbioru danych z dużą próbką oraz wspólne opracowanie scenariuszy z menedżerami zasobów wodnych, planistami i rolnikami11.