Ներողություն ենք խնդրում կայքի մեքենայական թարգմանության հնարավոր սխալների համար
Ջրի հոսքի սեզոնային կարգավորման ES
(InVEST սեզոնային ջրի բերքատվության մոդել)
GIS մոդելավորում և վերլուծություն, արդյունքների ներկայացում – Էդուարդ Կազակով (NextGIS OÜ, Էստոնիա);
Արդյունքների վերլուծություն և ներկայացում – Ելենա Բուկվարևա (BCC-Հայաստան);
Հայաստանի իրական ջրբաժանների հետ արդյունքների համապատասխանության գնահատում – Վարդան Ասատրյան (ՀՀ ԳԱԱ)
Մոդելավորման արդյունքների համապատասխանության գնահատում Հայաստանի համար. արտաքին փորձագետ Ալեքսանդր Առաքելյան (ՀՀ ԳԱԱ)
The ES of seasonal flow redistribution and baseflow provision is extremely important for Armenia, which has a seasonal climate with dry summers over a significant part of the territory. The ES was estimated and mapped using InVEST Seasonal Water Yield model. The model takes into account the monthly amount of precipitation, soil permeability, and the water balance of each pixel, including moisture that comes into it from the overlying pixels. The main resulting indicators for assessing the ecosystem service are quick flow (QF), that is, the generation of streamflow with watershed residence times of hours to days; and baseflow (B), that is the generation of streamflow with watershed residence times of months to years. Values of cumulative baseflow (Bsum) show the flow through a pixel, contributed by all upslope pixels (Figure 31A2-1). Baseflow ensures runoff maintaining during the dry season and possible droughts.
Figure 31A2-1. The InVEST Seasonal Water Yield model diagram (from Hamel et al., 2020)
Մոդելավորման համար օգտագործված գործակիցները կարող եք գտնել այստեղ
InVEST մոդելի արդյունքները միջնորդ փոփոխականներ են, որոնք պետք է մեկնաբանվեն հարաբերական առումով, այլ ոչ թե ֆիզիկական մեծությունների: Այնուամենայնիվ, բացահայտված արժեքներն օգտակար են երկրի տարածքում ծառայությունների տարածական բաշխվածության և դրանց հավասարակշռությունը բնակչության և տնտեսության կողմից ծառայությունների օգտագործման ցուցանիշների հետ վերլուծելու համար:.
Այս բաժինը ներկայացնում է էկոհամակարգային ծառայությունների գնահատման և քարտեզագրման մոդելների փորձարկման նախնական արդյունքները: Հետագայում, եթե որոշում կայացվի օգտագործել այս մոդելները, ապա դրանք պետք է չափորոշվեն Հայաստանում կատարված հիդրոլոգիական չափումների միջոցով։
Մեթոդաբանությունը և արդյունքները մանրամասն կներկայացվեն առաջիկա հրապարակման մեջ.
ES provided by terrestrial ecosystems
The baseflow values (B) for the actual land cover (ESRI 2023) are much higher, while the quick flow (QF) values, on the contrary, are lower compared to hypothetical cases where all natural ecosystems are replaced with bare ground or croplands (Table 31A2-1). The difference between runoff values for the actual land cover and bare ground can be interpreted as the ES provided by terrestrial vegetation: B = 47.8 mm in average, and QF = -22.2 mm in average. Thus, ecosystems provide 93% of baseflow and reduce quickflow by 22%.
ES maps show that under the bare ground and cropland scenarios, baseflow is almost absent (Figure 31A2-2), meaning that the existing baseflow is almost entirely provided by terrestrial ecosystems..
Table 31A2-1. Mean ES indicator values for Armenia under different scenarios
| Սցենար | Բազային հոսք, մմ (B) | Արագ հոսք, մմ (QF) | (B+QF) | B-ի մասնաբաժինը ընդհանուր հոսքում, % |
| Հողածածկ 2017թ | 51.97 | 97.01 | 148.97 | 34.88 |
| Հողածածկ 2023 թ | 51.28 | 98.04 | 149.32 | 34.34 |
| Բաց հողի սցենար | 3.43 | 120.22 | 123.65 | 2.78 |
| Բուսական հողերի սցենար | 3.58 | 124.96 | 128.54 | 2.78 |
Figure 31A2-2. Maps of ES indicators for different scenarios. For detailed maps see section “Seasonal Water Yield” in the project WebGIS
With the current landcover, the baseflow is on average 35% of the total flow (from 28 to 40% in different watersheds). With the bare ground scenario, the baseflow is only 3% (from 2 to 4%) (Table 31A2-2; Figure 31A2-3). The values of these indicators across the marzes largely mirrors them for the corresponding watersheds (Figure 31A2-4).
Figure 31A2-3. Baseflow and quick flow under the current land cover and the bare ground scenario across watersheds
Figure 31A2-4. Baseflow and quick flow under the current land cover across marzes
We estimated ES volume provided by terrestrial vegetation as difference between ES indicator values for the current land cover in 2023 and the bare ground scenario where all grasslands and trees were replaced by bare ground. Across watersheds, ecosystems provide 92%–95% of baseflow (Table 31A2-2; Figure 31A2-5) and reduce quick runoff by 13%–36% (Table 31A2-2; Figure 31A2-6).
Figure 31A2-5. Ecosystem effects on baseflow: (a) baseflow under current land cover and the component provided by ecosystems; (b) share of baseflow provided by ecosystems (%).
Figure 31A2-6. Ecosystem effects on quick flow: (a) quick flow under current land cover and its part reduced by ecosystems; (b) share of quick flow reduced by ecosystems (%).
Table 31A2-2. ES indicators under the current land cover and the bare ground scenario across watersheds
| Ցուցանիշներ | Voro- tan | Մեծամոր | Հրազդան | Դեբեդ | Արփա | Ախուրյանը | Աղստեւ | Հայաստան | ||
| Հողատարածք ESRI 2023 | Բազային հոսք, մմ | Բ2023 | 31.0 | 50.4 | 53.4 | 78.3 | 37.3 | 73.2 | 42.4 | 51.3 |
| Արագ հոսք, մմ | QF2023 | 79.2 | 77.5 | 132.0 | 115.6 | 59.3 | 120.7 | 87.6 | 98.0 | |
| B-ի մասնաբաժինը ընդհանուր հոսքում, % | 28.2 | 39.4 | 28.8 | 40.4 | 38.6 | 37.7 | 32.6 | 34.3 | ||
| Բաց հողի սցենար | Բազային հոսք, մմ | Բbg | 2.3 | 2.8 | 3.1 | 6.5 | 1.9 | 4.9 | 3.6 | 3.4 |
| Արագ հոսք, մմ | QFbg, | 104.6 | 90.7 | 149.9 | 148.2 | 80.6 | 138.1 | 117.3 | 120.2 | |
| B-ի մասնաբաժինը ընդհանուր հոսքում, % | 2.1 | 3.0 | 2.0 | 4.2 | 2.3 | 3.5 | 3.0 | 2.8 | ||
| Ecosystem effects | Baseflow provided by ecosystems, mm | Բeco= Բ2023-Բbg | 28.7 | 47.6 | 50.3 | 71.8 | 35.4 | 68.2 | 38.8 | 47.8 |
| Share of baseflow provided by ecosystems, % | Բeco*100/B2023 | 92.6 | 94.4 | 94.2 | 91.7 | 94.9 | 93.2 | 91.5 | 93.2 | |
| Reduction of quickflow by ecosystems, mm | QFeco = QF2023-QFbg, | -25.5 | -13.2 | -17.8 | -32.6 | -21.3 | -17.4 | -29.7 | -22.2 | |
| Share of quick flow reduced by ecosystems, % | QFeco* 100/ QF2023, | -32.2 | -17.0 | -13.5 | -28.2 | -35.9 | -14.4 | -33.9 | -22.7 |
Աղյուսակ 3. Ստորին գետի կուտակային բազային հոսքը (Bsum) հասնում է մինչև 2023թ. փաստացի հողածածկույթը և այն սցենարը, երբ բոլոր բնական էկոհամակարգերը փոխարինվում են մերկ հողով:
| Հողածածկ, 2023 թ | Բաց հողի սցենար | Տարբերությունը Բգումար 2023-Բմերկ գումարը | բաժնեմասը Բգումար ապահովված է ցամաքային էկոհամակարգերով, % | |
| Աղստեւ | 686063872 | 53702228 | 632361644 | 92.2 |
| Ախուրյանը | 1562785152 | 82415360 | 1480369792 | 94.7 |
| Արգիչի | 75929520 | 3704933 | 72224587 | 95.1 |
| Արփա | 313119360 | 11837469 | 301281891 | 96.2 |
| Ազատ | 349204672 | 11007498 | 338197174 | 96.8 |
| Դեբեդ | 2672169472 | 179894192 | 2492275280 | 93.3 |
| Գավառագետ | 422896480 | 13807621 | 409088859 | 96.7 |
| Հրազդան | 3543649536 | 154972240 | 3388677296 | 95.6 |
| Մասրիկ | 43279976 | 1879218.75 | 41400757.25 | 95.7 |
| Մեծամոր | 1676233984 | 72298784 | 1603935200 | 95.7 |
| Վեդի | 134872176 | 6922562 | 127949614 | 94.9 |
| Ողջի | 75920344 | 4065560.75 | 71854783.25 | 94.6 |
| Որոտան | 300472096 | 15817257 | 284654839 | 94.7 |
Table. Average baseflow and quick flow in provinces with land cover 2023
| Բազային հոսք, մմ | Արագ հոսք, մմ | |
| Արագածոտն | 61 | 91 |
| Արարատ | 28 | 50 |
| Արմավիր | 22 | 47 |
| Գեղարքունիք | 55 | 151 |
| Կոտայք | 67 | 90 |
| Լոռի | 74 | 116 |
| Շիրակ | 85 | 126 |
| Սյունիք | 31 | 79 |
| Տավուշ | 37 | 86 |
| Վայոց ձոր | 39 | 66 |
.
Figure. The share of baseflow (%) provided by terrestrial ecosystems in watersheds and in ստորին գետը հասնում է
Analysis of ES indicators across vegetation zones shows that the highest average baseflow values occur in alpine and subalpine zones, while woody vegetation exhibits levels similar to those of various grassland types. The range of values within individual polygons of certain vegetation zones is quite large — in alpine, subalpine, meadow-steppe, forest zones, and juniper woodlands (Figure 31A2-7a). Proportion of current baseflow contributed by ecosystems is lowest in forest, steppe, and marsh zones (88–89%), and highest (98%) in the desert zone due to the specific characteristics of the only small area where it is currently found in Armenia. In other vegetation zones, this indicator ranges from 93% to 96% (Figure 31A2-7b). 
Figure 31A2-7. Baseflow indicators across vegetation zones. Red dots indicate average values, while colored bars represent the range of values within individual polygon
The highest baseflow values (155 and 137 mm) are calculated for alpine and subalpine grasslands, while the forest zone showed a minimal baseflow of 29 mm, similar to that of the steppes (25 mm); both are lower than those in the semidesert and desert zones. The proportion of baseflow contributed by ecosystems is also minimal in the forest and steppe zones (89%). This counterintuitive result, in our view, is explained by the combined effects of multiple factors that determine baseflow—precipitation, terrain slope, and soil permeability. Very high absolute baseflow values in mountain grasslands result from the high precipitation in the mountains. The low baseflow values in the forest zone are most likely the result of forests occurring predominantly on the steep slopes of gorges and mountains. According to our assessment, the highest mean slope among the vegetation zones occurs in the forest and juniper zones —about 20°, whereas mountain grasslands and steppes occupy gentler slopes from 10° to 17°, and the semideserts and the single desert patch lie on plains with an average slope of about 6° (Figure 31A2-8). The moderate baseflow of 38 mm and the high proportion of it contributed by ecosystems (94%) in the semidesert zone are most likely due to its location in areas with the gentlest relief and a high proportion of highly permeable soils. The only small desert patch remaining in Armenia exhibits moderate baseflow of 34 mm an extremely high proportion of baseflow provided by ecosystems (98%), probably because it is entirely located on soils with the highest permeability.
Figure 31A2-8. Average slope for vegetation zones (the zones are arranged by elevation for illustration only; in reality, the elevation distribution may differ).
2. ES-ը փոխվում է 2017-ից մինչև 2023 թ
All changes identified are determined only by changes in the landcover (Section 2.2.B). Weather and climate changes are not taken into account. The maps show that the changes are sporadic and oppositely directed (Figure 31A2-9). Nevertheless, overall, they can be characterized as negative. In cases where ES indicators changed significantly, baseflow (B) decreased while quick flow (QF) increased. This means that the ability of terrestrial ecosystems to sustain baseflow during dry periods is declining. The only exception is the Arpa basin, where B has increased.The most significant negative changes occurred in the Shirak province and the Akhuryan watershed corresponding to land cover changes in the Shirak province. The reason is the expansion of the croplands at the expense of the grasslands in the Shirak province (տես այստեղ).
Figure 31A2-9. Changes in baseflow (a) and quickflow (b) from 2017 to 2023, mm. For detailed maps see sections “Seasonal Water Yield – Dynamics” in the project Web GIS
The most notable changes across watersheds and marzes are a decrease in baseflow and an increase in quickflow. These negative changes are most pronounced in the Akhuryan and Metsamor watersheds and in the marzes Shirak and Aragatsotn. A slight increase in baseflow is observed in the Arpa and Aghstev watersheds and in Ararat marz (Table 31A2-3; Figures 31A2-10 and 31A2-11).
Figure 31A2-10. Changes in baseflow and quick flow from 2017 to 2023 in absolute (a) and relative (b) indicators across watersheds
Figure 31A2-11. Changes in baseflow and quickflow from 2017 to 2023 in absolute (a) and relative (b) indicators across marzes
Table 31A2-3. Changes in baseflow and quick flow from 2017 to 2023 across watersheds and marzes
| Բազային հոսք, միջին | Արագ հոսք, նշանակում է | ||||||||
| B 2017, մմ | B 2023, մմ | B-ի փոփոխություն, մմ | B, %-ի փոփոխություն 2017թ | QF 2017, մմ | QF 2023, մմ | QF-ի փոփոխություն, մմ | QF-ի փոփոխություն, % 2017-ի համեմատ | ||
| Ջրբաժաններ | Աղստեւ | 42.1 | 42.4 | 0.2 | 0.6 | 87.3 | 87.6 | 0.3 | 0.3 |
| Ախուրյանը | 79.0 | 73.2 | -5.9 | -7.4 | 116.6 | 120.7 | 4.1 | 3.5 | |
| Արփա | 37.0 | 37.3 | 0.2 | 0.7 | 59.3 | 59.3 | 0.0 | 0.0 | |
| Դեբեդ | 78.6 | 78.3 | -0.3 | -0.4 | 114.3 | 115.6 | 1.3 | 1.1 | |
| Հրազդան | 53.5 | 53.4 | 0.0 | -0.1 | 131.4 | 132.0 | 0.6 | 0.5 | |
| Մեծամոր | 51.7 | 50.4 | -1.2 | -2.4 | 75.8 | 77.5 | 1.7 | 2.3 | |
| Որոտան | 31.0 | 31.0 | 0.0 | 0.1 | 78.6 | 79.2 | 0.5 | 0.7 | |
| գավառներ | Արագածոտն | 62.5 | 61.0 | -1.5 | -2.5 | 88.5 | 90.5 | 2.0 | 2.2 |
| Արարատ | 27.6 | 28.1 | 0.5 | 1.7 | 49.5 | 49.6 | 0.1 | 0.3 | |
| Արմավիր | 22.2 | 21.9 | -0.3 | -1.3 | 46.2 | 46.8 | 0.6 | 1.2 | |
| Գեղարքունիք | 54.6 | 54.5 | -0.1 | -0.1 | 150.5 | 151.1 | 0.6 | 0.4 | |
| Կոտայք | 66.4 | 66.7 | 0.3 | 0.4 | 89.8 | 90.3 | 0.5 | 0.6 | |
| Լոռի | 74.1 | 74.0 | -0.2 | -0.2 | 114.5 | 115.8 | 1.3 | 1.2 | |
| Շիրակ | 91.3 | 85.1 | -6.2 | -6.8 | 121.3 | 125.7 | 4.4 | 3.6 | |
| Սյունիք | 30.9 | 30.9 | 0.0 | 0.1 | 78.6 | 79.1 | 0.5 | 0.7 | |
| Տավուշ | 36.7 | 36.9 | 0.1 | 0.3 | 86.2 | 86.3 | 0.1 | 0.1 | |
| Վայոց ձոր | 39.3 | 39.4 | 0.1 | 0.3 | 66.0 | 66.0 | 0.0 | 0.0 | |
| Հայաստան | 52.0 | 51.3 | -0.7 | -1.3 | 97.0 | 98.0 | 1.0 | 1.1 | |
Ստորին հոսանքի կուտակային բազային հոսքը (Bsum) ամենաշատը նվազել է Ախուրյան գետում (-9.25%), և Մեծամոր գետում (-3.03%), մինչդեռ Ողջի գետում այն տեսանելիորեն աճել է (+5.28%)
Ստորին գետում կուտակային բազային հոսքի փոփոխությունները հասնում են 2017-ից 2023 թվականներին
| Գետեր | Բսում 2023, մմ | Bsum 2017, մմ | Փոփոխություններ, մմ | Փոփոխություններ, % համեմատ 2017 թ |
| Աղստեւ | 686,063,872 | 685,638,016 | 425,856 | 0.06 |
| Ախուրյանը | 1,562,785,152 | 1,722,085,760 | -159,300,608 | -9.25 |
| Արգիչի | 75,929,520 | 75,653,040 | 276,480 | 0.37 |
| Արփա | 313,119,360 | 309,875,328 | 3,244,032 | 1.05 |
| Ազատ | 349,204,672 | 345,485,152 | 3,719,520 | 1.08 |
| Դեբեդ | 2,672,169,472 | 2,684,654,848 | -12,485,376 | -0.47 |
| Գավառագետ | 422,896,480 | 425,799,968 | -2,903,488 | -0.68 |
| Հրազդան | 3,543,649,536 | 3,543,354,368 | 295,168 | 0.01 |
| Մասրիկ | 43,279,976 | 42,957,568 | 322,408 | 0.75 |
| Մեծամոր | 1,676,233,984 | 1,728,691,072 | -52,457,088 | -3.03 |
| Վեդի | 134,872,176 | 132,632,480 | 2,239,696 | 1.69 |
| Ողջի | 75,920,344 | 72,110,264 | 3,810,080 | 5.28 |
| Որոտան | 300,472,096 | 300,990,784 | -518,688 | -0.17 |
3. Supply-use balance
Since the baseflow and quick flow values computed by the SWY model are, in accordance with the InVEST User Guide, relative measures, the baseflow volume VԲ was calculated as follows: VԲ=BFI⋅Vtot where Vtot is the river discharge according to the data reported by the Government of Armenia and ARMSTAT (Table below); BFI is baseflow index, BFI=B/(B+QF), where B and QF are baseflow and quick flow values computed by SWY model.
This analysis is approximate, because we used river flow data by watersheds and water-use data by provinces. The boundaries of provinces and watersheds largely coincide, which allows such a rough estimate. For watersheds that include two provinces, the data for those provinces were summed. Obviously, the balance should be refined in the future using data from the same EAAs.
| Watershed | River flow, millions of m3 | Data source | Details |
| Աղստեւ | 260 | ArmStat Regional Statistics handbooks 2023 https://armstat.am/en/?nid=651 | – |
| Ախուրյանը | 700.7 | Basin management plan https://www.arlis.am/hy/acts/112336 | River flow was calculated as the sum of utilizable river flow (506.2) and environmental flow (194.5). |
| Արփա | 1177.0 | Basin management plan https://www.arlis.am/hy/acts/112336 | The sum of river flows of Azat, Vedi, and Arpa |
| Դեբեդ | 960 | ArmStat Regional Statistics handbooks 2023 https://armstat.am/en/?nid=651 | – |
| Հրազդան | 712.2 | Basin management plan https://www.arlis.am/hy/acts/171449 | – |
| Մեծամոր | 777.8+265.5=1043.3 | Basin management plans https://www.arlis.am/hy/acts/171449 https://www.arlis.am/hy/acts/112336 | Total river flow was calculated as the sum of Qasagh (265.5) and Metsamor (777.8) rivers. Metsamor river flow wascalculated as the sum of utilizable river flow (1,786.7) and environmental flow (106.0) excluding water available from the reach downstream of the confluence of the Araks and Akhuryan rivers (1,114.9). |
| Որոտան | 1319.6 | Basin management plan https://www.arlis.am/hy/acts/106124 | The sum of river flows of Southern BMA |
The overwhelming majority of water consumption is accounted for by the agriculture, fisheries, and forestry sector, which underscores the importance of assessing the ecosystem service of baseflow provision. In two watersheds —Metsamor and Arpa—agricultural water consumption exceeds the baseflow volume (Figure 31A2-12a). In the Metsamor and Arpa watersheds, baseflow provides 56% and 55% of agricultural water consumption, respectively; in the Hrazdan watershed baseflow slightly exceeds agricultural consumption—by 16%. In the other watersheds baseflow exceeds water consumption by many times (Figure 31A2-12c). The water-use data for Tavush marz used in the analysis pertain to an area larger than the Aghstev River watershed for which river flow data are available. Therefore, in reality the total river flow from the entire area exceeds water consumption even more than indicated by these results.
Figure 31A2-12. Supply–consumption balance of ES: (ա) river flow, baseflow and water consumption across watersheds (* data reported by ARMSTAT and the Government); (բ) Difference between baseflow volume and water consumption in agriculture, fishery and forestry; (c) Share of baseflow in agriculture, fisheries and forestry water consumption (%).
Table 31A2-5. Supply–use balance of ES
| Water use in agriculture, fish breeding, forestry*, millions of m3 | Total water consum-ption*, millions of m3 | Baseflow, millions of m3 | River flow*, millions of m3 | Difference between baseflow volume and water consumption in agriculture, fishery and forestry, millions of m3 | Share of baseflow in agriculture, fisheries and forestry water consumption (%) | |
| Aghstev (a part of the marz Tavush) | 11.10 | 20.90 | 84.78 | 260.00 | 73.68 | 763.78 |
| Akhuryan (marz Shirak) | 53.00 | 63.10 | 264.46 | 700.70 | 211.46 | 498.99 |
| Arpa (marzes Ararat and Vayots Dzor) | 820.30 | 835.90 | 454.50 | 1177.00 | -365.80 | 55.41 |
| Debed (marz Lori) | 3.00 | 18.00 | 387.69 | 960.00 | 384.69 | 12923.03 |
| Hrazdan (marzes Kotayk and Gegharkunik) | 176.10 | 218.60 | 205.13 | 712.20 | 29.03 | 116.49 |
| Metsamor (marzes Armavir and Aragatsotn) | 740.40 | 769.10 | 411.06 | 1043.30 | -329.34 | 55.52 |
| Vorotan (marz Syunik) | 6.00 | 89.50 | 371.51 | 1319.60 | 365.51 | 6191.75 |
*Data reported by ARMSTAT and the Government of Armenia
4. Կրթության փոքր կողմնակի արտադրանք
Սեզոնային ջրաբերություն մոդուլը պահանջում է տվյալներ Հայաստանի տարածքում ամսական տեղումների բաշխման վերաբերյալ: Այնուամենայնիվ, մոդուլը հաշվի չի առնում ձյան ժամանակահատվածը, երբ ձյունը չի հալվում, և, հետևաբար, հեղուկ տեղումներ չկան: Այս խնդրի լուծման համար մենք հաշվարկել ենք Հայաստանում հեղուկ տեղումների քանակը ամիսներով։ Ի լրումն մոդելում օգտագործելու, այս արդյունքները կարող են օգտակար լինել Հայաստանի կլիմայական բնութագրերը ավելի լայն լսարանի բացատրելու համար:
Հղումներ
Hamel, P., Valencia, J., Schmitt, R., Shrestha, M., Piman, T., Sharp, RP, Francesconi, W., & Guswa, AJ (2020): Լանդշաֆտի կառավարման համար սեզոնային ջրի բերքատվության մոդելավորում. կիրառություններ Պերուում և Մյանմարում: Բնապահպանական կառավարման ամսագիր, 270, 110792. https://doi.org/10.1016/j.jenvman.2020.110792
