Influence of land cover indices and surface temperature on the metals bioaccumulation by three Macrophytes in Lake Burullus, Egypt

Nowadays, the importance of using macrophytes in accumulation of heavy metals has gained great concerns. So, this study aimed at extracting the land use/cover types of three indices and surface temperatures in the habitats inside 100 m buffers from recent satellite images around three highly economic macrophytes namely; Phragmites australis, Typha domingensis and Potamogeton pectinatus species. In addition to land surface temperature (LST), three important indices expressing the land cover of habitats namely; normalized different vegetation index (NDVI), normalized different water index (NDWI), and normalized different moisture index (NDMI) were extracted to find out there influence on the efficiency of macrophytes in the accumulation of these metal ions; Fe, Cu, Zn, Cd and Pb. The Polynomial regression models were calculated to predict the accumulation factors of plants within the remotely sensed indices and LST. Results showed different accumulation values for individual or more metals in the below-ground and above-ground parts of macrophytes within different habitats. This study considers as an innovative approach using remote sensing technique and satellite images for the selecting of species that can accumulate more metals within different habitats. The obtained results will be useful for the optimal management of these macrophytes in Lake Burullus, a Ramsar site.


Introduction
Heavy metals in the environment have become a global issue due to the increase of human impact (Dar et al. 2020;Nabi 2021). They are severe contaminants in environments due to their persistence, toxicity, and bioaccumulation (Nabi and Dar 2022). Wetlands and coastal waters are ecosystems that are particularly vulnerable to heavy metal inputs (Mitsch and Gosselink 2007;Halpern et al. 2008). Natural processes typically do not remove heavy metals from these ecosystems (Bargagli 1998). As soon as heavy metals get accumulate in bottom sediments, they begin to move up the food chain, often biomagnifying at higher trophic levels and ultimately causing potential disorders in humans and animals (Barwick and Maher 2003;Roberts et al. 2008). Recently, there has been awareness towards using biological indicators such as plants for monitoring and quantifying different pollution types (air, water, and soil) (Al-Yemni et al. 2011). Plants have the ability to absorb all metals, especially those essential for their growth and development (Kabata-Pendias 2011). Macrophytes, in particular, play a fundamental role in wetland geochemistry because they are the principal living accumulators of heavy metals through active and passive absorption (Vodyanitskii and Shoba 2015;Dar et al. 2022a). They show significant variety in their accumulative capacities to heavy metals and transfer them to above-ground organs (Baldantoni et al. 2009).
Aquatic macrophytes are aquatic vascular plants that are broadly distributed in several wet environments, from freshwater to saltwater. The shallow areas of lakes, ponds, pools, marshes, streams, and rivers are where aquatic macrophytes are most frequently found (Dar et al. 2022b). They may be emergent, submerged or floating, rooted or unrooted in habit, with associated adaptations to the leaves, stems and/or roots matching the requirements of these aquatic environments (Bornete and Puijalon 2011;Peters and Lodge 2009;Rejmankova 2011). The spread of macrophytes affects the water circulation and subsequently may affect the quality of the water and fish in the lake. These plants are considered a huge source of raw material for industrial production of paper pulp, biofuel and natural therapeutics (Dar et al. 2021a, b). Besides, aquatic macrophytes have shown high efficiency to remove pollutants and recover nutrients from a wide variety of domestic, industrial and agricultural effluents, which validates their role in the bio-remediation of polluted water (Serag 1996;Haroon 2022).
Land surface temperature (LST) represents an important factor in global climatic change and is used in the applications of various fields such as meteorology, climatology and hydrology. Also, the change of land use/cover (LU/LC) has been long established to have an impact on the climate through variable features that modulate precipitation and LST. LULC, on the other hand, produces an impression that may affect the earth's energy balance, thereby altering the region's climate. One of the most common modifiers is LST (Solanky et al. 2018;El Garouani et al. 2021;Ovalle et al. 2021).
This study aimed to investigate the effect of land surface temperature and different land cover types in the surrounding plant habitats on the accumulation efficiency of metal ions in Lake Burullus, a Ramsar site.

