Next Article in Journal
Comparison of Supervised Learning and Changepoint Detection for Insect Detection in Lidar Data
Previous Article in Journal
SiameseNet Based Fine-Grained Semantic Change Detection for High Resolution Remote Sensing Images
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Remote Sensing Analysis for Vegetation Assessment of a Large-Scale Constructed Wetland Treating Produced Water Polluted with Oil Hydrocarbons

by
Khaled Al-Jabri
1,
Yaseen Al-Mulla
1,2,*,
Farid Melgani
3 and
Alexandros Stefanakis
4,5
1
Department of Soils, Water and Agricultural Engineering, Sultan Qaboos University, Al-Khoud 123, Muscat P.O. Box 34, Oman
2
Remote Sensing and GIS Research Center, Sultan Qaboos University, Al-Khoud 123, Muscat P.O. Box 33, Oman
3
Department of Information Engineering and Computer Science, University of Trento, Via Sommarive, 9, 38123 Trento, Italy
4
Laboratory of Environmental Engineering and Management, School of Chemical and Environmental Engineering, Technical University of Crete, 73100 Chania, Greece
5
The European University on Responsible Consumption and Production—EURECA-PRO, 73100 Chania, Greece
*
Author to whom correspondence should be addressed.
Remote Sens. 2023, 15(24), 5632; https://doi.org/10.3390/rs15245632
Submission received: 2 November 2023 / Revised: 25 November 2023 / Accepted: 28 November 2023 / Published: 5 December 2023
(This article belongs to the Section Ecological Remote Sensing)

Abstract

:
The identification and assessment of plant stress using wetland satellite images is a major task in remote sensing. In this study, one of the largest constructed wetlands (CWs) in the world, located in the Sultanate of Oman, was examined, assessed, and evaluated using remote sensor data from Sentinel-2. This CW system treats produced water generated during oil exploration activities in a desert environment; thus, CW vegetation is subjected to stress induced by oil hydrocarbons and water salinity. This study examined the plant stress and detected changes between the years of 2017 and 2019. Sentinel satellite images were evaluated for vegetation status extraction. The Normalized Difference Vegetation Index (NDVI), Modified Soil-Adjusted Vegetation Index (MSAVI), and Normalized Difference Salinity Index (NDSI) were used to evaluate the vegetation change. The results showed a comprehensive mapping identification of the plant stress and water flow parameter factors including oil in water contamination (OIW), dissolved oxygen (DO), water temperature (WT), and water conductivity (COND). Among the three indices, it was found that the NDVI showed a very good correlation with all parameters in both years with average R2 = 0.78, 0.67, 0.75, and 0.60 for OIW, DO, WT, and COND, respectively. The same trend was found for MSAVI but with R2 = 0.59, 0.48, 0.55, and 0.56 for OIW, DO, WT, and COND, respectively. This shows that the NDVI performed better than the MSAVI in evaluating the water flow parameters. On the other hand, the NDSI showed a strong correlation with one flow parameter, that is, water conductivity, especially at the outlet cells of the CW with R2 = 0.86 and 0.82 for winter time and summer time, respectively. The synchronization and correlation between the water flow parameters and remote sensing vegetation indices in this study lead to a new approach to large-scale landscape wetland monitoring that improves and helps predict any degradation or stress on vegetation growth. Furthermore, the results of this work can help decision makers potentially modify the wetland design and water flow path to improve future expansion phases. The mapping of such a critical and massive industrial CW should consider the use of high spatial resolution sensors where identifications and classifications are further improved. In summary, this research demonstrates that it is feasible to estimate vegetation stress within the constructed wetland using remote sensing techniques across extensive regions when an ample dataset comprising field data, satellite imagery, and supporting information is accessible.

1. Introduction

