1-Introduction and Background
The geography of ancient Egypt and its mythology are closely connected, with the annual flooding of the Nile shaping both the physical landscape and the cultural worldview. In ancient Egyptian cosmology, the world emerged from a primordial watery void (Nun) as the first land, or "primeval mound". This study focuses on Esna, Upper Egypt, home to the Temple of Khnum and several other unexcavated or demolished temple sites. Inscriptions from the temple, such as the Festival of Seizing the Staff, metaphorically describe the local terrain as marshlands and document processional routes. While traditional Egyptological methods have correlated some text-based place names with physical locations, they often lack quantitative spatial analysis, leaving the precise geomorphological context of these narratives largely undefined. Macro-topography remains stable for a long time, barring massive modern anthropogenic intervention. This research adopts a landscape-first approach. By leveraging open-source geographic information systems (GIS) and remote sensing, we mathematically translate qualitative ancient texts into a DEM-driven topographic model to reconstruct the historical landscape that inspired the legends.

2-Study Area and Open Data Sources
The study area is Esna, with five documented temple sites surrounding it: Temple of Khnum, House of God, the Temple of Isis, and Kom Mir. The construction of the Aswan High Dam has masked historical paleochannels and floodplains; insufficient data are available for hydrological reconstruction. Consequently, this study heavily prioritizes Digital Elevation Models (DEMs). We lean on declassified 1960s CORONA panoramic stereo pairs to extract a pristine bare-earth microtopography, capturing the landscape long before recent agricultural expansion and infrastructure projects leveled it. We then complement this historical topographic baseline with decades of multispectral data from Sentinel-2 and the Landsat program (Landsat 1–9) to monitor remaining vegetation patterns and water indices. Aligning with the core philosophy of FOSS4G, our primary computational platform (QGIS) and all incorporated remote sensing datasets are entirely open-access.

3-Methodology
To bridge the gap between mythological narratives and spatial reality, we developed a reproducible, DEM-centric workflow using QGIS and Python spatial libraries. The best method we can use is layer analysis, which can turn qualitative contents into different GIS layers. The layers can be easily linked and operated with different calculations. The methodology includes three parts:
First, textual descriptions are converted into topologically validated GIS layers. Layer 1 represents the "Primeval Mound" (elevated, unflooded zones), Layer 2 represents "Marshes & Lakes" (low-lying retention basins), and Layer 3 delineates the "Western Mountains" (the absolute safety boundary).
Second, we execute rigorous Topographic Surface Modeling. Based on terrain stability, we use QGIS terrain analysis algorithms to compute critical geomorphological variables from the historical DEM. By calculating slope, aspect, Topographic Position Index (TPI), and the Topographic Wetness Index (TWI), we quantitatively define the physical landscape characteristics, isolating potential ancient mounds from natural depressions.
Third, we perform Hydro-conditioning and Simulated Routing. Since modern hydrology is severely disrupted, we reverse-engineer the ancient floodpaths through the terrain. We apply open-source algorithms to hydro-condition the DEM—executing pit-filling, flow-direction, and flow-accumulation routing—to establish a hydrologically correct surface. A terrain-based inundation algorithm, such as the Height Above Nearest Drainage (HAND) model, is then deployed. By routing simulated water levels across this stable topography, the inundation results are meticulously calibrated to match the spatial descriptions in the ancient texts.

4-Preliminary Results and Archaeological Implications
Currently, the DEM-driven spatial analysis successfully bridges the text-to-terrain gap. Our analysis of Layer 3 (the "Western Mountains") exposes a significant divergence between mythological narratives and physical geography. While the inscriptions describe this western margin as an impassable, absolute safety zone, the DEM analysis reveals it to be a modest ridge with an average elevation of only 200 to 300 meters. Furthermore, extracting elevation profiles from our decadal remote-sensing time series confirms that this topography has remained geomorphologically static, ruling out historical degradation. This discrepancy suggests that the ancient characterization was a phenomenological exaggeration, likely stemming from limited mobility and the imposing visual perspective of looking westward from the low-lying Esna basin. Identifying this spatial hyperbole is archaeologically significant; it demonstrates a broader tendency for cognitive exaggeration within the temple texts, providing a critical, data-driven foundation for reinterpreting other geographic claims in the inscriptions.
Our geomorphological computations explicitly identified elevated landforms and topographical depressions that align flawlessly with the text-derived topological layers. The hydro-conditioned HAND model effectively simulates the historical flood recession zones, demonstrating how rising waters would inundate the natural basins (Layer 2) and expose the structurally sound primeval mounds for temple construction (Layer 1).