Study area
Lake Burullus extends along the Egyptian Deltaic Coast of Mediterranean sea from latitude 31˚ 15′ N to 31˚ 40′ N and from longitude 30˚ 20′ to 31 10′E (Fig. 1). It is located between Damietta and Rosetta branches of River Nile, and considers one of the largest natural lakes in Egypt after Lake Manzala, covering an area of 420 km 2 . It has a rectangle shape covers a distance of 47 km along the NE-SW axis with a width ranging between 4 and 14 km (Okbah 2005;Shaltout and Khalil 2005). The northeastern edge of the Lake has a short canal called Boughaz El-Burullus that connects it to the Mediterranean Sea, through which sea water comes into the lake easily during the periods of low Nile water inflows. The lake is known for its many islands, and the majority of these islands run from south to north (Abd el-Sadek et al. 2022). Other islands are either parallel or perpendicular to the current shore. Lake Burullus is home to about seventy five islands with different habitat types (Balah 2012). Many floating vegetation types, such as reeds and submerged in the lake, altering the water circulation. These plants serve as a vital role in preventing the collapse of the lake's interior shores (Abd el-Sadek, et al., 2022). The area of Lake Burullus has changed as a result of shrinking that led to the increase in islands size (Khedr 1999;Balah 2012).
The northern Mediterranean part of the Nile Delta belongs to the arid zone as shown by the world distribution map of arid regions. As a result, Lake Burullus is characterized by an arid climate with temperatures range between (20-30° C) and (10-20° C) in warm summers and mild winters, respectively (UNESCO 1977;UKMO, 2013).
Lake Burullus receives drainage water from seven drains along its southern edge and freshwater from the Brinbal Canal in its southwest corner (Okbah and Hussein 2006). The amount of drainage waters coming from agricultural lands to Lake Burullus equal to nearly 4 billion m 3 /year (El Shinnawy, 2002), which accounts for 97% of the water inflow (Shaltout and Khalil 2005;). Lake Burullus is affected by agricultural drainage water assorted with different types of drained waters from fish ponds, and industrial and municipal wastewater discharges through the drains.