Wetland ecosystems in nature are important for local communities as they provide fishery products, timber goods, water purification, recreational uses, and many other ecosystem services [1,2]. Coastal wetland ecosystems protect local communities from flooding and sea level rise [3]. Natural processes and components are currently exploited in human-made wetland ecosystems that replicate the multiple functions of natural wetlands under a controlled environment. These systems, known as Constructed Wetlands (CWs), are nature-based solutions that have been used for many years for flood protection, water storage, habitat creation, and water quality improvement [4,5]. The research and development advances in this field of ecological engineering enabled the widespread use of CWs across the world for the treatment of wastewater from many sources such as municipal [6,7], industrial [8,9,10], and agro-industrial sources [11,12,13]. Wetlands provide various benefits like cleaning water, preventing floods, shielding coastlines, saving soil and water, filtering mud, removing pollution, and offering beauty and recreation [14]. A specific and unique application of CW is for water containing fuel and oil hydrocarbons [15], especially for produced water treatment, i.e., water polluted with oil hydrocarbons from oil exploration activities, due to their low cost and sustainable character compared to the conventional technologies based on physical and chemical processes [16]. Natural wetlands are shallow water basins that are fully covered with wetland plants. They are often called earth’s kidneys due to their advantages in improving water quality, protecting shorelines, and providing habitats for many plants, birds, and animals.
Since the 1900s, around 50% of the earth’s wetlands have been lost [17]. Hence, different studies on the effects of climate change, seasonal changes, plant species, and environmental challenges explore these ecosystems’ existence and the current conditions as well as their geographical distributions to promote their preservation and restoration efforts. One of the recent tools used for wetland monitoring and restoration is remote sensing. This technique can be used to monitor the status of a CW, and it can be used as a useful tool that can provide valuable information, e.g., about plant health, in a relatively short time [1,2,16,18].
Remote sensing techniques have been widely used to classify vegetation monitoring as well as to evaluate the targeted plants, soil, and water changes [14]. These techniques have assisted in dealing with some impacts of natural and anthropogenic causes by using spectral data for decision making such as the vegetation identification of different regions using different satellites and methods of vegetation extraction [18,19]. Various vegetation indices, which are pivotal for evaluating vegetation and crop well-being, have emerged as standard tools in vegetation health assessments. Indices such as the Normalized Difference Vegetation Index (NDVI), Soil-Adjusted Vegetation Index (SAVI), modified SAVI (MSAVI), and Normalized Difference Salinity Index (NDSI) have been instrumental, utilized by the National Aeronautics and Space Administration (NASA) Earth Observatory for global vegetation health monitoring.
Abood et al. [20] investigated the use of remote sensing technology in analyzing and monitoring natural resources such as soil salinization in the Mesopotamian valley in Iraq. They showed the effectiveness of using different indices, including the NDSI, NDVI, and SAVI, to analyze and evaluate the extent of saline soils in this valley of Iraq. Guo et al. [17] provided a comprehensive review that showed how effective the use of different remote sensing techniques is in classifying, evaluating, and monitoring wetlands. Lu et al. [21] reviewed the applications of multispectral and hyperspectral remote sensing data in extracting biophysical and biochemical properties of wetland vegetation such as the water content, vegetation, and leaf area index. The challenges of wetland remote sensing applications in selecting the proper spatial and spectral resolutions were highlighted as well as in picking the appropriate method for the spectral information extraction of wetland vegetation health monitoring [22]. Dronova, I. [23] explored the use of the object-based image analysis (OBIA) method in characterizing wetlands and concluded that this method will play an important role in wetland remote sensing, especially with the advancement of the dataset’s very high resolutions. Valderrama-Landeros et al. [2] assessed mangrove species growing in wetlands in Mexico and found that WorldView-2 satellite images provided the highest classification accuracy for these species in comparison to the Landsat-8, SPOT-5, and Sentinel-2 satellite images. Singh et al. [24] studied the changes in land use and land cover of disturbed areas on forest wetlands from 1990 to 2014 in Goalpara and Bongaigaon, India using remote sensing and geographic information system (GIS) technologies. Their results showed that there was a degradation of forest cover and a loss of wetland area over the 24 years of study. Pan et al. [25] studied the climate change impact on the wetland of Dunhuang Yangguan National Nature Reserve in China using the NDVI on Landsat images for the years between 1988 and 2016. They related the annual precipitation reduction to wetland vegetation and soil surface water detention. On the other hand, Uddin et al. [26] mapped the change detection of the Koshi Basin wetlands between the years of 1990 and 2010 using an ERDAS Imagine software analysis.
One of the largest CW facilities in the world is in Oman at the Nimr oilfield [27]. This CW facility receives brackish produced water generated during oil exploration activities in the nearby oilfield. This massive CW in the desert of Oman provides a series of ecosystem services since it represents a green oasis that attracts thousands of migratory birds and improves and regulates the local microclimate [27]. The CW system has a zero-energy demand for the treatment, since natural treatment processes for the transformation and removal of pollutants occur naturally in this nature-based technology, while gravity flow is applied across the wetland cells. The cells have equal surfaces (approx. 10 hectares each) per the initial design of the facility to distribute the flow equally. The flow enters each cell at the upstream point (from the buffer pond) and flows with gravity along the cell series to the downstream evaporation ponds naturally without the use of an external power; it flows just by using the local topography. To optimize its effectiveness and minimize any negative impact on the environment, a human-made wetland was created based on the natural characteristics of the location. The shape of the wetland was determined by the pre-existing topography, geology, and the amount of available land. The number of cells required for the wetland would depend on factors such as the topography, hydrology, and water quality. In areas with a level surface, cells can be constructed with the use of dikes, while on a sloping terrain, terraced cells can be employed [28].
The CW is planted with different reed species, making the wetland a polyculture [16]. Over its 10-year operation, it has been observed that some of the planted reed species had a better tolerance to the water salinity levels and water quality in different seasons than other species that suffered from water salinity stresses at different times of the year and in different locations of the wetland. Furthermore, its indicated that plants thriving in contaminated water are expected to experience greater levels of stress compared to plants growing in unpolluted areas within the same ponds [29].
The CW facility is located under desert environmental conditions, treating an industrial effluent, while other smaller CWs in the area treat domestic wastewater [7], highlighting the capacity of this green technology to provide wastewater treatment services under the harshest environments. After more than 10 years of operation and three consecutive expansion phases, today, CW cells cover an area of 490 hectares, followed by 780 hectares of evaporation ponds [16].
The treatment capacity of the system is 175,000 m³/day, representing about 65% of the total oily produced water generated at that oilfield [16] or half of the daily water consumption of the Muscat governorate, the capital of Oman [27,28]. Various sustainable activities also take place in that facility such as the reuse of the treated effluent for the irrigation of commercial plants, compost production using the reed vegetation, etc.z [5,16,27]. Compared to the previous high-cost and energy-intensive practice of deep well injection for produced water management, the CW system provides not only an excellent treatment (effluent oil-in-water content is below 0.5 mg/L), but also a unique environmental performance, with more than 99% reduced greenhouse gas emissions [27].
Initially, phragmites australis reeds were planted as the main and only plant species. However, different reed species were later added in the system, such as Typha domingensis, Schoenoplectus littoralis, Juncus rigidus, and Cyperus spp., to enhance the vegetation production and the resilience and health of the ecosystem. The main biodegradation mechanism is based on cyanobacterial mats trapping oil and degrading hydrocarbons [30].
The data presented in Table 1, which was produced by the authors of [27,31], demonstrate the effectiveness of the CW system in removing a range of pollutants from the wastewater. The removal of OiW, suspended solids, BOD, and nutrients is attributed to the various physical, chemical, and biological processes that occur within the wetland. The results indicate that CW systems can provide an environmentally sustainable and cost-effective solution for treating wastewater and mitigating water pollution [27].
Given the large size of the CW, it is difficult to determine the species that suffered from salinity stress more than others and at what time and place this stress occurred. The use of remote sensing technology can address this difficulty and help in assessing the determination of such reed species temporally and spatially. Moreover, to our knowledge, no remote sensing studies have been conducted on large-scale CW systems or CWs treating oily produced water, especially under a hot and arid climatic region. This study focuses on assessing the health of plants in a massive constructed wetland in Oman. This wetland treats water from oil exploration in a challenging desert environment, exposing the plants to stress from the oil and high salinity. This research examines plant health and changes in the vegetation between 2017 and 2019 using data from Sentinel-2 satellites. Spectral indices including NDVI, MSAVI, and NDSI were used to analyze the plants’ behaviors with these stressors.
Therefore, the objectives of this study were to use remote sensing technology to (1) conduct a change detection analysis for the CW area in Oman, (2) evaluate the performance of the nature-based treatment system in this project, and (3) examine and determine the suffering level of the existing reed species over different seasons of the year.

2. Materials and Methods

2.1. Study Site and Weather Conditions

The CW treating oily produced water is located at the Nimr oilfield in the southern part of Oman (18°40′8.99″N and 55°46′39.55″E) and in a desert environment about 700 km southwest of Muscat (Figure 1). The temperature in the area can reach as high as 50–55 °C in the summer period (June-August), with a respective average monthly temperature between 30 and 32 °C, and average values for the colder period (December–February) between 20 and 21 °C [27]. The annual average precipitation is practically negligible (approx. 30 mm).
A constructed wetland typically consists of a shallow basin containing a substrate, commonly soil or gravel, and is planted with vegetation that can withstand saturated conditions [32]. The CWs are designed with an underlying mineral layer composed of locally sourced soil material. This design aims to mitigate environmental consequences and expenses associated with the high-density polyethylene liner positioned above it [28]. The mineral sealing layer was applied to all wetland cells, whereas the high-density polyethylene liner was exclusively employed within the inlet buffer pond. The inclusion of soil was imperative to offer a suitable medium for the growth of plant roots. A layer of 35 cm of finely screened red soil (ranging from 0 to 20 mm in size) was incorporated above the construction material, facilitating the cultivation of reed species [27,28]. Thin and variable soil layers can present challenges for traditional soil testing methods, making it less practical to carry out accurate and meaningful tests [33].
The CW started its operation in December 2010 [27], and the third expansion phase which was completed in 2019, increasing its treatment capacity to 175,000 m³/d of produced water. The oily produced wastewater is sent via a pipeline to the facility; in the first treatment stage, the separation and recovery of most of the oil content (>85%) takes place in passive hydro cyclones without the use of energy or chemicals. Then, the water is distributed into the surface flow CW cells via a long buffer channel without the use of pumps [16]. The water flows through the wetlands via gravity, and the treated effluent flows into a series of evaporation ponds, where water evaporation results in salt formation that can be processed into industrial-grade salt.
The total wetland area of 490 ha is divided into three phases named Phases A, B, and C, respectively (Figure 1). Phase A has 24 cells (A1–A24), Phase B has 12 cells (B1–B12), and Phase C has 10 cells (C1–C10). The water inlet from the first two cells in Phase A and Phase B is called terrace one, while terrace two starts from mid evaporation ponds. Phases A and B are built on two terrace stages with 9 cell lines of inlet water flow (i.e., A1–A19, A4–A22, B1–B9, and B5–B12) that allow the gravity flow of water from the higher terrace to the lower terrace. Similarly, in Phase C, each terrace contains two cells, and there are five lines of cells for the inlet water flow in each terrace (C1–C6 and C3–C8). The water in Phases A and B flows along 9 tracks, each with 4 cells, e.g., A1-A7-A13-A19, A6-A12-A18-A24, B1-B3-B6-B9, and B5-B8-B11-B12. In Phase C, there are 5 tracks, each with 2 cells.

2.2. Satellite Images

