Satellite Advanced Spaceborne Thermal Emission and Reflection Radiometer Mineral Maps of Australia Unmixed of Their Green and Dry Vegetation Components: Implications for Mapping (Paleo) Sediment Erosion–Transport–Deposition Processes

: The 2012 satellite ASTER geoscience maps of Australia were designed to provide public, web-accessible, and spatially comprehensive surface mineralogy for improved mapping and solutions to geoscience challenges. However, a number of the 2012 products were clearly compromised by variable green and/or dry vegetation cover. Here, we show a strategy to first estimate and then unmix the contributions of both these vegetation components to leave, as residual, the target surface mineralogy. The success of this unmixing process is validated by (i) visual suppression/removal of the regional climate and/or local fire-scar vegetation patterns; and (ii) pixel values more closely matching field sample data. In this process, we also found that the 2012 spectral indices used to gauge the AlOH content, AlOH composition, and water content can be improved. The updated (new indices and vegetation unmixed) maps reveal new geoscience information, including: (i) regional “wet” and “dry” zones that appear to express “deep” geological characters often expressed through thick regolith cover, with one zone over the Yilgarn Craton spatially anti-correlated with Archaean gold deposits; (ii) a ~1000 km wide circular feature over the Lake Eyre region defined by a rim of abundant “muscovite” that appears to coincide with opal deposits; (iii) a N–S zonation across the western half of the continent defined by abundant muscovite in the south and kaolinite in the north, which appears to reflect opposing E ↔ W aeolian sediment transport directions across the high-pressure belt; (iv) various paleo-drainage networks, including those over aeolian sand covered the “lowlands” of the Canning Basin, which are characterized by low AlOH content, as well as those over eroding “uplands”, such as the Yilgarn Craton, which have complicated compositional patterns; and (v) a chronological history of Miocene barrier shorelines, back-beach lagoons, and alluvial fans across the Eucla Basin, which, to date, had proved elusive to map using other techniques, with potential implications for heavy mineral sand exploration. Here, we explore the latter three issues.


