Application of artificially intelligent systems for the identification of discrete fossiliferous levels

The separation of discrete fossiliferous levels within an archaeological or paleontological site with no clear stratigraphic horizons has historically been carried out using qualitative approaches, relying on two-dimensional transversal and longitudinal projection planes. Analyses of this type, however, can often be conditioned by subjectivity based on the perspective of the analyst. This study presents a novel use of Machine Learning algorithms for pattern recognition techniques in the automated separation and identification of fossiliferous levels. This approach can be divided into three main steps including: (1) unsupervised Machine Learning for density based clustering (2) expert-in-the-loop Collaborative Intelligence Learning for the integration of geological data followed by (3) supervised learning for the final fine-tuning of fossiliferous level models. For evaluation of these techniques, this method was tested in two Late Miocene sites of the Batallones Butte paleontological complex (Madrid, Spain). Here we show Machine Learning analyses to be a valuable tool for the processing of spatial data in an efficient and quantitative manner, successfully identifying the presence of discrete fossiliferous levels in both Batallones-3 and Batallones-10. Three discrete fossiliferous levels have been identified in Batallones-3, whereas another three have been differentiated in Batallones-10.


INTRODUCTION
The Batallones Butte, located 30 km to the south of Madrid (Spain, Fig. 1A) and 1 km to the east of Valdemoro (Fig. 1B), is home to nine Late Miocene paleontological sites (named Batallones-1, Batallones-2 and so on). These sites hypothetically correspond to hourglass-shaped cavities with upper openings (Fig. 1C), formed as a consequence of pseudokarstic processes Calvo et al., 2013), where mammals became Palaeoecological, palaeoenvironmental and taphonomical studies carried out at Batallones Butte (Antón & Morales, 2000;Domingo et al., 2011Domingo et al., , 2013aDomingo et al., , 2013bDomingo et al., , 2016 have treated the upper and lower part of the hourglass-shaped cavities separately based on their different depth and their different taxonomical composition (herbivores vs. carnivores, respectively). Due to the challenges posed by the structure of the deposits that filled the Batallones Butte cavities, such as lateral facies changes, deformations, depositional asymmetries, collapse structures, slickensides and/or local tilting typical of cave deposits (Calvo et al., 2013), geological sub-levels have not been identified within each of these parts. However, the question still remains on whether or not discrete fossiliferous levels are found embedded in the apparently geologically continuous deposits.
AIAs, including those trained through Machine Learning and Deep Learning (ML & DL) techniques, present numerous possibilities for the processing of highly complex and noisy data sets (Bishop, 2006). Combined with their versatility, in many cases ML and DL approaches have outperformed human experts in a multitude of specialized tasks. Here we present a means of utilizing AIAs for the detection of patterns in 3D spatial data, establishing a method to objectively and efficiently detect discrete levels that human experts may have mistakenly overlooked. This study includes both unsupervised and supervised ML, joined through a hybrid Human-AI collaborative approach in order to establish a final model that can be used to describe the palaeontostratigraphic nature of each site.
The main goal of this study is to establish whether or not ML analysis can be used to quantitatively identify discrete fossiliferous levels in paleontological and archaeological sites based on the study of the spatial distribution of fossils. Specifically, this method has been applied and tested in two Batallones Butte sites (Batallones-3 and Batallones-10) in order to determine if the deposits are homogeneous or heterogeneous with discrete fossiliferous levels.

MATERIALS AND METHODS
The site of Batallones-3 has undergone systematic field excavations between the years 2001 and 2017, whereas excavations at Batallones-10 started in 2007 and are ongoing. Batallones-3 corresponds to the carnivoran-rich lower part of the hour-glass structure previously mentioned. It is not clear whether an upper part ever existed for this fossil site, but, if it existed, it was destroyed by slope-erosion processes (this site is located in the slope of the Batallones Butte). In turn, Batallones-10 corresponds to the herbivore-rich upper part of the hour-glass structure and, if present, will have a lower part, yet to be excavated. At both sites, standard fossil vertebrate excavation protocols were followed in the extraction of paleontological remains (Eberth, Rogers & Fiorillo, 2007). On-site documentation of excavated remains has primarily consisted of taxonomical and anatomical identification as well as in depth documentation of the spatial distribution of remains. Spatial data obtained thus include standard x, y and z coordinates according to their position within a grid, as well as the trend and plunge of elongated fossil remains. Additional data, for future taphonomical studies, were also collected concerning the degree of articulation between anatomical elements, the overall preservation of the remains and the restoration techniques used.
For the purpose of this study, the data used for both supervised and unsupervised ML applications consisted solely in the spatial distribution of remains, including their x, y and z coordinates for in-depth 3D analysis. Manageable slices were extracted, similar to the approach proposed by Canals, Vallverdú & Carbonell (2003), 2 for each site across the entire length of the excavated area. In the case of Batallones-10, one 2-m-wide slice was taken across the x and one across the y axis, searching for representative areas with the greatest object densities to capture a global vision of the entire site. In the case of Batallones-3, located in a domically-shaped cave with a debris cone in the middle under an inferred opening, two 2-m-wide slices were obtained along the site's x axis, one on each side of the debris cone.
Once slices were selected, the coordinates of each of the remains were imported into the free R software, x64 v.3.5.1 (www.rproject.org, R Development Core Team, 2018), for further analysis.

Unsupervised machine learning
For initial detection of hidden patterns among the levels of Batallones Butte site, an unsupervised density based clustering algorithm was used. Unsupervised ML algorithms are highly efficient AIAs that are exceptional for their use in tasks such as pattern recognition, anomaly detection, noise and dimensionality reduction, feature engineering as well as generative modeling (Patel, 2019). While unsupervised learning performs poorly in specifically defined tasks, their greatest advantage can be seen through their flexibility when applied to versatile datasets for general feature extraction. A popular application of unsupervised algorithms lies in clustering tasks, a component of pattern recognition that is, considered useful for processing and searching for hidden trends in highly noisy datasets.
A non-parametric Density-Based Spatial Clustering of Applications with Noise (DBSCAN) algorithm from the "fpc" R package was used for pattern recognition in 3D spatial data. Proposed as a means of overcoming clustering issues where groupings of samples are not straightforward (Ester et al., 1996), DBSCAN is highly efficient at detecting patterns in arbitrary and noisily distributed samples (Sander et al., 1998;Schubert et al., 2017).
DBSCAN performs through establishing areas of significant point densities, thus providing an efficient means of finding non-parametric patterns among noisy data sets while easily detecting outliers in low-density regions. This is especially useful when used for pattern recognition in spatial studies, where traditional clustering models such as partitioning and hierarchical algorithms tend to focus groupings around a centroid or subsamples in a convex manner. This frequently creates conflict in real-life circumstances.
DBSCAN has proven efficient when applied to 2D, 3D and/or higher dimensional feature spaces (Ester et al., 1996;Sander et al., 1998;Schubert et al., 2017). The algorithm works by separating points within the cluster (i.e., core points) from those on the border of the cluster (i.e., border points). The differentiation between the two is established mathematically depending on the reachability of neighboring points (q) from a core point (p). To define this, DBSCAN requires two main hyperparameters, the ε value defining the neighborhood (a.k.a. Eps) of a point and the minimum number of points (MinPts) that are required to form the cluster (i.e., the density of points). From here DBSCAN defines the neighborhood of p within dataset D as N Eps (p) = {q ∈ D | dist(p, q) ≤ ε}, adjusting for border points and the exclusion of noise through p ∈ N Eps (q) and |N Eps (q)| ≥ MinPts (Ester et al., 1996). The definition of density-reachable points for clustering is thus established if the points can be linked in a chain, whereby p 1 , …, p n = q, p n = p ensuring p i+1 is directly reachable through density from p i . Through this definition, border points may not be density reachable if they do not fulfill the core point condition |N Eps (q)| ≥ MinPts. The final components that define point density for clustering are established through density-connected points, whereby a point o is density-reachable from both p and q, essentially connecting both border and cluster points.
Once density has been established within the sample, the assigning of these points into cluster C requires the fulfillment of conditions; (1) ∀ p, q: if p ∈ C and q fulfill the definition of being density reachable, then q ∈ C, while (2) ∀ p, q ∈ C: p must be density-connected to q. These defined conditions are named Maximality and Connectivity, respectively (Ester et al., 1996). In accordance with the aforementioned conditions, the final definition of noise can be logically established through {p ∈ D | ∀ i: p ∉ C i }.
DBSCAN was thus trained on the entire dataset in an unsupervised manner, establishing an average MinPts value between 3 and 5. These values were intuitively chosen depending on slices with tightly packed densities (MinPts ≈ 5) and slices with a lower number of remains (MinPts ≈ 3). This abides by general rules recommended by numerous authors (Sander et al., 1998;Schubert et al., 2017), who define an approximate optimum MinPts = 2 × n dimensions in the dataset which is then adjusted depending on the complexity of point distributions and the analyst's knowledge of the domain under consideration (Sander et al., 1998;Schubert et al., 2017). ε values were established in accordance with both the correspondent MinPts values and the slices under study. For optimization of this hyperparameter, k-distance graphs were plotted according to the nearest neighbor (Ester et al., 1996), employing the "elbow" technique to find a ε value according to both the MinPts parameter and the actual dataset (Thorndike, 1953;Patel, 2019;Satopa et al., 2017). The use of this heuristic, as provided by the "dbscan" R package, can then be adjusted to find the smallest ε value possible that can define a final model for training. DBSCAN then establishes clusters within the dataset, separating noise into a separate "0" cluster and creating convex hulls through the use of border points q.

Collaborative intelligence learning
Considering the currently unquantifiable nature of numerous geological features that ML algorithms may be unable to detect, a human-in-the-loop collaborative strategy was proposed for this study (Kamar, 2016;Holzinger, 2016;Dellermann et al., 2019). The objectives of including this hybrid intelligence strategy meant that human interaction could complement the strengths of the pattern detection algorithms (Simard et al., 2017;Dellermann et al., 2019), whereby creating a bridge between both unsupervised and supervised techniques for model creation.
The Expert-in-the-Loop (EitL) approach adopted here consisted in the revision of the clustered dataset by considering possible underlying geological components that could be generating noise undetectable by the model. Components of this nature could easily include disconformities, geological faults, erosive surfaces, uneven sedimentation or even sedimentation of geological elements such as large boulders that could be separating point-densities that would essentially belong to the same layer.
The strategies employed consequently used EitL (in this case geologist-in-the-loop) interaction to manually correct clustered patterns, assessing which of the clustered groups were geologically separated and assigning these clusters to a new labeled layer. This was performed as objectively as possible through complimenting paleontological data routinely collected from the site during excavation and evaluating the nature of each of DBSCAN's groupings. EitL thus benefits through Machine Teaching (MT) in as much as the expert knowledge is used for troubleshooting and debugging (Dellermann et al., 2019; also known as a sense-making approach). Moreover, the amount of human input for models was monitored in a collective manner, using numerous domain experts to ensure accuracy (Dellermann et al., 2019). In cases where clusters seemed dubious or could not be clearly defined by either the EitL participants or the MT algorithm, these were stripped of their labels and included in the noise (cluster "0") group for objective classification by the trained ML algorithms in the supervised phase of the system. Once revised by the EitL participants, each cluster group was assigned a new label that could be used to train supervised algorithms in the following part of the workflow.

Supervised learning
Once spatial data points had been passed through the DBSCAN algorithm for initial pattern detection and then corrected by geologist-in-the-loop interventions, supervised algorithms were employed to define the final fossiliferous level models and classify those points in cluster group "0".
Considering the amount of data available, no bootstrapping procedures were deemed necessary prior to supervised training. Each algorithm was thus trained using 70:30% [train:test] splits using k-fold cross-validation (k = 10) in order to ensure the model could efficiently adjust its weights. While some authors propose the use of Spatial Cross-Validation (SCV) for studies regarding geographic data (Miller, 2004;Brenning, 2012;Pohjankukka et al., 2017;Lovelace, Nowosad & Muenchow, 2019), these algorithms were not seen to produce any significant change to the quality of results. Additionally, hyperparameter optimization was performed using a random search loop function programed in R; establishing the best hyperparameter values for each model using a combination of random values until finding the optimum settings (Bergstra & Bengio, 2012). These loop algorithms ran for 50 iterations and were then extrapolated and used for the final classification models.
Two primary supervised ML algorithms were used for the final fine tuning of the fossiliferous level models. Support Vector Machines (SVM). SVMs map out input vectors into a non-linear high dimensional feature space, using hyperplanes to calculate the degree of separation between samples (Cortes & Vapnik, 1995). In order to overcome traditional limitations imposed by linearity, a kernel function is used to define the feature space (Bishop, 2006). The constructed hyperplane can thus be used as a discriminant classifier decision surface which uses maximized margin or decision boundaries to reduce chances of overfitting (Cortes & Vapnik, 1995;Bishop, 2006). The consequent hyperplane can then be used to plot the maximized margin separations between samples and can provide a visual means of dividing strata, in an efficient and objectively computed way. For SVM applications the "e1071" R package was used. Random Forest (RF). RF can be described as a robust and highly complex form of decision tree which has proven useful in the past for the processing of spatio-temporal data (Hengl et al., 2018). The RF algorithms use small random numbers of data set variables rather than the whole dataset, constructing an independent decision tree with each subsampling (Breiman et al., 1984;Breiman, 2001). The random variable selection is additionally performed using bootstrap aggregation, using a technique frequently referred to as out-of-bag observations. RF is then able to calculate the number of iterations needed to minimize the out-of-bag error. Once the number of trees has been selected, the algorithm averages results to produce a robust classification model (Breiman, 2001). For RF applications the "caret" R package was used.
Model evaluations were performed following standardized ML protocol, evaluating the Kappa (κ), Sensitivity, Specificity and Balanced Accuracy values obtained through confusion matrices (Kuhn & Johnson, 2013;Lantz, 2013). Sensitivity, Specificity and Balanced Accuracy are values derived from evaluations of Type I and Type II statistical errors in proportion with the rest of the calculated confusion matrix. These values are presented as numbers between 0 (poor) and 1 (high performing classifiers). κ values are a further statistical adjustment of the accuracy metric that measure model agreement relative to what would be expected by chance (Kuhn, 2008). These values are similarly presented as numbers between 0 and 1 although negative values can (although rarely) occur. Further evaluation of this statistic uses 0.8 as the threshold between poor and powerful classification models. Finally, loss metrics were recorded to evaluate the predictive power of each model. For this, the Mean Squared Error (MSE) metric was employed: taking into consideration the target (t i ) and predicted (p i ) values for each classification. MSE values are further interpreted considering AIAs of this type are trained to reduce the error produced in CV sets, therefore the smaller the MSE value the more powerful the predictive model. Model evaluation was performed using the "caret" R package.

Fine tuning of fossiliferous levels
Once the optimal models were defined, these were used to classify any of the points separated as noise by DBSCAN as well as the clusters that the EitL participants could not objectively define as a part of either fossiliferous level (so as to avoid subjectivity). In order to empirically class any point as a member of a fossiliferous level, the trained SVM and RF algorithms were used to predict the class label of each point. A threshold of 80% security was established as acceptable. Any points that could not be classed by the model with over an 80% predictive decision boundary were thus rejected and therefore not classified into any level in the final fossiliferous levels models.

RESULTS
Unsupervised machine learning clustering DBSCAN clustering proved highly efficient at quickly processing each slice, grouping fossil remains into abundant clusters in both Batallones-3 (Fig. 2) and Batallones-10 (Fig. 3). DBSCAN results in most cases present very noisy profiles, with the detection of large numbers of tightly packed clusters. This is a result of the irregular density patterns detected in the data and probably results from numerous geological and paleontological disconformities that the algorithm is unable to quantify in some areas. Nevertheless, profiles become less noisy towards the extremities of each slice, helping the differentiation between levels, identifying at least three different fossiliferous levels.
On average, DBSCAN identified 166 points (5.6%) as noise in Batallones-10 and 62 (7.76%) in Batallones-3, with the highest levels of noise appearing in Batallones-10 x axis slice (n = 141) and Batallones-3 right slice (n = 56). For Batallones-3 right slice, this can be attributed to certain areas of the profile having little separation between the layers (Fig. 2C). Nevertheless, excluding noise, most clustered groups seem to follow a trend or pattern that highlights a separation between different fossiliferous levels (see "Collaborative Intelligence Learning").

Collaborative intelligence learning
While DBSCAN results produced a large number of clusters, detailed analysis through the EitL process identified numerous groupings that belonged to the same fossiliferous level. After careful evaluation, each of these separations could be recognized as products of minute or large disconformities. These can be attributed to the cavity's asymmetrical geomorphology at Batallones-10, visible in the x axis slice (Fig. 4A) or to a fallen large carbonate block in Batallones-10, as seen in the y axis slice (Fig. 4B). This highlights the importance of the geologist-in-the-loop hybrid intelligence system for clarification and unification of clusters.
In other cases, numerous examples can be seen where depositional processes have placed fossiliferous levels closer together such as the sedimentary onlap between levels II and III in Batallones-10 x axis slice (Fig. 4A) and Batallones-3 right slice (Fig. 5B). Nevertheless, considering the algorithms' effectiveness processing 3D data, DBSCAN is still able to detect density patterns across the slice's x, y and z axes, thus separating the different levels.

Supervised learning
Once processed and cleaned by both DBSCAN and the EitL specialists, databases were channeled into supervised algorithms for training. In all cases, each of the databases generated by unsupervised learning algorithms proved to be highly efficient for the training of their supervised counterparts. Both SVM and RF obtained exceptional results with over 90% accuracy when differentiating between groups and small MSE values on all accounts (Table 1). While SVM proved to have greater balanced accuracy on training sets, RF in general obtained optimal MSE values in testing, especially in the case of Batallones-10. Interestingly Batallones-3 (both left and right slices) seems to create the most amount of confusion, nevertheless on all accounts κ values remain above the acceptable threshold. Using AIAs to clean the final profiles and de-noise the datasets, both SVM and RF proved to be confident classifiers when making predictions. In most cases, over half the points denoted as noise by the previous phases were successfully classified and assigned into a palaeontostratigraphic level (Table 2). For Batallones-10 y axis slice and Batallones-3 left slice, RF outperformed SVM with greater MSE values and success in classifying indeterminable points. Nevertheless, in either case, RF tended to have greater decision making capabilities on some individual points than SVM, even though SVM obtained greater overall MSE in others.
Upon further evaluation of model performance, SVM maximized margin decision boundaries (dotted line in Fig. 6) are seen to expand or contract closer to the identified sedimentary onlap areas (in both Batallones-3 and Batallones-10) while an increase in MSE is observed for these particular areas as well (Table 2). Fine-tuned models Once the data had been processed by each of the phases of the workflow, a final fine-tuned model of the fossiliferous levels was created.
The fine-tuned model produced for Batallones-3 shows three discrete fossiliferous levels on either side (left and right) of the debris cone (Fig. 7). The left profile (Fig. 7A) is laterally less extensive (~3 m) than the right profile (Fig. 7B), which extends over 5 m. All identified levels in both profiles dip from the debris cone outwards, with dip angles decreasing towards the outermost parts of the cavity and towards the top of the infilling, with Level 3 dipping only slightly in Batallones-3 left profile (Fig. 7A) and nearly horizontal in the outer part of the Batallones-3 right profile (Fig. 7B).
Batallones-10 fine-tuned model shows three different discrete fossiliferous levels (Fig. 8). A first level (Level I), dips northward (Fig. 8A) and towards the East and West dips

DISCUSSION
The methodology proposed in this study, based on artificially intelligent systems, has been deemed effective in quantitatively identifying discrete fossiliferous levels in paleontological and archaeological sites based on the study of the spatial distribution of fossils. This new, more quantitative approach could be considered a valuable substitute for previous qualitative techniques that have been applied in the past. Evaluating variable importance through testing for mean decrease in node impurity for RF algorithms (Louppe et al., 2013) explains how a boost in the weight of the third dimension of each slice helps support the differentiation between levels (Fig. 9). Batallones-3, for example, shows a shift in the value of the x axis between the left-hand and  right-hand sides of the debris cone, having more weight in the final decision making capabilities of RF than the y dimension. Data of this nature thus supports the need for 3D over 2D analysis in order to ensure an efficient detection of patterns that 2D data may not be able to reveal. Fine-tuned models produced for Batallones-3, in the lower part of a hypothetical hourglass-shaped cavity, clearly show three discrete levels on either side of the debris cone (Fig. 7). Cave asymmetry, as explained in Calvo et al. (2013), is clearly visible in these profiles, with greater cave development towards the Northeast, in Batallones-3 right profile (Fig. 7B). Batallones-10 fine-tuned models, in the upper part of the hourglass-shaped cavities, show another three discrete fossiliferous levels (Fig. 8). Similar geological structures were observed by Calvo et al. (2013) in Batallones-9, located 50 m to the north of Batallones-10.
These newly discovered levels in both the lower levels (Batallones-3) and upper levels (Batallones-10) of the hourglass-shaped cavities of Batallones Butte site support the preliminary inferences made by Martín Escorza & Morales (2005) about the possibility of discrete fossiliferous levels within Batallones-1. Careful revision should therefore be carried out considering previous paleoecological, palaeoenvironmental and taphonomical analyses. As of this point, future studies at Batallones-3 and Batallones-10 should consider each of these individual discrete levels, instead of for the whole upper or lower part of the cavity.
As previously explained, Batallones-3 deposits and fossils are located in the lower part of the cavity and the fossil site is considered exhausted as the encasing rock was already found, for this reason, the three identified levels can be considered definitive. In the case of Batallones-10, the encasing rock has not yet been found and therefore fossil material recovered in future excavations may reveal additional discrete and separable fossiliferous layers. Thus, the increase in the number of levels at Batallones-10 cannot be discarded.
Future studies should analyze the embedded faunal assemblages separately for each of the described levels rather than the entire upper or lower part of the cavity, paying special attention to micromammal remains in order to establish chronological intra-site relationships, such as those inter-site relationships described by López-Antoñanzas et al. (2010). These detailed micromammal studies would greatly refine the geochronological framework of the Batallones Butte sites. The overall workflow of this hybrid intelligence system has proven effective in detecting the presence of 3 levels in both Batallones-10 and Batallones-3. While in some slices the detection of these levels is clearer than in others, the use of two representative, thick and potent slices for each site helps empirically support their identification. Through the implementation of unsupervised algorithms for the initial detection of patterns, expert (geologist)-in-the-loop interactions for sense making and the final tuning of profiles using supervised algorithms, the overall system can be used to find areas where identification of fossiliferous levels is clearer. This can then be extrapolated for the supervised classification of the rest of the site. Moreover, while some isolated points suspiciously appear to be classified into one layer and can be arguably considered out of place, future efforts should try to combine more data for the fine tuning of teaching and training processes. These could include taxonomic or taphonomic variables (Brain, 1981;Domínguez-Rodrigo et al., 2018), data regarding object density (Kreutzer, 1992;Lam et al., 1998Lam et al., , 2003Lam, Chen & Pearson, 1999), weight and size (Bunn et al., 1980;Brain, 1981;Bunn, 1987;Bunn & Pickering, 2010;Domínguez-Rodrigo et al., 2018), as well as orientation and plunge patterns (Woodcock, 1977;Woodcock & Naylor, 1983;Lenoble & Bertran, 2004;Domínguez-Rodrigo & García-Pérez, 2013). While models and systems without humans-in-the-loop would be an optimal solution for the future (Holzinger, 2016), numerous geological factors that need extensive investigation for modeling are required before this can be achieved. Likewise, a number of underlying depositional and post-depositional processes may be generating confusion that machines are unable to understand or process for the time being. Such is the case for the confusion generated amongst algorithms caused by sedimentary onlaps, causing fossiliferous levels to lie closer to each other, at both Batallones-3 and Batallones-10. This may, however, be solved by more complex AIAs and the use of Deep MT systems, such as clustering AIAs using auto-encoders (Xie, Girshick & Farhadi, 2016;Guo et al., 2017;Mrabah et al., 2019;Yang et al., 2019), or those used for reinforcement learning tasks (Lake et al., 2014;Mnih et al., 2015;Holzinger, 2016;Simard et al., 2017). Efforts should therefore be made to investigate the effects of these numerous geological components on pattern detection algorithms.

CONCLUSIONS
This study presents a novel use of Artificially Intelligent Systems for the quantitative identification of discrete fossiliferous levels in paleontological and archaeological sites based solely on the study of the spatial three-dimensional distribution of fossil remains. These results have been able to reveal new discoveries in the Batallones Butte site inner structure, including three discrete levels at both Batallones-3 and Batallones-10.
The two lowermost levels of Batallones-3, Level I and Level II, dip outwards from the debris cone towards the outer limits of the cave, becoming progressively more horizontal, whereas the uppermost level, Level III is sub-horizontal. On the other hand, Batallones-10's lowermost Level I dips northward into the cavity and adapts laterally to the cave limits whereas Levels II and III also show northern dips but are more influenced by cave limits and dip towards the center of the cavity in the outermost areas.
Through these discoveries, AIAs have been proved to be a highly efficient and objective means of detecting spatial patterns in paleontological sites. Unsupervised and unsupervised techniques can greatly reduce (but not eliminate) the amount of subjectivity in discrete fossiliferous level identification. The possibilities provided by the combination of these algorithms and expert-in-the-loop systems are multiple, however it is important to point out that in order for these methods to be widely used, a need for further experimentation and investigation is essential. This should include an increase in the data used to feed these AIAs. For example, the incorporation of taphonomic features of the remains in the analyzed fossil may provide an even finer tuning of these models.