This study used a total of twelve satellite images acquired from the Sentinel-2 satellite system for the respective periods of 2017 and 2019. These two years were selected because in 2018, the construction of the third expansion (expansion C) took place. Accordingly, we chose a study period comprising one year before the expansion and another year during the expansion based on this criterion. The gathered water parameters were synchronized with the collected satellite images on a bi-monthly basis for the years 2017 and 2019, ensuring a precise correlation that accurately reflected plant stress. Plants can be stressed due to any unfavorable condition or substance that affects or blocks the plant’s metabolism, growth, or development, which can be induced by various natural and anthropogenic stress factors. The accessibility of optical remote sensing information has advanced with the Sentinel-2A and Sentinel-2B satellites, which capture multispectral data comprising 13 bands across the visible and shortwave infrared (SWIR) wavelength ranges. This enhancement offers enhanced opportunities for monitoring agricultural fields [34].
The Sentinel-2 initiative, a program within the Copernicus project, captures optical imagery characterized by a high spatial resolution ranging from 10 to 60 m, systematically gathering data over land and coastal bodies of water for the purpose of Earth observation. The pre- and post-processing as well as the classification of the satellite images were conducted using the ERDAS Imagine software. The imagery consists of a set of twelve spectral bands that provide coverage over a strip of land measuring 290 km in width. Comprehensive information about the dates and specifications of the obtained satellite imagery is presented in Table 2.
Remote sensing techniques were used to identify the change detection over the years 2017 and 2019 by using the acquired satellite images with different spectral index information including the Normalized Difference Vegetation Index (NDVI) [35], Modified Soil-Adjusted Vegetation Index (MSAVI) [36], and Normalized Difference Salinity Index (NDSI) [37]. The equations used to calculate these indices are as follows:
N D V I = N I R R E D N I R + R E D
M S A V I = 2 N I R + 1 2 N I R + 1 2 8 ( N I R R E D ) 2
N D S I = R E D N I R R E D + N I R

2.3. Performance Evaluation

This study employs the NDVI as a numerical index of the red and near-infrared spectral bands that are highly correlated with the vegetation content. Higher NDVI values indicate areas with denser and healthier plants, as they reflect more in the near-infrared spectrum. To account for the influence of the soil background, especially in areas with low vegetative cover, the MSAVI is used to correct the NDVI outcomes. Furthermore, the NDSI is employed to delineate stressed reed plants due to the salinization effect. This research also shows the change detection during different seasons (winter and summer) between both years. This study involved the computation of the average values of each spectral index within the constructed wetland (CW), which are subsequently categorized into three distinct classes based on the index values: unhealthy (with values below 0.3), medium health (ranging between 0.3 and 0.5), and healthy (exceeding 0.5). The stress index values of plants within the specified ponds were computed at both the inlet and outlet of the water flow utilizing the index output derived from ERDAS Imagine. These values were averaged and represented as the Average Stress Index, aiming to establish a correlation with the water flow parameters.
To analyze the water quality parameters of the collected produced water samples from the CW inlet and outlet, the Central Analytical and Applied Research Unit, a certified laboratory at Sultan Qaboos University in Muscat, Oman, was used. The analyzed water parameters included the water temperature, pH, oil in water contamination (OiW), salinity, conductivity, dissolved oxygen, and oxidation–reduction potential. Multiparameter Water Quality Meters and Probes and ICP spectroscopy were the tools used to measure the collected monthly CW ponds’ water samples. The determination coefficient (R2) was computed to furnish a statistical measure of the regression goodness of fit between the water parameters and the vegetation index.
This research employed multivariate analysis techniques, specifically Multiple Linear Regression. The coefficient of determination, denoted as R2, was utilized to assess the correlation coefficient, providing an evaluation of the dispersion of data points in relation to the fitted regression line for the water parameters and index test. This study integrated field data, ancillary water flow information, and high-resolution satellite imagery to detect and highlight vegetation stress. Figure 2 provides a summary of the approach employed in this study.

3. Results and Discussion

3.1. Plant Stress Detection

In this study, an imaging index system was established to appraise plant reactions to stress induced by climatic changes and variations in water parameters across distinct seasons. Employing remote sensing and GIS processing tools along with Microsoft Excel, the imaging data underwent a thorough analysis. This involved the creation of both qualitative false color representations and quantitative graphs, facilitating the discernment of variations in the plant stress levels.
To create a false color image, an initial step involves the generation of a color index file derived from the raw index image. This file preserves the calculated index score for each pixel, denoted as “color index values”. The index scores span from −1.0 to 1.0, with each color index value category corresponding to an incremental range of these index scores. By utilizing this incremental system, an index score of −1.0 denotes an unhealthy status, while a score of 1.0 corresponds to a color index value representing a healthy state. The creation of false color images involves the application of a color lookup table (as illustrated in Figure 3 and Figure 4) to the color index file image. These images exhibit artificially generated color schemes that can be customized to highlight pixels with specific color index values in distinct colors, serving the user’s preferences. To produce false color images using the Photo Monitoring plugin, specific settings are applied. These settings involve selecting the option to stretch the near-infrared (NIR) band before generating the index. Additionally, the minimum index value for scaling the color index image is set to −1.00, while the maximum is set to 1.00.
The assessment of plant stress was conducted across various seasons within the specified years, involving a quantitative analysis of the targeted indices juxtaposed with the concurrent calculations of the water flow parameters. This evaluation discerns the impact of high or low scoring index values and their correlations with the water parameters, unveiling the status of healthy or unhealthy plant tissues. Statistical analyses were carried out by employing both Excel’s data analysis tool and the t-test function in R, yielding consistent results. At the conclusion of each year, subsequent to the exposure of targeted ponds, the shift in the peak color index values in the plants exhibited statistically significant differences in comparison to the alterations observed in the peak color index values associated with the water flow parameters.
The dense reed vegetation and the characteristics of the reed clusters make the monitoring and evaluation of each species more challenging with low-resolution mapping images. Therefore, the monitoring and evaluation of the wetland vegetation must be conducted at the scale of the entire reed cluster. In addition, since reed clusters often contain multiple reed species, the differentiation between the different species becomes even more difficult.
Remote sensing methodologies prove to be valuable for detecting alterations in vegetation cover and plant stress, with the spectral indices employed in this investigation (NDVI, MSAVI, and NDSI) demonstrating their effectiveness in capturing information pertaining to salinity, plant stress, and system degradation details induced by produced water. The differences observed between the VSWIR index and other indices indicate the potential of spectral indices to enhance and delineate vegetation stress cover details in an image. The findings of the multivariate analysis conducted on the targeted years indicate moderate saline stress in the wetlands, which is highlighted in the intensity of salinization in different areas depicted in the generated maps. The stressed areas with less vegetation cover can be identified through the remote sensing maps (Figure 3 and Figure 4).
The utilization of Sentinel-2 imagery and the spatial resolution of the data, in this study, demonstrated significant effectiveness in accurately delineating the distribution of wetland plant vegetation. To evaluate plant performance in the CW system during the years 2017 and 2019, the similarities and differences between the two years were investigated. The stress on the same cells appeared to be similar in both years, with the plants in the upstream cells close to the produced water inflow exhibiting more stress than those in the downstream cells, as illustrated in Figure 5 and Figure 6 for the year 2017.
The cells that showed a weakness in plants growing on their top parts compared to other cells in the CW system were Phase A cells (A1 to A6), Phase B cells (B1, B2, and B5), and Phase C cells (C1 to C5). These cells received inflow with high salinity and hydrocarbon contents, which likely contributed to the observed stress, as indicated by the Vegetation Near-Infrared and Short-Wave Infrared indices.
The research also showed changes in the NDVI and MSAVI indices. The results indicated that the NDVI and MSAVI indices can capture information about plant stress, salinity, and system degradation details. The spectral indices were found to have immense potential in enhancing and delineating vegetation stress cover details in an image. Thus, the Sentinel-2 imagery and the NDVI and MSAVI indices proved to be valuable tools in assessing the response of wetland plant vegetation to changes in the water quality parameters.

3.2. Numerical Analysis