5-Conclusion and Future Work
This study establishes a robust methodological blueprint for landscape archaeology. By pivoting from a purely hydrological focus to a terrain-driven analysis, we demonstrate how open-source GIS tools can bypass modern infrastructural disruptions to reconstruct ancient environments. The calibrated, DEM-based inundation model not only contextualizes the past but also serves as a predictive tool.
Future work will focus on:
First, we will broaden the temporal depth of our spatial database by incorporating and georeferencing early historical cartography (e.g., 18th- and 19th-century expedition maps). This process will cross-validate our findings and further corroborate the millennial topographic stability of the Esna geomorphology.
Second, we will expand our textual dataset—currently focused on the Festival of Seizing the Staff—to include inscriptions from other major local events, such as the Festival of Raising the Sky. Cross-referencing these distinct mythological narratives will allow us to generate additional topological layers and conduct rigorous statistical evaluations to assess their spatial validity and accuracy.

References:
[1] Sauneron, S. (1962). Le temple d'Esna. Tome V: Les fêtes religieuses d'Esna aux derniers siècles du paganisme. Le Caire: IFAO.
[2] Abdel-Raham, A. M. (2009). The Lost Temples of Esna. BIFAO, 109, 1-8.
[3] Assmann, J. (1996). The Mind of Egypt: History and Meaning in the Time of the Pharaohs. Metropolitan Books.

1-Introduction and Study Area
Cladophora is a filamentous green alga native to the North American Great Lakes. Its excessive proliferation not only causes foul odors and impairs public beach recreation but also triggers severe ecological issues, including avian botulism outbreaks. Since the 1990s, the filtering effect of invasive species such as dreissenid mussels has significantly increased water clarity, allowing sunlight to penetrate to greater depths. This has led to massive Cladophora blooms even under relatively low nutrient concentrations. The study area of this research focuses on the nearshore waters along the southern shore of Lake Ontario (the United States side). To achieve precise calibration of remote sensing observations, the spatial scope of the study is strictly defined as two independent 6 km × 6 km square regions, centered respectively around two key hydrological and biological monitoring stations established by the United States Geological Survey (USGS): the OIR station (Irondequoit, near Rochester) and the OOL station (Olcott).
These two core USGS stations provide substantial, highly valuable ground-truth data for this study. These comprehensive datasets encompass multi-depth water flow velocities, water turbidity, and various critical chemical constituents in the water column (such as nutrient concentrations). More importantly, the stations provide net weight data of Cladophora samples collected in situ across different depth gradients. These multi-dimensional, high-precision ground truth indicators not only serve as an irreplaceable validation foundation for evaluating and calibrating various spectral remote sensing indices within our open-source computational architecture, but also enable us to deeply investigate the complex mechanisms underlying the relationships between micro-environmental physicochemical variables and nearshore benthic algal outbreaks.

