Targeting Paleovalley-Related Ferricrete Units in Yilgarn Craton Using High-Resolution Aeromagnetic Data and Spatial Machine Learning

: The ferricrete units (Fe oxide cemented colluvial-alluvial sediment) of the Yilgarn Craton in Western Australia formed during the humid tropical and sub-tropical climates of the Cenozoic. Ferricretes are generally developed on long-lived paleodrainage systems and are products of the ferruginisation of detritus provided by the continuous erosion of upslopes. These iron-rich accumulations can become Au-enriched, as is the case in several locations previously discovered in the Yilgarn Craton; many of these host economic secondary gold deposits (e.g., Moolart Well, Mt Gibson, and Bulchina), typically occurring downslope of low saprolite hills and near paleovalleys (i.e., inset-valleys). Inset-valleys are a common paleotopographic feature buried under Quaternary alluvial and colluvial sedimentary cover. Maps of these ancient channel networks can be used as a proxy for targeting ferricrete gold deposits. These inset-valley systems generally form dendritic and noisy patterns in high-resolution aeromagnetic data due to the presence of maghemite-rich nodules and detrital magnetic pisoliths on their ﬂanks. The main aim of this study was to use high-resolution aeromagnetic data to target ferricrete units related to inset-valleys systems across the Yilgarn Craton. A spatial predictive model was used to learn and predict the geological units of interest from pre-processed aeromagnetic data. The predicted inset-valleys systems were able to conﬁne the exploration space and deﬁne a new exploration frontier for ferricrete gold deposits.