Through the application of laboratory testing on produced water during both the summer and winter seasons, the numerical values of the indices were effectively merged with the corresponding values of the water quality parameters. This integration facilitated a comprehensive understanding of the vegetation conditions prevailing at the CW, thereby providing valuable insights into the complex interactions between the water quality and ecological systems. The outcomes of this approach were instrumental in identifying and assessing the impact of seasonal variations on water quality, as well as enabling the development of informed strategies for sustainable water resource management in the context of the CW ecosystem.
The Modified Soil-Adjusted Vegetation Index (MSAVI) was utilized to rectify the outcomes of the NDVI in the presence of a low vegetative cover and high soil background. By calculating the average maximum MSAVI values of each cell in the CW during both 2017 and 2019, the studied area was classified into three categories based on vegetation health, namely unhealthy for an MSAVI value of <0.3, medium for 0.3 ≤ MSAVI ≤ 0.5, and healthy for an MSAVI value of >0.5. The distribution of the MSAVI served as an example of the classification of plants’ health conditions, as depicted in Figure 7. Notably, most cells near the outlet exhibited the presence of unhealthy plants under stress due to elevated water salinity levels. Conversely, cells close to the inlet, such as A2, A3, A5, A6, B1, B2, and B5, displayed healthier plants, as inferred from the MSAVI values. This trend was consistent over the two-year period from 2017 to 2019. However, during 2019, some cells malfunctioned, with A3 being out of operation for maintenance reasons, which led to it being dried out. A similar situation was also observed in cells A4 and A9 during the late summer of 2019.
Figure 8, Figure 9 and Figure 10 depict the outcomes of the average water flow parameters in the wetland cells of both Phase A and Phase B. The water flow parameters were measured from the first CW cells receiving the produced water inflow to the last cells that discharge water to the evaporation ponds during the winter and summer times of 2017 and 2019. The findings revealed that the water temperature and oil in water (OiW) levels were higher in the inlet CW cells of the first cells, and a subsequent reduction was observed as the produced water flowed towards the last cells in terrace two. Furthermore, the treated produced water showed a further reduction in the temperature and OiW levels upon discharge from the last CW cells to the evaporation ponds. In contrast, the water conductivity levels exhibited an increasing trend as the produced water flowed from the first cells at the first terrace to the last cells at terrace two and subsequently towards the evaporation pond.
The impact of different seasons (winter and summer) on the water flow parameters was assessed for the years 2017 and 2019. The findings revealed that the water temperature levels were higher in the inlet CW cells of the first cells (A1 and B5), whereas a decrease of 3–9 °C was observed as the produced water flowed towards the last cells (A19 and B12) across the reed beds in both evaluated years (Figure 8). The level of oil contamination in the produced water was similarly reduced by 5–23 mg/L (Figure 9), and it nearly approached zero as the produced water reached the last cells (A19 and B12), indicating the effectiveness of the CW system. Conversely, the water conductivity levels showed an increase of 2–3 mS/cm as the produced water moved from the inlet at the first cells towards the last cells (Figure 10).
The data presented in Figure 11, Figure 12, Figure 13, Figure 14, Figure 15 and Figure 16 and in Table 3 and Table 4 demonstrate the influences of various water quality parameters on two seasons in 2017 and 2019, namely the water temperature, oil contamination, water conductivity, dissolved oxygen, and vegetation indices, including the NDVI, NDSI, and MSAVI. The data were collected from the CW tracks of the first cells in Phase A and the last cells in Phase B. The results indicate a positive correlation between the vegetation indices, specifically the NDVI, and water flow parameters. Strong relationships were found between the NDVI and oil contamination in the first cells (A1 to A19) during both seasons of 2017, with an R2 of 0.87 in the winter and 0.72 in the summer, and 2019, with an R2 of 0.94 in the winter and 0.66 in the summer. Furthermore, there was a correlation between the NDVI and water temperature in the first cells during the season of 2019, with an R2 of 0.915 in the winter and 0.704 in the summer (Figure 12 and Table 3).
In contrast, the last cells in Phase B (B5 to B12) showed a strong correlation between the NDVI and water conductivity in 2017, with an R2 of 0.98 in the winter and 0.56 in the summer (Figure 11 and Table 4). However, no significant associations were observed in 2019 during any of the seasons for any of the water flow parameters (Figure 12 and Table 4). Additionally, the NDSI showed a strong correlation with water conductivity in both seasons of both targeted years, 2017 and 2019, with an R2 of 0.95 in the winter and 0.87 in the summer for 2017 and an R2 of 0.704 in the winter and 0.766 in the summer for 2019 (Figure 13 and Figure 14 and Table 4).
Overall, these findings suggest that the NDVI and NDSI are useful indicators for water quality assessment and monitoring, particularly in relation to oil contamination and water conductivity, respectively. The results also indicate that the correlation between the vegetation indices and water flow parameters varies depending on the season and location within the CW tracks, emphasizing the importance of considering the temporal and spatial variability when conducting water quality assessments.
The analysis of the vegetation indices (NDSI) in both the winter and summer seasons of the targeted years of 2017 and 2019, along the cells of the CW tracks of the first cells in Phases A did not show significant correlations with the water flow parameters of the first cells (A1 to A19). However, the last cells in Phase B (B5 to B12) exhibited a strong correlation between the vegetation indices (NDVI) and water conductivity in the year of 2017 (R2:0.95/winter and 0.87/summer) and with the dissolved oxygen in the same year (R2:0.818/winter and 0.482/summer), as indicated in Figure 11 and Table 4. Moreover, the water conductivity demonstrated a correlation between the last cells in Phase B (B5 to B12) and the NDSI indices in the year of 2019 (R2:0.70/winter and 0.766/summer), as presented in Figure 14 and Table 4. These findings suggest that the NDVI and water conductivity are useful indicators for water quality assessment and monitoring in the last cells of Phase B. However, the NDSI did not show any significant correlations with the water flow parameters in the first cells of Phase A.
The Modified Soil-Adjusted Vegetation Index (MSAVI) was evaluated for both seasons (winter and summer) of the years 2017 and 2019 along the cells of the CW in Phase A and the last cells in Phase B. The results showed that the MSAVI is correlated with the OIW flow parameters of the first cells (A1 to A19) in both years in Phase A with the R2 values of 0.608 (winter) and 0.626 (summer) in the year of 2017, and R2 values of 0.947 (winter) and 0.567 (summer) in the year of 2019, as shown in Figure 15 and Figure 16 and Table 3 and Table 4.
Moreover, the last cells in Phase B (B5 to B12) showed a good correlation between the MSAVI and water conductivity in the year of 2017 with R2 values of 0.977 (winter) and 0.5866 (summer), and in the year of 2019 with R2 values of 0.568 (winter) and 0.218 (summer), as shown in Figure 15 and Figure 16 and in Table 3 and Table 4.
The findings of the study suggest that both the MSAVI and NDSI exhibited similar reflectance correlations with the water flow parameters. The MSAVI was found to be useful in correcting the outcomes of the NDVI by accounting for the influence of the soil background, especially in areas with a low vegetation cover. Additionally, the NDSI was effective in identifying the salinity stress on reed plants and was able to capture this stress through the correlation between the NDSI and water conductivity in the last cells of Phase B (B5 to B12). These results are presented in Figure 13 and Figure 14 as well as in Table 3 and Table 4.
Figure 17 and Figure 18 and Table 3 and Table 4 provide an overview of the R2 values obtained for the water collection parameters in both 2017 and 2019 for the winter and summer seasons. The results highlight the usefulness of the correlation indicators in water quality assessment and monitoring. The figure indicates that the vegetation index NDVI showed higher average R2 values compared to the other parameters (Figure 18).

3.3. Discussion