2-Evaluation of Traditional Indices and Experimental Derivation of a Novel Index
In the preliminary remote sensing analysis phase, we developed a Python-based workflow to extract Sentinel-2 image bands and automatically calculated various traditional spectral indices, including NDVI, FAI, NDAVI, and SABI. Statistical analysis of multi-temporal imagery (from May to August 2023) revealed that the mean and median values of these indices were frequently negative or extremely low, accompanied by disproportionately large standard deviations. For instance, across multiple summer observation dates, the median values for NDVI and FAI consistently hovered near zero (ranging from -0.012 to 0.025). At the same time, NDAVI and SABI exhibited even deeper negative medians (often between -0.05 and -0.09). Furthermore, the high standard deviations—frequently exceeding 0.25 for NDVI and 0.50 for SABI—demonstrated massive signal noise. This statistical analysis demonstrates that vegetation indices based on the Near-Infrared (NIR) band exhibit severe absorption failures in aquatic environments, rendering them inadequate for precise mapping of submerged benthic Cladophora.
To address this optical challenge and identify the optimal spectral response, we designed a controlled physical experiment. A 3m × 3m water tank was used, with an incandescent light source simulating solar irradiance. A receiver simulated the satellite sensor to capture reflectance from a green surrogate representing benthic algae. Strikingly, the experimental results revealed that the strongest reflectance signals emerged in the Blue and Short-Wave Infrared (SWIR) bands, significantly diverging from the band selections of traditional vegetation indices. Based on these empirical findings, we are currently conducting rigorous mathematical derivations utilizing the Blue and SWIR bands to formulate a novel, water-penetrating spectral index specifically optimized for Cladophora detection.

3-Automated Open-Source Cloud-Masking Algorithm to Bypass API Limitations
To achieve high-frequency monitoring of Cladophora, we aimed to build a fully open-source, automated data acquisition architecture. However, querying the Copernicus Data Space API inevitably encounters strict request frequency limits and download volume quotas. Furthermore, the official API only provides the average cloud cover percentage at the full-scene level. For our small 6 km × 6 km Region of Interest (ROI), this macroscopic cloud assessment is highly inaccurate. A scene with a low average cloud percentage might still have dense clouds completely obscuring our study area, leading to massive invalid downloads and wasted bandwidth. Additionally, a single remote sensing image rarely covers the target area perfectly without clouds, necessitating the seamless mosaicking of multiple images and stricter screening for high-quality data.
To overcome this core bottleneck, we designed and implemented a regional cloud-masking algorithm based on image Quicklooks (previews) within our workflow. Since Quicklook files are extremely small and consume negligible download bandwidth, the program automatically prioritizes retrieving them. Given that Quicklooks do not inherently contain geographic coordinates, the algorithm first extracts the boundary coordinates of the scene's footprint polygon from the metadata. Subsequently, it correlates and standardizes the ROI's geographic coordinates against this footprint boundary. Based on this geometric translation, the system can precisely reverse-engineer the specific pixel rectangle corresponding to the study area on the unreferenced Quicklook image. Ultimately, the algorithm computes the proportion of white pixels exclusively within this localized bounding box to accurately assess the true cloud cover rate within the ROI. Only when the ROI's cloud cover meets strict clear-sky thresholds does the system automatically trigger the API to download the heavy, high-resolution original imagery. This algorithm successfully achieves precise "on-demand downloading," effectively circumventing API bandwidth restrictions while dramatically improving the efficiency of acquiring the cloud-free data required for subsequent image mosaicking.

4- Conclusion and Future Works
This study successfully established a highly efficient, Python-based open-source remote sensing download architecture that practically circumvents API limitations. It also highlighted the severe shortcomings of traditional vegetation indices through both satellite data statistics and controlled physical experiments. Future research will focus on advancing two primary tasks:
First, further refining the Quicklook-based cloud-masking algorithm to automate the acquisition of extensive multi-temporal imagery for seamless spatial mosaicking. To ensure complete reproducibility, this process will be integrated into an end-to-end Python pipeline, with the full source code made freely available on GitHub.
Second, finalizing the mathematical formulation of our novel Blue-SWIR spectral index based on the water tank experiment, and deploying it within our open-source pipeline to precisely map the spatial distribution and evolutionary dynamics of Cladophora during peak summer blooms.

References:
[1] Howell, E. T. (2018). A decadal-scale perspective on the occurrence of Cladophora on the north shore of Lake Ontario. Environmental Monitoring and Assessment.
[2] Wright, N., et al. (2024). CloudS2Mask: A novel deep learning approach for improved cloud and cloud shadow masking in Sentinel-2 imagery. Remote Sensing of Environment, 306, 114122.
[3] Copernicus Data Space Ecosystem. (2024). Quotas and Limitations Documentation.