Start - Overview

Welcome to the 360 degree pond bathymetry study tour
Along the red line, we explain how pond depth can be determined from aerial photo data using photogrammetry. Every red cone marks one step on your tour. Use the arrows to navigate from pond to pond or select the processing step from the list in the top left corner.
Use this overview to briefly familiarise yourself with the location, the control in the 360 images and to perceive the dimensions based on the people located below you.
R/V Polarstern campaign PS106/PASCAL ice station
You are now standing on an ice floe north of Svalbard in June 2017. Tethered to the floe with the ship, we can both measure airborne and simultaneously collect ground-truth data around and in the ponds.
The study site is located in between a rafted ice zone, where deformation has caused ice to submerge. As a result, numerous ponds of different colours have formed. They are neither visually nor geometrically different from melt ponds and yield a perfect study site for aerial remote sensing of pond depth as shown with a spectral method in König et al. (2020).

The study site from above:

Aerial images
Aerial images were recorded with a consumer grade CANON camera (14mm lens) mounted nadir on a helicopter.
The typical ground resolution of the setup ranges between 5 and 10cm.

The helicopter above the pond site:

Melt ponds
Ponds play a key role in the Arctic energy budget and ecosystem. Vertical geometry is largely underrepresented in studies and documented only by relatively few in situ measurement points. Remote sensing of pond bathymetry raises the possibility to better implement pond depth and its effect on albedo and meltwater distribution into simulation models.
Photogrammetric reconstruction
A parallel flight grid provides overlap between images to realistically reconstruct the surface using photogrammetric triangulation (ideally 80% forward, 60% lateral overlap). The full-colour images are stitched together and projected onto the surface as an orthomosaic. We use Agisoft Metashape for this processing and correct for offsets caused by the continuous ice drift. Ground control points were not available for the model, but a laid out 2x2m pattern was used in the processing to scale the surface correctly (±1%).
Surface classification
A surface type classification of ponds based on a Random Forest Classifier is used to detect ponds and derive pond polygons.
Pond water level
Framing the pond polygon outline in the elevation model yields the pond water level for each pond individually.

The blue area in this pond shows the derived imaginary water surface. It differs partially from the pond outline as can be seen below this info button due to height differences in the elevation model. A statistical analysis of the vertical deviation however shows, that for all ponds in this study except one, half of the raster points along the outline vary by less than ±2cm from the waterlevel. 99.3 % of the points are usually within ±5cm.
Pond bathymetry
Cutting edge photogrammetry tools can not only reconstruct over-water topography, but are by now powerful enough to trace intersection points even at the bottom of clear shallow waters. However, refraction at the water surface is not taken into account in such automated calculations and must be considered manually. We have developed a 2-step correction method as a convenient extension to the largely automated workflow (see next box). The pond bathymetry eventually results from the vertical distance between corrected topography of the pond bottom and previously retrieved pond water level.
From this step onwards, all ponds are displayed with their depth layer, the corresponding colormap is linked further to the right.
Colorbar

Refraction correction I
Photogrammetry relies on colinearity of intersection lines. Refraction at the air-water interface bends the beam lines and causes horizontal and vertical translation of intersection points. We developed a two step correction method to minimize the impact on the pond depth estimation.
I: Images are masked to a maximum incident angle of 40° at the ground to keep horizontal mismatch between intersection lines below ground resolution.
The impact can be expressed by the relation between deviation and pond depth, we call it k. k depends on the incident angle of the different camera positions to the intersection point:

Keeping k below 0.028 means the horizontal deviation becomes max. 2.8cm for 1m deep ponds and remains therefore below the ground resolution.
Refraction correction II
II: Obtained pond depth is multiplied by the refractive index of water (n=1.335) to correct the vertical error:

The correction factor for pond depth rises exponential with the incident angle. Note the variant y-axis. We neglect the slow rise to 1.527 towards the maximum opening angle to keep the algorithm simple.
In situ measurements
We collected in total 52 manual depth measurements (shown by the spheres) in 14 ponds on two days, one clear sky day shown here and one cloudy day with only diffuse irradiance. The exact localization of the ground truth data is challenging. To compensate for inaccuracies, we average the data within a radius of 30cm. Pond depth ranged from 5cm to 40cm.
Accuracy

Pond bathymetry estimation compared to all in situ observations from a clear sky and a cloudy day. Point size is proportional to pond size. The greatest relative error is found in the three tiny ponds at 10 cm.
All in situ points were located towards the center parts of the ponds, at the edge the accuracy is probably slightly lower. The reason is that the elevation model is an estimated and smoothed reconstruction of the ice surface in which sharp edges, such as those that sometimes occur at the edges of ponds, are unlikely to be resolved.
Colorbar

Takk for turen - thanks for the trip!
Contact niels.fuchs@uni-hamburg.de for more information.

The data was recorded during RV Polarstern cruise PS106/1 PASCAL (Grant-No. AWI_PS106/1_2-00)

Acknowledgements:

Many thanks to the wonderful open source communities behind QGIS, Blender and MarziPano, which were used to create this tour.
Furthermore, special kudos to TomDee (MELLOWMESHER studio) for the great advice on 3D processing and illumination.
Start again?
Good luck! ;)
Outlook
Airborne remote sensing of ponds using photogrammetry yields unprecented opportunities to study the hydraulic pond system using drones, helicopters or airplanes.
With regular grid surveys as flown over the MOSAiC expedition Central Observatory, the whole pond system of the floe can be resolved in four dimensions: horizontally, vertically and with time.

Pond size, pond water level and pond depth measured on the MOSAiC floe on 30 June 2020 using photogrammetry.

For more MOSAiC data and a comparison to other methods please visit EGU21-12860 (Linhard et al.) on Friday 30 April.
A publication about the study presented here is in preparation.
Conclusion
Pros:

- photogrammetry allows for a high resolution and accurate pond bathymetry estimation (BIAS<1cm, RMSE=3.8cm)
- the method is based solely on aerial images taken with a consumer grade camera and thus, depending on the carrier platform (e.g. drones), relatively economical
- estimation of the whole pond bathymetry instead of single pond depth measurements
- coverage of large surfaces becomes possible to increase statistical robustness of pond geometry analysis
- accuracy independend of pond color, size and irradiance conditions (with few exceptions listed below)

Cons:

- requires dedicated survey flight pattern with great forward and lateral overlap of images
- water depth must not be too great in comparison to the image resolution (e.g. 2m limit for 5.6cm ground resolution)
- reconstructed surface inaccurate towards pond edges
- local depression of very small ponds (approx. <10x ground resolution) are not resolved in the elevation model
- pond bottom needs to have some recognisable structure (fulfilled in all our test studies)
- transparent pond surface required without sun glint