Analyses of Recent Sediment Surface Dynamic of a Namibian Kalahari Salt Pan Based on Multitemporal Landsat and Hyperspectral Hyperion Data

This study combines spaceborne multitemporal and hyperspectral data to analyze the spatial distribution of surface evaporite minerals and changes in a semi-arid depositional environment associated with episodic flooding events, the Omongwa salt pan (Kalahari, Namibia). The dynamic of the surface crust is evaluated by a change-detection approach using the Iterative-reweighted Multivariate Alteration Detection (IR-MAD) based on the Landsat archive imagery from 1984 to 2015. The results show that the salt pan is a highly dynamic and heterogeneous landform. A change gradient is observed from very stable pan border to a highly dynamic central pan. On the basis of hyperspectral EO-1 Hyperion images, the current distribution of surface evaporite minerals is characterized using Spectral Mixture Analysis (SMA). Assessment of field and image endmembers revealed that the pan surface can be categorized into three major crust types based on diagnostic absorption features and mineralogical ground truth data. The mineralogical crust types are related to different zones of surface change as well as pan morphology that influences brine flow during the pan inundation and desiccation cycles. These combined information are used to spatially map depositional environments where the more dynamic halite crust concentrates in lower areas although stable gypsum and calcite/sepiolite crusts appear in higher elevated areas.