Ecological and botanical description of macrophytes
Aquatic plants including Phragmites australis, Typha domingensis and Potamogeton pectinatus; are macrophytes that are distributed all over the world. These species are rhizomatous perennial herbaceous plants that grow in condensed mono-specific stands in natural lakes with stagnant, shallow water and, sediments of muddy nature (Pignatti 1982). These plants have been employed to detect, monitor, and remediate water contamination (Wolverton and McDonald 1978;Peng et al. 2008;Bonanno and Giudice 2010;).
i Phragmites australis: It is believed to be one of the most widely distributed species in the world (Holm et al. 1977). In Egypt, Phragmites australis occurs in all phytogeographical areas (Täckholm 1974;Zahran and Willis 2009;Boulos 2005). It has found in the main habitats of the Lake Burullus area including: salt marshes, sand sheets, lands that have been cut off from the lake, terraces, slopes, and water edges as well as open water zones of drains, lake shores and the open water of the lake (Shaltout and Al-Sodany 2008). It is an emergent plant and one of the important macrophytes. It is a perennial reed that grows from elongated rhizomes or stolons. It is 1-6 m tall, and forms dense stands. Stems are erect, hollow, reed-like, simple, 150-600 cm long, 5-15 mm thick and hollow internodes. Leaves are linear, flat, drooping, and leaf blades are deciduous at the ligule; 20-60 cm long; 8-32 mm wide, with pointed tips. Flowers happen in August and September and form bushy panicles (Clayton et al. 2006;Klein 2011). Phragmites australis forms a dense network of roots and rhizomes that can go down up to two meters in depth to reach deep groundwater (MA DCR 2002). It grows in marshes and swamps, along streams, lakes, ponds, ditches, and wet wastelands, often weedy and very difficult to eradicate. It grows best in firm mineral clays and tolerates moderate salinity, where the water level varies from 15 cm below the soil surface to 15 cm above. Ranging from cool temperate steppe to wet through the tropical desert to moist forest life zones, the reed is reported to tolerate annual precipitation of 3.1 to 24.1 dm, the annual temperature of 6.6 to 26.6 ºC and pH of 4.8 to 8.2 (Duke 1978(Duke , 1979. Besides, the reed tolerates soil conductivity up to 12 mS cm-1 and pH 7.0 to 9.3 (Serag, 1996). Phragmites australis has been used as a bio-indicator for heavy metals, and to store heavy metals to some extent (Bonanno 2011;Salem et al. 2014;Morari et al. 2015). ii Typha domingensis is an emergent plant native to warm temperate and tropical climates that grows in ditches and marshy areas all across Egypt. (Täckholm 1974;Boulos 2005). ). It is an erect, perennial, freshwater macrophyte that can grow 3 or more meters in height. The linear cattail leaves are thick, ribbon-like structures with a spongy cross-section exhibiting air channels. The subterranean stem arises from thick creeping rhizomes (Smith 1962(Smith , 1967. Typha domingensis is one of the main components of vegetation that stands along the shores of Lake Burullus close to the Deltaic Mediterranean coast (Shaltout and Al-Sodany 2008). Typha domingensis is used in constructed wetlands for the enhancement of water quality (Abdel-Ghani et al. 2009) due to its high growth rate and great capacity for heavy metal accumulation (Newman et al. 1996;Lorenzen et al. 2001). iii Potamogeton pectinatus is a submerged perennial aquatic macrophyte (Boulos 2005) with a parvopotamid growth form (Hogeweg and Brenkert 1969). It is characterized by slender round shoots up to 3 m long with narrow linear leaves (Kłosowski and Kłosowski 2007). It also occurs in almost all climatic areas and has a widespread distribution (Pilon et al. 2002), and occurs in a variety of habitats including water of different trophic levels, standing and running water, alkaline, fresh and brackish waters. (van Wijk 1988). It can survive in environments with high salinity and pollution (Casagranda and Boudouresque 2007).

Sampling protocol
Fifteen sampling sites were selected along Lake Burullus. Samples considered were water, sediments and plant. The samples of water, sediments, and macrophytes were collected from the same sites. The sampling program was conducted during the summer season of the years 2020 and 2021, respectively.

Water sampling
Water samples were collected from various locations along the shoreline and islets of Burullus Lake (Table 1; Fig. 1). The collected water samples were kept in an ice box, then, it was transferred to the laboratory for heavy metals determination.

Sediment sampling
Surface sediment samples were collected using a Van-Veen grab coated with polyethylene (Amini Ranjbar 1998), and analyzed for Fe, Cu, Zn, Cd, and Pb. The samples were kept in plastic bags and transported to the laboratory, air-dried at room temperature and stored in plastic bags until analysis.

Plant sampling
At each sampling point, 4-6 samples of Phragmites australis, Typha domingensis, and Potamogeton pectinatus were collected from the lake shore and islets within a 5 m x 2 m plot. To remove sediments, roots and rhizomes were washed in the lake water and kept in plastic bags for transferring to the laboratory. from fishfarms that are distributed at the northern side of the Lake, c a canal which is located at the northern side to the Lake, d East El-Burullus area, e and f West El-Burullus Drain, g Phragmites australis species which is distributed nearby the connected area between El-Kashaa drain and the Lake, h and i Typha domingensis and Eicchornia crassipes species which were distributed along the shoreline and drain

Water analysis
For the determination of heavy metals in water samples, the EPA digestion method was used according to Gregg (1989). A 100 ml of the representative water samples was put into Pyrex beakers containing 10 ml of concentrated HNO 3 . The samples were slowly heated and then evaporated on a hot plate to the lowest possible volume (about 20 ml). The beakers were allowed to cool and another 5 ml of Conc. HNO 3 was added. The heating was continued with the addition of Conc. HNO 3 as necessary until digestion was completed. The samples were evaporated again to dryness (but not baked) and the beakers were cooled, followed by the addition of 5 ml of HCl solution (1:1 v/v). After warming the solutions, 5 ml of 5 M NaOH was added and then were filtered. The filtrates were transferred to 100 ml volumetric flasks and diluted to the mark with distilled water. These solutions were then used for the elemental analysis. A total of five metallic elements namely; Fe, Cu, Zn, Cd, and Pb were determined in the pre-treated samples of water using Atomic Absorption Spectrophotometry as described by Gregg (1989).

Sediment analysis
The concentration of heavy metals in sediments was determined according to EPA-ROC (1994). Where, conventional aqua regia digestion was prepared in glass beakers with volumes of 250 ml covered with watch glasses. A 0.5 g of sample was digested in 12 ml of aqua regia on a hotplate for three hours at 108 °C. After the evaporation process near drying, the samples were diluted with 20 ml HNO 3 of 2% (v/v with H 2 O), then it was transferred into a 100 ml volumetric flask after filtering through Whatman No. 42 paper for dilution to 100 ml with double distilled water (DDW). Heavy metals (Fe, Cu, Zn, Cd and Pb) were analyzed using the atomic absorption spectrophotometer, and the results were expressed as microgram per gram (µg/g). Accuracy and precision were checked by using reference material (SD-M-2/IM).

Plant analysis
Plant samples were first dissected and divided into two parts above-ground and below-ground parts of plants; the first part (A) includes leaves and stems, and the second part (B) includes roots and rhizomes; to identify the various bioaccumulation capacities in the above-ground and belowground parts. Plant organs were cut off using stainless steel scissors, and kept at 2 °C for further analysis. The sampled plants were rinsed thoroughly with distilled water, separated as mentioned before, dried at 70ºC for 72 h, crushed and digested using an acid mixture of concentrated H 2 SO 4 and HClO 4 (Grimshaw 1987). The studied heavy metals (Fe, Cu, Zn, Cd and Pb) were determined using the flame atomic absorption spectrophotometer (FAAS, GBC-932), whereas data were expressed as µg/g.

Statistical analyses
Linear correlations between measured metals in the aboveground and below-ground parts of plants, water and sediments were tested through Pearson's r coefficient. The analysis of regression was also conducted. All statistical calculation and Polynomial regression analyses were performed using the software of PAST program. After the analyses, the values were determined for bioconcentration or bio-accumulation and translocation factors to assess element mobility in the study species. The values obtained were based on the following: Where, C plant parts and C sediment were the concentrations of metal ions in plant parts and medium (water and/or sediments) in µg/l and µg/g, respectively. BAF expresses the efficiency of plant species to accumulate metal ions from the surrounded medium. Higher BAF values imply a greater bioaccumulation capability (EPA 2007).

Downloading satellite images and preprocessing
Downloading of a Landsat image from this site; https:// earth explo rer. usgs. gov/ was done. The acquisition date is at 31/07/2021 close to the time of sampling collections. The radiometric corrections were occurred to convert it from digital numbers into reflectance using QGIS 3.16 program.

Analysis of Land Surface Temperature (LST) from Landsat images
Scientific theory of obtaining LST To obtain the LST, different steps should be considered;

I) Conversion of digital number into radiance
The use of Band 10 to recover the LST using ArcGIS 10.5 for digitizing the spectral radiance of B10 using the top of atmosphere (TOA) according to Barsi et al. (2014): Where; AL is the band-specific additive rescaling factor; ML is the band-specific multiplicative rescaling factor; Q cal is the band 10 image; O i is the band 10 correction.

II) Conversion to Maximum brightness temperature
The data of used band could be transformed into brightness temperature (BT) using the metadata file's thermal constant according to this equation; Where: K1, K2 = Bands Specific thermal conversion from the metadata; TB: Temperature of satellite brightness (Celsius); Emissivity correction is crucial to decreasing these inaccuracies and it was occurred to lastly obtain the LST from BT.

III) Land surface emissivity of NDVI (LSE)
The emissivity of the ground surface can be calculated using three equations based on the effectiveness of transporting thermal energy through the surface to the atmosphere. LSE is a proportionality factor that scales blackbody radiance (Planck's law) to predict emitted radiance according to the following equations; LSE is a proportionality coefficient that adjusts the radiance of the black body (Planck's law) to predict the emitted radiation.
Where ε v : the vegetation emissivity, ε s is the soil emissivity, P v is the proportion of vegetation, F is a shape factor whose mean value is equal to 0.55, after assuming several geometrical distributions.