Introduction
The ferricrete units of the Yilgarn Craton are Fe oxide-cemented colluvial-alluvial sediments [1] that are generally developed on long-lived drainage systems of paleovalleys. These units are formed during the humid tropical and sub-tropical climates of the early to mid-Cenozoic [2]. Ferricretes are generated through the ferruginisation of nodules and pisoliths supplied by the continuous erosion of older duricrust and lithic fragments of less-weathered materials from upslope saprock and saprolite [3,4]. Pisoliths and nodules are mainly hematite-, goethite-and kaolinite-rich. Ferricretes have been overprinted continuously through the re-weathering of sediments by groundwater in the presence of organic matter [3,4].
Au-enriched ferricretes have been discovered in multiple locations of the Yilgarn Craton and some of them host economic secondary gold mineralisation (e.g., Moolart Well, Mt Gibson, and Bulchina). Ferricrete gold deposits are typically small (<500,000 ounces) and low grade (<1-5 g/t Au). This type of deposit generally consists of flat-lying zones (a few hundred metres to several kilometres long, less than a kilometre wide, and 2-10 m thick) of Au enrichment within the upper ferruginous zone of lateritic profiles [3,4]. They can be efficiently mined because they are generally located near the surface [4][5][6][7]. The physical and chemical behaviour of Au during the formation and evolution of ferricretes are not fully known [3]. For instance, significant unweathered mineralisation beneath Au-enriched ferricretes cannot be located easily. This is because the ferricrete deposits are mechanically transported from the upper ferruginous zone of lateritic profiles and form an extensive blanket compared to unweathered mineralisation. Anand et al. [3] thoroughly investigated the ferricrete Au deposits of the Yilgarn Craton and identified the possible mechanisms of dispersion and enrichment of Au during their formation. The evolution of the ferricrete Au deposits generally consists of multiple stages of mechanical, chemical, and biological mobilisation and reprecipitation of Au in a paleodrainage environment. Au enrichment in ferricretes is likely sourced from underlying or near upslope-located primary Au-enriched rocks. The Fe oxide-rich paleodrainage systems of the Yilgarn Craton indicate that significant subterranean lateral influx and groundwater migration have transported dissolved Fe and Au along and above an unconformity during an extended period of biological and chemical weathering [8].
Au-enriched ferricretes typically occur downslope of the low saprolite hills and proximal to inset-valleys (i.e., buried paleovalleys). The predicted map of inset-valleys [9] can be used as a proxy for targeting ferricrete gold deposits as they are developed along the flanks of inset-valleys. Ferricretes are generally covered by 2-10 m of younger transported sediments (Figure 1). Inset-valleys are a common feature of regolith in the Yilgarn Craton [2,[10][11][12][13]. Large-scale structures such as faults and shears may have influenced the position of paleovalleys; however, these structures may not be associated with Au mineralisation [14]. Well, Mt Gibson, and Bulchina). Ferricrete gold deposits are typically small (<500,000 ounces) and low grade (<1-5 g/t Au). This type of deposit generally consists of flat-lying zones (a few hundred metres to several kilometres long, less than a kilometre wide, and 2-10 m thick) of Au enrichment within the upper ferruginous zone of lateritic profiles [3,4]. They can be efficiently mined because they are generally located near the surface [4][5][6][7]. The physical and chemical behaviour of Au during the formation and evolution of ferricretes are not fully known [3]. For instance, significant unweathered mineralisation beneath Au-enriched ferricretes cannot be located easily. This is because the ferricrete deposits are mechanically transported from the upper ferruginous zone of lateritic profiles and form an extensive blanket compared to unweathered mineralisation. Anand et al. [3] thoroughly investigated the ferricrete Au deposits of the Yilgarn Craton and identified the possible mechanisms of dispersion and enrichment of Au during their formation. The evolution of the ferricrete Au deposits generally consists of multiple stages of mechanical, chemical, and biological mobilisation and reprecipitation of Au in a paleodrainage environment. Au enrichment in ferricretes is likely sourced from underlying or near upslopelocated primary Au-enriched rocks. The Fe oxide-rich paleodrainage systems of the Yilgarn Craton indicate that significant subterranean lateral influx and groundwater migration have transported dissolved Fe and Au along and above an unconformity during an extended period of biological and chemical weathering [8].
Au-enriched ferricretes typically occur downslope of the low saprolite hills and proximal to inset-valleys (i.e., buried paleovalleys). The predicted map of inset-valleys [9] can be used as a proxy for targeting ferricrete gold deposits as they are developed along the flanks of inset-valleys. Ferricretes are generally covered by 2-10 m of younger transported sediments ( Figure 1). Inset-valleys are a common feature of regolith in the Yilgarn Craton [2,[10][11][12][13]. Large-scale structures such as faults and shears may have influenced the position of paleovalleys; however, these structures may not be associated with Au mineralisation [14].  Maghemite-rich nodules and pisoliths in inset-valley systems form dendritic and noisy patterns in high-resolution aeromagnetic data. Aeromagnetic data are often used to map paleovalleys [15,16], and resolution permitting, are accompanied with radiometric or different types of electromagnetic data [17,18]. Maghemite is a common mineral in ferricretes overlying saprolite units; however, its levels decrease sharply down the weathering profile [2]. The heating of goethite during bush fires may be a reason for maghemite formation [19]. High-resolution aeromagnetic data can be used efficiently to identify maghemite-rich ferruginous nodules and pisoliths in paleovalley systems. These maghemite-rich units form noisy and dendritic patterns visible in imaged aeromagnetic data [2,14]. Periodic accumulation or the reworking of maghemite-rich ferruginous nodules and pisoliths from higher elevations in the paleodrainage system likely have generated more than one magnetic edge feature in inset-valleys. Inset-valleys with magnetic features vary in width from 100 m to 1.5 km and are up to 30 km in length.
Applied machine learning is considered an established tool in the geosciences, including solid Earth geoscience [20,21], geochemistry [22], mineral potential mapping [23][24][25], remote sensing [26], and seismic interpretation [27]. However, it is important to acknowledge the shortcomings of machine learning in geoscience applications [28,29]. The particularities of geosystems and geoscience data (e.g., multi-source, multi-scale, spatial and temporal heterogeneity, auto-and cross-correlations) must be fully understood before any implementation of machine learning models. Suitable machine learning algorithms for geoscience applications should be accurate, interpretable, and physically realistic.
The geological features of interest in this study are the paleovalley-related ferricrete units with unique magnetic characteristics. These units form noisy and dendritic patterns in gridded aeromagnetic data. Capturing such complex spatial patterns is not possible via classical machine learning algorithms because these techniques only use the information available at single points. Spatially-aware machine learning algorithms, such as convolutional neural networks [30,31] and spatial random forests [32], can be used to predict such complex spatial patterns automatically.
The main aim of this study is to use high-resolution aeromagnetic data and spatial machine learning to automatically target paleovalley-related ferricrete units with magnetic characteristics. High-resolution aeromagnetic data are pre-processed to enhance the shortwavelength features produced by shallow sources and suppress long-wavelength magnetic features. A spatial predictive model (i.e., spatial random forests, Talebi, et al. [32]) is subsequently used to learn paleovalley-related ferricrete units with magnetic characteristics from pre-processed aeromagnetic data. The trained spatial random forests model is finally used to predict ferricrete units in other areas across the Yilgarn craton where high-resolution aeromagnetic data are available. The predicted ferricrete units can be used later to confine the exploration space and define new exploration targets for ferricrete gold deposits.

Study Areas
Three areas in the Yilgarn craton ( Figure 2) were selected for the proposed experiment. There are publicly available high-resolution aeromagnetic data for Tarmoola, Bronzewing, and Moolart Well in Western Australia, collected by the various exploration companies over the last few decades. The following sub-sections discuss the geological settings of these three areas.

Figure 2.
Location of three case study areas (black rectangles crosshatched) within the Yilgarn Craton in grey colour. The regional extent of the known mapped paleovalleys in blue after Bell, et al. [33].

Tarmoola
The Tarmoola gold deposit is located 30 km to the north-west of Leonora in the Eastern Goldfields and is 265 km north of Kalgoorlie ( Figure 2). It is the largest granitoidhosted gold deposit in the Yilgarn Craton and produces ~28 Mt of ore with an average grade of 1.8 g/t Au; secondary mineralisation occurs in multiple styles including shallow ferricrete gold and supergene-enriched saprolite.

Bronzewing
The Bronzewing ferricrete Au deposit is located in the Yandal greenstone belt, 375 km north of the Kalgoorlie gold mine. No primary mineralisation has been reported beneath the deposit [5]. In total, 19,000 ounces of gold (average ore grade of 1.49 g/t Au) have been mined at this location. This deposit is located on alluvial systems that are underlain by 30-50 m deep inset-valleys. Variably indurated yellowish-brown ferricrete at the flanks of the inset-valley hosts the Au mineralisation. The ferricrete units are overlain by about 7 m of younger transported cover. The lower part of the younger transported cover is Figure 2. Location of three case study areas (black rectangles crosshatched) within the Yilgarn Craton in grey colour. The regional extent of the known mapped paleovalleys in blue after Bell, et al. [33].

Tarmoola
The Tarmoola gold deposit is located 30 km to the north-west of Leonora in the Eastern Goldfields and is 265 km north of Kalgoorlie ( Figure 2). It is the largest granitoid-hosted gold deposit in the Yilgarn Craton and produces~28 Mt of ore with an average grade of 1.8 g/t Au; secondary mineralisation occurs in multiple styles including shallow ferricrete gold and supergene-enriched saprolite.

Bronzewing
The Bronzewing ferricrete Au deposit is located in the Yandal greenstone belt, 375 km north of the Kalgoorlie gold mine. No primary mineralisation has been reported beneath the deposit [5]. In total, 19,000 ounces of gold (average ore grade of 1.49 g/t Au) have been mined at this location. This deposit is located on alluvial systems that are underlain by 30-50 m deep inset-valleys. Variably indurated yellowish-brown ferricrete at the flanks of the inset-valley hosts the Au mineralisation. The ferricrete units are overlain by about 7 m of younger transported cover. The lower part of the younger transported cover is significantly Au-enriched (up to 600 ppb). The Au contained in the ferricrete is probably derived from nearby deposits [5].

Moolart Well
The Moolart Well is the largest ferricrete Au deposit of the Yilgarn Craton, located 400 km NE of Kalgoorlie in the Duketon greenstone belt. Intermediate and tholeiitic dolerite sills (up to 200 m thick) have intruded the interflow sedimentary horizons [4,34] and formed the bedrock of the area. Low-grade primary mineralisation (0.3-0.5 g/t Au; 103,000 oz Au) is predominantly hosted within diorite intrusives and a north-trending shear zone. The economic Au mineralisation is of a secondary nature and confined to the upper flat-laying ferricrete units (i.e., yellowish-brown and reddish-brown ferricretes extend over 4 km N-S and up to 1 km E-W) within the paleovalley sediments. The ferricrete deposit has an average thickness of 4 m. In addition, the underlying Au-bearing saprolite is derived from the weathering of primary mineralisation [4]. The Moolart Well Au deposit is situated in a paleovalley system bounded by low hills of saprolite and saprock.

Materials and Methods
The high-resolution aeromagnetic data were pre-processed to enhance the characteristic patterns of the inset-valley's ferricrete units and to prepare the input raster layers for analysis through the spatial random forest predictive model workflow. The implemented pre-processing steps and the spatial random forests algorithm are discussed in the following subsections.

High-Resolution Aeromagnetic Data
The high-resolution aeromagnetic data used in this study were sourced from the Western Australia Department of Mining, Industrial Regulation and Safety (DMIRS) data portal (dasc.dmirs.wa.gov.au, accessed on 10 June 2021). Data was extracted from various company surveys recorded between 1996 and 2003. Details of these surveys are available in Table 1. The first step in the preparation of a consistent aeromagnetic input layer for all three areas involved generating total magnetic intensity (TMI) grids of consistent resolution, using line data from the surveys listed in Table 1. The upward and downward continuation transformations [35] were used to reproduce the aeromagnetic data at the designated height (25 m).
The aeromagnetic dataset for the Moolart Well area was firstly upward-continued for 5 m and then gridded by the use of a 25 m grid cell size. The aeromagnetic data for the Tarmoola were downward-continued for 25 m and then gridded with a 25 m grid cell size. The individual surveys in the Bronzewing area were upward/downwardcontinued to reproduce a signal at the same height of 25 m prior to gridding. The same gridding algorithm (i.e., minimum curvature [36]) was used for all areas. These TMI grids of consistent resolution were used as a basis for the reduction to pole transformation, as described below.
The texture and amplitude variation of the magnetic signal strongly depend on various structural elements (e.g., faults, folds, and lithological factors) and the variable distribution of the magnetic minerals among different units. In order to enhance features of interest, a few transformation techniques were applied to the aeromagnetic data. The main goal was to emphasise shallow inset-valley-related ferricrete units. In addition, long-wavelength magnetic features should have been suppressed. Since the available data is of high resolution with considerable high-frequency content, applying almost any type of high and low pass filters would cause Gibb's ringing, a common artifact with this type of processing.
Firstly, the reduction to pole transformation (RTP) was applied to the TMI data. To avoid processing artifacts, the upward-continued version of the RTP data were generated at the height of an additional 50 m on top of the flying height (25 m). This upward-continued data were subtracted from the original RTP data, suppressing long-wavelength magnetic content and producing a data set with high frequency (HF) features. The short-wavelength features of shallow sources were further enhanced by applying the first vertical derivative transformation (1VD) to the HF content data. This final transformed data set was used as an input for SRF. All these transformations are further described in Table 2. The processing of the aeromagnetic data was undertaken in Seequent Oasis Montaj. Figure 3 shows the pre-processing sequence for the high-resolution aeromagnetic data in the Tarmoola area.

Spatial Machine Learning
The spatial random forests technique [32] is a spatially-aware machine learning algorithm [28] on the basis of higher-order spatial statistics and is used for the analysis and modelling of geospatial data. Unlike non-spatial machine learning techniques that use the spectral information in a single pixel as input predictors, the SRF technique uses local spatial-spectral information to learn heterogeneity, auto-and cross-correlations, and complex spatial patterns. This study used an SRF classifier to target inset-valley-related ferricrete units.
The most important parameters for the SRF technique are the size of the local patterns, the number of spatial trees in the forest, and the number of randomly selected

Spatial Machine Learning
The spatial random forests technique [32] is a spatially-aware machine learning algorithm [28] on the basis of higher-order spatial statistics and is used for the analysis and modelling of geospatial data. Unlike non-spatial machine learning techniques that use the spectral information in a single pixel as input predictors, the SRF technique uses local spatial-spectral information to learn heterogeneity, auto-and cross-correlations, and complex spatial patterns. This study used an SRF classifier to target inset-valley-related ferricrete units.
The most important parameters for the SRF technique are the size of the local patterns, the number of spatial trees in the forest, and the number of randomly selected predictors at each split. More information regarding the details of the SRF technique and suitable strategies used for tuning hyperparameters can be found in Probst, et al. [39] and Talebi, et al. [32].

Model Training
The processed (i.e., RTP, HF, and 1VD-filtered) high-resolution aeromagnetic layer for the Tarmoola area was used to train the SRF classifier. A balanced dataset was selected manually from this area, including 850 samples from the dendritic and noisy patterns associated with the inset-valley-related ferricrete units (blue dots in Figure 4) and 850 randomly selected samples from the rest of the area (red dots in Figure 4). In addition, validation datasets (500 samples from each class) were selected from the three study areas to assess the overall performance and generalisation capability of the SRF classifier.

Model Training
The processed (i.e., RTP, HF, and 1VD-filtered) high-resolution aeromagnetic layer for the Tarmoola area was used to train the SRF classifier. A balanced dataset was selected manually from this area, including 850 samples from the dendritic and noisy patterns associated with the inset-valley-related ferricrete units (blue dots in Figure 4) and 850 randomly selected samples from the rest of the area (red dots in Figure 4). In addition, validation datasets (500 samples from each class) were selected from the three study areas to assess the overall performance and generalisation capability of the SRF classifier. A grid-search was used to find the optimal hyperparameters of the SRF model, which resulted in 400 pixels (i.e., 20 pixels × 20 pixels) as the optimum size of local patterns, 1000 spatial trees in the forest, and the default value, √400 20, for the number of randomly selected predictors at each split. This trained SRF classifier was used to predict the dendritic and noisy patterns associated with the inset-valley-related ferricrete units in the Tarmoola, Bronzewing, and Moolart Well areas.

Tarmoola
The processed high-resolution aeromagnetic layer for the Tarmoola area is shown in Figure 5. The predicted results of the SRF classifier were draped over the high-resolution aeromagnetic data. The blue colour represents areas predicted as the positive class- A grid-search was used to find the optimal hyperparameters of the SRF model, which resulted in 400 pixels (i.e., 20 pixels × 20 pixels) as the optimum size of local patterns, 1000 spatial trees in the forest, and the default value, ceiling ( √ 400) = 20, for the number of randomly selected predictors at each split. This trained SRF classifier was used to predict the dendritic and noisy patterns associated with the inset-valley-related ferricrete units in the Tarmoola, Bronzewing, and Moolart Well areas.

Tarmoola
The processed high-resolution aeromagnetic layer for the Tarmoola area is shown in Figure 5. The predicted results of the SRF classifier were draped over the high-resolution aeromagnetic data. The blue colour represents areas predicted as the positive classmapping the extent of the magnetic gravels-and the red colour represents the areas predicted as the negative class, devoid of these gravels. The majority of the dendritic and noisy patterns related to the inset-valley ferricretes units are accurately predicted. Table 3 shows the confusion matrix for the SRF classifier on the basis of the validation dataset from the Tarmoola area. The overall accuracy of the SRF model in this area was 91.9%.
Minerals 2022, 12, x FOR PEER REVIEW mapping the extent of the magnetic gravels-and the red colour represents the are dicted as the negative class, devoid of these gravels. The majority of the dendr noisy patterns related to the inset-valley ferricretes units are accurately predicted. shows the confusion matrix for the SRF classifier on the basis of the validation from the Tarmoola area. The overall accuracy of the SRF model in this area was 91 Figure 5. Processed high-resolution aeromagnetic layer for the Tarmoola area overlain by dicted results via the SRF classifier. The blue colour represents the areas predicted as the class whereas the red colour represents the areas predicted as the negative class. The yellow shows the point location of the main deposit, with a yellow outline giving the approximate e of the mineralised ferricrete unit.

Bronzewing
The processed high-resolution aeromagnetic layer for the Bronzewing area is in Figure 6. Although the majority of the dendritic and noisy patterns related to th Figure 5. Processed high-resolution aeromagnetic layer for the Tarmoola area overlain by the predicted results via the SRF classifier. The blue colour represents the areas predicted as the positive class whereas the red colour represents the areas predicted as the negative class. The yellow marker shows the point location of the main deposit, with a yellow outline giving the approximate envelope of the mineralised ferricrete unit.

Bronzewing
The processed high-resolution aeromagnetic layer for the Bronzewing area is shown in Figure 6. Although the majority of the dendritic and noisy patterns related to the insetvalley ferricrete units were accurately predicted, the proportion of misclassified areas was higher than the Tarmoola area. This was expected because the SRF model was trained via the use of a dataset from the Tarmoola area. Table 4 shows the confusion matrix for the SRF classifier on the basis of the validation dataset from the Bronzewing area. The overall accuracy of the SRF model in this area was 82.5%.
Minerals 2022, 12, x FOR PEER REVIEW Figure 6. Processed high-resolution aeromagnetic layer for the Bronzewing area overlain predicted results via the SRF classifier. The blue colour represents the areas predicted as the class whereas the red colour represents the areas predicted as the negative class. The yellow shows the point location of the main deposit, with a yellow outline giving the approximate e of the mineralised ferricrete unit.

Moolart Well
Even though the SRF classifier was trained on the basis of the dataset from t moola area, the extent of the ferricrete units hosting the Moolart Well Au depo accurately predicted (yellow dot in Figure 7). This is partly related to the impleme Figure 6. Processed high-resolution aeromagnetic layer for the Bronzewing area overlain by the predicted results via the SRF classifier. The blue colour represents the areas predicted as the positive class whereas the red colour represents the areas predicted as the negative class. The yellow marker shows the point location of the main deposit, with a yellow outline giving the approximate envelope of the mineralised ferricrete unit.

Moolart Well
Even though the SRF classifier was trained on the basis of the dataset from the Tarmoola area, the extent of the ferricrete units hosting the Moolart Well Au deposit was accurately predicted (yellow dot in Figure 7). This is partly related to the implementation of a consistent pre-processing step for the high-resolution aeromagnetic data and the generalisation capability of the SRF model. Table 5 shows the confusion matrix for the SRF classifier based on the validation dataset from the Moolart Well area. The overall accuracy of the SRF model in this area was 92%.

Discussion
Aeromagnetic data are often used to map paleovalley systems. However, the mapping process is a time-consuming manual task [33]. Manually delineated potential paleovalley system maps are subjective and depend on the operator's experience and familiarity with the study area. The use of spatial machine learning to detect paleovalley systems is more objective and faster, potentially allowing for the large-scale analysis of extensive areas-provided the appropriate resolution of the geophysical data.
This study used high-resolution aeromagnetic data and spatial machine learning to predict inset-valley-related ferricrete units. Implementing a consistent approach to the pre-processing of the aeromagnetic data was necessary to improve the generalisation capability of the predictive model. A sequence of transformations (e.g., RTP, multiple upward continuations and 1VD) was used to enhance the dendritic and noisy patterns associated with the inset-valley-related ferricrete units and suppress the long-wavelength magnetic content. The resolution of the magnetic data was crucial; while the presented method worked well with detailed surveys flown between 20-50 m above ground at 50-100 m line spacing (Table 1), it will fail to detect patterns in regional-scale data. The patterns became non-apparent in the regional data due to higher ground separation, which

Discussion
Aeromagnetic data are often used to map paleovalley systems. However, the mapping process is a time-consuming manual task [33]. Manually delineated potential paleovalley system maps are subjective and depend on the operator's experience and familiarity with the study area. The use of spatial machine learning to detect paleovalley systems is more objective and faster, potentially allowing for the large-scale analysis of extensive areasprovided the appropriate resolution of the geophysical data.
This study used high-resolution aeromagnetic data and spatial machine learning to predict inset-valley-related ferricrete units. Implementing a consistent approach to the pre-processing of the aeromagnetic data was necessary to improve the generalisation capability of the predictive model. A sequence of transformations (e.g., RTP, multiple upward continuations and 1VD) was used to enhance the dendritic and noisy patterns associated with the inset-valley-related ferricrete units and suppress the long-wavelength magnetic content. The resolution of the magnetic data was crucial; while the presented method worked well with detailed surveys flown between 20-50 m above ground at 50-100 m line spacing (Table 1), it will fail to detect patterns in regional-scale data. The patterns became non-apparent in the regional data due to higher ground separation, which attenuated the high-frequency signal and the additional along-and across-flight line signal smoothing inherently present in most gridding algorithms. The overall accuracy of the SRF classifier in the Bronzewing area was lower than in the other two study areas. The lower prediction accuracy might be related to the fact that the Bronzewing dataset was composed of multiple surveys and the slightly different near-surface geological setting in this area.
Ferricrete units form complex spatial patterns (i.e., dendritic and noisy patterns) in processed aeromagnetic layers. These patterns cannot be detected via classical non-spatial machine learning algorithms. A spatially aware machine learning algorithm such as SRF or convolutional neural networks (CNN) is needed to detect such complex spatial patterns. The SRF classifier predicted the location of the host regolith units, i.e., the magnetic ferricretes and gravels, successfully. The model could be extended to target the ferricrete Au deposits directly. For this purpose, more pre-processed geophysical, geochemical, geological, and paleotopographic layers could be added to provide additional information that would further constrain the search space. However, the SRF model is computationally demanding, and for the prediction of larger areas such as the entire Yilgarn craton or using many input layers, access to fast computers would be necessary.
Aeromagnetic mapping of the ferricrete deposits in inset-valleys is relevant to gold exploration in the Yilgarn Craton. These ferricretes are dominated by ferruginous pisoliths and nodules, which represent the coarser fraction (>2 mm in diameter)-compared to the finer fractions (<2 mm)-of the ferricrete. Maghemite only occurs in the coarse fraction, which displays the highest content of Fe 2 O 3 , Au, and As, as compared to the fine fractions. Anand et al. (2019) [3] separated the magnetic (hematite-maghemite-rich) from the nonmagnetic (hematite-goethite-rich) ferruginous nodules and pisoliths of the coarse fraction and found that both magnetic and non-magnetic fractions occur in the ferricrete at Moolart Well; the non-magnetic fraction is more abundant than the magnetic fraction, and both are highly anomalous in Au. However, Au is much more abundant in the magnetic fraction compared to the non-magnetic fraction. Therefore, mapping inset-valley systems can confine the exploration space and define the exploration frontier for ferricrete gold deposits. The proposed predictive model for mapping inset-valley systems uses high-resolution aeromagnetic data which are generally available for brownfield exploration.

Conclusions
This study used high-resolution aeromagnetic data and spatial machine learning to target ferricrete units that are related to inset-valley systems across the Yilgarn Craton. The high-resolution aeromagnetic data were pre-processed by implementing a sequence of transformations (e.g., RTP, multiple upward continuations, and 1VD) to enhance the dendritic and noisy patterns associated with the inset-valley-related ferricrete units and to suppress the long wavelength magnetic content. A spatial random forest predictive model was trained by the use of the pre-processed aeromagnetic layer in the Tarmoola area. The trained SRF model was used to efficiently predict the ferricrete units in the Tarmoola, Bronzewing, and Moolart Well areas in Western Australia. The predicted map of the inset-valley-related ferricrete units can be used as a proxy for targeting ferricrete gold deposits.