Introduction
Natural salt pan environments, also known as salt flats or playas, are some of the most geomorphological dynamic environments on Earth [1,2].According to [3], there are approximately 50,000 salt pans on Earth.Although the area of individual salt pans is small, mostly less than a few square kilometers, except for prominent examples such as Salar de Uyuni (Bolivia) or Etosha pan (Namibia), the study of salt pan is important because of increasing population, agricultural demand and the hydrology of these dry regions becoming greatly stressed [4].Salt pans are also one of the biggest sources of mineral dust emission in arid regions [5,6] especially in the southern Africa Kalahari desert (Botswana, Namibia, and Republic of South Africa) [7] and the specific salt pan mineralogy influences the susceptibility for dust emission [8][9][10].For continued advances in understanding dust sources a spatially and temporally detailed assessment of their surface sediments is necessary [11].Salt pans are also a potential significant carbon stock in the Kalahari environment [12] and the dynamic of salt pan surfaces is of great importance for the understanding of the carbon cycle and carbon fluxes of a pan systems.Furthermore, the sensitive interrelation between regional changes in climate patterns and tectonic activity within continents is often record in the pan sedimentary deposit, and insights on current surface processes are needed for a proper interpretation of paleoenvironmental conditions [13].
The harsh nature [14] and limited accessibility of salt pan environments makes them hard to study by conventional field mapping.This problem is exacerbated by ephemeral nature of salt pan evaporite minerals that is caused by sedimentological processes occurring at the surface as well as changes of the brine composition and the associated sequence of minerals one can expect to form in episodic dry-wet-dry cycles.The saline pan cycle, described by, e.g., Lowenstein and Hardie [15], provides a conceptual model of pan process dynamic and its effect on the pan deposits and associated evaporite mineralogy.In this model, the stages of salt pan cycles are described beginning at the desiccation stage (dry saline pan), followed by the flooding stage (brackish pan), than by the evaporative concentration stage (saline pan) and finally the return to the first status, the desiccation stage (dry saline pan).The ephemeral nature of salt pan evaporite minerals also emphasizes the need for multitemporal observations [2,14].Over the last decades a limited number of studies have applied remote sensing techniques to study salt pan evaporite minerals and were able to provide spatial information on pan surfaces in terms of evaporite mineral discrimination supported by field observations.Most of them rely on traditional optical multi-spectral satellite sensors [1,13,[16][17][18][19] that often do not provide the necessary spectral resolution to cover the spectral complexity of evaporite mineral assemblages on salt pan surfaces [20] and only provide a mono-temporal assessment.In comparison, analysis based on advanced hyperspectral imagery, also named imaging spectroscopy, are more sensitive for the identification and quantification of a wider range of evaporite minerals [13,21], and usually achieve higher modeling results [20].EO-1 Hyperion is the first hyperspectral satellite that operates across the full solar-reflective spectrum with nominal spectral coverage form 0.4-2.5 µm and 10 µnm spectral response functions [22].Although the applications of Hyperion imagery is limited by the sensor's relatively low signal-to-noise ratio [16,23], spectral mapping results with Hyperion data compare reasonably well with results obtained from airborne hyperspectral sensors, including Hyperspectral Mapper (HyMap) [24] and Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) [25,26].Nevertheless, the availability of hyperspectral data, spaceborne or airborne, is very limited, which prevents time-series analyses of salt pan surfaces.
In this context, this paper is looking at the potential of modern remote sensing analyses for the study of the surface processes and dynamics of salt pans.This study combines both the potential of hyperspectral EO-1 Hyperion satellite data for the identification and mapping of major salt crust types and the potential of freely available 30 years Landsat time-series archives evaluating the temporal stability of the pan surface, linked with accurate elevation information from LIDAR data.The Omongwa salt pan in Namibian Kalahari, one of the largest and most studied salt pan in its region, is known for episodic flooding events and ongoing formation of evaporates minerals at the surface [27,28].The Omongwa salt pan is taken as the object of investigation in this paper, as representing other extended saline pan bodies in arid regions of the world that could be similarly studied with newly available remote sensing data.The pan's spatial extend allows it to be complete covered by current spaceborne hyperspectral datasets (Hyperion) and therefore presents a suitable opportunity to test the potential of such datasets on a regionally important site with coherent coverage.The focus of the paper is placed on the development of the relationships between evaporite mineral mapping, time-series based surface dynamic and pan morphology to differentiate depositional environment of the investigated salt pan.State-of-the-art methodological tools in optical remote sensing for change detection analyses based on multitemporal Landsat data and for surface spectral mapping based on hyperspectral imagery are used in order to deepen our knowledge and understanding of the surface processes in the Omongwa pan.

Study Area
The Omongwa pan is located in eastern Namibia near Aminuis, ~260 km southeast of Windhoek (Figure 1).With a size of ~20 km 2 , it is the largest one in a series of pans that appear to be aligned along a former branch of the Nossob River system, with a southward palaeodrainage direction [27,29,30].
The climate in the Aminuis region is hot and semi-arid.It classifies as BWh according to the Köppen scheme [31] with average annual temperature and yearly precipitation of 19.8 • C and 200-250 mm respectively.The rainfall is highly seasonal with 90% of the total precipitation occurring in the wet season from December to April.The interannual precipitation variation is very high, resulting in occasional draught years with down to ~40% of the average amount [32].The salt is very flat with a mean elevation of ~1200 m above sea level.The surrounding Kalahari landscape is characterized by an undulating linear dune system typically in a NW-SE direction (as shaped by the prevailing winds) with elevation magnitudes of ~1-3 m between dune crest and interdune valley.South of the pan, a lunette dune rises up to ~50 m above the pan floor level.
Remote Sens. 2017, 9, 170 3 of 24 in the wet season from December to April.The interannual precipitation variation is very high, resulting in occasional draught years with down to ~40% of the average amount [32].The salt is very flat with a mean elevation of ~1200 m above sea level.The surrounding Kalahari landscape is characterized by an undulating linear dune system typically in a NW-SE direction (as shaped by the prevailing winds) with elevation magnitudes of ~1-3 m between dune crest and interdune valley.South of the pan, a lunette dune rises up to ~50 m above the pan floor level.The pan evolved in calcareous deposits of Kalahari Group which are exposed in the paleo river valley and in interdune space of shallow dunes [27,33,34].Along these paleo river courses, Early Kalahari deposits were cemented by pedogenic and groundwater calcretes driven by re-occurring changes in moisture and sedimentary regime during the Tertiary until the Late Pleistocene [35].Prominent outcrops of calcareous material appear along the pan margins [36], e.g., ~1 m high scarp at the northeastern margin of Omongwa pan.Detailed mineralogical information of the pan deposits are known from a diagonal sampling survey from the SW to NE of Omongwa pan conducted by Mees [27] in 1991 and 1992.The upper lithological unit (top ~50 cm) of the pan deposits is mainly sand-sized detrital grains with varying contents in halite (>30%), gypsum (5%-40%) and calcite (10%-20%) mixed with the clay and mica minerals sepiolite, muscovite, and smectite.The lower units have The pan evolved in calcareous deposits of Kalahari Group which are exposed in the paleo river valley and in interdune space of shallow dunes [27,33,34].Along these paleo river courses, Early Kalahari deposits were cemented by pedogenic and groundwater calcretes driven by re-occurring changes in moisture and sedimentary regime during the Tertiary until the Late Pleistocene [35].Prominent outcrops of calcareous material appear along the pan margins [36], e.g., ~1 m high scarp at the northeastern margin of Omongwa pan.Detailed mineralogical information of the pan deposits are known from a diagonal sampling survey from the SW to NE of Omongwa pan conducted by Mees [27] in 1991 and 1992.The upper lithological unit (top ~50 cm) of the pan deposits is mainly sand-sized detrital grains with varying contents in halite (>30%), gypsum (5%-40%) and calcite (10%-20%) mixed with the clay and mica minerals sepiolite, muscovite, and smectite.The lower units have evaporite-filled cracks in the upper part (vadose zone) and generally increase in calcite content with depth until the calcareous mudstone bedrock is met in up to 3.5 m below the surface [27].It is further noted that the mineralogy of the pan is sensitive to changes of the hydrological conditions.The point based mineralogical description together with the analysis of the general pan environment is very helpful for the evaluation of the remote sensing based mapping of the surface mineralogy as it provides a baseline for the interpretation of this study.In general, the Kalahari Group is relatively thin in this region with a thickness of <10 m [37], which has been recognized from the cuttings of several water boreholes in the Aminuis region [35,38].One borehole (WW26165) directly north of the Omongwa pan revealed only 1 m of Kalahari Group deposits [38].They overlie more than 50 m of Pre-Kalahari basalts of the Kalkrand Formation.The surroundings landscape of the pans are generally covered by the reddish sands of the Kalahari Group [30,35], which are classified as Rubi-Feralic Arenosols, ARo (red, iron rich, sandy soils) according to the WRB system [39], whereas the calcareous deposits along the paleo riverbed form greyish soils that classify as Petric Calcisols, CLp (soils with secondary lime accumulation and strongly cemented layers) [34].
The pan lacks significant surface inflow, as surface runoff is minimized by the high infiltration rates of the Kalahari sands and topographically limited to by distance of the next dune crest.The hydrogeological situation of the study region is very complex due to the occurrence of several pre-Kalahari normal faults as well as rapid facies changes in the Kalahari Group deposits [35,40].Borehole data near Aminuis suggest that the groundwater level is reached below the Kalahari Group deposits between 13 and 43 m in the Pre-Kalahari basalt layers [38].However the calcareous mudstone of the Omongwa pan is likely to support a locally perched groundwater tables above the regional aquifer due to its low permeability.This is supported by the observations of fluctuating groundwater levels of 25-150 cm between the wet and the dry season of 1991-1992 [27] as well as temporally formation of shallow open waterholes fed by groundwater discharge [30].

Field and Laboratory Analysis
Two field campaigns for surface characterization and sampling took place in October 2014 and June 2015, which represent the end and beginning of the dry season respectively.Prior to both visits small precipitation events of <5 mm led to slight moist conditions in the central part of the Omongwa pan [41], which limited the accessibility to the most outer ~500 m of the pan, and lead to a sample size of 17 surface samples.
Fieldwork consisted of acquisition of surface samples and field reflectance spectroscopy measurements along transects into the pan and for sites of interest representative for the variability of surface conditions and presumed sediment mixtures of the pan (see Figure 1).Sampling was performed in homogeneous areas collecting the top surface crust (<2 cm).These samples were composed of 5 to 10 sub-samples collected at random locations within a 5 m wide square around the center point for each site of interest.Prior to sampling, field spectra of the untouched surfaces were collected at the selected locations using an FieldSpec 3 point spectroradiometer (ASD Inc., Boulder, CO, USA), covering the VNIR-SWIR spectral range [42].Per site of interest, 5 spectra were acquired at nadir configuration from 1 m height with a ~25 cm target radius from the center of each location.The measurements were calibrated to a barium-sulfate white reference standard to determine relative reflectance, and were subsequently corrected for detector offset and averaged for higher signal-to-noise ratio.All spectra and center of sampling location were geolocated using the Global Positioning System (GPS).
Geochemical and VNIR-SWIR spectral analyses of surface samples were conducted in the laboratory including the parameters mineralogy, grain size, pH, and conductivity.Initial sample preparation comprised of sub-sampling and drying at 105 • C.After the removal of organic material and carbonates using 10% H 2 O 2 and 17% HCl, respectively, sample texture was determined by a combination of wet-sieving (particles > 63 µm) and sedimentation (particles < 63 µm) following [43].Electrical conductivity was measured in a 1:5 ratio with bi-distilled water after [44] and pH in a 0.01 M CaCl 2 suspension after [43].The sub-samples used for powder X-ray diffraction (XRD) were crushed in a vibrating disc mill to a particle size <63 µm.Powder samples were measured using Cu-Kα radiation.Diffractograms were recorded from 5 • to 85 • 2θ using a PANalytical Empyrean diffractometer (PANalytical GmbH, Kassel, Hesse, Germany).Results were analyzed using the software Bruker EVA (Bruker Corporation, Billerica, MA, USA) for mineral identification and AUTOQUAN (Agfa NDT Pantak Seifert GmbH & Co. KG, Ahrensburg, Schleswig-Holstein, Germany) for quantification.Laboratory spectral characterisation was performed with the ASD FieldSpec 3 point spectroradiometer under controlled light illumination simulating spaceborne observations (sensor nadir viewing, light source azimuth 35 • ).Finally, two spectral libraries of Omongwa sediment surfaces were created associated with field and laboratory optical signatures.

Remote Sensing Change Detection Analyses
The multitemporal analysis uses all available, cloud free Landsat images of the study area at World Reference System-2 path/row 176/76 acquired through the USGS EarthExplorer data portal [45].The time series covers the last 3 decades (1984-2015) and is composed of imagery from three different satellite and sensor generations of the Landsat Mission.From a total of 146 cloud-free scenes, the majority (98) originate from Landsat 5 (TM), 16 from Landsat 7 (ETM+) and 32 from Landsat 8 (OLI).The three sensors have the same spatial resolution of 30 m, but slightly different spectral coverage or even additional bands [46].Only the six reflective bands (Blue, Green, Red, NIR, MIR, and SWIR) common to all sensors are used.On average, the data series is composed of five images per year, but the distribution is highly irregular.Data gaps occur during the mid-1980s and from late 2011 to early 2013 when Landsat 5 data were discontinued and Landsat 8 was not yet operative.Between 2003 and 2013, data availability further decreased to 1-4 scenes/year after the failure of the Scan Line Corrector (SLC) of the Landsat 7 ETM+ sensor, which makes the data unusable for change-detection approaches such as in this study.
The Landsat data were pre-processed to top of atmosphere (TOA) reflectance by the Landsat Ecosystem Disturbance Adaptive Processing System (LEDAPS) [47] including the conversion of calibrated Digital Numbers (DNs) to absolute units of at-sensor spectral radiance and reduction of in scene-to-scene variability due to different solar/sensor geometry and Sun/Earth distance.Top of atmosphere (TOA) reflectance is used to normalize scene radiances differences due to variations in solar illumination, sensor viewing geometry, and seasonality (Earth-Sun distance).Surface reflectance products for all optical bands are also available [48,49].However the atmospheric correction applied by LEDAP follows the popular dense vegetation (DDV) approach [50], which is known to provide limited accuracy in arid environments [51], as it depends on finding sufficient and evenly distributed dark objects such as dense vegetation or water bodies that are completely absent in the study area.Because of the lack of such surfaces in the study area, it is unclear on which basis the estimation of the aerosol properties is performed and it is likely that an artificial bias is introduced to the data instead of a correction.Furthermore the provisional AOT estimation for the Landsat 8 OLI images is based on coarse (0.05 • ) MODIS data [52], which can lead to tilling artifacts in the processed image [49] that were also identified on images of the Omongwa pan.Due to the current limitations in corrections to surface reflectance especially applied in arid areas, TOA reflectance is used in this paper for the change detection analyses.An adapted procedure for change detection analysis is used that limits the atmospheric influence.
In a first step, to restrict the change detection analyses to the dynamic of the pan surface deposits in dry state, Landsat images related to pan flooded state have to be excluded from the dataset.For this purpose, Xu's Normalized Differenced Water Index (NDWI) [53], a band ratio using the green and mid-infrared (MIR) band, is derived from the Landsat record.This adaption of the NDWI is sensitive even to shallow water bodies [54].A conservative threshold of 0.5 was selected to filter scenes featuring surface water.This restriction robustly excludes images that show an inundated pan surface, as well as scenes which have significant moisture in the top layer of pan deposits.Exactly how the selected NDWI threshold relates to the water content of the surface sediments is scene dependent and cannot be profoundly estimated from the remote sensing image.A total of 30 scenes were water covered or in very wet surface conditions and discarded for further analyses.
On the remaining 116 Landsat images a change detection analysis is applied.The described discontinuities and gaps of the Landsat record restrict popular approaches of time series analyses based on temporal trajectories such as the Breaks For Additive Season and Trend (BFAST) [55] and Reversed Ordered Cusum (ROC) [56] procedures.Discontinuous datasets do not meet the need of many curve-fitting time series methods to establish a "stable" reference model from historical data [57], e.g., 16-day time series over 2 years [58].Accordingly, a temporally more robust approach of bi-temporal change detection, the Iterative-reweighted Multivariate Alteration Detection (IR-MAD) developed by A.A. Nielsen [59] is selected.The IR-MAD has been successfully used for different applications and sensors, e.g., to detect changes caused by floods from coarse AVHRR imagery [60], mining activity [61,62] and forest disturbance from Landsat data [63], as well as earthquakes from high resolution imagery of GF-1/PMS and GeoEye-1 [64,65].The core idea of MAD is to find maximum differences between two images by removing correlations between them as much as possible [66].The MAD algorithm transforms two sets of multivariate observations (in this case two multispectral satellite images covering the same geographical area acquired at different points in time) into a difference between two linear combinations of the original variables explaining maximal change (i.e., the difference explaining maximal variance) in all variables (wavelength) simultaneously [67].IR-MAD is an improvement of the MAD approach.In order to increase the detection of observations, whose status over time is uncertain, weights are put on difficult observations that exhibit small change according to their change probability in the previous run.In each iteration a better no-change background is identified, achieving a better separability between the change and no-change classes [59].
The IR-MAD can be classified as a bi-temporal linear data transformation type of change detection [68].This change-detection approach is most suitable to observe the type of changes that can be expected in the salt pan environment, due to the following properties: (1) It can detect modifications (changes in crust composition) and not only discrete conversions (crust type A to crust type B); (2) it is sensible to progressive as well as abrupt changes (slow deflation of particles, loss and build-up of new crust after inundation); (3) it can work with the temporal sampling rates of the observations (irregular, data gaps) [68].A further advantage over most other change detection methods is that IR-MAD is invariant to linear transformations of the original image intensities [59,69].The method is robust to differences in atmospheric conditions or sensor calibrations at two acquisition times, which greatly reduces the requirements on radiometric correction or normalization [66].This is an important advantage for this study because it limits the atmospheric influence on the change detection result when using TOA reflectance data.A general constraint to bi-temporal change detection methods such as the IR-MAD is the loss of data due to a number of contaminations or errors such as SLC-off data gaps [57] and cloud cover [70].Although well-performing correction approaches have been proposed for data gaps with information derived from neighboring pixels [71] or even other sensors [72], interpolation or the inclusion of extraneous data can introduce additional errors [73,74].For this, SLC-off acquisitions as well as clouded scenes are discarded.
In this study, the python implementation of IR-MAD by [75] is used to which some modifications were made that allow for batch processing of the Landsat image stack.For the Omongwa pan subset a total of 115 bi-temporal change maps are produced.The IR-MAD processing always converged to a no-change background before hitting the defined maximum of 50 iterations.Besides providing the basis for the no-change probability calculation, the sum of squared, standardized MAD variates, also called χ 2 image [76] is used as a general measure of change of the pan surface similar to [65,77,78].From each bi-temporal IR-MAD analysis of the Landsat dataset, the final χ 2 image is used to derive aggregated change over the time-series to differentiate dynamic pan areas from radiometrically more stable regions.

EO-1 Hyperion Mineralogical Mapping
One EO-1 Hyperion hyperspectral image is used in this study.It has been acquired on 7 September 2014 at 07:34 UTC with a sun azimuth of 63 • and a sun elevation of 36 • .The Hyperion sensor covers the VNIR and SWIR spectral region (400-2500 nm) with 198 calibrated bands of ~10 nm width at a spatial resolution of 30 m [79] similar to the Landsat imagery.The spatial extent is 7.5 km wide 105 km long.The signal-to-noise ratio (SNR) is specified as ~140:1 to 190:1 for the VNIR and 96:1-38:1 for the SWIR detector [79].
An automated pre-processing chain for EO-1 Hyperion data developed by [80] is used that provides georeferenced surface reflectance data.The L1 radiance data (L1R), L1T radiance terrain corrected data, acquisition parameters from the metadata and an ASTER 30 m DEM, freely distributed through the USGS EarthExplorer data portal [45] are used.Data processing includes rescaling to radiance, removal of the spectrally overlapping bands, bad band detection, reductions for dead pixel and erroneous detector columns, as well as corrections for intra-band spatial shifts, keystone, erroneous co-registration, radiometric miscalibration.Part of the processing chain is the atmospheric correction as proposed in [81,82] originally designed for the EnMAP Box [83].It includes a complex Radiative Transfer Modeling (RTM), Aerosol Optical Thickness (AOT), Columnar Water Vapor (CWV) retrieval as well as spectral Smile and adjacency correction [80].Georeferencing is performed via GCP (Ground Control Point) selection and first-order transformation using ENVI 5.3 [84] between the corrected image and the orthorectified L1T Hyperion image.The geometrical transformation is applied at the end on the results to prevent resampling effects on the spectral analyses.The pre-processing chain was used previously for geological applications, e.g., by Mielke et al. [85,86] proving that it can deal with the sensor limited quality and still delivers reflectance data that are adequate for mineralogical identification analyses.
Spectral mixture analyses (SMA) [87,88] implemented in ENVI 5.3 [84] is used to determine the mineralogy of the pan surface based on the EO-1 Hyperion image.SMA is a physically-based model that transforms radiance or reflectance values to physical variables that are linked to the subpixel abundances of surface components within each pixel [89,90].SMA proved useful specifically for the mapping of salt and gypsum surfaces [13,91], but also in a wide range of other applications, such as vegetation (e.g., [92]), mineralogical and lithological mapping (e.g., [93,94]), and soil degradation assessment [95].The SMA analysis assumes a linear mixing of the scene constituents in the sensor field-of-view and is particularly adapted to arid landscapes with aerial mixing.In SMA, the reflectance of a pixel is determined by the sum of the reflectance of each material within a pixel multiplied by its fractional cover.Most crucial for the quality of an unmixing model is the selection of endmembers that account for the spectral variability of the data, while remaining physically meaningful [89,96].Instead of focusing on specific absorption features, SMA uses the whole reflectance of the VNIR-SWIR selected endmember to map subpixel abundances.This is especially helpful in the study of salt affected surfaces that often do not show diagnostic absorption features outside of the known water (vapor) absorption region [97].The inherent dimensionality of the EO-1 Hyperion data is examined via minimum noise fraction transformation [98] and image endmembers are extracted from the edges of the data cloud.Following [99], candidate endmembers extracted from the image are further redefined in minimum noise fraction transformed space to assign the purest and least noisy pixel spectra as final endmembers.In order to focus on the fractional cover of different mineralogical crust types and to reduce the spectral variability of the data set, the areas outside the Omongwa pan are excluded from the image for SMA analyses using an albedo threshold.

LiDAR Digital Surface Model
Airborne Light Detection and Ranging (LiDAR) data were acquired with a Riegl LMS-Q780 at an altitude of 2850 m above ground level over the Omongwa pan on 6 June 2015 during a GFZ/DIMAP airborne campaign [100].The pan surface was partly moist, but not flooded at the time of overflight.A total of nine flightlines were recorded with a swath width of ~3 km and an overlap of approximately 60% and an average point density of 2.1 pts/m 2 .Fundamental accuracy assessed through calibration flights showed ±10 cm vertical and ±15 cm horizontal accuracy within a of 95% confidence level.For the topography of the Omongwa pan a 1 m grid Digital Surface Model (DSM) was generated from the first return LiDAR mass point data by converting the points to a Triangulated Irregular Network (TIN) surface mesh and the TIN to a regular raster grid with Natural Nearest Neighbor interpolation using the program Terrasolid TerraScan.Values are given in meters above sea level (a.s.l.)

Field and Laboratory Analysis
Field observations showed highly heterogeneous pan environments that could be grouped into three different crust types (Figure 2a-c): (a) halite (NaCl) is present over the whole pan surface in variable amounts, and is identifiable as bright white surfaces, with puffy structure and blisters; (b) the sulfate mineral gypsum (CaSO 4 •2H 2 O) is a major constituted of the crust at the pan margin and appear as a smoother and darker surface; and (c) carbonate crust appears in the form of some sediment accumulation mixed with salt composition, linked with carbonate outcrops at the border of the pan (Figure 2e).Surface roughness as well as the general brightness are very variable between the crust types and are also effected by anthropogenic disturbances (Figure 2f).In general, the pan surface is bare of any vegetation with the exception of the most northern part that is sparsely covered with patchy bushes that increase in density in the direction of the pan margin (Figure 2d).

Field and Laboratory Analysis
Field observations showed highly heterogeneous pan environments that could be grouped into three different crust types (Figure 2a-c): (a) halite (NaCl) is present over the whole pan surface in variable amounts, and is identifiable as bright white surfaces, with puffy structure and blisters; (b) the sulfate mineral gypsum (CaSO4•2H2O) is a major constituted of the crust at the pan margin and appear as a smoother and darker surface; and (c) carbonate crust appears in the form of some sediment accumulation mixed with salt composition, linked with carbonate outcrops at the border of the pan (Figure 2e).Surface roughness as well as the general brightness are very variable between the crust types and are also effected by anthropogenic disturbances (Figure 2f).In general, the pan surface is bare of any vegetation with the exception of the most northern part that is sparsely covered with patchy bushes that increase in density in the direction of the pan margin (Figure 2d).From field spectral measurements as well as laboratory geochemical analysis, the major mineralogical components in the salt pan sediments could be identified and related to the different crust types.Figure 3 shows the field spectra for the three representative crusts types and Table 1 the geochemical results of all samples.According to the XRD analysis, Omongwa halite field endmember (sample 141, Figure 2a) is almost pure halite (NaCl) (94%) mixed with small amounts of gypsum From field spectral measurements as well as laboratory geochemical analysis, the major mineralogical components in the salt pan sediments could be identified and related to the different crust types.Figure 3 shows the field spectra for the three representative crusts types and Table 1 the geochemical results of all samples.According to the XRD analysis, Omongwa halite field endmember (sample 141, Figure 2a) is almost pure halite (NaCl) (94%) mixed with small amounts of gypsum (CaSO 4 •2H 2 O) (3%).The associated field spectrum shows a constant high reflectance level in the VNIR with almost no slope, which is typical for pure halite [101] or halite rich salt pan sediments [24].The absorption and general lower reflectance level of the Omongwa halite field spectrum compared to standard laboratory halite powder spectral reflectance from the USGS spectral library [102] originates from the partially moist conditions at the time of field spectral acquisition as well as from the intermixture with some gypsum at that location.Laboratory spectral measurements of this sample after air drying for several days result in reflectance level up to 70% in the VNIR, close to reference halite spectral characteristics.The field spectrum of the second crust type (sample P65, Figure 2b) is characterized by strong absorption features in the SWIR at ~1750 nm and ~2200 nm as well as the unique triplet absorption around 1500 nm, which are characteristic for gypsum spectral properties [101,103,104].The XRD analysis also confirms the high gypsum content (83%) mixed with some halite.The third surface crust type (sample P68, Figure 2c) is dominated by high calcite (CaCO 3 ) content (45%), mixed with the clay mineral sepiolite (Mg 4 Si 6 O 15 (OH) 2 •6H 2 O) (16%).The associated field spectrum shows a strong absorption feature at ~2340 µm.Both calcite and sepiolite have their characteristic absorption features at this wavelength [101].
Remote Sens. 2017, 9, 170 9 of 24 absorption and general lower reflectance level of the Omongwa halite field spectrum compared to standard laboratory halite powder spectral reflectance from the USGS spectral library [102] originates from the partially moist conditions at the time of field spectral acquisition as well as from the intermixture with some gypsum at that location.Laboratory spectral measurements of this sample after air drying for several days result in reflectance level up to 70% in the VNIR, close to reference halite spectral characteristics.The field spectrum of the second crust type (sample P65, Figure 2b) is characterized by strong absorption features in the SWIR at ~1750 nm and ~2200 nm as well as the unique triplet absorption around 1500 nm, which are characteristic for gypsum spectral properties [101,103,104].The XRD analysis also confirms the high gypsum content (83%) mixed with some halite.The third surface crust type (sample P68, Figure 2c) is dominated by high calcite (CaCO3) content (45%), mixed with the clay mineral sepiolite (Mg4Si6O15(OH)2•6H2O) (16%).The associated field spectrum shows a strong absorption feature at ~2340 µm.Both calcite and sepiolite have their characteristic absorption features at this wavelength [101].

EO-1 Hyperion Analyses
Four images endmembers are extracted from the Hyperion imagery and shown in Figure 4. Overall the main spectral properties of the first three extracted endmembers are similar to the field spectra of the different crust types.However some differences are observed due to variable sensor characteristics and surface conditions.The extracted Hyperion endmember spectrum (EM 1 in Figure 4) for halite crust has a high reflectance level around 70% and significant spectral slope in the shorter wavelength compared to the field halite spectrum.The higher general albedo can be attributed to the dry conditions of the central pan surface at the time of Hyperion image acquisition, whereas the differences in the shorter wavelength could result from a higher proportion of silicates exposed due to the larger Hyperion pixel or that are present in the central part of the pan not accessible during the field campaigns.In the image endmember that corresponds to the gypsum crust (EM 2 in Figure 4), the triplet absorption around 1500 nm visible in the gypsum field spectrum is within the water vapor atmospheric absorption bands and cannot be observed, but the other absorption features at ~1750 and ~2200 nm are well defined.The characteristic absorption feature of the calcite and sepiolite crust at ~2340 nm exceeds the usable spectral range of the Hyperion scene.However the left slope and shoulder of the absorption feature can still be identified at ~2300 nm in the image endmember (EM 3 in Figure 4).Because the observed absorption at ~2300 nm can originate from both minerals, we refer this endmember as calcite/sepiolite.In the surface reference samples analyzed by XRD the minerals mostly occur together with ~2-3× higher share of calcite.Besides the spectrally and mineralogically well-defined endmembers of halite, gypsum, and calcite/sepiolite, a fourth image endmember had to be added to account for the spectral variability within the scene.This endmember is associated with disturbed halite crust of lower reflectance (EM 4 in Figure 4).In general, all endmembers have their abundances well mapped with no component exceeding 1.The accuracy of the unmixing model is assessed through the root mean squared error (Figure 5e).The RMSE image shows that the overall pan surface is well explained with the linear mixture of the four selected endmembers except for higher RMSE areas at: (1) the northeast corner and western pan border associated with heavy cattle disturbances; (2) the mainly diagonal SW-NE lines crossing the pan, which are used as driving shortcuts across the pan during the dry season as well as cattle path ways; and (3) the parallel, vertical stripes (~11°), associated with remaining sensor noise that were not removed through the pre-processing procedure.In general the performance of the unmixing model is very good with an overall mean RMSE of ~1% comparable to equivalent studies, e.g., [105] achieved an mean RMSE of ~2% using linear unmixing on Hyperion data to map salt effect soils or [106] getting the same accuracy range for the abundances of vegetation, soil and limestone bedrock.The SMA mapping shows that more than half of the pan surface is dominated by halite crust in variable amounts (Figure 5a).The highest abundances are mapped in the central part and in a circular shape in the southern half of the pan, whereas the border region and the anthropogenic pathways have minimal fractions of halite crust.In most regions the distribution of gypsum crust appears to be inversely related to the halite crust.Highest gypsum composition is found at the southern pan border, in the south-central circular shape and in a ~1.5 km long stripe in the northern pan area (Figure 5b).The occurrence of calcite and sepiolite crust is confined to the northern part of the pan and to the pan border (Figure 5c) close to outcrops of the carbonatic bedrock that were observed in the field (Figure 2e).The disturbed, dark crust endmember (Figure 5d) is mainly associated with the cattle and human disturbance of the salt crust through the pathways crossing the pan and at the border regions (e.g., Figure 2f).In these areas the top crust is at partly removed and the subsurface material is exposed.The subsurface is usually less concentrated in evaporite minerals and more mixed with the clastic, quartz-rich material of the pan.This is reflected by the geochemical analyses of subsurface field samples 142 (halite: 9%, quartz: 63%) and 172 (halite: 17%, quartz: 41%), directly sampled under the surface crust of samples 141 (halite: 94%, quartz: 1%) and 171 (halite: 50%, quartz: 27%), respectively at a depth at a depth of ~10 cm.The sub-surface samples show decreased halite and increased quartz content compared to surface samples (Table 1).In addition to mineralogical differences, the physical properties of the salt pan surface are altered by cattle trampling, increasing the surface roughness [107], which may also cause the reflectance to decrease even at similar salt content by casting micro  The accuracy of the unmixing model is assessed through the root mean squared error (Figure 5e).The RMSE image shows that the overall pan surface is well explained with the linear mixture of the four selected endmembers except for higher RMSE areas at: (1) the northeast corner and western pan border associated with heavy cattle disturbances; (2) the mainly diagonal SW-NE lines crossing the pan, which are used as driving shortcuts across the pan during the dry season as well as cattle path ways; and (3) the parallel, vertical stripes (~11 • ), associated with remaining sensor noise that were not removed through the pre-processing procedure.In general the performance of the unmixing model is very good with an overall mean RMSE of ~1% comparable to equivalent studies, e.g., [105] achieved an mean RMSE of ~2% using linear unmixing on Hyperion data to map salt effect soils or [106] getting the same accuracy range for the abundances of vegetation, soil and limestone bedrock.The SMA mapping shows that more than half of the pan surface is dominated by halite crust in variable amounts (Figure 5a).The highest abundances are mapped in the central part and in a circular shape in the southern half of the pan, whereas the border region and the anthropogenic pathways have minimal fractions of halite crust.In most regions the distribution of gypsum crust appears to be inversely related to the halite crust.Highest gypsum composition is found at the southern pan border, in the south-central circular shape and in a ~1.5 km long stripe in the northern pan area (Figure 5b).The occurrence of calcite and sepiolite crust is confined to the northern part of the pan and to the pan border (Figure 5c) close to outcrops of the carbonatic bedrock that were observed in the field (Figure 2e).The disturbed, dark crust endmember (Figure 5d) is mainly associated with the cattle and human disturbance of the salt crust through the pathways crossing the pan and at the border regions (e.g., Figure 2f).In these areas the top crust is at partly removed and the subsurface material is exposed.The subsurface is usually less concentrated in evaporite minerals and more mixed with the clastic, quartz-rich material of the pan.This is reflected by the geochemical analyses of subsurface field samples 142 (halite: 9%, quartz: 63%) and 172 (halite: 17%, quartz: 41%), directly sampled under the surface crust of samples 141 (halite: 94%, quartz: 1%) and 171 (halite: 50%, quartz: 27%), respectively at a depth at a depth of ~10 cm.The sub-surface samples show decreased halite and increased quartz content compared to surface samples (Table 1).In addition to mineralogical differences, the physical properties of the salt pan surface are altered by cattle trampling, increasing the surface roughness [107], which may also cause the reflectance to decrease even at similar salt content by casting micro shadows [108].The rougher surface crust is also likely to be more moist [109] adding to the decrease in reflectance.

Remote Sensing Change Detection Analyses
Figure 6 shows the result of the IR-MAD analysis.It represents the mean of the 115 IR-MAD χ2 images derived from the Landsat scene dry surface conditions of the period 1984-2015.The change magnitude map reveals that the salt pan is a highly dynamic and also heterogeneous landform.Compared to the changes of pan surface, the magnitude of spectral change in the surrounding Kalahari savannah is insignificant, which include the seasonal dynamic of the vegetation.Other change signals outside of the pan are of anthropogenic origin associated with housing activity in the southwest of the pan (see few pixels south of the pan in Figure 6) and a gravel road in the north of the pan.Within the pan, a gradient of change is determined from a very stable pan border to a highly dynamic central pan.The small northern sub-pan, where P68 and P69 are located (Figure 1), is very stable compared to the main pan body.The changes in the central part of the pan are very heterogeneous.Two regions of highest change (red colour in Figure 6) are identified located in the north of a circular area in the central southern part of the pan and a corridor in the northern part of the pan.In relation to the highest change in the pan center, several distinctive regions of lower change can be identified.For example straight lines crossing the pan as well as the mentioned circular shape have lower change magnitude compared to their surroundings.The largest of these lines crosses the

Remote Sensing Change Detection Analyses
Figure 6 shows the result of the IR-MAD analysis.It represents the mean of the 115 IR-MAD χ 2 images derived from the Landsat scene dry surface conditions of the period 1984-2015.The change magnitude map reveals that the salt pan is a highly dynamic and also heterogeneous landform.Compared to the changes of pan surface, the magnitude of spectral change in the surrounding Kalahari savannah is insignificant, which include the seasonal dynamic of the vegetation.Other change signals outside of the pan are of anthropogenic origin associated with housing activity in the southwest of the pan (see few pixels south of the pan in Figure 6) and a gravel road in the north of the pan.Within the pan, a gradient of change is determined from a very stable pan border to a highly dynamic central pan.The small northern sub-pan, where P68 and P69 are located (Figure 1), is very stable compared to the main pan body.The changes in the central part of the pan are very heterogeneous.Two regions of highest change (red colour in Figure 6) are identified located in the north of a circular area in the central southern part of the pan and a corridor in the northern part of the pan.In relation to the highest change in the pan center, several distinctive regions of lower change can be identified.For example straight lines crossing the pan as well as the mentioned circular shape have lower change magnitude compared to their surroundings.The largest of these lines crosses the pan center and is associated with a pathway that connects settlements in the west and northeast of the pan.
Remote Sens. 2017, 9, 170 13 of 24 pan center and is associated with a pathway that connects settlements in the west and northeast of the pan.

Pan Surface Characterisation
The combined remote sensing results over the Omongwa pan based on Landsat change detection analyses, Hyperion mineralogical identification and mapping and the LIDAR elevation model are summarized in Figure 7. Hyperion mineralogical identification agrees with the findings of previous studies [27] based on conventional point-based field sampling.The main mineralogical constituents of the top layer of the pan deposits described in [27], namely: quartz, halite, gypsum, calcite and sepiolite could be identified in our XRD analyses and furthermore Hyperion satellite imagery was able to provide spatial distribution of the major crust types and abundance maps over the whole pan.However, no indication of muscovite or smectite (described to occur in traces in [27]) were found in our analyzed field samples nor were these minerals recognized in the spectral properties of in-situ or laboratory measurements.This discrepancy may be the result of different sampling locations or depth.In this study the top <2 cm of the pan deposits were sampled, unfortunately [27] does not provide spatial differentiation of the mineralogy or state the specific sampling depth, which prevents further evaluation.Another possible explanation is a change in exposed surface mineralogy.The mineralogical data published by [27] were sampled in 1991 and 1992.During this >20 year time frame the surface alteration, e.g., by redistribution of surface particles due to surface flooding [15] and or wind erosion [9] are possible scenarios.The formation of fresh halite crust after significant rainfall and flooding events may superimpose underlying less soluble (e.g., clastic) minerals [110] and hide these (minor) components from surface observations, especially in the central halite dominated, dynamic part of the salt pan.The detection of minor crust components may be also limited by the low signal-to-noise ratio of the Hyperion data [16,23].However, the spectral properties of the most abundant crust components are represented in the selected

Pan Surface Characterisation
The combined remote sensing results over the Omongwa pan based on Landsat change detection analyses, Hyperion mineralogical identification and mapping and the LIDAR elevation model are summarized in Figure 7. Hyperion mineralogical identification agrees with the findings of previous studies [27] based on conventional point-based field sampling.The main mineralogical constituents of the top layer of the pan deposits described in [27], namely: quartz, halite, gypsum, calcite and sepiolite could be identified in our XRD analyses and furthermore Hyperion satellite imagery was able to provide spatial distribution of the major crust types and abundance maps over the whole pan.However, no indication of muscovite or smectite (described to occur in traces in [27]) were found in our analyzed field samples nor were these minerals recognized in the spectral properties of in-situ or laboratory measurements.This discrepancy may be the result of different sampling locations or depth.In this study the top <2 cm of the pan deposits were sampled, unfortunately [27] does not provide spatial differentiation of the mineralogy or state the specific sampling depth, which prevents further evaluation.Another possible explanation is a change in exposed surface mineralogy.The mineralogical data published by [27] were sampled in 1991 and 1992.During this >20 years time frame the surface alteration, e.g., by redistribution of surface particles due to surface flooding [15] and or wind erosion [9] are possible scenarios.The formation of fresh halite crust after significant rainfall and flooding events may superimpose underlying less soluble (e.g., clastic) minerals [110] and hide these (minor) components from surface observations, especially in the central halite dominated, dynamic part of the salt pan.The detection of minor crust components may be also limited by the low signal-to-noise ratio of the Hyperion data [16,23].However, the spectral properties of the most abundant crust components are represented in the selected endmember spectra of the generalized major crust types.Figure 7 shows that the different areas of surface change appear to be well correlated to the crust mineralogical mapping as well as pan morphology.Areas associated with high change magnitudes are located in the central part of the pan, which is dominated by halite mineralogy and lower topographical areas.Gypsum crust is the second most abundant crust type of the Omongwa pan (Figure 7b).Unlike the halite crust, the gypsum area mapped with Hyperion is associated with more stable or intermediately dynamic surface regions in the multitemporal analysis (Figure 7a) and slightly higher topography (Figure 7c).Although the difference in pan height of 10-15 cm between halite dominated and gypsum mixed crust is rather small, it seems to be enough to affect the surface mineralogy.
Remote Sens. 2017, 9, 170 14 of 24 which is dominated by halite mineralogy and lower topographical areas.Gypsum crust is the second most abundant crust type of the Omongwa pan (Figure 7b).Unlike the halite crust, the gypsum area mapped with Hyperion is associated with more stable or intermediately dynamic surface regions in the multitemporal analysis (Figure 7a) and slightly higher topography (Figure 7c).Although the difference in pan height of 10-15 cm between halite dominated and gypsum mixed crust is rather small, it seems to be enough to affect the surface mineralogy.Figure 8 presents the relationship between endmember abundance and digital elevation model resampled to Hyperion spatial scale.For clarity not all Hyperion pixels are shown, rather a Random Stratified Sampling (RSS) pixel selection is represented.RSS is used to respect the pixel distribution over all elevations and put emphasis on the extreme elevations that represent ~33% of the total pan area.The figure shows that pan depressions (areas lower than 1227.25 m, ~18% of the total pan area) are dominated by halite crust (box A in Figure 8) and higher topographical positions (areas higher than 1227.4m, ~15% of the total pan area) are dominated by gypsum crust (box B in Figure 8), whereas the majority of pixels with more mixed crust abundances are present at the medium elevations.Figure 8 presents the relationship between endmember abundance and digital elevation model resampled to Hyperion spatial scale.For clarity not all Hyperion pixels are shown, rather a Random Stratified Sampling (RSS) pixel selection is represented.RSS is used to respect the pixel distribution over all elevations and put emphasis on the extreme elevations that represent ~33% of the total pan area.The figure shows that pan depressions (areas lower than 1227.25 m, ~18% of the total pan area) are dominated by halite crust (box A in Figure 8) and higher topographical positions (areas higher than 1227.4m, ~15% of the total pan area) are dominated by gypsum crust (box B in Figure 8), whereas the majority of pixels with more mixed crust abundances are present at the medium elevations.over all elevations and put emphasis on the extreme elevations that represent ~33% of the total pan area.The figure shows that pan depressions (areas lower than 1227.25 m, ~18% of the total pan area) are dominated by halite crust (box A in Figure 8) and higher topographical positions (areas higher than 1227.4m, ~15% of the total pan area) are dominated by gypsum crust (box B in Figure 8), whereas the majority of pixels with more mixed crust abundances are present at the medium elevations.

Interpretation of Pan Surface Processes and Depositional Environment
The interpretation of the observed differences in pan surface dynamic, mineralogy and morphology based on the remote sensing findings confirms that Omongwa pan follows the conceptual model of the depositional cycle of salt pan environments given by Lowenstein and Hardie [15], and that this sediment dynamic can be spatially mapped with advanced remote sensing analyses such as in this paper.The conceptual model describes the reaction of pan deposits to a precipitation event.The stages of these pan cycle are: (1) desiccation stage (dry saline pan); (2) flooding stage (brackish pan); (3) evaporative concentration stage (saline pan); and the return to (1) desiccation stage (dry saline pan).Starting from the dry pan state (stage 1) ephemeral inundation events occur that cover parts or rarely even total surface area of the pan.Evaporite minerals will be dissolved by undersaturated meteoric water and redistributed over the pan mainly following the topography (stage 2).Halite is one of the less stable (more soluble) evaporite minerals and will be quickly leached into the brine and removed from the uppermost surface layer [111], whereas gypsum and calcite are more stable (less soluble) [112] and therefore are more likely to remain in the surface crust in higher abundancies, which is what we observe in the Hyperion mapping where high gypsum and calcite areas are more concentrated on the borders of the pan, which are also the most stable areas as identified in the Landsat time-series change analyses.After some time of continued evaporation the increasingly concentrated brine will accumulate in the local depressions (stage 3).In these saline brine pools the dissolved solutes will (re-)precipitate in reversed order of solubility.Gypsum crystals will form first in the brine pool and halite will be the last evaporite mineral that forms on top of the other deposits [111] resulting in a most pure efflorescent salt crust surface layer.In summary, the increasing concentration of brine leads to a horizontal zonation of evaporite minerals ordered by solubility.Such a "bulls eye" effect of lateral zonation of evaporite minerals from carbonates at the edge, through sulfates to chlorides in the topographically lower, central part are also observed in other salt pans [2,113].Influences that complicate this idealized model of evaporite mineral deposition and zonation that lead to a more complex distribution are (1) high groundwater levels during flooding events that affect the surface brine chemical composition and its saturation stage in respective to the ions of evaporite minerals [112]; (2) strong winds can move the standing waterbody on the very low-angle salt pan [1] during the brine concentration stage (stage 3) ignoring the minimal elevation differences.Subsurface brines of the Omongwa pan that were sampled in 1991 at the end of the raining season by [27] were Na-K-Cl-SO 4 dominated with a total dissolved solids concentration of about 260 g/L, which basically reflect the elemental composition of the mapped surface evaporite minerals.When the diluted meteoric water mix with these groundwater brines it lowers the potential for dissolving the surface crust evaporite minerals.
The highest change magnitude in the halite crust region detected with the multitemporal remote sensing analyses correlates very well with the regular build-up of bright efflorescent halite crust described by the model of the pan depositional cycle.During each flooding-and-desiccation cycle, the surface albedo will significantly lower after a rainfall event due to (partly) dissolved halite curst and potentially from some remaining moist sediments at the surface and will increase again, when the surface has dried and the halite crust is reformed.This means, episodic flooding events are likely to be the main driver for the observed Landsat change in surface reflectance linked to crust dynamic.The reformation of a surface salt crust during the desiccation stage is typical for the wet playas type with shallow (<5 m), fluctuating ground water table described by [9].In these types of salt pans capillary action in the sediments allows for continuous evaporation from shallow ground water that results in the formation of evaporite minerals in the capillary fringe zone as well as directly on the surface [9].This categorization is supported by the observed evaporite filled cracks and fissures in the first 50 cm of sediments reported by [27] throughout vertical profiles of the Omongwa pan.These efflorescent salt crusts that form directly from evaporation at the surface are commonly very soft and fluffy and have a high potential for dust emission even under moderate winds [10].The soft consistency of the upper ~10 cm of halite rich sediments has been also observed for the Omongwa pan [27].Over time theses soft surfaces have the potential to become crusted surfaces sufficiently thick to shut down all dust emission, through receding groundwater levels during the dry season or extended periods without rainfall [9].However, the high change magnitudes in the Landsat multitemporal analysis indicate that the halite crust is periodically dissolved and reprecipitated, thus do not support that in Omongwa long-remaining hard layered halite crusts are regularly formed.
Following the idea of [114], proportion maps of the mineralogical crust type can be used to classify the depositional environment of a salt pan.In this study the approach is adapted to the different geomorphological setting of the Omongwa pan (bordering calcite bedrock instead of clastic alluvial facies) and complemented by the topographic information and the change detection result.Figure 9 shows an exemplary South-North profile of Hyperion mineralogical mapping compared with Landsat ~30 years dynamic mapping and LiDAR topography that crosses the different depositional zones of the Omongwa pan (location of the profile is giving in Figure 7c).Six zones with three different levels of activity are recognized.The most southern surface is dominated by high proportions of gypsum crust (40%-90%) with some calcite/sepiolite crust (up to 30%) and lower gypsum content at the very beginning.This very stable gypsiferous surface represents the southern pan margin with increasing topography in direction of the lunette dune.As the profile moves north an increase in halite crust abundance (up to 70%) and a decrease in gypsum abundance (around 30%) is observed.This mixture dominated saline mud flat [114] is more active compared to the pan margin, but still more stable than the lower central depression of the salt pan.Further north the gypsum fraction in the surface sediment quickly declines and halite increases further.This very dynamic region dominated by high halite proportions (between 50% and 95%) can be interpreted as the central area of saline pan sediments [13,15].This is the aforementioned region where build-up of bright efflorescent halite crust alternates with more moist phases and lower surface albedo in each recurring inundation cycle.To the north the saline pan zone is bordered by a slight increase in topography (~10 cm) that results in conditions similar to the southern mud flat with lower surface dynamic and mineralogical mixture of halite, gypsum, but also northward increase of calcite/sepiolite.North of this elevated ridge position the halite fraction increases again and comparable condition to the central saline pan are met.This linear shaped depression in the north of the pan originates from a short, poorly developed channel that feeds into the depression at the northeastern pan margin and is active during significant rainfall events (Figure 7c).The channel was observed to hold water for several days after a rainfall in October 2014 and June 2015.The transect ends at the stable pan margin dominated by the calcite/sepiolite mudstone bedrock that seems to deposit into the adjacent halite dominated depression.Additionally, anthropogenic disturbance lead to modification of the original salt crust especially at the pan border, where the fresh efflorescent halite crust is removed and mineralogically more mixed material emerges from the subsurface that is less cohesive and can more easily contribute to dust emission.

Conclusions
In this paper, combined multitemporal and hyperspectral remote sensing analyses over the Omongwa salt pan in the Namibian Kalahari region are used to spatially differentiate and map depositional environments within the pan that deepen our knowledge of the local surface processes.Additional field and laboratory measurements and LiDAR topographical survey support the interpretation.The results show that the major salt pan mineralogical crust types (halite, gypsum, sepiolite/calcite, disturbed dark crust) could be spatially and spectrally well differentiated by Spectral Mixture Analysis on the basis of to EO-1 Hyperion imagery.The spatial distribution patterns of the mineralogical crust types correlate well with the Landsat ~30 years change map of the pan surface derived with IR-MAD analysis.A Highly dynamic pan center is identified that is dominated by halite crust, whereas at the more stable pan margins and the southern parts of the pan gypsum and calcite/sepiolite crusts are more common.The surface change map relates well to different susceptibility to pan inundation events according to the solubility of the evaporite minerals as described in conceptual models of the depositional cycle within salt pan environments.Whereas halite is more easily dissolved even by small precipitation events, gypsum and calcite are more likely to remain in the surface layer.Our analyses provide for the first time a deep understanding of the spatial distribution of the fresh efflorescent salt crust formation and the influence of pan morphology.

Conclusions
In this paper, combined multitemporal and hyperspectral remote sensing analyses over the Omongwa salt pan in the Namibian Kalahari region are used to spatially differentiate and map depositional environments within the pan that deepen our knowledge of the local surface processes.Additional field and laboratory measurements and LiDAR topographical survey support the interpretation.The results show that the major salt pan mineralogical crust types (halite, gypsum, sepiolite/calcite, disturbed dark crust) could be spatially and spectrally well differentiated by Spectral Mixture Analysis on the basis of to EO-1 Hyperion imagery.The spatial distribution patterns of the mineralogical crust types correlate well with the Landsat ~30 years change map of the pan surface derived with IR-MAD analysis.A Highly dynamic pan center is identified that is dominated by halite crust, whereas at the more stable pan margins and the southern parts of the pan gypsum and calcite/sepiolite crusts are more common.The surface change map relates well to different susceptibility to pan inundation events according to the solubility of the evaporite minerals as described in conceptual models of the depositional cycle within salt pan environments.Whereas halite is more easily dissolved even by small precipitation events, gypsum and calcite are more likely to remain in the surface layer.Our analyses provide for the first time a deep understanding of the spatial distribution of the fresh efflorescent salt crust formation and the influence of pan morphology.The change detection map also does not support the presence of long-lasting stable halite crust but rather confirm and detail the dynamic aspect of the distribution of halite crust, due to climatic and anthropogenic influences.The study also shows that the detailed morphological information provided by the LiDAR DSM could complement the analysis of surface mineralogy and dynamic and could contribute to the classification of salt pan's depositional environment related to topographical position.While dynamic halite crust dominates most of the surface area and especially the lower regions, more stable gypsum and calcite/sepiolite crust appear to concentrate in higher areas.
The results demonstrate the potential of current advanced remote sensing optical methodologies for the improvement of our knowledge on pan surface processes.To further deepen the understanding of the pan environment and its responses to climatic events, future work should focus on high temporal monitoring of hydrological parameters such as groundwater flow, precipitation and flooding events that could lead to insights into major driving factors for the observed surface change and mineralogical patterns both in the dry and in the wet seasons.Further improvement on the classification of the pan surface and its properties lies in the combination of the presented approach with RADAR based multitemporal assessments of the surface roughness.Nevertheless, in this paper, with newly available Earth observation techniques the full complexity of a saline pan dynamic could be spatially described, combining mineralogical and topographical mapping with ~30 years change detection analyses.The Omongwa pan is presented as an exemplarily object of study to test the potential of new remote sensing techniques.The approach used in this study is not pan-specific and could be well transferred and applied for the study of arid landscapes and salt pan characterization in the region and in other arid regions of the world.The proposed framework of methods is also relevant regarding upcoming remote sensing technology.For example, regarding the planned and upcoming Earth Observation missions equipped with high quality imaging spectrometers, such as the German Environmental Mapping and Analysis Program (EnMAP) mission [115], hyperspectral analyses, such as the proposed mineralogical pan surface characterization, will become available even for extensive sites on a more regular temporal basis.A better, wider characterization of pan mineralogy would also support further studies on pan susceptibility for wind erosion and contribute to assessment of potential dust emissions.

Figure 1 .
Figure 1.Location and overview of the Omongwa salt pan showing surface sample location, position of lunette dunes and extent of EO-1 Hyperion image (Basemap source: DigitalGlobe RGB image, September 2013, provided by Google Inc. (Mountain View, CA, USA)).

Figure 1 .
Figure 1.Location and overview of the Omongwa salt pan showing surface sample location, position of lunette dunes and extent of EO-1 Hyperion image (Basemap source: DigitalGlobe RGB image, September 2013, provided by Google Inc. (Mountain View, CA, USA)).

Figure 2 .
Figure 2. Different types of pan surfaces composed of detrital and evaporite mineral mixtures and environments of the Omongwa pan in October 2014 and June 2015: (a) bright efflorescent halite crust, puffy salt blisters; (b) smooth structure gypsum surface; (c) rough carbonate/sepiolite crust; (d) vegetated patches of the northeastern pan area; (e) ~50 cm high outcrop of the carbonatic mudstone bedrock, northern pan margin; and (f) surface disturbed by cattle treading at the southern pan margin; length of pen used for scale ~13 cm.

Figure 2 .
Figure 2. Different types of pan surfaces composed of detrital and evaporite mineral mixtures and environments of the Omongwa pan in October 2014 and June 2015: (a) bright efflorescent halite crust, puffy salt blisters; (b) smooth structure gypsum surface; (c) rough carbonate/sepiolite crust; (d) vegetated patches of the northeastern pan area; (e) ~50 cm high outcrop of the carbonatic mudstone bedrock, northern pan margin; and (f) surface disturbed by cattle treading at the southern pan margin; length of pen used for scale ~13 cm.

Figure 3 .
Figure 3. Field spectra (straight lines) of the three characteristic pan surface types of the Omongwa pan associated with halite crust (sample 141), gypsum crust (sample P65), and sepiolite/calcite crust (sample P68) and USGS laboratory spectra [102] (dotted lines) of the main minerals.The grey overlay indicates wavelength regions of atmospheric absorption and sensor lower signal-to-noise ratio not useable for Hyperion analyses.

Figure 3 .
Figure 3. Field spectra (straight lines) of the three characteristic pan surface types of the Omongwa pan associated with halite crust (sample 141), gypsum crust (sample P65), and sepiolite/calcite crust (sample P68) and USGS laboratory spectra [102] (dotted lines) of the main minerals.The grey overlay indicates wavelength regions of atmospheric absorption and sensor lower signal-to-noise ratio not useable for Hyperion analyses.

Figure 4 .
Figure 4. Hyperion endmember image spectra used for the spectral mixing analysis.Arrows indicate the position of diagnostic absorption features for gypsum and calcite/sepiolite crust.

Figure 5
Figure 5 show the SMA result of the identified pan crust components based on the Hyperion image.Each abundance map represents the fraction of the surface type on the Omongwa salt pan in: (a) bright halite crust; (b) gypsum crust; (c) calcite and sepiolite crust; and (d) disturbed, dark crust.In general, all endmembers have their abundances well mapped with no component exceeding 1.The accuracy of the unmixing model is assessed through the root mean squared error (Figure5e).The RMSE image shows that the overall pan surface is well explained with the linear mixture of the four selected endmembers except for higher RMSE areas at: (1) the northeast corner and western pan border associated with heavy cattle disturbances; (2) the mainly diagonal SW-NE lines crossing the pan, which are used as driving shortcuts across the pan during the dry season as well as cattle path ways; and (3) the parallel, vertical stripes (~11°), associated with remaining sensor noise that were not removed through the pre-processing procedure.In general the performance of the unmixing model is very good with an overall mean RMSE of ~1% comparable to equivalent studies, e.g.,[105] achieved an mean RMSE of ~2% using linear unmixing on Hyperion data to map salt effect soils or[106] getting the same accuracy range for the abundances of vegetation, soil and limestone bedrock.The SMA mapping shows that more than half of the pan surface is dominated by halite crust in variable amounts (Figure5a).The highest abundances are mapped in the central part and in a circular shape in the southern half of the pan, whereas the border region and the anthropogenic pathways have minimal fractions of halite crust.In most regions the distribution of gypsum crust appears to be inversely related to the halite crust.Highest gypsum composition is found at the southern pan border, in the south-central circular shape and in a ~1.5 km long stripe in the northern pan area (Figure5b).The occurrence of calcite and sepiolite crust is confined to the northern part of the pan and to the pan border (Figure5c) close to outcrops of the carbonatic bedrock that were observed in the field (Figure2e).The disturbed, dark crust endmember (Figure5d) is mainly associated with the cattle and human disturbance of the salt crust through the pathways crossing the pan and at the border regions (e.g., Figure2f).In these areas the top crust is at partly removed and the subsurface material is exposed.The subsurface is usually less concentrated in evaporite minerals and more mixed with the clastic, quartz-rich material of the pan.This is reflected by the geochemical analyses of subsurface field samples 142 (halite: 9%, quartz: 63%) and 172 (halite: 17%, quartz: 41%), directly sampled under the surface crust of samples 141 (halite: 94%, quartz: 1%) and 171 (halite: 50%, quartz: 27%), respectively at a depth at a depth of ~10 cm.The sub-surface samples show decreased halite and increased quartz content compared to surface samples (Table1).In addition to mineralogical differences, the physical properties of the salt pan surface are altered by cattle trampling, increasing the surface roughness[107], which may also cause the reflectance to decrease even at similar salt content by casting micro

Figure 4 .
Figure 4. Hyperion endmember image spectra used for the spectral mixing analysis.Arrows indicate the position of diagnostic absorption features for gypsum and calcite/sepiolite crust.

Figure 5
Figure 5 show the SMA result of the identified pan crust components based on the Hyperion image.Each abundance map represents the fraction of the surface type on the Omongwa salt pan in: (a) bright halite crust; (b) gypsum crust; (c) calcite and sepiolite crust; and (d) disturbed, dark crust.In general, all endmembers have their abundances well mapped with no component exceeding 1.The accuracy of the unmixing model is assessed through the root mean squared error (Figure5e).The RMSE image shows that the overall pan surface is well explained with the linear mixture of the four selected endmembers except for higher RMSE areas at: (1) the northeast corner and western pan border associated with heavy cattle disturbances; (2) the mainly diagonal SW-NE lines crossing the pan, which are used as driving shortcuts across the pan during the dry season as well as cattle path ways; and(3) the parallel, vertical stripes (~11 • ), associated with remaining sensor noise that were not removed through the pre-processing procedure.In general the performance of the unmixing model is very good with an overall mean RMSE of ~1% comparable to equivalent studies, e.g.,[105] achieved an mean RMSE of ~2% using linear unmixing on Hyperion data to map salt effect soils or[106] getting the same accuracy range for the abundances of vegetation, soil and limestone bedrock.The SMA mapping shows that more than half of the pan surface is dominated by halite crust in variable amounts (Figure5a).The highest abundances are mapped in the central part and in a circular shape in the southern half of the pan, whereas the border region and the anthropogenic pathways have minimal fractions of halite crust.In most regions the distribution of gypsum crust appears to be inversely related to the halite crust.Highest gypsum composition is found at the southern pan border, in the south-central circular shape and in a ~1.5 km long stripe in the northern pan area (Figure5b).The occurrence of calcite and sepiolite crust is confined to the northern part of the pan and to the pan border (Figure5c) close to outcrops of the carbonatic bedrock that were observed in the field (Figure2e).The disturbed, dark crust endmember (Figure5d) is mainly associated with the cattle and human disturbance of the salt crust through the pathways crossing the pan and at the border regions (e.g., Figure2f).In these areas the top crust is at partly removed and the subsurface material is exposed.The subsurface is usually less concentrated in evaporite minerals and more mixed with the clastic, quartz-rich material of the pan.This is reflected by the geochemical analyses of subsurface field samples 142 (halite: 9%, quartz: 63%) and 172 (halite: 17%, quartz: 41%), directly sampled under the surface crust of samples 141 (halite: 94%, quartz: 1%) and 171 (halite: 50%, quartz: 27%), respectively at a depth at a depth of ~10 cm.The sub-surface samples show decreased halite and increased quartz content compared to surface

Figure 8 .
Figure 8. Relationship of pan surface elevation and crust endmember abundance.Box A marks high halite content at low elevations and box B marks high gypsum content at high elevations.

Figure 8 .
Figure 8. Relationship of pan surface elevation and crust endmember abundance.Box A marks high halite content at low elevations and box B marks high gypsum content at high elevations.

Table 1 .
Properties of Omongwa surface samples.Locations of field samples are mapped in Figures1 and 7b.