IV) Calculation of LST°C
Once, obtaining the emissivity images, the LST can be calculated using the following equation: Where: BT is the brightness temperature in Celsius ( • C ). LST is expressed in Celsius ( • C ). λ (11.5 μm) is the wavelength of emitted radiance: = h * c∕ = 1.438 * 10 −2 mK, σ is the constant of Stefan-Boltzmann, ε is the land surface emissivity, c is the velocity of light, and h is Planck's constant, and (LSE) as described by Avdan and Jovanovska (2016). In the research, the methods are based on the online methodology for an image of summer season at 15 July 2021, respectively according to produced models of Parastatidis et al. (2017).

Calculations of remotely sensing indices
The remotely sensed indices which were represented the land cover were calculated as illustrated in Table 2.

Results
The extracted values of LST s and land cover types' indices (NDVI, NDWI and NDMI) in a buffer of 100 m around the sampling plants were illustrated in Table 3 and Fig. 3. The average concentrations of heavy metals in both habitats of lake shores and islets of Lake Burullus were as follow; for the water of lake shores were; Fe > Pb > Zn > Cu > Cd; for the water around islets were; Fe > Zn > Pb > Cu > Cd. These orders in sediments take the following sequences; Fe > Zn > Cu > Pb > Cd in both the habitats. The average concentrations of metal ions in the above-ground parts of Phragmites australis were Fe > Zn > Pb > Cu > Cd in both. In the below-ground parts of Phragmites australis, this order was; Fe > Zn > Cu > Pb > Cd from the lake shore and Fe > Zn > Pb > Cu > Cd in the islets habitat. The concentrations of metal ions in the above parts of Typha domingensis in the lake shore were taken this sequence; Fe > Cu > Zn > Pb > Cd, in islets, this sequence was; Fe > Zn > Pb > Cu > Cd. Also these sequences were; Fe > Zn > Cu > Pb > Cd and Fe > Zn > Cu > Pb > Cd in the below-ground parts of Typha domingensis in both the habitats, respectively (Table 4). The raw concentrations of metal ions in water, sediments and plant species were illustrated in Appendix I. The orders of average BAF s in the studied area were as follow; Cd > Fe > Zn > Pb > Cu; Cd > Cu > Zn > Fe > Pb and Fe > Zn > Pb > Cu > Cd for Typha domingensis, Phragmites australis and Potamogeton pectinatus species, respectively as illustrated in Appendix II. The  (Leprieur et al. 2000) Normalized difference moisture index NDMI = NIR−SWIR NIR+SWIR (Gao 1995) Fig. 4, it is obvious that the intake of Fe, Zn and Pb by Potamogeton pectinatus may be associated with the values of temperature as there is high significant correlation (R 2 = 0.98, 0.91 and 0.90), respectively. While, there was a highly negative significant correlation between low sparse vegetation (low vegetation cover) and the intake or accumulation of Pb. Other species of Phragmites australis and Typha domingensis showed low alterations with extracted LST s values. There is also positive correlation with NDWI.
Results showed that some metals could be accumulated more in high vegetative habitats as Cd and Pb in the above parts of Typha domingensis. Whereas, other metals as Cu and Pb were accumulated in below parts of Phragmites australis in moderate vegetative areas. It was obtained that there were high positive significant correlation between the accumulation of Cu and Pb in below and above parts of Phragmites australis in moderate vegetative areas and mid-low canopy cover. In the above-ground parts of Typha domingensis species there were significant positive correlation between accumulation factor of Fe in moderate vegetation habitat and negatively with average canopy cover. While the accumulation of Cd and Pb increases in high vegetative locations. For below parts of Typha domingensis, results give indication that Zn accumulation decreases within low or sparse vegetation.
The regression analysis role was applied to obtain unknown information based on field studies (Austin 1971). The determination coefficient R 2 indicated that the most fitted regression models is for Fe with average canopy cover ACC and mid-high canopy cover MHCC habitats (R 2 = 0.94 and 0.92, respectively). For Pb, the regression analysis R 2 equal 0.77. In the above parts of Phragmites australis, the fitted polynomial regression was for Pb in conditions of high vegetative and mid to low canopy cover with R 2 equal to 0.77 and 0.93, respectively. Other models showed different significant ranged between low to moderate as illustrated in Fig. 5, 6, and 7; Table 5. Appendix III explained survey data of different vegetation of Burullus ecosystem.

