The precepetation that occurred in the north of Iran on October 14, 2012 caused heavy rainfall in the three provinces of Mazandaran, Golestan, and East-Azarbaijan. The Arasbaran region in the East-Azarbaijan province, as one of the affected areas, has a Mediterranean climate, but the Caspian and Caucasian climates had a major impact on this event (Ramezani et al. 2019). Because of the high altitude of the Arasbaran region, the area has a variety of climates that naturally affect the vegetation of the area. Kaleybar city in the Arasbaran region was affected by the 2012 flood, causing serious damage to many buildings and cars. Fortunately, one of the Iranian permanent GPS stations was located in Kaleybar (named KLBR) and was used to help us in the tropospheric WV study. In this research, the tropospheric delay and the PWV have been evaluated using five existing stations in the Azerbaijan region during Day of Year (DOY) 286 to 291 in 2012. The Arasbaran flood occurred on DOY 289. The distribution of the five GPS stations is shown in Fig. 2, where KLBR and SKOH stations are located at a higher elevation than other stations.
Among the five GPS stations, AHAR, KLBR, and SKOH are equipped with meteorological sensors that can record air pressure, humidity, and temperature. Based on the meteorological data, KLBR station had the highest air pressure as well as the maximum humidity during the flooding days of DOY 287 to 291. The temperature range also decreases during the duration DOY 287 to 291 (Fig. 3).
GPS data of the five stations was processed using Bernese software (Dach et al. 2015) with a protocol of double-difference equations at a sampling rate of 30 seconds, and an elevation cut off angle of 5˚. The global mapping function (GMF) model was used to estimate ZHD (Rohm et al. 2014). The ZTD estimate for GPS stations varies from 1.80 to 2.20m (Fig. 4-a). The maximum amount of ZTD is observed at station AHAR and the minimum amount at station TABZ. The GPS-PWV estimates were evaluated in 5 phases for all five GPS stations with a focus on station KLBR (Fig. 4-b): during the first phase, a jump of 2 mm in a few hours before the flood, during the second phase, an increase of 3 mm, during the third phase, a jump of 3.5 mm, during the fourth phase (end of precipitation), a decrease of 3 mm, and during the fifth phase (the day after), a decrease of 4 mm was observed. The PWV jumps are marked with red squares (first, second, and third) and the PWV falls (fourth and fifth) with blue squares (Fig. 4-b). Stations that are not in flood areas have fewer variations than others.
ERA5 is a climate reanalysis dataset including hourly estimates of a large number of climatic, terrestrial, oceanic, and atmospheric variables from 2010 to 2017 with a horizontal resolution of 31 km grid spacing) Dong and Jin, 2018). In reduced spatial and temporal resolutions, the ERA5 has information about uncertainty for all variables in which temporal resolution reaches one hour (Jiang et al. 2020).
Based on the ERA5 model, the PWV value can be calculated based on specific humidity and air pressure at a horizontal resolution of 0.25°⋅0.25° and a vertical resolution of 37 pressure levels (Zhang et al. 2019). By integrating GPS-derived PWV and ERA5-PWV (Wang and Lio, 2019), tomography and regularization were performed using Tikhonov to achieve 4 Dimension (4D) PWV. The PWV increases 7–8 mm at altitude 3 km on the day before the flood, reaches 5-8mm at altitude 3–4 km on the flood day, and decreases 4–5 mm at altitude 3–4 km on the day after the flood. Also, the PWV reaches 10 mm at an altitude of 1–2 km, and is less than 4 mm at an altitude of more than 5 km on the day before the flood (Fig. 5).
On the day of the flood, an increase of PWV is observed from altitude 2 km to 4 km over KLBR and AHAR stations, and on the day after the flood, a significant decrease of PWV is estimated from altitude 2 km to 4 km. PWV changes in higher than altitude 3 km can be observed in which PWV increase is abnormal. By simultaneously examining the PWV time series (Fig. 4) and PWV tomography (Fig. 5), we can see an increase of PWV in intervals of less than 3 hours in phase 2 before the flood and a reduction of PWV in intervals of less than 3 hours in phase 4 and 5 after the flood. The increase in PWV at intervals of less than 3 hours can cause flooding in the region. It can be a feature of climate, topography, or other influential factors that require future flood studies to estimate the increase in PWV. Examining several floods and the rate of increase of PWV in the study area can be developed as an indicator for a flood forecast.
Fig. PWV tomography on a) DOY 288, b) DOY 289 and c) DOY 290
Tomographic validation is essential to verify the correctness of the proposed method and confirm the results. For this purpose, in this research, two sensors: Radiosonde data at TABRIZ station and Moderate Resolution Imaging Spectroradiometer (MODIS) satellite have been used. Radiosonde data at TABRIZ station (TABZ) have been applied for the validation of the tomography results (Fig. 6). The obtained results from the tomography and the radiosonde on the day before the flood (DOY 288), the day of the flood (DOY 289), and the day after the flood (DOY 290) are in good agreement (Fig. 6). The correlation between radiosonde and tomography with a difference of -1.5 to 1.5 mm and a significant correlation of 0.95 is also shown in Fig. 6-d. The difference on the flood day was more than 1mm at an altitude of ~ 1.5 km, but it was lower at other altitudes (Fig. 6-c).
MODIS is capable of recording phenomena information in 36 bands ranging from 4.4 to 4.14 microns with variable resolution (2 bands 250 m, five bands 500 m, and 29 bands 1 km). In this study, MOD05-L2 data with a resolution of 1˚×1˚ was used for the final validation of PWV for an 8-day duration (from 8 October 2012 to 15 October 2012). Details of the MODIS-Near-Infrared (NIR) Water Vapor algorithm are discussed in Gui et al. (2017). Figure 7 shows the PWV values obtained from the GPS and MODIS-NIR in the study area over the 8 days. The differences between GPS-PWV and MODIS-PWV for the 5 GPS stations are given in Table 1.
Table 1
Differences of GPS-PWV and MODIS-PWV for the 5 GPS stations
GPS station | GPS-PWV (mm) | MODIS-PWV (mm) | Difference (mm) |
TABZ | 14.62 | 13.87 | 0.75 |
SKOH | 13.81 | 13.10 | 0.71 |
KHJE | 22.56 | 20.79 | 1.77 |
AHAR | 20.37 | 22.77 | -2.4 |
KLBR | 23.55 | 28.72 | -5.17 |
According to the obtained results, stations KLBR, AHAR, TABZ, and SKOH have a difference of less than 3 mm, and KHJE station has a difference of about 5 mm. The absolute error between GPS-PWV and MODIS-NIR-PWV is not more than 4.4 mm. In addition, MODIS-NIR-PWV data is known to have larger uncertainty than typical GPS data even for strictly co-located observations. This is probably another factor contributing to the large difference (up to 5 mm) between the GPS and MODIS data. GPS spatial resolution is low due to various constraints, while MODIS offers a much better spatial resolution. In contrast, MODIS like radiosonde has a lower temporal resolution than GPS with a sampling rate of 30 seconds. This is a prominent advantage of GPS in meteorological studies.
High ZTD/PWV does not create floods. It is high precipitation together with other factors that can generate floods. High PWV does not necessarily mean there will be precipitation. On the other hand, there is no precipitation without high PWV. Unlike other studies, in this study, an increase in PWV was observed at several stations before the flood (Supartaet al. 2012). Many studies have tried to link high PWV to precipitation, indicating many times a high correlation between PWV and precipitation, which have been used in attempts of predicting flood generating precipitation events. The most convincing ones stem from regions in the tropics where events are extremely severe. There the analysis of extremely high PWV, or the change in it, could be linked to the onset of heavy rain events (Huelsinget al. 2017; Yao et al. 2017; Chen, 2017; Dong and Jin, 2018; Biet al. 2006).
What will increase the potential of GPS for meteorology is the determination of the real time PWV. This can be used in artificial intelligence algorithms to estimate flash floods and heavy rainfall (Carreau al. 2021; Dtissibe al. 2020). These high-accuracy algorithms can help detect flood signals with PWV GPS.