Introduction
The launch of the ASTER (Advanced Spaceborne Thermal Emission and Reflection Radiometer) satellite imaging sensor [1] in 1999 and its subsequent successful program of global mapping the Earth's entire land surface (~six repeats) <70 • latitude using its "mineral-tuned", multi-spectral VNIR-SWIR-TIR (visible to near infrared; shortwave infrared; thermal infrared) sensors presents an opportunity to map a diverse suite of processcritical minerals (e.g., rock type, metamorphic alteration, weathering, and transported cover).This "mineral mapping" is based on the spectroscopic measurement of electromagnetic energy and its interaction with minerals [2], which have diagnostic, wavelengthdependent absorption and/or reststrahlen features.For example, absorption at 2200 nm is caused by the (polarized) vibration of di-octahedrally-coordinated Al-OH bonds in phyllosilicates like muscovite [3] and kaolinite [4].This 2200 nm absorption is specifically targeted by ASTER band 6.
The first continental-scale ASTER mineral maps were publicly released in 2012 for Australia [5].These geoscience maps are based on spectral band ratios, which target a range of mineral-diagnostic features across ASTER's three wavelength modules.These modules include: (i) three bands in VNIR (15 m pixel resolution); (ii) six bands in the SWIR (30 m pixel resolution); and (iii) five bands in the TIR (90 m pixel resolution).The 2012 ASTER maps of Australia have since been investigated for a range of geoscience applications [6][7][8][9][10][11][12][13][14][15], though their value has arguably been limited by both the use of a heavy green vegetation mask [5], which culled >40% of all pixels (100% in more wet, vegetated regions), with a number of products especially over arid regions showing the effects of variable dry vegetation cover [6].Therefore, a new strategy for correcting vegetation is required to improve the geoscience value of ASTER imagery [7].
Estimating the green vegetation content is relatively simple using ASTER as its VNIR bands target specific plant properties, like chlorophyll absorption (ASTER bands 1 and 2) and the NIR photosynthetically-active reflectance plateau (ASTER band 3) [16].However, ASTER was not designed to measure dry vegetation features, like cellulose absorption at 2080 nm or lignin absorption at 2270 nm [17].This deficiency in sensor design is common to all satellite multi-spectral systems launched since ERTS-1 in 1972 [16].This despite dry vegetation being the dominant vegetation cover type across much of the Earth's drier regions as well as during seasonal senescence in more temperate latitudes.Thus, accurate estimations of the dry vegetation component from multispectral data are (and will remain) problematic [7,[18][19][20][21][22].
Here, we test dry vegetation indices that use different wavelength modules of ASTER, namely: (i) an NIR-SWIR index [22], which targets the shoulder of the 2080 nm cellulose and 2270 nm lignin features, though these also overlap with mineral absorptions (e.g., chlorite); and (ii) a TIR index [22], which targets a broad cellulose feature centered at 10,500 nm [17,23], though this also overlaps with mineral reststrahlen features (e.g., quartz).That is, both indices are likely to yield residual mineral-related artifacts in certain geological environments.These ASTER dry vegetation indices, as well as a green vegetation index, are then used in an unmixing strategy [7,21,[24][25][26][27][28][29][30][31][32][33] to unmix the Australian ASTER reflectance/emissivity mosaic [34], which comprises >3000 cross-calibrated ASTER images that were collected between 2000 and 2007 [35].We also reassess a number of indices from the 2012 Australian ASTER geoscience product suite [5] to determine whether they are optimized for mapping mineral-sensitive (paleo) erosion-transport-deposition processes.We validate these results by (i) observing how well vegetation effects evident as fire scars and regional rainfall patterns are suppressed/removed at different scales; and (ii) comparing these ASTER indices with laboratory spectral and physiochemical measurements of surface samples collected as part of the National Geochemical Survey of Australia (NGSA) [36][37][38].
Along this journey, new geoscience information is revealed using these updated ASTER mineral maps about the geology and mineralized systems across Australia, including: (i) Archaean hydrothermal gold across the Yilgarn Craton; (ii) opal (silica) development across central Australia; and (iii) natural, REDOX-related hydrogen degassing [39].These are to be tackled, ideally, in future work.But here we focus on the networks of (paleo) rivers and marine basins that were active during the Cenozoic, when sea levels and rainfall were higher [40][41][42][43][44][45][46][47][48][49][50].Some of these paleo-systems are now the focus of heavy mineral sand and valley-calcrete-uranium exploration [43,44,48].
To date, paleo-valleys have been primarily mapped using digital elevation models (DEMs) [51] and satellite thermal infrared imagery (e.g., NOAA-AVRR and ASTER) [41,50], with the latter sensing surface kinetic temperature patterns best observed during stable (no wind) night-time conditions [41,50].Other "deeper" geophysical sensing techniques have also been tested, such as airborne electromagnetics, targeting saline groundwater [52].However, surface sensing techniques, like airborne gamma radiometrics [53] and the 2012 ASTER products [5], have been less attractive because of the perception that the surfacedetected composition is too complex to unravel [54].We instead argue that accurate surface mineral maps that are diagnostic of erosion-transport-deposition processes will not only assist geoscientists in unravelling the surface regolith complexity but also reveal new and valuable geoscience information/knowledge not obtainable from any other mapping data type.
Here, we focus on improving the accuracy/value of a selection of the 2012 ASTER mineral products [5] using this vegetation unmixing strategy.These "regolith" products include: (i) the "silica index" [55,56], which is largely driven by the content of quartz sand; (ii) "AlOH content", which is a gauge of the amount of "weathering-related", AlOH-bearing "clay" minerals [57]; (ii) "AlOH composition", which is sensitive to the species of "clay" minerals, including illite (~muscovite), montmorillonite, and kaolinite; and (iii) water content or "soil moisture" [6], which can also be driven by the type of clay mineral, like montmorillonite = "wet" while kaolinite = "dry" [6].Not included here but also likely to assist in "unravelling" the complexity of the regolith are other ASTER products [5] such as REDOX-related ferric oxyhydroxides and ferrous oxides/silicates/carbonates, especially given their potential importance for hydrogen exploration [39].
Based on the updated ASTER mineral maps, we therefore asked the following: (i) are (paleo-) rivers, shorelines, and fans/deltas better expressed and/or is new geoscience information revealed; (ii) does this compositional information provide clues regarding understanding Cenozoic sediment transport systems, especially sediment provenance, depositional chronology, the paleo-climate, flow direction, energy of depositional environments, and tectonics/the climate; and (iii) does this information create new opportunities for mineral exploration?
In summary, the primary objective of this paper is (i) to describe a vegetation unmixing method that compensates for both the green components and, especially, the dry vegetation components using the ASTER's limited spectral band configuration to generate more accurate maps of surface mineral information (compared with [5]); and (ii) to demonstrate the efficacy of the updated mineral maps for mapping paleo-geomorphic systems, especially those across more arid regions of Australia obscured by dry vegetation cover.

ASTER Data Pre-Processing
The 2012 Australian cross-calibrated ASTER reflectance (VNIR and SWIR bands at 30 m pixel resolution) and emissivity (TIR bands at 90 m pixel resolution) mosaics [34] were processed using the software package, Environment for Visualizing Images (ENVI) (Version 3.5), to generate a single mosaic combining all 14 VNIR, SWIR, and TIR bands resampled to a 100 m pixel size (nearest neighbor sampling) to allow for vegetation unmixing across all three ASTER wavelength modules.Essential to the task of isolating and enhancing compositional information is the process of pixel normalization, such as band ratios [58], which cancels the often-dominant, multiplicative, wavelength-independent effects, such as topographic illumination, surface scattering, and albedo [59,60].This step is only effective when all of the additive effects, like atmospheric scattering (in the VNIR) and instrument error (e.g., the ASTER SWIR cross-talk effect) [61], have first been corrected.Given this, then, the resultant normalized images should appear "flat" over variable topographic illumination, with areas in deep shade (and water in the SWIR) appearing as "salt-andpepper" (shot) noise.The additive error assessment and correction step was implemented as part of the 2012 Australian ASTER reflectance mosaic generation [5].

Vegetation Cover Issue
The vegetation cover challenge [62] is demonstrated in Figure 1, where the published paleo-geomorphic features [47], like rivers (cyan lines), shorelines and lagoons (yellow lines), are obscured by a patchwork of variable vegetation cover.The 2012 ASTER "AlOH" content (Figure 1b) and "silica index" (Figure 1c) products [5] are especially obscured by the dry vegetation component, which is highlighted by the edges of the fire scars (red lines), while the 2012 ASTER "iron oxide content" product (Figure 1d) is instead affected by green vegetation (green lines).That is, the green and dry vegetation cover not only compromises the mapping of paleo-geomorphic elements (Figure 1b-d) but also disrupts the integrity of the pixel spectral mineral signatures.

Vegetation Cover Issue
The vegetation cover challenge [62] is demonstrated in Figure 1, where the publishe paleo-geomorphic features [47], like rivers (cyan lines), shorelines and lagoons (yellow lines), are obscured by a patchwork of variable vegetation cover.The 2012 ASTER "AlOH content (Figure 1b) and "silica index" (Figure 1c) products [5] are especially obscured b the dry vegetation component, which is highlighted by the edges of the fire scars (re lines), while the 2012 ASTER "iron oxide content" product (Figure 1d) is instead affecte by green vegetation (green lines).That is, the green and dry vegetation cover not onl compromises the mapping of paleo-geomorphic elements (Figure 1b-d) but also disrupt the integrity of the pixel spectral mineral signatures.

ASTER Vegetation Unmixing Process
Three ASTER vegetation indices were tested using ENVI (Version 3.5), including: (i) A green vegetation index: This index is based on the normalized difference vegetation index [16], which target photosynthetic activity and cancels cross-track surface scattering variation.

ASTER Vegetation Unmixing Process
Three ASTER vegetation indices were tested using ENVI (Version 3.5), including: (i) A green vegetation index: This index is based on the normalized difference vegetation index [16], which targets photosynthetic activity and cancels cross-track surface scattering variation.
(ii) A SWIR-derived dry-vegetation index [21]: This dry vegetation index targets the cellulose and lignin absorption rates (ASTER bands 5 and 7, respectively), as well as a lack of the photosynthesis-active NIR plateau [17] (ASTER bands 3 and 4).It is also more comprehensive than the ASTER dry vegetation index of [18], which used a ratio of ASTER bands 6 and 7 and, thus, did not consider the impact of "clay" mineral-related AlOH absorption.
(iii) A TIR-derived dry-vegetation index [22]: This TIR-based index (Equation (3) targets a broad reflectance peak (emissivity low) for dry vegetation at wavelengths > 9200 nm (centered at 10,500 nm) (ASTER bands 13 and 14) [17,22,23,63].ASTER band 12, centered at 9100 nm, was omitted because it is where clay minerals typically have their reflectance peak (i.e., minimize an inverse clay silicate mineral response).Band 10 is centered at 8250 nm, positioned over a quartz reflectance peak, while band 11 (8625 nm) in the denominator helps to minimize the contribution of quartz.
The vegetation unmixing process begins by assessing for suitable mixing behaviors.This can be achieved using a 2-dimensional (2D) scattergram, with the vegetation index (Equations ( 1)-( 3)) plotted against the target mineral index.An "ideal" mixing pattern is where the data cloud approximates a right triangle, with its perpendicular corner positioned at zero-zero.Given this geometry for both green vegetation (V green ) and dry vegetation (V dry ) with the target mineral index (T), the process of unmixing for a given pixel "i" can be subsequently described as follows: The gains "a" and "b" were iteratively determined using the image data by gauging how well they suppress vegetation-related "edge" effects, such as fire scars (Figure 1), with the assumption being that the surface mineral signature across these "edges" is more slowly changing.A suitable gain (i) suppresses these edge effects across a wide range of (all) geological surface compositions for all areas, and (ii) generates a dimensionless 2D data cloud.Note that we first scaled all input spectral indices between 0 and 1 for simplicity.Though not included in Equation (4), there is also wavelength-dependency involved in the process.For example, green vegetation can affect minerals with diagnostic VNIR spectral features, such as iron oxides (Figure 1d), whereas dry vegetation can affect SWIR-and TIR-active minerals, like the 2200 nm AlOH absorbers (Figure 1b) and silicates (Figure 1c), respectively; that is, those wavelengths where vegetation is "aspectral" (transparent and with no spectral features) require no correction.

NGSA Validation Data
The NGSA samples used for the validation of the ASTER products comprise >1000 unconsolidated, fluvial-channel (including paleo-channels) overbank samples taken from the top 10 cm of the regolith of (paleo) river catchments across Australia [36].A range of laboratory measurements were taken from these samples, including geochemistry, particle size distribution, and VNIR-SWIR-TIR reflectance measurements.The reflectance data were collected using an Analytical Spectral Device (ASD) FieldspecPro-3 bi-directional reflectance spectrometer and a Bruker Vertex 70 directional-hemispherical spectrometer [38].These spectral data were then convolved to ASTER bandpass responses using ENVI image processing software (Version 5.3).
Three high spectral resolution, NGSA absorption band-depth indices [58] were used for the assessment/validation of ASTER indices/products, namely: Note here the use of the lower-case "b" mnemonic for the ASD indices in contrast to the upper-case "B" used for the ASTER indices (Equations ( 1)-( 3)).Statistical assessment (Table 1) for the same 165 NGSA samples used for the 2016 study [6].Correlations here are deemed significant at the 90% confidence level if the regression coefficient of correlation (R 2 ) exceeds 0.35.All relationships were tested using linear functions, though some R 2 values would have been improved if non-linear functions were fitted instead.Given that at ASTER spectral resolution only two 2165-2200 nm absorbing mineral components are measurable (degrees of freedom), we described these endmembers as being either "muscovite" or "kaolinite".The "muscovite" endmember can also include minerals like paragonite, phengite, lepidolite, brammallite, illite, and montmorillonite, while the "kaolinite" endmember can also include minerals like dickite, nacrite, halloysite, pyrophyllite, and alunite.

NGSA-ASTER Mineral Indices
Table 1 presents the R 2 results of the spectral (ASD and ASD-ASTER convolved), chemical, and particle-size parameters for the 165 randomly selected NGSA samples.The seven "target" parameters include: 1.

2.
The geochemical indices: • Al/(Al + Si), which is an indicator of the proportion of Al-bearing clay minerals to coarser quartz grains; • Al/(Al + Ca + Mg + Na + K), i.e., based on the chemical index of alteration called the CIA [65]; • Loss-on-ignition (LOI), which is a measure of the volatiles released after heating, including water, organic compounds, and carbonates [66].

NGSA-Water Content
Unlike the ASTER's VNIR bands, all of the ASTER's SWIR bands show inverse relationships (better modelled using an exponential decay function) with the ASD 1900 nm water absorption depth (Equation ( 7)), with correlations increasing with the ASTER's SWIR band number.That is, ASTER band 9 has the highest correlation (R 2 = 0.39), followed by ASTER band 8 (R 2 = 32).This is consistent with the "shoulder" effect to the major, fundamental stretching vibration of molecular water at ~2700 nm [64,[68][69][70].Thus, ratios with ASTER band 9 in the denominator generate higher correlations; for example: . The latter two also yield significant correlations with the 2165 and 2207 nm absorption depths (Equations ( 5) and ( 6)), respectively).Replacing ASTER band 9 with band 8, to simulate the WV-3 satellite sensor data, also yields a significant correlation with the 1900 nm water absorption depth (Table 1).

NGSA-AlOH Clay Content
The 2012 ASTER's AlOH content ratio (B 5 + B 7 )/B 6 is not correlated (weakly anticorrelated) with the ASD 2200 nm absorption depth (Equation ( 5)) (Table 1).ASTER indices that do show significant correlations include (i) B 7 /B 6 (R 2 = 0.69) and (ii) B 7 /B 8 (R 2 = 0.62).There is no correlation with respect to the B 7 /B 5 ratio (R 2 = 0.02), which means that ASTER band 5 should not be used for gauging the 2200 nm AlOH absorption depth.Using ASTER band 8 for this purpose is also problematic, as a range of minerals (for example, chlorite, jarosite, and gibbsite), which are not well represented in the NGSA sample suite, also absorb in the ~2250 nm wavelength region.Thus, we recommend using the ASTER B 7 /B 6 ratio for gauging the 2200 nm AlOH absorption depth.

NGSA-AlOH Clay Composition
No ASTER or ASD spectral measure was found to correlate with the CIA geochemical index (Al 2 O 3 /(Al 2 O 3 + CaO* + Na 2 O + K 2 O) × 100), including the ASD 2165 nm absorption depth (Equation ( 6)) (R 2 = 0.08; Table 1).Given that the 2165 nm absorption in these NGSA samples is caused by kaolin, we thereby conclude that the CIA geochemical index is not a reliable gauge of kaolin content.Furthermore, common regolith minerals, like calcite (CaCO 3 ) and muscovite (KAl 2 (AlSi 3 O 10 )(F,OH) 2 ) [71], would adversely affect any clay mineral interpretation using the CIA.
absorb in the ~2250 nm wavelength region.Thus, we recommend using the ASTER B7/ ratio for gauging the 2200 nm AlOH absorption depth.
Finally, the inverse relationship between the 2165 nm absorption depth (Equation ( and the 1900 nm water absorption depth (Equation ( 7)) (R 2 = 0.37, Table 1) is consiste with kaolinite having low affinity with water, unlike illite and, to a greater extent, mon morillonite, which absorb water molecules in their interlayers [73].This relationship h previously been used to help separate illite (and muscovite) from montmorillonite [6].Directed principal components [72] could help isolate the minor AlOH composition component from the dominant AlOH content component, though we also suggest using a color image (R:G:B) combination of red: B 7 /B 6 , green: B 7 /B 5 , and blue: B 7 /B 6 , with the intensity depicting the AlOH content and the hue (color) expressing the target AlOH composition, with kaolinite-rich pixels generating green tones and muscovite-rich pixels generating magenta tones.
Finally, the inverse relationship between the 2165 nm absorption depth (Equation ( 6)) and the 1900 nm water absorption depth (Equation ( 7)) (R 2 = 0.37, Table 1) is consistent with kaolinite having low affinity with water, unlike illite and, to a greater extent, mont-morillonite, which absorb water molecules in their interlayers [73].This relationship has previously been used to help separate illite (and muscovite) from montmorillonite [6].

Vegetation Unmixing Process
Figure 3 presents a suite of 2D scattergrams of various indices collected from the entire Australian ASTER dataset.The spatial coherency of these data clouds is a testament to the success of the L1T calibration [74] and cross-calibration [35] of >3000 ASTER images.These 2D scattergrams also show a variety of geometries.One type approaches an amorphous "blob" (Figure 3a,b,g,h), which indicates no apparent relationship between the input variables.The first two involve the water content ((B 5 + B 8 )/(B 7 + B 9 )), with the lack of a relationship with vegetation cover also evident in the associated imagery (Figure 4i).Thus, no vegetation unmixing is required for the ASTER's water content product.Having said this, the 2D scattergram of water content versus TIR-derived dry vegetation (Figure 3c) shows a "T" pattern involving two overlapping populations, neither of which is apparently correlated with the water content.
Even more complex 2D geometries are expressed by the AlOH content versus the SWIR-derived and TIR-derived dry vegetation indices (Figure 3e and Figure 3f, respectively).Nevertheless, using Equation ( 4) to unmix the AlOH content using both of these dry vegetation indices results in (i) the SWIR-derived dry vegetation data cloud collapsing to a dimensionless blob (Figure 3g), while the TIR-derived dry vegetation data cloud increases its spread (in the direction of the white arrows in Figure 3i), and (ii) the associated unmixed images show suppression of fire scars (when comparing Figure 4d with Figure 4e,f).This indicates the "dry vegetation" indices using either the SWIR and TIR wavelength regions (Equations ( 2) and (3), respectively) are modeling essentially the same information.Thus, either dry vegetation algorithm (Equations ( 2) and ( 3)) can be used for the unmixing process (Equation ( 4)), though we suggest using the SWIR-derived index (Equation ( 2)) simply because it can be implemented at a 30 m pixel resolution in comparison to the 90 m pixel resolution of the TIR index.Note also the narrow zones of elevated response (black/white arrows in Figure 4), which are likely caused by residual mineral information in the dry vegetation indices (Figure 4b,c).These result in a neutral-colored (~white) response in the AlOH composition map (Figure 4h).
A comparison of Figure 4g versus Figure 4h shows the impact of vegetation unmixing on the AlOH composition product; that is, a coherent geological pattern is revealed after unmixing (Figure 4h), including most/all of the published paleo-geomorphic elements as well as clear demarcation of muscovite-rich "uplands" in the center and northwest and kaolinite-rich "lowlands" in the southeast, none of which is apparent in the nonunmixed products (Figure 4g).However, we emphasize the importance of gauging the appropriate gain levels (a and b in Equation ( 4)), which, in this case, were determined using visual criteria (suppression of fire scars) from only the area shown in Figure 4 and not at the continental-scale, which was the case for the other "unmixed" images shown in this paper; that is, the magnitude of these user-defined gains will impact the information clarity/accuracy of the unmixed products, so careful iteration and validation are critical.
a relationship with vegetation cover also evident in the associated imagery (Figure 4i).Thus, no vegetation unmixing is required for the ASTER's water content product.Having said this, the 2D scattergram of water content versus TIR-derived dry vegetation (Figure 3c) shows a "T" pattern involving two overlapping populations, neither of which is apparently correlated with the water content.The AlOH content (Figures 3d,e and 4d) shows very different behavior to the w content (Figures 3a,b and 4i).The 2D scattergram of AlOH content versus green vegeta content (Figure 3c) approximates a right triangle (highlighted by a yellow line in Fig 3d ), with its three apices representing (i) 100% green vegetation (0,1); (ii) 100% AlOH m eral content (1,0); and (iii) 0% of both (0,0).Given this "ideal" mixing geometry, the step in the unmixing process (Equation ( 4)) involves transforming each pixel's estima of its green vegetation content (examples shown by the white vectors) onto the AlOH tent axis (Equation ( 4)).This transformation yields, as expected, an amorphous "b

The Australian Continent ASTER Results
The impact of unmixing at the continental-scale is also pronounced.For example, unmixing the AlOH content (compare Figure 5a,b) shows overall increase across the >250 mm (magenta line) annual rainfall zone, which results in better agreement with the associated NGSA data (colored dots), e.g., areas like "B" and "C".The unmixing of the silica index (Figure 5e) yields similar increases across the >250 mm zone, including, for example, "D", which is associated with deep sandy soils and sand dunes of the Little Desert, which was not apparent in the original ASTER geoscience maps [5,6].In arid areas (<250 mm rainfall), the AlOH content is reduced along the tract of the Darling River (highlighted by small yellow arrows) and for large parts of the Lake Eyre ("A") and Canning Basins (Figure 5b) but is elevated in the southern Yilgarn Craton (cyan polygon) and the rim of the Lake Eyre Basin ("A"), which coincides with opal deposits (red triangles).
The unmixed AlOH composition (Figure 5c) generates a different pattern compared with the original 2012 product [5,6]; for example (i) the northeastern Yilgarn Craton (dashed cyan polygon), Arunta-Musgrave complex (dashed white polygon), and Gawler Block (dashed yellow polygon), as well as the regolith between them, were originally mapped as dominantly kaolinite-rich but are now relatively muscovite-rich; (ii) the Eucla Basin, which was originally mapped as muscovite-rich [5], is now kaolinite-rich; (iii) the central part of the Lake Eyre basin ("A") is more kaolinite-rich; and (iv) the <250 mm rainfall zone is now largely muscovite-rich, while the >250 mm rainfall zone remains kaolinite-rich.The latter may reflect the high solubility/mobility of potassium [75].Importantly, the NGSA kaolinite-muscovite spectral index (color-coded dots in Figure 5) better matches the updated results.
After unmixing, a N-S zonation becomes apparent especially in the western half of the continent across ~25 • S latitude, with kaolinite developed in the north and abundant muscovite in the south (Figure 5b,d).This zonation is consistent with (i) the climatedriven high-pressure belt and associated prevailing wind pattern (yellow open arrows) [76]; (ii) relatively thin crust under the Canning and Lake Eyre Basins [77]; and (iii) the proposed Neogene "tilt axis" (black dashed line in Figure 5b), which is based on Cenozoic marine sediments being 250 m above the present-day sea level in the south while they are at present-day sea level in the north [78].The vegetation unmixing also appears to better define the "Menzies Line" (yellow line in Figure 5b,d,f), which marks a transition from alkaline, fresh groundwater in the north (muscovite-rich) to acid, saline groundwater in the south (kaolinite-rich) [75,79,80].
The water content product (Figure 5f), which did not require vegetation unmixing, shows similarity with the NGSA 1900 nm absorption depth (similar color-coded dots), except for "wet" latitudes north of ~19 • S, with some individual satellite tracks standing out, e.g., north of the Canning Basin.This seasonally "wet" region is likely caused by the timing of ASTER images being acquired during the wet-tropical summer to maximize the high sun angle.Across the arid zone (<250 mm), the ASTER water content largely follows the ASTER AlOH products (when comparing Figure 5f with Figure 5b,d), including a "wet" kaolinite-rich Eucla Basin and "dry" muscovite-rich regions of the Yilgarn, Gawler, and Arunta-Musgrave terranes.The Lake Eyre region ("A"), however, shows departure, e.g., "wet" muscovite east of "A", which likely reflects the development of a hygroscopic clay like montmorillonite instead of "muscovite" [6].
the "Menzies Line" (yellow line in Figure 5b,d,f), which marks a transition from alkaline, fresh groundwater in the north (muscovite-rich) to acid, saline groundwater in the south (kaolinite-rich) [75,79,80].  1) masked for green vegetation only.All maps have been masked to remove to the top ~10% of pixels with abundant green.NGSA validation data embedded as colored dots.The magenta line is the annual 250 mm rainfall isohyet [48].The dashed, cyan, white, and yellow polygons show the Yilgarn, Arunta-Musgrave, and Gawler Blocks, respectively.Yellow arrows highlight the Darling River.The yellow line marks the "Menzies Line" [79].Red triangles show the major opal deposit fields.Red dots (with labels) are sites of the ocean drilling program  1) masked for green vegetation only.All maps have been masked to remove to the top ~10% of pixels with abundant green.NGSA validation data embedded as colored dots.The magenta line is the annual 250 mm rainfall isohyet [48].The dashed, cyan, white, and yellow polygons show the Yilgarn, Arunta-Musgrave, and Gawler Blocks, respectively.Yellow arrows highlight the Darling River.The yellow line marks the "Menzies Line" [79].Red triangles show the major opal deposit fields.Red dots (with labels) are sites of the ocean drilling program (ODP) [76].The black dashed line is a proposed Neogene "approximate tilt axis" [80].A to D are locations referred to in the body text.

Paleo-Valleys of the Canning Basin Region
The Paleozoic Canning Basin [81,82] in northern Western Australia (Figures 5f and 6) was superimposed by a network of paleo-valleys [41,47,83] (Figure 6a) before superposition by aeolian quartz sands of the Great Sandy, Little Sandy, Tanami, and Gibson Deserts.This quartz sand cover is well-mapped using the (vegetation unmixed) ASTER silica index (warmer colors in Figure 6b).The present land surface is also obscured by a patina of fire scars, which is evident in both the ASTER false color composite (Figure 6a) and nonunmixed AlOH content (Figure 6c) products.After vegetation unmixing, the AlOH content (Figure 6d) clearly reveals this paleo-river network (Figure 6a), including fine-scale paleovalleys, such as those near "E".Many of these paleo-valleys are mapped "under" aeolian sand cover (e.g., "G" in Figure 6b), being characterized with a low AlOH content, similar to the present-day Darling River (yellow arrows in Figure 5b).The possible causes of this low AlOH content associated with river channels (Figures 5b and 6d) include (i) mobile, "platy" AlOH-bearing minerals (e.g., muscovite) being preferentially removed by (paleo) flood events; and/or (ii) AlOH-poor sand and/or heavy mineral grains preferentially accumulating in valley lows and effecting the optical depth of the spectral measurement [67,84,85].The reason why paleo-valleys are mappable using the ASTER across the entire Canning Basin, irrespective of the aeolian sand cover, may be related to the role of the swales between the sand dune ridges exposing underlying AlOH mineralogy.

Paleo-Valleys of the Eastern Yilgarn Craton Region
The Cenozoic paleo-valleys in the southeast region of Western Australia (Figures 5f and 7b) once transported sediment that was eroded from a variety of bedrock geologies (Figure 7a) [86], including (i) Archaean granite-greenstone belts of the Yilgarn Block (cyan polygon in Figure 6); (ii) Archaean-to-Proterozoic high-grade metamorphic rocks of the Albany-Fraser Orogen; (iii) a series of Proterozoic to Cambrian basins of largely siliciclastic sediments; (iv) Carboniferous to Permian glaciogenic sediments of the Canning Basin; and (v) Cretaceous siltstones and sandstones of the Gunbarrel Basin (Figure 7a).The Neoproterozoic to Late Devonian Officer Basin underlies the Carboniferous to Cretaceous sequences in this area [87].Today, the paleo-rivers are expressed over low relief (~600 m vertical change over ~600 km distance) (Figure 7b) [88,89] and include: (i) eastward flowing systems like the Carey, Raeside, Yindarlgooda, and Lefroy Paleorivers, with catchments spanning Archaean rocks of the Yilgarn Block and Proterozoic rocks of the Albany-Fraser Orogen; and (ii) southward flowing systems like the Throssel and Baker Paleorivers, with catchments spanning Phanerozoic rocks to the north.That is, this region contrasts with the Canning Basin area (Figure 6) as it comprises diverse, "upland" (exposed and erosional) geological environments.
Similar to the Canning Basin results (Figure 6c,d), vegetation unmixing is effective in suppressing fire scars, as shown by a comparison of the before/after AlOH content products (Figure 7d,e).However, this unmixing reveals a change from higher AlOH contents and muscovite-rich in the south to lower AlOH contents and kaolinite-rich in the north (Figure 7e,g), consistent with the continental-scale results (Figure 5b,d).This change occurs over a narrow (~100 km wide) zone (white dashed line in Figure 7d,g) that crosscuts the Phanerozoic geology (Figure 7a) and landform (Figure 7b) and is not apparent in any other ASTER (Figure 7f,h) or geophysical products (Figure 7c); that is, it is specifically an AlOH mineral signature.
The ASTER AlOH composition map (Figure 7g) shows the surface of the Yilgarn Block (Figure 7a) to be muscovite-rich, which is consistent with the high K response in the airborne gamma radiometric image (Figure 7c).However, this compositional association is not universal.For example, a pronounced high K radiometric response of the Eucla Basin, as well as smaller K-rich areas over "upland" rocks (highlighted by dashed-yellow polygons in Figure 7c), are relatively kaolinitic in the ASTER AlOH composition product (Figure 7c,g).This apparent inconsistency between sensing techniques can be explained by (i) the K-bearing minerals being feldspars, such as microcline (KAlSi 3 O 8 ) or orthoclase (KAlSi 3 O 8 ), and not an AlOH-bearing phyllosilicate, like muscovite; or (ii) the different sensing penetration depths.That is, gamma rays can persist through a meter of dry, porous, bare ground [90], while "optical" sensing penetrates the top few microns [84,91].For example, a thin kaolinite-rich weathering rind developed on a microcline-bearing rock.
In contrast with the Canning Basin (Figure 5d), the composition of paleo-drainage across this diverse, (paleo-)erosional geological environment (Figure 7) is expressed as low AlOH content for only some, intermittent segments of paleo-valleys and only for some types of exposed bedrock, such as those parts arrowed off "K" over the Yilgarn Craton (Figure 7e).However, even over this bedrock type, the surface composition of paleo-valleys alternates (10-100 km) between (i) low AlOH content, kaolinite-rich, low quartz sand content, high water content, high opaque content (Figure 7e, 7f, 7g, 7h, and 7i, respectively), and high ASTER gypsum content [86], for example, "K" and "L", and (ii) high AlOH content, high quartz sand content, muscovite-rich, low water content, low ASTER opaque content [86] (Figure 7e, 7f, 7g, 7h, and 7i, respectively), and low ASTER gypsum content [86]; for example, "H".The latter type seamlessly merges with the surrounding "upslope" geology/regolith.This alternating pattern indicates that a number of factors have driven the surface composition of paleo-valleys.One possibility is groundwater condition, especially as this area represents the transition across the "Menzies Line" [75,79]; that is, (i) above is related to shallow, acidic groundwater, while (ii) above is related to deeper and/or neutral groundwater.Unfortunately, there are too few groundwater sample points [92,93] to corroborate an association between groundwater chemistry and surface mineral composition.
Over the Carboniferous-Permian and Cretaceous rocks (Figure 7a), paleo-valleys are more consistently richer in quartz sand content (Figure 7f), though this signature merges/broadens downslope into similarly high responses over the aeolian sand plains/dunes of the Gibson and Great Victoria Deserts (Figure 7b,f).The pattern of quartz-sand-filled paleo-valleys over the Phanerozoic rocks both mirrors and helps explain the low gamma-radiometric response (Figure 7c).
The alternating compositional pattern of paleo-valleys across the Yilgarn Craton (Figure 7) terminates against a ~150 km wide, 600 km long, NW-SE-trending "dry" zone (red polygon in Figure 7g).Note that this "dry" zone is not an atmospheric effect nor caused by recent rainfall because it is measurable over dozens of ASTER scenes collected over many years, and most importantly, it aligns with the surface mineralogy, including a lack of opaque minerals (Figure 7i) and gypsum [86], as well as a dearth of hydrothermal gold deposits (Figure 7i).Down-stream from this "dry" zone, the paleo-valleys return to being "wet", AlOH-poor, kaolinitic, gypsiferous, and rich in opaque minerals (e.g., "N").
This "dry" zone (red polygon in Figure 7h) also coincides with a zone of higher deformation (i.e., tectonic thrust stacking) and increased granite intrusions, as interpreted from the 2001 seismic (01AGS-NY1 black dashed line in Figure 7i) and gravity data [94].There is also a related paucity of gold-associated greenstone keels [95] (green units in Figure 7a).There is another, albeit smaller/weaker, "dry" zone further to the east (red dashed line in Figure 7h), which, from the seismic and (11GA-YO1) drill-hole (DH) data [95] (DH) [95,96] (magenta triangles in Figure 7i), corresponds to a rapid, west-to-east thickening of the Phanerozoic cover stratigraphy from the Empress DH (~500-1500 m) to Yowalga DH (>4100 m depth).The thinner western part of this zone is developed over a strongly reflective middle crust and a gravity high [95].Thus, the surface "dry" zones mapped using the ASTER have deeper crustal associations that appear to persist through sedimentary/regolith cover.Ideally, these ASTER results should be combined with drill-core data [97] to better map/understand the 3D architecture.
The updated ASTER mineral maps (Figure 7) also provide new information about the Eucla Basin, much of which appears to cross-cut the published Eocene depositional features (blue dotted, dashed, and solid lines in Figure 7b) [47,50] that host heavy mineral sand occurrences (green dots) [54].These new features include (i) a possible barrier shoreline ("M", which is highlighted by a black dotted line) characterized by a linear edge/zone of elevated quartz sand (Figure 7f) and AlOH content (Figure 7e) rich in muscovite (Figure 7g); (ii) a possible back-beach lagoon ("J") characterized by low quartz sand and AlOH contents and is relatively rich in kaolinite as well as comprising an apparent meandering (paleo) channel ("O" in Figure 7b) that extends southward from the mouth of Throssel Paleoriver and terminating at the barrier shoreline "M"; and (iii) a multi-lobate, deltaic-like fan ("I") outboard of the mouth of the Carey Paleoriver that comprises elevated quartz sand and AlOH contents and is muscovite-rich.These, and other depositional features across the Eucla Basin, are described in more detail below.

Miocene Littoral Deposits across the Eucla Basin
The published geology of the Eucla Basin [50,81,82,86,100-106] includes (i) Eo barrier shoreline/lagoon systems (Figure 8a); (ii) Oligocene-to-Miocene limestones sandstones; and (iii) quaternary sands, with the major marine regression occurring i middle-to-late Miocene as the climate cooled and aridity increased.The early-to-mid ocene, ~30 m thick Nullarbor Limestone grades northward into the fossiliferous Co

Miocene Littoral Deposits across the Eucla Basin
The published geology of the Eucla Basin [50,81,82,86,[100][101][102][103][104][105][106] includes (i) Eocene barrier shoreline/lagoon systems (Figure 8a); (ii) Oligocene-to-Miocene limestones and sandstones; and (iii) quaternary sands, with the major marine regression occurring in the middle-to-late Miocene as the climate cooled and aridity increased.The early-to-mid-Miocene, ~30 m thick Nullarbor Limestone grades northward into the fossiliferous Colville Sandstone, which is <23 m thick and up to 120 km wide (the published "approximate" boundary [105] is shown by the orange line in Figure 8b-d).The Colville Sandstone comprises thin-interbedded, micaceous sandstones, sandy calcarenites, shales, siltstones, and conglomerates, which were deposited on the margin a marine shelf (30-45 m depth) characterized by good circulation.In the eastern part of the Eucla Basin, the thinner Yarle Sandstone (<10 km wide) is a silicified and ferruginous sandstone unit containing rare clasts of Nullarbor Limestone (as well as fossils) that were deposited along a moderateenergy beach strandline.These early-mid-Miocene Eucla-group rocks became emergent (and remain so to this day) by the Pliocene, driven by global cooling (falling sea level) and tectonic uplift (at 8-9 Ma) across southern Australia [50,78].
The mapping of the mid-and late-Eocene shorelines across the Eucla Basin is wellestablished [47] (Figure 8a), though mapping the Miocene shoreline/s, especially the Colville Sandstone, has been more of a challenge [50].In the field, the Colville Sandstone is distinguishable by a surface pediment rich in quartz sand variably covered by mulga woodlands.This changes southward to myall scrub developed over calcareous soils associated with the Nullarbor Limestone [100].Another challenge is unravelling the relationship between the contemporaneous (?) Miocene Colville and Yarle Sandstones.
These issues are targeted here using the ASTER products: (i) water content (Figure 8a); (ii) silica index (Figure 8b); (iii) AlOH content (Figure 8c); (iv) AlOH composition (Figure 8d); and (v) the carbonate index (B 13 /B 14 ) [107].Note that the ASTER carbonate index image is not presented in Figure 8 because it is visually compromised by image-to-image miscalibration errors largely caused by uncorrected water vapor [5].It nonetheless has sufficient "mineralogical" contrast to define the approximate boundary between regions of high versus low carbonate contents, with this boundary shown as a white dashed line in Figure 8b-d.Importantly, all of the ASTER products presented in Figure 8 show close similarity with their respective field validation NGSA measurements (colored dots), underlying the reliability of the geoscience information available from these new ASTER products but not present in the original 2012 ASTER products [5].
The water content image (Figure 8a) shows large (~500 km wide), coherent domains, including: (i) a relatively "wet" (=warmer tones) Eucla Basin; (ii) moderate water contents across the mid-Eocene barrier-beach lagoon zone, except for southwest of the Carey Paleoriver, where it is "dry"; (iii) moderate levels over much of the Phanerozoic cover rocks, to the north, including over the Officer Basin; and (iii) "dry" (=cooler tones) zones over parts of the Yilgarn Block and Albany-Fraser Mobile Zone in the west (discussed previously) and the Gawler Craton in the east (red lines).Anomalously "wet" satellite overpasses produce step effects, e.g., green arrows, which were likely caused by rainfall occurring shortly before image acquisition.Ideally, the affected images should be replaced in the National ASTER mosaic with "drier" images.
The uniformly "wet" ASTER signature of the Miocene Eucla Basin (Figure 8a) is cross-cut by a pattern of quartz sand in the north versus carbonate in the south (Figure 8b), with the boundary between these largely coincident, except near (i) "P", where there is a ~100 km wide area of elevated carbonate and silica contents, and (ii) "S", which is an alluvial fan/lobe-like feature off a paleo-river mouth to the west.This lobe also shows an apparent image/track calibration error (yellow arrows in Figure 8b) but is evident in the ASTER TIR module only.That is, this TIR mis-calibration also likely impacts the carbonate index but not the SWIR mineral information (Figure 8c,d).The pattern of AlOH content (Figure 8c) cross-cuts the "wet" AlOH signature of the Eucla Basin (Figure 8a) and broadly coincides with high silica content, though it also shows opposite patterns, e.g., "Q" and "R".Based on the ASTER maps (Figure 8b-d), the Yarle Sandstone ("T") is characterized by (i) a high silica index, i.e., high quartz sand content; (ii) a low carbonate content; (iii) a low AlOH content; and (iv) is kaolinitic.The Yarle Sandstone ("T") exhibits a pronounced north-south gradient in the silica index (Figure 8b) and has a geometry that is relatively smooth, curvilinear, and westward-widening, terminating close to the carbonate index transition (white dashed line) between U and T.These characteristics contrast with the Colville Sandstone (spans "Q", "P" and "R") which is characterized by (i) moderate levels of quartz sand content; (ii) variable carbonate content; (iii) generally high AlOH content; (iv) is generally muscovite-rich; and (v) cross-cuts the carbonate index transition.The Colville Sandstone also shows a different pattern of relatively diffuse N-S and E-W compositional gradients/boundaries, as well as a series of E-W linear, muscovite-rich zones (purple arrows in Figure 8c).
Based on this new ASTER information (Figure 8), we argue that the Miocene Yarle and Colville Sandstone units formed at different times and under different geomorphic/tectonic/climatic conditions (discussed in more detail below).

Discussion
The results above show that more accurate geoscience information is generated from the Australian ASTER reflectance/emissivity mosaic when vegetation unmixing is applied to affected products (Figures 1 and 4-8).These "affected" products are driven by the wavelength regions used for the mineral indices and whether these overlap with either green or dry vegetation spectral features; that is, where the vegetation is effectively "aspectral", it has no apparent effect on the mineral index.Thus, the SWIR and TIR mineral indices (e.g., the ASTER AlOH and silica indices, respectively) are prone to dry vegetation effects, while the VNIR mineral indices (e.g., iron oxyhydroxides) are prone to green vegetation effects.
Although estimation of the dry vegetation content using ASTER remains problematic, which is a consequence of the lack of diagnostic cellulose/lignin spectral bands, the SWIRand TIR-based indices tested here (Equations ( 2) and ( 3)) were nevertheless effective in suppressing regional rainfall (Figure 5) and local fire scar (Figures 4, 6 and 7) vegetation patterns.Having said this, the method tested here is subjective, requiring the user to determine the "optimum" gains in Equation ( 4).Arguably, a more objective "Model Optimization Process", as part of an automated "Machine Learning" strategy, that corrects for both regional, climate-driven and local, fire-scare-driven vegetation patterns is required.
The few updated ASTER compositional maps presented here show that new and more accurate surface mineralogy is generated and that much of this information is not apparent from other "traditional" mapping data.This new mineral information includes (i) the extent/nature of the Colville Sandstone (Figure 7), with implications for heavy mineral sand exploration; (ii) a continental-scale N-S change in AlOH content/composition, which likely reflects a change in aeolian sediment transport controlled by the high-pressure belt at ~25 • S (Figure 5b,d); (iii) "dry" zones (Figures 7h and 8a) that are anti-correlated with Archaean gold deposits and consistent with "buried" tectonic-structural zones; and (iv) a 1000 km wide, circular, muscovite-rich feature spatially associated with opal deposits (Figure 5b,d).Much of this new geological information should be examined in greater detail in future work, ideally in combination with spectral measurements (ASTER-simulated) of drill cores [98] as well as other (updated) ASTER mineral maps to generate a 3D perspective of the mineral/geologic architecture of specific mineral systems.
The results presented here show that the surface mineralogical expression of paleovalleys and their catchments is complex, being dependent on terrane, geological exposure, groundwater, and aeolian processes.The apparent contrast in complexity between different upland/lowland terranes (Figure 6 versus Figure 7, respectively) also helps explain the "cold" signature of paleo-valleys often observed in night-time thermal infrared imagery [41,50,89].That is, the detected cold signature in TIR kinetic temperature imagery is likely a cold air drainage "landform" effect rather than evaporative cooling of groundwater or a volumetric, thermal inertia effect, though either could be critical where the landform is clearly not the driver, e.g., ref. [108].
Here, we discuss in more detail how the updated ASTER maps can be interpreted to unravel the Miocene depositional history across the Eucla Basin [50].The first stage is defined by the northern limit of carbonate development associated with the deposition of the Nullarbor Limestone, which we suggest is mappable using the ASTER carbonate index (black dashed line in Figure 9a).At this time, the bulk of the terrigenous sediment entering the marine basin was sourced from the northeast, presumably driven by localized tectonic uplifts [78]), and carried by the Meramangye Paleoriver (thick red arrow), as well as other rivers from the north (thin red arrows).Even though most of the eroding catchments were rich in quartz and muscovite (Figure 8b-d), the transported AlOH mineral component was either deposited in a low-energy environment before reaching the shoreline, such as a back-beach lagoon (e.g., X in Figure 8c), or separated from the quartz sand component in the high-energy shoreline environment before being transported by easterly longshore drift to be deposited further to the east.This resulted in the Yarle Sandstone being quartz sand-rich but AlOH-poor, with the bulk of this coarser sediment preferentially deposited near the mouth of the Meramangye Paleoriver ("T" in Figure 8).
The second stage (Figure 9b) involved a westward shift in the area of upland tectonic uplift [50,78,89] during a continuing marine regression, possibly also involving a pronounced "step-down" event in sea level caused by the tectonic activity.This resulted in the deactivation of the Yarle shoreline system and the initiation (expansion) of the Colville shoreline system, which operated both westward and southward of the previous axis of Yarle long-shore drift (Figure 9a).This also marked a change from a higher-energy, quartz sand-dominated shoreline to a lower-energy, AlOH-mineral-dominated shoreline, with the latter characterized by (i) a west-to-east gradient from quartz sand-rich (eastern limit shown by the orange line) to AlOH-rich (purple and green shading), which presumably reflects distance from the primary sediment source, direction of longshore drift and the particle size/shape of transported mineral grains; and (ii) a southward progradation of the shoreline, characterized by linear, muscovite-rich zones that intersect/start in the west and spread and fade eastward (e.g., purple arrows in Figure 8d).The Carey Paleoriver and, to a lesser degree, the Baker Paleoriver (magenta arrows in Figure 9) delivered much of this muscovite-rich sediment, while the Throssel, Wanna, and Noorina Paleorivers delivered kaolinite-rich sediments (green arrows in Figure 9b).
The third and final stage (Figure 9c) was associated with a continuing westward shift in tectonic uplift, now centered over the Yilgarn Block (Figure 5b,d), as well as more localized uplift, such as that over the basement block bounded by the Rodona and Mundrabilla Shear zones [109][110][111][112] (black dotted line in Figure 9d), which resulted in the emergence of the Bunda Plateau (blue dot-dashed polygon in Figure 9c).A suite of faults (yellow lines in Figure 9d) and rotated blocks [111,112] are also evidence of this uplift.
This localized uplift is expressed in the gamma radiometric signal (Figure 9d).That is, lower gamma radiometric signal generally corresponds to higher ASTER AlOH content (compare Figures 8c and 9d) associated with alluvial fans (e.g., "S" and "T"), which extend for ~200 km to the edge of the uplifted "basement" block (black dotted line), with the AlOH composition of each fan reflecting its sedimentary provenance.For example, the Throssel Paleoriver ("Q") delivered kaolinite-rich sediments (green arrow and green shading) from Phanerozoic cover rocks (Figure 7g), whereas paleo-rivers to the south (e.g., Carey Paleoriver "T") delivered muscovite-rich sediments (magenta arrows and magenta shading) from the Yilgarn Craton (Figure 7g).This third phase of sedimentation was likely contemporaneous with a pulse of late-Miocene K-rich terrigenous sediments identified [76] in ocean drill core sediments (DH U1459 in Figure 5d) on the western side of the Yilgarn Craton.This K-rich (muscoviterich) ocean sediment was proposed by the authors of [76] to be caused by increased rainfall over the Yilgarn Craton driven by a brief northward migration in the Miocene of the high-pressure belt.However, we suggest that the Miocene K-rich ocean sediment recorded at U1459 (Figure 5d) is the result of the activation of fluvial systems across the Yilgarn Craton caused by westward-intensifying, tectonic uplift across a continentalscale "hinge-line" (black dashed line in Figure 5) [78].At any rate, the high-pressure belt is (again?)today located at ~25 • S (yellow dot-dashed line in Figure 5), generating opposing E-W wind patterns (yellow arrows in Figure 5d) that, respectively, transport aeolian AlOH-rich sediments across Australia (Figure 5d) and into the neighboring marine environments [6,113,114].That is, a flow of kaolin-rich aeolian sediment westwards in the north and muscovite aeolian sediment eastwards in the south (large open yellow arrows in Figure 5d), which are clearly mappable using the updated ASTER AlOH mineral maps (Figures 5d and 7e,g).
Given that the "Colville" shoreline mapped using the ASTER (Figures 8 and 9b) is correct, this presents a targeting opportunity for heavy mineral sand exploration across the central part of the Eucla Basin, especially given that Miocene age deposits (reworkings of Eocene deposits) have been identified [97,115] on the margins of the basin (red triangles in Figure 9b).The Yilgarn Craton source rocks of these interpreted Miocene shorelines are particularly attractive, given that world-class heavy mineral sand deposits on the west coast of Australia also sourced their economic endowment from the Yilgarn Craton.However, the interpreted low-energy nature of the "Colville" shoreline potentially limits the significance of this opportunity.
Finally, the updated Australian ASTER products presented here clearly show new geoscience information not apparent in other public mapping datasets.Even though this unique mineral information is sensed from only the top few microns, its complexity speaks of a multiplicity of processes that, in theory, can be untangled to generate potentially unexpected information about our Earth system, such as that shown here for the Miocene geomorphic history of the Eucla Basin.Importantly, the methods and maps presented in this paper are a significant improvement of those published in 2012 for the public release of the Australian ASTER (version 1) geoscience maps [5,6], not just because ~90% of all pixels now provide compositional information (<60% in 2012) but also because these pixels provide more accurate mineral information.Whether this improvement assists in the journey of capturing the value of the Australian ASTER reflectance/emissivity mosaic [34], and ideally the entire ASTER Global Data Archive [116], remains to be seen.However, we also recognize that even though this opportunity can deliver significantly new and unique geoscience information, the level of mineralogical detail is constrained by the ASTER's (multi) spectral resolution and spatial (30-90 m pixel) sampling.The current suite of operational, hyperspectral VNIR-SWIR imaging systems, including (a) China's GF-5; (b) Italy's PRISMA; (c) Japan's HISUI; (d) Germany's EnMap; and the (e) USA's EMIT, should, in theory, improve the specificity of this mineral information.Indeed, NASA has just published a preliminary, composite hematite-goethite-kaolinite mineral map of the world's desert regions [117], including a large part of the Australian continent.From these data, a more precise map of AlOH-bearing minerals can, in theory, be generated, as well as a (much) better estimation of the dry vegetation cover, e.g., using the 2080 nm cellulose absorption rate.At any rate, ASTER's multispectral VNIR-SWIR-TIR global data archive provides a historical baseline of the Earth's land surface composition of the early 2000s, which, following appropriate processing, has the potential to still surprise us with new and valuable geoscience information.
data archive provides a historical baseline of the Earth's land surface composition of the early 2000s, which, following appropriate processing, has the potential to still surprise us with new and valuable geoscience information.(a-c) use a DEM [51] as the base map (warmer colors are higher elevation), while (d) uses a ternary gamma radiometric map (R:G:B; K:Th:U) [53].The boundary of a basement rock gravity and magnetic high feature [118] is shown as a dotted line in Figure (c,d).Yellow lines in (d) are faults with surface expression [112].Q, S and T are locations referred to in the body text.

Conclusions
This study aimed to improve the geoscience information content extractable from the publicly available Australian ASTER reflectance/emissivity mosaic [5,6,34] through a reassessment of selected spectral indices used to generate selected mineral products suitable for mapping (paleo)geomorphic features across the arid parts of Australia, including the silica index (~quartz sand content), AlOH content and composition (e.g., kaolinite and muscovite/illite), and water content.In particular, these ASTER products were tested using an unmixing strategy to suppress/remove the effects of both the green and dry vegetation components to leave, as residual, the target surface mineral information.These results were validated using the NGSA field sample data suite, as well as by visual improvements in the derived imagery at different spatial scales.The major results include: 1.
The AlOH content is more accurately gauged using B 7 /B 6 ; 2.
The AlOH composition is more accurately gauged using B 7 /B 5 , though in color imagery, we recommend a R:G:B combination of B 7 /B 6 : B 7 /B 5 : B 7 /B 6 ; 3.
Despite the ASTER's lack of spectral bands covering diagnostic dry vegetation features, both N-SWIR (Equation ( 2)) and TIR (Equation ( 3)) band combinations yield useful approximations, albeit with some residual mineral information remaining and the N-SWIR index showing lesser artefacts (and better pixel resolution); 5.
The ASTER vegetation unmixing strategy requires a stepwise, iterative approach to assess whether vegetation is compromising the pixel mineral information.We recommend using both 2D scattergrams of the target mineral and vegetation indices, a selection of appropriate scale factors (gains in Equation ( 4)), and a visual assessment of the resultant unmixed imagery.Note that the size of the gains used in the unmixing algorithm can have pronounced, potentially erroneous effects on the derived product if not properly gauged; 6.
The AlOH content/composition and silica index products were compromised by vegetation, but the water content product was not; 7.
At a continental-scale, the vegetation unmixing suppresses/removes the vegetation patterns driven by climate-driven rainfall, while at the regional-level, the vegetation unmixing removes the effect of fire scars; 8.
The new methodology resulted in 90% of the pixels being mapped for mineralogy, in contrast with the 60% provided in the 2012 Australian ASTER mineral products; 9.
Vegetation unmixing revealed a regional N-S zonation across the ~25 • S latitude, especially in the western half of the continent, characterized by abundant AlOH minerals, especially muscovite, in the south versus less-abundant AlOH minerals, which are relatively kaolinite-rich in the north.This zonation is likely driven by the high-pressure belt and its effect on the prevailing wind patterns and associated transport of aeolian sediment; 10.Vegetation unmixing revealed the paleo-drainage network across the aeolian sandcovered Canning Basin to comprise low AlOH content, similar to the current-day Darling River.In more complex, "upland" geological environments, like the Yilgarn Craton, paleo-rivers have a more complicated, discontinuous surface mineralogical pattern driven by a variety of processes, including catchment rock type and groundwater condition; 11.Vegetation unmixing better maps the composition of paleo-rivers (and their source catchments), Miocene shorelines, and alluvial fans across the Eucla Basin, including the Colville Sandstone, which was distinguished by linear patterns in quartz sand and AlOH minerals that present a new opportunity for heavy mineral sand exploration; 12.We interpreted the updated ASTER maps of the Eucla Basin to show three stages of basin sedimentation, including: (i) stage 1-high-energy shoreline, Yarle Sandstone development related to tectonic uplift in the NE and relatively higher sea level and removal of the AlOH mineral component; (ii) stage 2-lower-energy shoreline, Colville Sandstone development related to westward shift in tectonic uplift and lowering sea level and the deposition of AlOH minerals from different source regions; and (iii) Stage 3-cessation of the shoreline system caused by a regional westward shift in tectonic uplift coupled with a localized Eucla Basin basement block uplift and lowering sea level resulting in the development of alluvial fans with different AlOH mineral compositions depending on paleo-river catchment rocks; 13.The non-vegetation unmixed water content product revealed a 100 km wide, 600 km long "dry" zone located in the eastern margin of the Yilgarn Craton, apparently associated with a deep tectonic structure (from published seismic data and drilling) and an apparent lack of Archaean gold deposits; 14.A 1000 km wide circular pattern over the Lake Eyre region, characterized by a broad rim of abundant "muscovite" (and other mineral products), is spatially coincident with opal fields.An explanation for this association is to be followed up in future work.
Author Contributions: T.C.: conception, data curation, methodology, analysis, writing; L.C.: data curation.All authors have read and agreed to the published version of the manuscript.
Funding: This research received no external funding.
Acknowledgments: I would like to thank my ASTER colleagues from around the world who helped achieve the vision of creating a complete suite of ASTER geoscience maps of Australia (see [6] for details).This paper commemorates the short life of William (Liam) Johannes Cudahy (5 April 1996 to 30 May 2021).Liam was beginning his career as an environmental scientist when the COVID-19 pandemic and its associated lockdowns caused him to postpone his university honors research project at the University of Western Australia on "Soil properties of variably-performing rehabilitated lands associated with the Greenbushes mine, Western Australia".Liam's honors thesis was to be based on processed multi-temporal satellite imagery using the unmixing techniques presented in this paper.
tion; 12.We interpreted the updated ASTER maps of the Eucla Basin to show three stages of basin sedimentation, including: (i) stage 1-high-energy shoreline, Yarle Sandstone development related to tectonic uplift in the NE and relatively higher sea level and removal of the AlOH mineral component; (ii) stage 2-lower-energy shoreline, Colville Sandstone development related to westward shift in tectonic uplift and lowering sea level and the deposition of AlOH minerals from different source regions; and (iii) Stage 3-cessation of the shoreline system caused by a regional westward shift in tectonic uplift coupled with a localized Eucla Basin basement block uplift and lowering sea level resulting in the development of alluvial fans with different AlOH mineral compositions depending on paleo-river catchment rocks; 13.The non-vegetation unmixed water content product revealed a 100 km wide, 600 km long "dry" zone located in the eastern margin of the Yilgarn Craton, apparently associated with a deep tectonic structure (from published seismic data and drilling) and an apparent lack of Archaean gold deposits; 14.A 1000 km wide circular pattern over the Lake Eyre region, characterized by a broad rim of abundant "muscovite" (and other mineral products), is spatially coincident with opal fields.An explanation for this association is to be followed up in future work.

Figure 2 .
Figure 2. Scattergrams of spectral parameters targeting the relative depths of the 2165 and 2207 n absorption rates, both of which are sensitive to AlOH clay minerals.(a) All NGSA ASD data at a nm spectral resolution.(b) All NGSA ASD data at a 40 nm ASTER/WV3 spectral resolution.(c) A TER dry vegetation unmixed clay composition ratios (obtained from the Australia-wide dataset)

Figure 2 .
Figure 2. Scattergrams of spectral parameters targeting the relative depths of the 2165 and 2207 nm absorption rates, both of which are sensitive to AlOH clay minerals.(a) All NGSA ASD data at a 10 nm spectral resolution.(b) All NGSA ASD data at a 40 nm ASTER/WV3 spectral resolution.(c) ASTER dry vegetation unmixed clay composition ratios (obtained from the Australia-wide dataset).

Figure 5 .
Figure 5. ASTER composition maps of Australia with comparable NGSA ASD validation data (colored dots).(a) ASTER AlOH content (B7/B6) map before green and dry (TIR) vegetation unmixing.(b) ASTER AlOH content (B7/B6) map after green and dry vegetation mixing.(c) ASTER clay composition (B7/B5) map before green and dry vegetation unmixing.Magenta tones indicate muscovite, while green tones indicate kaolinite.Intensity is related to the 2200 nm absorption depth, with brighter shades related to abundant dioctahedral clay mineral content (see Figure 2b).(d) ASTER clay composition (B7/B5) map of Australia after green and dry (TIR) vegetation unmixing.(e) The ASTER silica index (B13/(B10 + B11 + B12)) after green and dry vegetation unmixing.(f) ASTER water content (B5 + B8)/(B7 + B9) product (Table1) masked for green vegetation only.All maps have been masked to remove to the top ~10% of pixels with abundant green.NGSA validation data embedded as colored dots.The magenta line is the annual 250 mm rainfall isohyet[48].The dashed, cyan, white, and yellow polygons show the Yilgarn, Arunta-Musgrave, and Gawler Blocks, respectively.Yellow arrows highlight the Darling River.The yellow line marks the "Menzies Line"[79].Red triangles show the major opal deposit fields.Red dots (with labels) are sites of the ocean drilling program

Figure 5 .
Figure 5. ASTER composition maps of Australia with comparable NGSA ASD validation data (colored dots).(a) ASTER AlOH content (B 7 /B 6 ) map before green and dry (TIR) vegetation unmixing.(b) ASTER AlOH content (B 7 /B 6 ) map after green and dry vegetation mixing.(c) ASTER clay composition (B 7 /B 5 ) map before green and dry vegetation unmixing.Magenta tones indicate muscovite, while green tones indicate kaolinite.Intensity is related to the 2200 nm absorption depth, with brighter shades related to abundant dioctahedral clay mineral content (see Figure 2b).(d) ASTER clay composition (B 7 /B 5 ) map of Australia after green and dry (TIR) vegetation unmixing.(e) The ASTER silica index (B 13 /(B 10 + B 11 + B 12 )) after green and dry vegetation unmixing.(f) ASTER water content (B 5 + B 8 )/(B 7 + B 9 ) product (Table1) masked for green vegetation only.All maps have been masked to remove to the top ~10% of pixels with abundant green.NGSA validation data embedded as colored dots.The magenta line is the annual 250 mm rainfall isohyet[48].The dashed, cyan, white, and yellow polygons show the Yilgarn, Arunta-Musgrave, and Gawler Blocks, respectively.Yellow arrows highlight the Darling River.The yellow line marks the "Menzies Line"[79].Red triangles show the major opal deposit fields.Red dots (with labels) are sites of the ocean drilling program (ODP)[76].The black dashed line is a proposed Neogene "approximate tilt axis"[80].A to D are locations referred to in the body text.

Figure 8 .
Figure 8. Eucla Basin ASTER composition maps.(a) ASTER water content index ((B 5 + B 8 )/(B 7 + B 9 )) image with the validation NGSA ASD 1900 nm absorption depth as colored dots.Green arrows show satellite tracks collected over "wet" (recent rainfall) ground.Red lines are (stable) low water content

Figure 9 .
Figure 9. Geomorphic reconstruction of the Miocene Eucla Basin based on the ASTER mineral maps (Figure 8) for chronological geomorphic stages 1 (a); 2 (b); and 3 (c).(a-c)use a DEM[51] as the base map (warmer colors are higher elevation), while (d) uses a ternary gamma radiometric map (R:G:B; K:Th:U)[53].The boundary of a basement rock gravity and magnetic high feature[118] is shown as a dotted line in Figure(c,d).Yellow lines in (d) are faults with surface expression[112].Q, S and T are locations referred to in the body text.