Discussion
In the aquatic environments, the metal ions have been considered as a result of their toxicity, tendency to bio-concentrate and persistence. Nowadays, aquatic macrophytes were being used as functional intent for phytoremediation purposes. Recently, removal of metal ions from aqueous surface waters using accumulating roots and rhizomes (Pillai 2010;Xing et al. 2013).
The bioaccumulation factors are important to understand the availability of trace metals to plant species (Cheng 2003). Minerals content and compositions of plants differed significantly based on various species types (Kibar and Temel 2015). The ability of macrophytes to translocate and accumulate metal ions differs according to habitat species, tissues, prevailing climatic conditions, redox potential and pH (Eid et al. 2021). It was observed that the most accumulated element was Cd in Typha domingensis and Phragmites australis, while the most accumulated element ion in Potamogeton  Majid et al. (2014). Aquatic ecosystems were subjected to numerous stress factors; one of these is the increase in temperature due to climatic changes and metal disposal. So, thermal stress can magnify the impacts of metal ions on the aquatic macrophyes (Nin and Rodgher 2021). The LST is a significant indicator for the ecological and environmental changes in coastal wetlands and showed substantial spatiotemporal changes under severe sealand interactions and different anthropogenic activities (Chi et al. 2020). The increase in sediment temperatures made influence on the vegetation mechanisms, metabolisms and the sediment characteristics itself, so it may enhance the sediment-plant translocations of metallic ions. So, the global temperature may effect on these pollutants bio-accumulation (Cornu et al. 2016;Lee and Kim 2022). The soil temperatures affected the capacity of different plant to accumulate metal ions such as; Cd, Zn, Cu and Pb. Principle component analysis proved that temperature, physiological and photosynthetic factors play role in the metal translocation properties of plant-soil (Yu et al. 2013).
The differences in LST s of LULC types within different locations around or inside the lake area were clear. It may be interpreted as, saline basins and water areas caused lower LST, built-up areas and roads caused high LST, vegetative lands and islets caused medium LST, and barren areas possessed high LST (Chi et al. 2020). The correlation of BAF of Potamogeton pectinatus within temperature indicated the efficiency to accumulate metal ions. It is efficient to accumulate metals discharged in rainy waters at high temperature. Potamogeton pectinatus can accumulate more Zn (Fritioff et al. 2005).
Temperature showed low to no significance with BAF in Phragmites australis parts. There is low negative Fig. 4 A-E The relation between BAF of metal ions in A Potamogeton pectinatus, B below-ground parts of Phragmites australis, C above-ground parts of Phragmites australis, D above-ground parts of Typha domingensis and E below-ground parts of Typha domingensis with land surface temperature and cover types' indices ▸ significance correlation between above and below-ground parts of Typha domingensis and surface temperature for Cd. Typha domingensis has the ability to overcome the Cd toxicity and characterized by its potentiality of phytoremediation.
Vegetation Indices (VI s ) were insensitive indicators for monitoring the effects of metal in vegetation. As the spectrum alterations of leaf within different seasons may be caused by metal pollution (Zhou et al. 2018). The normalized difference moisture index (NDMI) explains the Fig. 5 Polynomial regression analysis between metal accumulation in the above-ground parts of Typha domingensis with LST and LC type in a buffer 100 m at Lake Burullus vegetation content of water. It is suggested for monitoring moisture of vegetation using remote sensing data (Gao 1996). The vegetation water content considers one of the vital biophysical characteristics of the healthy vegetation. The NDWI can aid in the evaluation of dryness stress on the aquatic vegetation of Mediterranean type ecosystems through plant available water (Serrano et al. 2019). The moisture content was also known as a function of plant sample's water content (Makarius et al. 2013). It could be used for monitoring the water stress in vegetation . It is obvious that canopy water content is a comprehensive indicator reflecting the vigor and health of vegetation growth. Vegetation water content is one of the significant biophysical of vegetation health features, and its remote valuation can be exploited to real-timely monitor vegetation water stress. It can be used to different factors as Fig. 6 Polynomial regression analysis between metal accumulation in the below-ground parts of Typha domingensis with LST and LC type in a buffer 100 m at Lake Burullus canopy water content for further water treatments (Zhang and Zhou 2019). It was may be interpreted that it had the ability to adsorb more metals and efficiency of water stress treatment increase. This was especially for Fe in the aboveground parts of Typha domingensis and Zn in the aboveground parts of Phragmites australis. Also, the water content of the plant shoots in below-ground parts in water habitats of low minerals content than those of high nutrients (Drew 1967). The vigor NDVI not only represented by the dense of habitat; as some metals may cause change in NDVI. For example, high concentrations of Zn in plant species may induce significant decreases in NDVI (Chen et al. 2009). It was clear from the correlations in this study.
It is clear that the abundance of watered habitats was vital for the submerged vegetation as Potamogeton pectinatus. It is obvious that the accumulation of Fe, Zn and Pb increases in water conditions, whereas the accumulation of Cu and Cd can occur with low water conditions. While NDWI represented low to no significant correlation with other two species.
Using multiple regressions in prediction process is a tool used to forecast dependent factors from a group of Fig. 7 Polynomial regression analysis between metal accumulation in the above and below ground parts of Phragmites australis with LST and LC type in a buffer 100 m at Lake Burullus independent factors (Eid et al. 2010a, b). Also regression models were suitable to predict metal accumulation in plants and more useful when being compared with numerous validation tools (Kumar et al. 2020). The coefficient factor (R 2 ) is important in the prediction models. Abd El-Hamid et al. (2022) stated that R 2 and root mean squares errors RMSE were fit to evaluate the accuracy of these models. It is observed from results that the most fitted models are for Fe with ACC and MHCC; and Pb with MLCC in the above-ground parts of Typha domingensis and Phragmites australis, respectively.

Conclusion
The study involved land use/cover of fifteen habitats within and around Lake Burullus. The integration of different remote sensing indices was aided more in identifying the differences in plant efficiencies of metal accumulation. Each plant species can accumulate different metals within different LST, NDVI, NDWI and NDMI values. Polynomial regression models were aided in the prediction of metal accumulators within different habitats in the case of integration of remotely sensed indices. Extracted surface temperatures showed more influences with the submerged species as Potamogeton pectinatus more than those of Phragmites australis and Typha domingensis. NDVI values may be interpreted in different two trends; one is the vegetation health and or the dense presence and distribution in the ecosystem. Potamogeton showed high accumulation within water conditions. The two species of Phragmites australis and Typha domingensis has the ability to accumulate metal ions in different habitats more than Potamogeton. This is an indication to the availability to grow and endure in polluted habitats. So these macrophytes species were highly recommended to be used for the remediation of polluted waters by metals in similar habitats' conditions.