The CW system is responsible for a significant increase in the total dissolved solids (TDS) concentration in the effluent. The high TDS concentration is attributed to the large surface area of the CW and the dense reed vegetations that result in high water losses through evapotranspiration [27]. The inflow TDS concentration is around 7–8000 mg/L, while the effluent TDS concentration exceeds 12,000 mg/L [28]. As the water flows with gravity from the upstream to the downstream wetland cells, a gradual increase in salinity is detected. Boron, which is relatively high in concentration in the inflow (Table 1), might also have an impact on plant growth, along with salinity. However, the boron levels generally remain lower than the reported level, which might induce phytotoxicity (i.e., 7 mg/L) [27]. Despite these conditions, the reed plants manage to survive and grow well, establishing, in time, a rich vegetation with a high density.
Figure 3 and Figure 4 present a bi-monthly estimation of plant stress change detection in the CW for the years 2017 and 2019. The analysis was conducted using the Vegetation Short-Wave Infrared (VSWIR) index and spectral indices (NDVI, MSAVI, and NDSI) obtained from Sentinel resolution imaging. The outcomes of the change detection analysis revealed noteworthy alterations in the vegetation cover within specific cells and more subtle variations in the others, alongside changes in the land cover during the two designated years. Concretely, cells A3, A4, A9, A10, A13, A15, A16, A19, A21, and A22 in Phase A demonstrated disparities between the two years, while the downstream cells of Phase B, specifically cells B8, B10, B11, and B12, displayed detectable changes.
The produced water samples were taken at the inlet of the cells and at the outlet of the last cells, with the first cells starting from A1, A7, A13, and A19 and leaving to the evaporation ponds, while the inlet of the last cells of phase B starts at B5, B8, B11, and B12 and leaves to the evaporation ponds.
The average content of oil in water (OiW) in the inflow was found to be approximately 280 mg/L, with occasional spikes that exceeded 500 mg/L, as reported in earlier studies [27,28]. However, the treatment process in the CW system resulted in a significant reduction in the OiW, with levels dropping to below 0.5 mg/L in the treated effluent. The average concentration of suspended solids was reduced from 28 mg/L in the inflow to 10 mg/L in the treated effluent, while the average BOD concentration dropped from 16 mg/L to less than 1 mg/L. This reduction in the BOD is particularly noteworthy, as it indicates the removal of biodegradable organic matter from the wastewater.
The CW system also showed an efficient removal of nutrients, such as nitrogen and phosphorus, which are typically found in wastewater. However, the concentration of these nutrients in the inflow was already low, with the levels not exceeding 2–3 mg/L. Consequently, the treated effluent did not have measurable concentrations of nitrogen and phosphorus.
During the operational lifespan of the wetlands, a discernible observation emerged, indicating that certain reed species, when compared to others, exhibited heightened resilience to varying water salinity levels and fluctuations in water quality across different seasons and diverse locations within the wetland. Conversely, several other species encountered challenges due to water salinity stresses, displaying susceptibility at different periods throughout the year and at various sections within the wetland [16]. The analysis of change detection in the constructed wetland (CW) in Oman provided valuable insights into the variations in wetland stress across the specific seasons and years under examination. The wetland showed that plants could be grown economically with irrigation-treated water of up to EC 16 mS cm (Figure 8). The greater application of produced water positively affected plant growth by accumulating a toxic level of salts in the lower soil horizons. The use of diluted seawater for barley irrigation is only possible if the leaching of excess salt from the root zone is implemented. However, the soil and plant data found in this study were in the same line of many published data, where the plants underwent stress as the salinity, temperature, and oil in water increased [15].
In addition to employing indices pertaining to the chlorophyll content, the integration of water flow parameters played a pivotal role in the estimation of various botanical attributes, encompassing photosynthetic pigments, water, nitrogen, and leaf mass per area (e.g., NDSI) [34]. Consequently, the inclusion of vegetation indices exerted a more pronounced impact on the classification results compared to the original short-wave infrared reflectance. Nevertheless, discernible distinctions in relative importance emerged among the three tested indices (NDVI, NDSI, and MSAVI). The computed importance values, denoted by R2, served as valuable indicators when correlating the indices with water flow, leading to a substantial enhancement in the accuracy of identifying stressed croplands. The research utilized cost-free satellite images to evaluate stress in reed grass plants induced by hydrocarbon water flow contamination within the wetland. Moreover, the study area experiences an extremely harsh climate, making it uncommon to find comparable research references for the assessment of spectral indices. A higher spectral resolution sensor image could potentially provide a more accurate estimation of the stress content in the plants [29].
In the third phase (C) of the study, which involved the cultivation of new plants, it was observed that the same stress challenges associated with inflowing water were present. This was evidenced by the appearance of stress on the new CW cells of Phase C within four months of planting, as was also observed in Phases A and B, regardless of the plant species. These findings highlight the presence of high levels of salinity and hydrocarbon stress during the first stage of produced water inflow. In addition, the NDVI and MSAVI revealed variations in plant health across Phase A and Phase B cells, with some cells exhibiting higher levels of vegetation stress than others, as illustrated in Figure 5 and Figure 6 for the year 2017. However, when new reed species were planted in Phase C, it was difficult to determine definitively whether the impact of produced water inflow stress on the plants was high or low within a year’s time. As such, further investigation is necessary in Phase C to detect stress more comprehensively on the plants.
This study represents a measure towards the potential of remote sensing in the treatment of hydrocarbon-contaminated water in wetlands via using correlation indicators in water quality assessment and monitoring. The vegetation indices showed their usefulness as tools for assessing water quality and monitoring the impact of various water quality parameters on the environment. By analyzing the R2 values, it is possible to determine the extent to which changes in the water quality parameters affect the vegetation index, which can help identify areas that require further investigation and remediation efforts. These findings have significant implications for environmental monitoring and management, as they provide a valuable means of assessing the impact of water quality on the environment and identifying areas that require targeted interventions to improve the water quality.
In addition, the use of correlation indicators such as the NDVI can help identify the specific water quality parameters that are affecting vegetation growth and health. This information can be used to inform targeted interventions that address the specific water quality issues in each area. Furthermore, the use of correlation indicators can help prioritize areas for further investigation and remediation efforts, which can be particularly useful in resource-limited settings.
The findings of this study suggest that correlation indicators such as the NDVI can be a valuable tool in water quality assessment and monitoring. By identifying areas that require targeted interventions, these indicators can help improve environmental management and protect ecosystems from the negative impacts of poor water quality.

4. Conclusions

The results of the multivariate analysis conducted on the targeted years indicate the presence of moderate saline stress in the wetlands. This study represents a novel attempt to investigate various vegetation spectral indices in a large constructed wetland (CW) system that treats oily produced water in a hot and arid climate using remote sensing techniques. The comparison of vegetation spectral indices between 2017 and 2019 revealed consistent plant stress along the water flow path, despite slight differences in the land cover between the two years. These findings suggest that the wetland plants were under stress during the study period.
The results of this study demonstrate the feasibility of estimating CW vegetation stress over large areas using remote sensing techniques, provided that sufficient field data, satellite imagery, and ancillary data are available. Given the dynamic nature of the CW system investigated in this study, the use of satellite-based remote sensors and low-cost GIS tools is necessary for effective management and monitoring. In recent years, satellite imagery has provided increasingly precise data in advanced forms, which can assist in risk assessment and prediction in various environmental contexts and study areas. The combination of the NDVI with the MSAVI and NDSI can reduce uncertainty in differentiating between plant stress and plant density, thereby improving the accuracy of vegetation stress estimates.
Overall, this study highlights the potential of remote sensing techniques for monitoring vegetation stress in constructed wetland systems. By providing valuable insights into the health of wetland plants, these techniques can aid in effective environmental management and promote the sustainable use of natural resources.

Author Contributions

Conceptualization and methodology, Y.A.-M. and K.A.-J.; writing—original draft preparation, K.A.-J. and Y.A.-M.; writing—review and editing, K.A.-J., Y.A.-M., F.M. and A.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Millennium Ecosystem Assessment. Ecosystems and Human Well-Being; Island Press: Washington, DC, USA, 2005; Volume 5. [Google Scholar]
  2. Valderrama-Landeros, L.; Flores-de-Santiago, F.; Kovacs, J.M.; Flores-Verdugo, F. An assessment of commonly employed satellite-based remote sensors for mapping mangrove species in Mexico using an NDVI-based classification scheme. Environ. Monit. Assess. 2018, 190, 23. [Google Scholar] [CrossRef]
  3. Hussain, S.A.; Badola, R. Valuing mangrove benefits: Contribution of mangrove forests to local livelihoods in Bhitarkanika Conservation Area, East Coast of India. Wetl. Ecol. Manag. 2010, 18, 321–331. [Google Scholar] [CrossRef]
  4. Wu, Q. GIS and remote sensing applications in wetland mapping and monitoring. In Comprehensive Geographic Information Systems; Elsevier: Oxford, UK, 2017. [Google Scholar]
  5. Oral, H.V.; Radinja, M.; Rizzo, A.; Kearney, K.; Andersen, T.R.; Krzeminski, P.; Buttiglieri, G.; Ayral-Cinar, D.; Comas, J.; Gajewska, M.; et al. Management of urban waters with nature-based solutions in circular cities—Exemplified through seven urban circularity challenges. Water 2021, 13, 3334. [Google Scholar] [CrossRef]
  6. Abou-Elela, S.I.; Golinielli, G.; Abou-Taleb, E.M.; Hellal, M.S. Municipal wastewater treatment in horizontal and vertical flows constructed wetlands. Ecol. Eng. 2013, 61, 460–468. [Google Scholar] [CrossRef]
  7. Al-Wahaibi, B.M.; Jafary, T.; Al-Mamun, A.; Baawain, M.S.; Aghbashlo, M.; Tabatabaei, M.; Stefanakis, A.I. Operational modifications of a full-scale experimental vertical flow constructed wetland with effluent recirculation to optimize total nitrogen removal. J. Clean. Prod. 2021, 296, 126558. [Google Scholar] [CrossRef]
  8. Schultze-Nobre, L.; Wiessner, A.; Bartsch, C.; Paschke, H.; Stefanakis, A.I.; Aylward, L.A.; Kuschk, P. Removal of dimethylphenols and ammonium in laboratory-scale horizontal subsurface flow constructed wetlands. Eng. Life Sci. 2017, 17, 1224–1233. [Google Scholar] [CrossRef]
  9. Stefanakis, A.I. Introduction to Constructed Wetland Technology; John Wiley & Sons, Ltd.: Chichester, UK, 2018. [Google Scholar]
  10. Ramírez, S.; Torrealba, G.; Lameda-Cuicas, E.; Molina-Quintero, L.; Stefanakis, A.I.; Pire-Sierra, M.C. Investigation of pilot-scale constructed wetlands treating simulated pre-treated tannery wastewater under tropical climate. Chemosphere 2019, 234, 496–504. [Google Scholar] [CrossRef]
  11. Sultana, M.-Y.; Akratos, C.S.; Vayenas, D.V.; Pavlou, S. Constructed wetlands in the treatment of agro-industrial wastewater: A review. Hem. Ind. 2015, 69, 127–142. [Google Scholar] [CrossRef]
  12. Tatoulis, T.; Stefanakis, A.; Frontistis, Z.; Akratos, C.S.; Tekerlekopoulou, A.G.; Mantzavinos, D.; Vayenas, D.V. Treatment of table olive washing water using trickling filters, constructed wetlands and electrooxidation. Environ. Sci. Pollut. Res. 2017, 24, 1085–1092. [Google Scholar] [CrossRef]
  13. Stefanakis, A.I. Constructed Wetlands for Wastewater Treatment in Hot and Arid Climates, 1st ed.; Springer Nature: Cham, Switzerland, 2022. [Google Scholar]
  14. Mahdavi, S.; Salehi, B.; Granger, J.; Amani, M.; Brisco, B.; Huang, W. Remote sensing for wetland classification: A comprehensive review. GIScience Remote Sens. 2018, 55, 623–658. [Google Scholar] [CrossRef]
  15. Mozaffari, M.-H.; Shafiepour, E.; Mirbagheri, S.A.; Rakhshandehroo, G.; Wallace, S.; Stefanakis, A.I. Hydraulic characterization and removal of metals and nutrients in an aerated horizontal subsurface flow “racetrack” wetland treating primary-treated oil industry effluent. Water Res. 2021, 200, 117220. [Google Scholar] [CrossRef] [PubMed]
  16. Stefanakis, A.I. Constructed wetlands for sustainable wastewater treatment in hot and arid climates: Opportunities, challenges and case studies in the Middle East. Water 2020, 12, 1665. [Google Scholar] [CrossRef]
  17. Guo, M.; Li, J.; Sheng, C.; Xu, J.; Wu, L. A review of wetland remote sensing. Sensors 2017, 17, 777. [Google Scholar] [CrossRef] [PubMed]
  18. Al-Mulla, Y.; Al-Ruheili, A.; Al-Lawati, A.; Parimi, K.; Ali, A.; Al-Sadi, N.; Al-Harrasi, F. Assessment of urban expansion’s impact on changes in vegetation patterns in Dhofar, Oman, using remote sensing and GIS techniques. IEEE Access 2022, 10, 86782–86792. [Google Scholar] [CrossRef]
  19. Al-Awadhi, T.; Al-Shukili, A.; Al-Amri, Q. The Use of Remote Sensing & Geographical Information Systems to Identify Vegetation: The Case of Dhofar Governorate (Oman). Available online: http://www.isprs.org/proceedings/2011/ISRSE-34/211104015Final00239.pdf (accessed on 10 September 2011).
  20. Abood, S.; Maclean, A.; Falkowski, M. Soil Salinity Detection in the Mesopotamian Agricultural Plain Utilizing WorldView-2 Imagery; Michigan Technological University: Houghton, MI, USA, 2011. [Google Scholar]
  21. Lu, B.; He, Y.; Dao, P.D. Comparing the performance of multispectral and hyperspectral images for estimating vegetation properties. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2019, 12, 1784–1797. [Google Scholar] [CrossRef]
  22. Beisel, N.S.; Callaham, J.B.; Sng, N.J.; Taylor, D.J.; Paul, A.L.; Ferl, R.J. Utilization of single-image normalized difference vegetation index (SI-NDVI) for early plant stress detection. Appl. Plant Sci. 2018, 6, e01186. [Google Scholar] [CrossRef]
  23. Dronova, I. Object-based image analysis in wetland research: A review. Remote Sens. 2015, 7, 6380–6413. [Google Scholar] [CrossRef]
  24. Singh, R.P.; Singh, N.; Singh, S.; Mukherjee, S. Normalized difference vegetation index (NDVI) based classification to assess the change in land use/land cover (LULC) in Lower Assam, India. Int. J. Adv. Remote Sens. GIS 2016, 5, 1963–1970. [Google Scholar] [CrossRef]
  25. Pan, F.; Xie, J.; Lin, J.; Zhao, T.; Ji, Y.; Hu, Q.; Pan, X.; Wang, C.; Xi, X. Evaluation of climate change impacts on wetland vegetation in the Dunhuang Yangguan National Nature Reserve in Northwest China using Landsat derived NDVI. Remote Sens. 2018, 10, 735. [Google Scholar] [CrossRef]
  26. Uddin, K.; Wahid, S.M.; Murthy, M.S.R. Mapping of Koshi Basin wetlands using remote sensing. In Proceedings of the 5th International Conference on Water & Flood Management, Dhaka, Bangladesh, 6–8 March 2015. [Google Scholar]
  27. Stefanakis, A.I.; Prigent, S.; Breuer, R. Integrated produced water management in a desert oilfield using wetland technology and innovative reuse practices. In Constructed Wetlands for Industrial Wastewater Treatment; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2018; pp. 23–42. [Google Scholar]
  28. PDO. 2019. Nimr Reed Beds. Technical Report. Available online: https://www.pdo.co.om/en/technical-expertise/nimr-reed-beds/Pages/default.aspx (accessed on 6 October 2022).
  29. Batkhuyag, O. Spectral Indicators for Assessing the Effect of Hydrocarbon Leakage on Vegetation; International Institute for Geoinformation Science & Earth Observation: Enschede, The Netherlands, 2008. [Google Scholar]
  30. Abed, R.M.M.; Al-Kharusi, S.; Prigent, S.; Headley, T. Diversity, distribution and hydrocarbon biodegradation capabilities of microbial communities in oil-contaminated cyanobacterial mats from a constructed wetland. PLoS ONE 2014, 9, e114570. [Google Scholar] [CrossRef]
  31. Alexandros, S.; Prigent, S.; Breuer, R. 6.5 CASE STUDY 4–NIMR WATER TREATMENT PLANT (OMAN); IWA Publishing: London, UK, 2019. [Google Scholar]
  32. Davis, L. A Handbook of Constructed Wetlands: A Guide to Creating Wetlands for: Agricultural Wastewater, Domestic Wastewater, Coal Mine Drainage, Stormwater. In the Mid-Atlantic Region. Volume 1: General Considerations; USDA-Natural Resources Conservation Service: Washington, DC, USA, 1995. [Google Scholar]
  33. Egbueri, J.C.; Igwe, O.; Omeka, M.E.; Agbasi, J.C. Development of MLR and variedly optimized ANN models for forecasting the detachability and liquefaction potential index of erodible soils. Geosyst. Geoenviron. 2023, 2, 100104. [Google Scholar] [CrossRef]
  34. Sonobe, R.; Yamaya, Y.; Tani, H.; Wang, X.; Kobayashi, N.; Mochizuki, K.-I. Crop classification from Sentinel-2-derived vegetation indices using ensemble learning. J. Appl. Remote Sens. 2018, 12, 026019. [Google Scholar] [CrossRef]
  35. Rouse, J.W., Jr. Monitoring the Vernal Advancement and Retrogradation (Green Wave Effect) of Natural Vegetation; NASA/GSFC Type III Final Report; NASA/GSFC: Greeenbelt, MD, USA, 1974; p. 371. [Google Scholar]
  36. Qi, J.; Chehbouni, A.; Huete, A.R.; Kerr, Y.H.; Sorooshian, S. A modified soil adjusted vegetation index. Remote Sens. Environ. 1994, 48, 119–126. [Google Scholar] [CrossRef]
  37. Odeh, I.O.A.; Onus, A. Spatial analysis of soil salinity and soil structural stability in a semiarid region of New South Wales, Australia. Environ. Manag. 2008, 42, 265–278. [Google Scholar] [CrossRef]
Figure 1. Illustration of the study area denoted with letter (X), which includes the CW site and its three phases: A, B, and C. Arrows illustrate direction of outflow from the three phases.
Figure 1. Illustration of the study area denoted with letter (X), which includes the CW site and its three phases: A, B, and C. Arrows illustrate direction of outflow from the three phases.
Remotesensing 15 05632 g001
Figure 2. Illustrates an overview of the study’s methodology and analytical procedures.
Figure 2. Illustrates an overview of the study’s methodology and analytical procedures.
Remotesensing 15 05632 g002
Figure 3. Sentinel images portraying the Oman CW over the course of the year 2017, featuring a series of six maps taken bi-monthly.
Figure 3. Sentinel images portraying the Oman CW over the course of the year 2017, featuring a series of six maps taken bi-monthly.
Remotesensing 15 05632 g003
Figure 4. Sentinel images portraying the Oman CW over the course of the year 2019, featuring a series of six maps taken bi-monthly.
Figure 4. Sentinel images portraying the Oman CW over the course of the year 2019, featuring a series of six maps taken bi-monthly.
Remotesensing 15 05632 g004
Figure 5. Plant stress signs at the inlet cells (inside the red circle) that receive the produced water.
Figure 5. Plant stress signs at the inlet cells (inside the red circle) that receive the produced water.
Remotesensing 15 05632 g005
Figure 6. (a) Depicts plant stress in the middle and downstream cells of phase (A & B) in 2017, and (b) illustrates the establishment of new reed plants in Phase C in 2019.
Figure 6. (a) Depicts plant stress in the middle and downstream cells of phase (A & B) in 2017, and (b) illustrates the establishment of new reed plants in Phase C in 2019.
Remotesensing 15 05632 g006
Figure 7. Demonstration of the distribution of MSAVI indices along the cells of the CW tracks in Phases A and B during the winter of the years (a) 2017 and (b) 2019.
Figure 7. Demonstration of the distribution of MSAVI indices along the cells of the CW tracks in Phases A and B during the winter of the years (a) 2017 and (b) 2019.
Remotesensing 15 05632 g007
Figure 8. The average water temperature (°C) of produced water from the initial cell to the last outlet cell during both the 2017 and 2019 seasons along the cells of the CW tracks, starting from the 1st cell in Phase A to the last cells in Phase B.
Figure 8. The average water temperature (°C) of produced water from the initial cell to the last outlet cell during both the 2017 and 2019 seasons along the cells of the CW tracks, starting from the 1st cell in Phase A to the last cells in Phase B.
Remotesensing 15 05632 g008
Figure 9. The average oil contamination (mg/L) in produced water from the initial cell to the final outlet cell during the 2017 and 2019 seasons along the cells of the CW tracks, spanning from the 1st cell in Phase A to the last cells in Phase B.
Figure 9. The average oil contamination (mg/L) in produced water from the initial cell to the final outlet cell during the 2017 and 2019 seasons along the cells of the CW tracks, spanning from the 1st cell in Phase A to the last cells in Phase B.
Remotesensing 15 05632 g009
Figure 10. The average water conductivity (mS/cm) in produced water from the initial cell to the final outlet cell during the 2017 and 2019 seasons along the cells of the CW tracks, ranging from the 1st cell in Phase A to the last cells in Phase B.
Figure 10. The average water conductivity (mS/cm) in produced water from the initial cell to the final outlet cell during the 2017 and 2019 seasons along the cells of the CW tracks, ranging from the 1st cell in Phase A to the last cells in Phase B.
Remotesensing 15 05632 g010
Figure 11. Comparison of water parameter data (WT-°C; OIW-mg/L; COND-mS/cm; and DO-mg/L) with NDVI along the cells of the CW tracks for (a) winter 2017 (first cell in Phase A), (b) winter 2017 (last cell in Phase B), (c) summer 2017 (first cell in Phase A), and (d) summer 2017 (last cell in Phase B).
Figure 11. Comparison of water parameter data (WT-°C; OIW-mg/L; COND-mS/cm; and DO-mg/L) with NDVI along the cells of the CW tracks for (a) winter 2017 (first cell in Phase A), (b) winter 2017 (last cell in Phase B), (c) summer 2017 (first cell in Phase A), and (d) summer 2017 (last cell in Phase B).
Remotesensing 15 05632 g011
Figure 12. Comparison of water parameter data (WT-°C; OIW-mg/L; COND-mS/cm; and DO-mg/L) with NDVI along the cells of the CW tracks for (a) winter 2019 (first cell in Phase A), (b) winter 2019 (last cell in Phase B), (c) summer 2019 (first cell in Phase A), and (d) summer 2019 (last cell in Phase B).
Figure 12. Comparison of water parameter data (WT-°C; OIW-mg/L; COND-mS/cm; and DO-mg/L) with NDVI along the cells of the CW tracks for (a) winter 2019 (first cell in Phase A), (b) winter 2019 (last cell in Phase B), (c) summer 2019 (first cell in Phase A), and (d) summer 2019 (last cell in Phase B).
Remotesensing 15 05632 g012
Figure 13. Contrast of water parameter data (WT-°C; OIW-mg/L; COND-mS/cm; and DO-mg/L) with NDSI along the cells of the CW tracks for (a) winter 2017 (first cell in Phase A), (b) winter 2017 (last cell in Phase B), (c) summer 2017 (first cell in Phase A), and (d) summer 2017 (last cell in Phase B).
Figure 13. Contrast of water parameter data (WT-°C; OIW-mg/L; COND-mS/cm; and DO-mg/L) with NDSI along the cells of the CW tracks for (a) winter 2017 (first cell in Phase A), (b) winter 2017 (last cell in Phase B), (c) summer 2017 (first cell in Phase A), and (d) summer 2017 (last cell in Phase B).
Remotesensing 15 05632 g013
Figure 14. Comparison of water parameter data (WT-°C; OIW-mg/L; COND-mS/cm; and DO-mg/L) with NDSI along the cells of the CW tracks for (a) winter 2019 (first cell in Phase A), (b) winter 2019 (last cell in Phase B), (c) summer 2019 (first cell in Phase A), and (d) summer 2019 (last cell in Phase B).
Figure 14. Comparison of water parameter data (WT-°C; OIW-mg/L; COND-mS/cm; and DO-mg/L) with NDSI along the cells of the CW tracks for (a) winter 2019 (first cell in Phase A), (b) winter 2019 (last cell in Phase B), (c) summer 2019 (first cell in Phase A), and (d) summer 2019 (last cell in Phase B).
Remotesensing 15 05632 g014
Figure 15. Contrast of water parameter data (WT-°C; OIW-mg/L; COND-mS/cm; and DO-mg/L) with MSAVI along the cells of the CW tracks for (a) winter 2017 (first cell in Phase A), (b) winter 2017 (last cell in Phase B), (c) summer 2017 (first cell in Phase A), and (d) summer 2017 (last cell in Phase B).
Figure 15. Contrast of water parameter data (WT-°C; OIW-mg/L; COND-mS/cm; and DO-mg/L) with MSAVI along the cells of the CW tracks for (a) winter 2017 (first cell in Phase A), (b) winter 2017 (last cell in Phase B), (c) summer 2017 (first cell in Phase A), and (d) summer 2017 (last cell in Phase B).
Remotesensing 15 05632 g015
Figure 16. Comparison of water parameter data (WT-°C; OIW-mg/L; COND-Ms/cm; and DO-mg/L) with MSAVI along the cells of the CW tracks for (a) winter 2019 (first cell in Phase A), (b) winter 2019 (last cell in Phase B), (c) summer 2019 (first cell in Phase A), and (d) summer 2019 (last cell in Phase B).
Figure 16. Comparison of water parameter data (WT-°C; OIW-mg/L; COND-Ms/cm; and DO-mg/L) with MSAVI along the cells of the CW tracks for (a) winter 2019 (first cell in Phase A), (b) winter 2019 (last cell in Phase B), (c) summer 2019 (first cell in Phase A), and (d) summer 2019 (last cell in Phase B).
Remotesensing 15 05632 g016
Figure 17. A comparison of average R2 values for the initial cells in both phases against water parameters for the years 2017 and 2019 during winter and summer.
Figure 17. A comparison of average R2 values for the initial cells in both phases against water parameters for the years 2017 and 2019 during winter and summer.
Remotesensing 15 05632 g017
Figure 18. Illustration of the average R2 comparison for the initial cells in both phases against water parameters, COND-Ms/cm (dashed line) in comparison to NDVI along the cells of the CW tracks, specifically for winter 2017; last cell in Phase B (dotted line) during the years 2017 and 2019.
Figure 18. Illustration of the average R2 comparison for the initial cells in both phases against water parameters, COND-Ms/cm (dashed line) in comparison to NDVI along the cells of the CW tracks, specifically for winter 2017; last cell in Phase B (dotted line) during the years 2017 and 2019.
Remotesensing 15 05632 g018
Table 1. Water quality data for samples collected at the inflow and outflow of the CW system in Oman [27,31].
Table 1. Water quality data for samples collected at the inflow and outflow of the CW system in Oman [27,31].
Parameter Inflow Inflow (mg/L)Outflow (mg/L)
Total dissolved solids700012,000
Suspended solids2810
Oil in water280<0.5
BOD15.7<1
Total nitrogen2.5<0.5
Ammonia nitrogen1.3<0.5
Total phosphorus0.03<0.5
Boron4.55.6
Table 2. Details of the acquired satellite images (* https://www.satimagingcorp.com/satellite-sensors/other-satellite-sensors/sentinel-2a/ (accessed on 6 October 2022)).
Table 2. Details of the acquired satellite images (* https://www.satimagingcorp.com/satellite-sensors/other-satellite-sensors/sentinel-2a/ (accessed on 6 October 2022)).
Detection Image TimeSentinel-2A and 2B Specifications
Date (M/D/Y)SatellitesBands* Wavelength (µm)* Resolution (m)
02/10/2017Sentinel-2A1—Coastal Aerosol
2—Blue
3—Green
4—RED
5—Vegetation Red Edge
6—Vegetation Red Edge
7—Vegetation Red Edge
8—NIR
8A—Vegetation Red Edge
9—Water Vapor
10—SWIR—Cirrus
11—SWIR
12—SWIR
0.443
0.49
0.56
0.665
0.705
0.74
0.783
0.842
0.865
0.945
1.375
1.61
2.19
60
10
10
10
20
20
20
10
20
60
60
20
20
04/21/2017Sentinel-2A
06/30/2017Sentinel-2A
08/29/2017Sentinel-2A
10/28/2017Sentinel-2A
12/27/2017Sentinel-2A
02/01/2019Sentinel-2B
04/16/2019Sentinel-2B
06/20/2019Sentinel-2A
08/19/2019Sentinel-2A
10/28/2019Sentinel-2A
01/01/2020Sentinel-2A
Table 3. R2 of 1st cells in Phase A verses water collection parameter data of both 2017 and 2019 in winter and summer seasons.
Table 3. R2 of 1st cells in Phase A verses water collection parameter data of both 2017 and 2019 in winter and summer seasons.
Year/IndicesSeasonWTOIWCONDDO
2017 NDVIWinter0.45430.87540.28120.2572
Summer0.40560.72170.12720.3786
2019 NDVIWinter0.91540.94090.63650.9189
Summer0.70420.66980.13940.2284
2017 NDSIWinter0.29070.00150.45760.4775
Summer0.02180.03730.19210.0293
2019 NDSIWinter0.00460.01160.09990.0505
Summer0.28120.027340.00340.0277
2017 MSAVIWinter0.07810.6080.01110.0114
Summer0.28730.62680.04980.2656
2019 MSAVIWinter0.930.94770.64110.9581
Summer0.58310.5670.04890.1822
Table 4. R2 of 1st cells in Phase B verses water collection parameter data of both 2017 and 2019 in winter and summer seasons.
Table 4. R2 of 1st cells in Phase B verses water collection parameter data of both 2017 and 2019 in winter and summer seasons.
Year/IndicesSeasonWTOIWCONDDO
2017 NDVIWinter0.8860.35760.95710.6942
Summer0.4740.1660.56520.1975
2019 NDVIWinter0.7550.94950.49020.6842
Summer0.01230.02440.26670.0382
2017 NDSIWinter0.92960.48950.95920.8188
Summer0.7180.47910.87150.4823
2019 NDSIWinter0.01130.42140.70470.226
Summer0.14120.13180.76610.3935
2017 MSAVIWinter0.82410.2730.97740.6258
Summer0.46680.16620.58660.176
2019 MSAVIWinter0.1510.52810.5680.0759
Summer0.04650.05610.2180.0091
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Al-Jabri, K.; Al-Mulla, Y.; Melgani, F.; Stefanakis, A. Remote Sensing Analysis for Vegetation Assessment of a Large-Scale Constructed Wetland Treating Produced Water Polluted with Oil Hydrocarbons. Remote Sens. 2023, 15, 5632. https://doi.org/10.3390/rs15245632

AMA Style

Al-Jabri K, Al-Mulla Y, Melgani F, Stefanakis A. Remote Sensing Analysis for Vegetation Assessment of a Large-Scale Constructed Wetland Treating Produced Water Polluted with Oil Hydrocarbons. Remote Sensing. 2023; 15(24):5632. https://doi.org/10.3390/rs15245632

Chicago/Turabian Style

Al-Jabri, Khaled, Yaseen Al-Mulla, Farid Melgani, and Alexandros Stefanakis. 2023. "Remote Sensing Analysis for Vegetation Assessment of a Large-Scale Constructed Wetland Treating Produced Water Polluted with Oil Hydrocarbons" Remote Sensing 15, no. 24: 5632. https://doi.org/10.3390/rs15245632

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop