Modern River-Sand Geochemical Mapping in the Manufahi Municipality and Its Surroundings, Timor-Leste: Implications for Provenance

: A geochemical mapping of regional modern river-sand is performed to clarify geological information in the study area of Timor-Leste. Several areas of Timor-Leste including the study area in particular have limited geological information due to limited accessibility and dense vegetation coverage, and deformed, weathered, and erosion-covered materials. A total of 53 modern river sand samples were collected and analyzed. Ten major elements were determined by using wavelength-dispersive X-ray fluorescence. Areas characterized by clastic sedimentary rocks are recognized clearly by elevated concentrations of SiO 2 , Al 2 O 3 , and K 2 O. Meanwhile, areas covered by carbonate sedimentary rocks are detected by significant concentrations of CaO and MnO. The occurrences of the altered clastic and carbonate sedimentary rocks of the Wailuli and Aitutu Formations due to metamorphic, silicification and other alteration processes were responsible for the elevated concentrations and positive correlation between SiO 2 , CaO, K 2 O, and MnO, and CaO, TiO 2 , and MnO in the midstream and near the downstream areas of the Clerec and Sahe River catchments. The positive correlation observed between TiO 2 , CaO and MnO may be ascribed to the presence of carbonate components associated with secondary Ti-bearing minerals, which are potentially formed through hydrothermal alteration processes.


Introduction
Timor-Leste is a relatively new independent country.It officially gained its independence in 2002.Timor-Leste is located between Southeast Asia and the Pacific region, surrounded by Indonesia's islands, with the western part of the island of Timor shared with Indonesian territory and Australia as its southern neighbor across the Timor Sea (Figure 1A).
Since Portuguese colonization, the geology of Timor-Leste has been studied and mapped.The majority of geological research studies conducted in the territory of Timor-Leste have been concerned with prospecting and discovering oil and natural gas, paleontology, tectonic settings, and stratigraphy [1][2][3][4][5][6][7][8][9][10][11][12].However, geochemical mapping studies have not yet been conducted in the majority of the Timor-Leste territory, including the study area.
Timor island is located in a zone where the Australian and the Eurasian plates collide, and it is characterized by active tectonics.This zone has led to the uplifting of land, formation of the landscape, mountain-building processes, installation of overthrust sheets, complex lithological formations, and other geological phenomena and processes within Timor Island.The terrain of the island is characterized by a rugged topography with steep slopes, large mountains, high peaks, and deep valleys.External geodynamic processes such as erosion and weathering also contribute to landscape shaping.The lithologic formations installed and deposited on Timor Island are largely composed of non-volcanic strata that range in age from Permian to Quaternary.Tectonically, these formations can be classified into allochthonous, para-autochthonous, and autochthonous units, and are similar to the Asian and Australian successions.Some of these formations have been involved in various geological events and tectonic processes, such as deformation, metamorphism, volcanic activity, mineralization, and hydrothermal alterations, before and after the orogenic phase, which occurred from the Late Jurassic to the Middle Miocene [7,9,[13][14][15][16][17][18][19][20].
In this work, the studied sediment samples are referred to as river sands [26][27][28][29][30][31].River sand samples were collected from the active channel below the water level and were dominated by fine to coarse sand fractions.Modern river sands are composite products of the weathering and erosion of rocks, minerals, and soils from upstream of the sampling site.These sediments are transported and deposited by stream water flow within a river catchment area [32][33][34][35][36].
The majority of Timor-Leste's land surface and study area in particular are covered by light to dense vegetation and weathered and eroded materials.Direct field geological mapping in these regions may also influence lithological identification capability, especially in characterizing more than one type of source or rock unit [61].This study demonstrated the potential of river sand geochemical mapping as a lithological mapping tool.It started from the research area; this study must lead to a future national project establishing a geochemical map of the country that can be utilized for environmental studies, mineral resource discovery, and lithological mapping all over the country.This study aimed to establish a geochemical database for a limited area of Timor-Leste and to discuss its characteristics and the relationship with the geology of the study area.

Geographical and Geological Setting of Timor-Leste and the Study Area
Officially known as Timor-Leste, the eastern portion of Timor encompasses the island of Atauro to the north, the islet of Jaco to the east, and the enclave of Oecusse, which is bordered by Indonesian West Timor, to the west (Figure 1A,B).
The study area is located in the municipality of Manufahi and its surroundings on the southern coast of Timor-Leste.The Timor Sea forms the southern border, while Manatuto, Aileu, and Ainaro Municipalities form the eastern, northern, and western borders, respectively (Figure 1B).The study area comprises a region of rugged mountainous topography characterized by varying slopes, ranging from almost flat to very steep.Here, it is possible to observe the activities of external geodynamic processes, such as erosion and weathering.The highest elevations are found in the eastern and northern parts of the territory; the highest peak is at Cablac Mountain, with a height of up to 2085 m.
Timor Island, the islet of Jaco, and several other islands in the southeastern region of the Indonesian archipelago are part of the Outer Banda Arc System because of their position in a complex tectonic setting, where two major plates, the Australian and Eurasian plates, collide.Timor Island is characterized by tectonic and geological features, including the predominance of Permian to Cretaceous sedimentary rocks with affinities to Australian and Asian geological sequences that have undergone significant deformation, resulting in complex structures, faulting, and folding.In contrast, the small island of Atauro is located in the Banda Sea between Timor and Wetar, is linked to the Inner Banda Arc, an area of geological interest due to its volcanic and tectonic characteristics, which are associated with subduction-related processes [2,6,10,11,19,[62][63][64][65][66][67].
There are five tectonostratigraphic units on the island of Timor [9,[68][69][70], which can be summarized as follows: (1) the Gondwana megasequence is largely composed of allochthonous and para-autochthonous sedimentary rocks from pre-rift sequences that were deposited on the northeastern shelf of the Gondwana continent, and its age ranges from the Permian to the Middle Jurassic; (2) the Australian-margin megasequence is mostly composed of para-autochthonous sedimentary rocks of the marine environment from post-rift sequences that were deposited on the plateau near the Australian mainland.Its age ranges from the Middle Jurassic to the Middle Miocene, and the initiation of the post-rift sequence deposition at 155 Ma is associated with a significant tectonic event, which is indicated by a break-up unconformity; (3) the Banda Terrane units are characterized by Cenozoic igneous rocks and their cover sediments as well as Mesozoic metamorphic basements that are exposed as thrust sheets and emplaced onto Timor-Leste during the collision phase.They are believed to be derived from the Banda Arc and its forearc arc, as well as components of the crustal blocks of Australia; (4) the Synorogenic Mélange is primarily composed of a matrix rich in scaly clay that is derived from mudrock-rich pre-and post-rift successions.Tectonic processes have transported and reassembled this sequence into a chaotic mixture; and (5) the Synorogenic megasequence is characterized by autochthonous sedimentary successions installed through uplifting processes and deposited from eroded sediments of existing rocks in the surrounding basin and depressions during the latest Miocene to the Holocene.
The five tectonostratigraphic units are well-preserved from north to south in the study area (Figure 1C), and the contact between these units is a fault and/or unconformity [1,9,14,16,63,[68][69][70][71].In the northern to central parts of the study area, the Banda Terrane and its sedimentary cover units were emplaced above the Synorogenic Mélange and Gondwana Megasequences.In the southern part, Synorogenic Megasequences were deposited above the Australian-Margin Megasequences during the latest Miocene to Holocene.

(a) Pleistocene-Holocene Suai Formation
It is mostly composed of rudites and arenites, with minor amounts of mud and marls.The particles of the formation were mostly sourced from the Viqueque Formation, Dilor Conglomerate, and Lolotoi Complex.This formation belongs to the Synorogenic Megasequence.(b) Pleistocene-Holocene Ainaro Formation (Ainaro Gravels) This formation is part of the Synorogenic Megasequence and is primarily composed of matrix-supported conglomerates that are believed to be sediments from an ancient river terrace.Occasionally, calcite lateritic cements fill these sediments, and the irregular surfaces of the river terrace sediments are frequently covered by ferruginous horizons.(c) Lower Pleistocene-Holocene Baucau Limestone (Baucau Formation) It is predominantly composed of coral reef limestones along with a minor proportion of calcarenites, calcirudites, and conglomerates (submature graywackes).This formation is expected to provide evidence of the uplift of the island.This formation belongs to the Synorogenic Megasequence.(d) Pliocene Dilor Conglomerate (Dilor Formation) This formation is part of the Synorogenic Megasequence and consists of conglomerates and sandstones with a significant contribution of detritus from the Lolotoi Complex, particularly quartzite.(e) Upper Miocene-Lower Pliocene Viqueque Formation (Synorogenic Viqueque Megasequence) Lithologically, this formation is separated into lower (described as "more clayey and silty") and upper (characterized as "more silty and sandy") sections.Large amounts of silty marls, marly siltstones, silty claystones, siltstones, and sandstones, along with minor proportions of calcilutites and biocalcarenites are present in the upper section.

Sampling Method
The sampling and sample preparation procedures were implemented in this study according to the general methodology described by Fletcher [32], Yamamoto et al. [40], Darnley et al. [80], Hale and Plant [81], Ohta et al. [82], and Tanaka et al. [83].A total of 53 modern river sand samples were collected during the dry season from ten permanent rivers or catchment areas (i.e., the Aiasa, Turon, Ermeti, Holarua, Karau Ulun, Soin, Merek, Laclo do Sul, Clerec, and Sahe Rivers) (Figure 1C).Coordinates and other field-related information, such as geology, geomorphology, and land use, were documented at each sampling point.To avoid sampling errors and heterogeneity of the analytical results, approximately 4 kg of four subsamples were collected over a distance of 20 m at each sampling point.Access conditions, topographic relief, and availability time influenced the sampling point distribution (Figure 2).River sand samples were transported to the laboratory under wet conditions for preparation and chemical analysis.To obtain grain-size fractions of 180-150 µm and <150 µm, the samples were dried in an oven at 105 °C and then sieved using a vibrating sieve shaker machine.Manual crushing using an agate mortar and pestle was followed by machine crushing using a disk mill.Subsequently, the organic matter and volatile substances were removed using the loss-on-ignition method, which was verified by weighing the weight loss after heating at 950 °C.

Analytical Method
According to Yamamoto and Morishita [84], glass beads are commonly used for sample preparation in XRF spectrometry analysis.Typically, glass beads are prepared by mixing a flux with a powdered sample.
Glass beads were prepared by fusing mixtures of powdered samples and alkali flux at a weight ratio of 1:10.An 8:2 mixture of lithium tetraborate (Li 2 B 4 O 7 ) and lithium metaborate (LiBO 2 ) was used as the alkali flux.The mixture was processed to prepare glass beads using a bead sampler machine (TK-4100 model, Amena Tech Co., Yokohama, Japan) after placing the mixed powder in a platinum crucible.At the Division of Instrumental Analysis, Gifu University, modern river sand samples were examined by wavelength-dispersive X-ray fluorescence (WDXRF).Ten major elements (SiO 2 , TiO 2 , Fe 2 O 3 , Al 2 O 3 , MnO, MgO, Na 2 O, K 2 O, P 2 O 5 , and CaO) were determined and reference samples of igneous and sedimentary rocks provided by the Geological Survey of Japan were used for calibration.

Statistical Analysis
The major element compositions to be statistically analyzed and discussed were normalized based on the estimated average concentration of elements in the Upper Continental Crust (UCC) [85] and Post-Archean Australian Shale (PAAS) [45] (Table 1).
In order to identify the enrichment or abundance of major elements in the river sand samples from the study area, the "enrichment factors" (EF) were calculated [86][87][88][89].Enrichment factors (EF) are calculated by comparing the concentrations of major elements in river sand samples with the average concentrations of Post-Archean Australian Shale (PAAS) [45].The following formula is used to determine these factors: XEF=[(X/Al)SAMPLE/(X/Al)PAAS], where X represents the element in question, and Al represents aluminum, with both expressed as their respective weight concentration values.
In this study, several ratios of selected major elements, the Pearson correlation coefficient matrix, and the index of compositional variability (ICV) were used to analyze and interpret provenance [50][51][52][53][90][91][92][93][94].The ICV was calculated as weight percent (wt%) [95].To avoid mistakes in the analysis and interpretation of ICV values due to carbonate materials, which are abundant in the majority of sediments and sedimentary rocks [96], as well as the carbonate rocks that contribute to the CaO enrichment in river sand samples, the CaO * value represents Ca, which is only present in silicate minerals.According to the concept that the molecular CaO/Na 2 O ratio of silicate is not greater than one, the CaO * value is considered to be equivalent to the Na 2 O content [46,78].
To simplify and evaluate the results of the geochemical analysis by reducing the dimensionality of the datasets, exploring elemental associations, and allowing the interpretation of variance within the dataset based on major controlling factors, a multivariate statistical technique, principal component analysis (PCA) [97], was used.To analyze and interpret PCA results in a biplot of joint graphical representations of variables and samples, only principal components with eigenvalues greater than or near 1.0 were used because they explain the majority of the variance in the data [98,99].
Chemical and statistical analyses were performed using Microsoft Office Excel version 2405 (Build 17628.20144)and JMP Pro 14.0.The maps presented in this work were created using ArcGIS 10.4.
Table 1.WDXRF geochemical results of modern river sands from the study area, including some selected ratios, ICV, and descriptive statistics of average and standard deviation for each basin.Post-Archean Australian Shales (PAAS) [45] and upper continental crust (UCC) [85] values are also presented for comparison.

Geochemical Features of River Sand Samples
The chemical analysis and statistical summary results of river sand samples from ten drainage basins in the study area are shown in (Table 1).The major elements of the river sand samples from the study area reveal considerable variations.The SiO respectively.The majority of these highest concentration values were observed near the midstream areas (also shown in Figure 3).In the Laclo do Sul River catchment, the highest concentrations of Al 2 O 3 and MgO were reported (can also be seen in Figure 3), the measured average values were 17.23 wt% (ranging from 16.40 to 17.90 wt%) and 4.41 wt% (ranging from 4.30 to 4.70 wt%), respectively.The highest concentrations of K 2 O and P 2 O 5 were reported in the upstream areas of the Aiasa and Soin River catchments (see also Figure 3), with average values measured at 1.94 wt% (ranging from 1.64 to 2.15 wt%) and 0.22 wt% (ranging from 0.20 to 0.27 wt%), respectively.In addition, the concentrations of CaO in the Turon, Sahe, and Holarua River catchments were measured to be the highest, while the concentrations of all other major elements were low (see also Figure 3).In these catchment regions, the average CaO concentrations were measured to be 20.22 wt% (ranging from 16.19 to 24.25 wt%), 13.02 wt% (ranging from 8.20 to 23.49 wt%), and 18.11 wt% (ranging from 13.43 to 23.21 wt%), respectively.

Comparison with UCC and PAAS
The results of the UCC-and PAAS-normalized patterns for the major elements of the river sand samples from the study area are shown in Figure 4.The bulk composition of the UCC is granodioritic and representative of the continental crust [85].In contrast, PAAS represents a crustal shale sedimentary rock bulk composition [45].In comparison with the average major element compositions of the UCC and PAAS (Figure 4A,B), the composition of the river sand samples from the ten drainage basins revealed a significant depletion in K 2 O, a slight depletion in SiO 2 , a slight depletion to considerable enrichment in MgO, a slight enrichment in P 2 O 5 , slight to considerable enrichments in TiO 2 , Fe 2 O 3 , and MnO, and considerable to significant enrichments in CaO.However, Turon and Holarua had more abundance and Laclo do Sul and Clerec had lower abundance than in the other river catchments.The UCC-normalized patterns (Figure 4A) showed considerable depletion of Na 2 O, and the concentration of Al 2 O 3 was enriched in the Aiasa, Laclo do Sul, and Clerec drainage basins and depleted in the other river catchments.In contrast, compared to the average composition of the PAAS (Figure 4B), Na 2 O was enriched and Al 2 O 3 was depleted.

Discussion
The results of the bulk geochemical composition analysis of the river sand samples from the study area showed an abundance of CaO, a slight depletion in SiO 2 and Al 2 O 3 contents, as well as depletion of K 2 O, and variations in Na 2 O concentrations.These findings suggested that geochemical composition of the river sands from study area were mostly controlled by the destruction of carbonate, clay and plagioclase minerals, which were related to the major contributions from the carbonate sedimentary rocks of the northern edge of the Australian continental shelf and from carbonate components that are integrated into clastic sedimentary rocks in different proportions [96,100,101].

The Enrichment of the Elements
River sand samples were collected from the ten drainage basins in the study area, which are dominated by fine to coarse sand fractions.In order to evaluate elemental mobility during weathering and transportation, the average concentrations of major elements in the analyzed samples were compared to the mean compositions of a sandstone reported by Turekian and Wedepohl [102].The results of the comparison showed that SiO 2 was depleted, while other major elements were enriched.
The enrichment of Al 2 O 3 concentration in river sand samples can be attributed to significant contributions from clay minerals, particularly those associated with fine-grained sedimentary rocks such as clay and silt [95,103,104].The results of the Pearson correlation coefficient matrix for major elements (Table 2) showed a positive correlation between most of the major elements and Al 2 O 3 in most of the study area's river basins.This suggests that clay minerals are also responsible for the distribution and enrichment of most major elements in the river sands from the study area.It is compatible with the results and conclusions of Hossain et al. [54], Biswas et al. [90], and Cai et al. [105].
The enrichment factor (EF) values have been established to evaluate the leaching and enrichment degrees of major elements in river sand samples.The values higher than 1 indicate enrichment, while values lower than 1 disclose depletion [88,89].Moreover, enrichment may be detected with EF values more than 3, while moderate to significant enrichment is revealed by EF values greater than 10 [86,87].The average calculations of the EF with respect to the mean concentrations of the PAAS (Table 3 and Figure 5) indicated that CaO exhibited detectable to moderate enrichment, K 2 O showed depletion, and TiO 2 , Fe 2 O 3 , MnO, MgO, Na 2 O, and P 2 O 5 demonstrated enrichment across all drainage basins in the study area.In addition, SiO 2 displayed slight depletion in the Aiasa and Laclo do Sul catchment areas, while showed slight enrichment in the other river basins.
CaO, Na 2 O, MgO and K 2 O are widely known to be mobile major elements due to their mobility and leachability during weathering processes, and they are frequently depleted in river sands [106][107][108][109].The lithology of the source rocks has a considerable influence on the enrichment of these major elements in the river sand samples from the study area.SiO 2 exhibited high concentrations in the samples were collected in the areas drained from clastic sedimentary bedrock lithologies, which are associated with Wailuli, Viqueque, Ainaro and Suai Formations (Figure 3).The Na 2 O/TiO 2 ratio (proposed by Yamamoto et al. [40]) values in the study area ranged from 0.79 to 2.47, with an average of 1.77, which is slightly higher than PAAS but lower than UCC (Table 1).The lowest Na 2 O/TiO 2 ratio values, which were slightly lower than PAAS (less than 1.20), and the highest TiO 2 values were mostly recorded in samples drained from clastic sedimentary rock source areas affiliated with the Wailuli Formation (Figure 3).On the other hand, the highest Na 2 O values were mostly observed in samples derived from metamorphic and igneous rocks affiliated with the Lolotoi Complex, Dilor Conglomerate, and Bobonaro Complex (as exotic blocks) that lie upstream.Additionally, the Pearson correlation coefficient matrix for most river basins of the study area showed that TiO 2 and Na 2 O were positively correlated with SiO 2 , Al 2 O 3 , Fe 2 O 3 , MnO, MgO, and P 2 O 5 .The enrichment of SiO 2 and TiO 2 in the study area could be attributed to significant contributions from quartz and clay minerals.The notable contributions from chlorite, biotite, amphibole, pyroxene, and garnet were likely responsible for the enrichment of Na 2 O in the study area.These minerals also contribute to the abundance of TiO 2 .These findings are consistent with the result and discussion of Yamamoto et al. [40] from a study of geochemical mapping in Aichi prefecture.
The results of the Pearson correlation for major elements (Table 2) showed a positive association between most major elements (particularly SiO 2 and Al 2 O 3 ) and CaO in most of the study area's drainage basins.Moreover, CaO revealed high concentrations in the samples mainly drained from carbonate sedimentary rock source areas, which are affiliated with Aitutu, Wailuli, and Viqueque Formations, as well as Dartollu Limestone (Figure 3).The abundance of CaO in the study area are possibly ascribed to the significant presence of carbonate minerals associated with carbonate and clastic sedimentary rocks in varying proportions.The highest concentrations of K 2 O were mostly observed in the samples derived from clastic sedimentary bedrock source areas, particularly those associated with the Wailuli and Viqueque Formations (Figure 3).In most of the river basins in the study area, K 2 O showed a positive association with SiO 2 , Al 2 O 3 , and Na 2 O, and a negative correlation with CaO and MnO (Table 2).These findings suggest that the distribution of K 2 O in the study area may have been controlled by the destruction of feldspar, mica, and clay minerals.2).This suggests that these elements are mostly accumulated in finer grain sizes, like those of clay minerals.These findings are compatible with the results and discussion of Lim et al. [110] and von Eynatten et al. [111].However, the distribution patterns of enrichment areas also showed some modest similarities between TiO 2 , Fe 2 O 3 , MnO, P 2 O 5 , and MgO (Figure 3).These elements also exhibited positive correlations with one another across the majority of the drainage basins (Table 2).This suggests that the enrichment of these elements could be related to the contribution of mafic, heavy, and accessory minerals (e.g., amphibole, pyroxene, biotite, ilmenite, hematite, sphene, rutile, and garnet), which are mostly associated with igneous and metamorphic rocks.These findings are consistent with the results and discussions of Biswas et al. [90], Armstrong-Altrin [112], and Armstrong-Altrin et al. [113].

Geochemical Characteristics of the Drainage Basins (1) Aiasa River catchment
The Aiasa River catchment is covered by the Lolotoi and Bobonaro Complexes, as well as the Viqueque and Ainaro Formations (Figures 1C and 3).This catchment area showed the highest concentrations of K 2 O and P 2 O 5 , with average values of 1.94 wt% (ranging from 1.64 to 2.15 wt%) and 0.22 wt% (ranging from 0.20 to 0.26 wt%), respectively (Table 1).The highest values were registered in the area near the upstream (as shown in Figure 3).12 wt% (ranging from 0.09 to 0.15 wt%), 1.17 wt% (ranging from 0.99 to 1.36 wt%), and 0.17 wt% (ranging from 0.15 to 0.18 wt%), respectively (Table 1).The highest and lowest values were recorded in the midstream area (Figure 3).The abundance of CaO might be sourced from carbonate components, which contributed much more than others, while other major elements could be associated with clay and silicate minerals.

(3) Ermeti River catchment
The Ermeti River catchment is mostly covered by the Ainaro Formation (Figures 1C and 3).This catchment area showed high content of SiO , along with low concentrations of CaO and MnO, were mostly recorded near the upstream areas (Figure 3).The midstream regions had low SiO 2 content and high concentrations of CaO, MnO, and Fe 2 O 3 .In addition, the downstream regions had high SiO 2 , CaO, and MnO concentrations, and low content values of P 2 O 5 , MgO, TiO 2 , Fe 2 O 3 , and Al 2 O 3 .The upstream areas could be influenced by clay and mica content as well as quartz, plagioclase, amphibole, chlorite, hematite, ilmenite, rutile, garnet, sphene, and apatite minerals, which are associated with clastic sedimentary and metamorphic rocks.Carbonate components and Mn minerals affiliated with carbonate sedimentary rocks, along with the clay content associated with clastic sedimentary rocks (e.g., shales and mudstones), made major contributions to the midstream regions.The downstream regions were mostly controlled by silicate minerals, such as quartz associated with clastic sedimentary rocks; carbonate components and manganese minerals affiliated with carbonate sedimentary rocks; along with clay content related to clastic sedimentary rocks.(7) Laclo do Sul River catchment The Laclo do Sul River catchment is covered by the Lolotoi Complex, Wailuli Formation, and Bobonaro Complex, as well as the Viqueque, Ainaro, and Suai Formations (Figures 1C and 3).The lowest content of CaO and highest values of Al 2 O 3 , MgO, Fe 2 O 3 , Na 2 O, and P 2 O 5 were registered near the upstream and downstream areas (Figure 3).The average values of the lowest and highest concentrations were measured to be 6.3).The destruction of clay and mica, along with quartz, plagioclase, amphibole, chlorite, hematite, ilmenite, rutile, garnet, sphene, and apatite minerals associated with clastic sedimentary and metamorphic rocks, contributed significantly to the composition of the river sands near the upstream regions.Near the downstream and continuing to the midstream regions, quartz and clay contents affiliated with clastic sedimentary rocks (e.g., shales and mudstones) made major contributions; however, in the sampling location SS045, there was significant inputs from mafic and heavy minerals related to the presence of igneous and metamorphic rocks that associated with the Bobonaro Complex, as exotic blocks that were incorporated into the clastic sedimentary rocks (such as shales and mudstones) [9,13,76].Silicate minerals and clay related to clastic sedimentary rocks contributed significantly to the downstream area.( 8) Merek River catchment The Merek River catchment is mostly covered by the Wailuli Formation and Lolotoi Complex (Figures 1C and 3).The highest and lowest concentrations of SiO 2 and TiO 2 were also recorded in this river catchment, with the average values reported at 55.76 wt% (ranging from 53.33 to 58.19 wt%) and 0.86 wt% (ranging from 0.78 to 0.94 wt%), respectively (Table 1).CaO and K 2 O appeared to increase toward downstream, while the other major element contents showed an decreasing trend toward downstream (Figure 3).Silicate (such as quartz, plagioclase, amphibole, and pyroxene), and clay minerals might significantly contribute to the geochemical composition of the river sand samples from this river catchment.The remaining contributions could be attributed to the presence of carbonate components.(9) Clerec River catchment The Clerec River catchment was registered as having the highest concentrations of SiO 2 , TiO 2 , Fe 2 O 3 , MnO, Na 2 O, and MgO (Table 1).This catchment area is mostly covered by the Lolotoi and Bobonaro Complexes, as well as the Aitutu, Wailuli, Viqueque, and Suai Formations (Figures 1C and 3).As demonstrated by the Pearson correlation of the Clerec River catchment (Table 2C), SiO 2 showed a positive association with CaO and K 2 O, Al 2 O 3 had a positive correlation with TiO 2 , Fe 2 O 3 , Na 2 O, P 2 O 5 , and MgO, and CaO showed a positive correlation with SiO 2 , MnO, and K 2 O.There was a moderate-to-very strong positive relationship between TiO 2 , Fe 2 O 3 , MgO, Na 2 O, and P 2 O 5 , indicating that the destruction of clay, mica, amphibole, pyroxene, biotite, ilmenite, hematite, sphene, rutile, garnet, and apatite, along with calcium carbonate, manganese, and alteration minerals, contributed significantly to the composition of the river sand from this river catchment.The positive correlations between SiO 2 , CaO, and K 2 O indicated that they may have been derived from the same source.The elevated K 2 O level may be related to the presence of secondary K-bearing minerals.The CaO content tended to increase downstream; however, the highest measured concentration in this river catchment was recorded at the sample site SS049.In the midstream regions, high contents of Al This indicate that the composition of river sand sample near the upstream area was mostly controlled by the destruction of quartz, clay, and mica minerals, which were affiliated with clastic sedimentary rocks; however, silicate and accessory minerals (such as quartz, muscovite and garnet minerals), which are associated with their altered rocks due to metamorphic processes, contributed to the sample location SS013.Sample sites SS014 and SS015 had major contributions from calcium-carbonate minerals associated with carbonate sedimentary rocks of the Aitutu Formation and their altered rocks due to silicification or certain alteration processes [13,72,114].There were notable inputs from quartz and secondary K-bearing minerals affiliated with altered carbonate sedimentary rocks intercalated with clastic sedimentary (e.g., shales) strata due to silicification processes in the sample location SS014, and considerable contributions to SS015 came from calcium carbonate, manganese, and Ti-bearing minerals, which were also associated with altered carbonate sedimentary rocks intercalated with clastic sedimentary (e.g., shales) strata due to certain alteration processes.The downstream regions were characterized by contributions from quartz, clay, carbonate, mica, amphibole, chlorite, ilmenite, garnet, and other minerals related to clastic sedimentary rocks.

Relationship between River Sand Geochemistry and Provenance Geology
The ratios of SiO 2 /Al 2 O 3 and K 2 O/Na 2 O, classification diagram of sediment using SiO 2 /Al 2 O 3 versus K 2 O/Na 2 O, and the geochemical index of ICV, are also frequently used to evaluate the abundance of quartz, clay, K-feldspar, plagioclase, and other less resistant minerals to identify provenance [50,52,90,92,94,[115][116][117][118].The SiO 2 /Al 2 O 3 ratio has different average values for determining source rock compositions and sediment maturities in clastic sediments; in basic igneous rocks, the ratio is approximately 3 and the ratio is approximately 5 in acidic igneous rocks.Hereafter, in clastic sediments, values >5 indicate compositionally mature sediments and values >10 indicate high maturity and potential recycling [53,92,94,117,119].
In the study area, the values of SiO 2 /Al 2 O 3 ratios ranged from 2.93 to 4.58, with an average of 3.59, slightly higher than PAAS but slightly lower than UCC (Table 1), suggesting that river sand compositions were not as silica and aluminum-rich as the average composition of typical crustal shale and granodiorite.The negative associations between the SiO 2 /Al 2 O 3 ratio and elements such as TiO 2 , MgO, Fe 2 O 3 , Na 2 O, and P 2 O 5 (Table 2A-D) in most of the selected drainage basins, suggesting that these elements were not highly associated with silicate minerals but were incorporated into clay minerals.On the other hand, these elements showed a positive correlation with the SiO 2 /Al 2 O 3 ratio in the Sahe River catchment.This finding might also suggest that the geochemical composition of the river sand samples from this river catchment is significantly influenced by mineral contributions related to the presence of the Aitutu Formation's clastic sedimentary rocks and their altered rocks as a result of silicification or other alteration processes.
The K 2 O/Na 2 O ratio also has variational values in defining the abundance of K-feldspar relative to plagioclase; with high values (>1) suggesting potential contributions from alkali feldspar or mica and/or possibly from secondary K-bearing minerals, whereas for the abundance of sodium-bearing minerals or plagioclase, the ratios were <1 [53,90,93,94].The K 2 O/Na 2 O ratios in the study area ranged from 0.41 to 1.74, with an average of 0.73, which was much lower than those of UCC and PAAS (Table 1).Compared with the UCC, the K 2 Oand Na 2 O concentrations showed greater depletion, which could be related to the low contribution of feldspathic minerals [53,100,101].These results indicated that the river sand content in the study area was also characterized by an abundance of aluminosilicate compared to feldspar.The highest values of the K 2 O/Na 2 O ratio (>1) were identified at several sampling points (Figure 6), which were mostly characterized by notable inputs from aluminosilicate and accessory minerals such as clays and muscovites.These minerals were mainly sourced from clastic sedimentary and metamorphic rocks rich in accessory minerals (e.g., muscovites), as well as shales, micaceous shales, micaceous sandstones, and their altered forms owing to metamorphic processes, which were affiliated with the Viqueque Formation, Dilor Conglomerate, and Wailuli Formation, respectively.However, sample locations SS014, and SS015 revealed features indicative of the presence of secondary K-bearing minerals and secondary Ti-bearing minerals associated with carbonate sedimentary intercalated with clastic sedimentary rocks that have undergone silicification and other alteration processes.These altered rocks are associated with the Aitutu Formations.The geochemical classification binary diagram using SiO 2 /Al 2 O 3 against K 2 O/Na 2 O (Figure 7), proposed by Wimmenauer [118], showed that most of the river sand samples from the study area fall within the "pelitic greywacke" field, and only 4 samples stray into the "greywacke" field.This finding indicated that clastic sedimentary rocks (pelitic greywackes) also contributed significantly to the geochemical composition of most of the river sands in the study area.Cox et al. [95] also determined ICV values greater than 1 for the sediments rich in non-clay minerals (chlorites, feldspars, amphiboles, and pyroxenes), whereas values less than 1 were found for sediments with significant concentrations of clay minerals (kaolinite, illite and montmorillonite) [50,53,90,92,94].The ICV values in the study area ranged from 0.93 to 1.55, with an average of 1.26, which was much higher than those of PAAS and UCC (Table 1).The enrichments of TiO 2 , Fe 2 O 3 , MnO, MgO, and P 2 O 5 (Table 3 and Figure 5), as well as the positive association between these elements in most drainage basins (Table 2), corroborate the presence of mafic sources containing heavy and accessory minerals (such as amphibole, pyroxene, biotite, ilmenite, hematite, sphene, rutile, and garnet) [90,112].These results revealed the abundance of chlorite minerals along with mafic, heavy, and accessory minerals in river sands from the study area.The occurrence of basic igneous and metamorphic rocks is mainly responsible for the contribution of these minerals to the research area.These rocks are affiliated with the Lolotoi Complex; however, they can be integrated into sedimentary rock fragments and matrices associated with the Aitutu, Wailuli, and Viqueque Formations and/or incorporated into mudstones and shales within the Bobonaro Complex as exotic blocks.

Principal Component Analysis
PCA was conducted on the major elemental compositions of stream sediments from the ten river catchments in the study area to determine their provenance.The first three principal components of the samples and variables with eigenvalues greater than or near 1.0 were processed for analysis and interpretation (Figure 8) [98,99].The findings of the PCA were consistent with the results of the SiO 2 /Al 2 O 3 and K 2 O/Na 2 O ratios as well as those of the ICV.PC2 was strongly and positively associated with TiO 2 (0.40), Fe 2 O 3 (0.37), CaO (0.28), and MgO (0.25), particularly in the following samples: SS04 (3.64), SS050 (2.64), SS045 (2.35), SS015 (2.06), SS019 (1.60), and SS025 (1.58).Input from SiO 2 -poor rocks may explain the positive values of PC2.The elemental association of the positive values of PC2, possibly related to the contributions from mafic, heavy, and accessory mineral (such as amphiboles, pyroxenes, biotites, ilmenites, hematites, titanites, rutile, garnet, apatite, and others) inputs, is mainly attributed to the occurrence of igneous and metamorphic rocks in the study area.These rocks are associated with the Lolotoi Complex, but they may also be integrated into clastic sedimentary rock fragments and matrices associated with the Aitutu, Wailuli, Viqueque, Ainaro, and Suai Formations as well as the Dilor Conglomerate.These rocks may have also been incorporated into clastic sedimentary rocks (e.g., mudstones and shales) as exotic blocks within the Bobonaro Complex.Component 2 also had strong negative loadings for K 2 O (−0.56) and SiO 2 (−0.47).Negative PC2 scores were observed for samples SS049 (−3.84),SS012 (−3.65),SS043 (−2.27),SS039 (−2.25),SS014 (−2.06), and SS09 (−1.98).SiO 2 -rich rocks may have been responsible for the negative PC2 scores.The inputs of mica (muscovite) and quartz minerals may have contributed to the higher K 2 O and SiO 2 concentrations, which consequently influenced the negative values of PC2.In the study area, these minerals are mostly ascribed to the existence of clastic sedimentary and metamorphic rocks (e.g., shales, micaceous shales and sandstones, conglomerates, sandstones, siltstones, claystones, mudstones, quartzites, and mica schists) of the Aitutu, Wailuli, Viqueque, and Ainaro Formations, Lolotoi Complex, and Dilor Conglomerate, as well as exotic blocks of these rocks from the Lolotoi Complex, Aitutu Formation, and Wailuli Formation within the Bobonaro Complex.

Conclusions
This study showed the characteristics of the geochemical composition analysis and their relationship to the underlying lithology in the study area with limited geological information because of limited accessibility and dense vegetation coverage, as well as deformed, weathered, and erosion-covered materials.The areas covered by clastic sedimentary rocks are recognized higher concentrations of SiO 2 , Al 2 O 3 , and K 2 O.Meanwhile, the area covered by carbonate sedimentary rocks are detected by higher CaO and MnO.This study also contributed significantly to the discovery of altered clastic and carbonate sedimentary rocks within the Wailuli and Aitutu Formation in the midstream and near the downstream areas of the Clerec and Sahe drainage basins.These findings were identified by the increased contents and positive correlation between SiO 2 , CaO, K 2 O, and MnO, and CaO, TiO 2 , and MnO in the areas.The contribution of quartz, accessory (e.g., muscovite and garnet), secondary K-bearing, and manganese minerals associated with clastic and carbonate sedimentary rocks, which may have undergone metamorphic and silicification processes, could account for the presence of SiO 2 , CaO, K 2 O, and MnO.The positive correlation observed between TiO 2 , CaO and MnO may be ascribed to the presence of carbonate components associated with secondary Ti-bearing minerals, which are potentially formed through hydrothermal alteration processes.The occurrence of altered rocks and potential mineral resources associated with the Wailuli and Aitutu formations has never been reported or discovered in the surrounding areas of the Clerec and Sahe River catchments.
However, in another areas of Timor-Leste and the Indonesian west Timor region, Audley-Charles [13] and Barkham [114] have reported the existence of altered carbonate sedimentary rocks of the Aitutu Formation by silicification processes.Meanwhile, Vicente et al. [120] have presented the discovery of hydrothermal mineralization associated with Gondwana Megasequence in the Maquelab area (Oecusse -Timor-Leste).Further regional geochemical investigations in the study area are warranted, and the findings of this study could have implications for future lithological mapping research.Future geochemical mapping investigations are advised to focus on the Laclo do Sul, Clerec, and Sahe River catchments, which are covered by Wailuli and Aitutu Formations.In order to provide valuable information in mapping the distribution patterns of altered rocks and identifying potential mineral resources.

Figure 3 .
Figure 3. Major element distribution in the study area on the geological map [24,25] with river sand sampling points.

Figure 5 .
Figure 5. Spider plots of the enrichment factor calculation with respect to the average concentrations of Post-Archean Australian Shale (PAAS)[45] for major elements in river sands from the study area.The scattered areas of enrichment observed in the distribution patterns of Al 2 O 3 , Fe 2 O 3 , and MgO are quite similar to each other (Figure3).Additionally, a positive relationship was observed between Al 2 O 3 , Fe 2 O 3 , and MgO in most of the study area's drainage basins (Table2).This suggests that these elements are mostly accumulated in finer grain sizes, like those of clay minerals.These findings are compatible with the results and discussion of Lim et al.[110] and von Eynatten et al.[111].However, the distribution patterns of enrichment areas also showed some modest similarities between TiO 2 , Fe 2 O 3 , MnO, P 2 O 5 , and MgO (Figure3).These elements also exhibited positive correlations with one another across the majority of the drainage basins (Table2).This suggests that the enrichment of these elements could be related to the contribution of mafic, heavy, and accessory minerals (e.g., amphibole, pyroxene, biotite, ilmenite, hematite, sphene, rutile, and garnet), which are mostly associated with igneous and metamorphic rocks.These findings are consistent with the results and discussions of Biswas et al.[90], Armstrong-Altrin[112], and Armstrong-Altrin et al.[113].

Figure 6 .
Figure 6.Geological map [24,25] and spatial distribution of K 2 O/Na 2 O ratio in the study area.

Figure 7 .
Figure 7. Geochemical classification binary diagram of SiO 2 /Al 2 O 3 vs.K 2 O/Na 2 O for river sand samples from the study area (after [118]).

Figure 9 .
Figure 9.Samples (A) and variables (B) plotting on the first two axes of the PCA.

Figure 10 .
Figure 10.Samples (A) and variables (B) plotting on the first and third axes of the PCA.
The lower section is mostly composed of marls, clayey marls, silty marls, claystones, silty claystones, calcilutites, and tuffs, with minor amounts of basal conglomerates and mottled marls.These rocks are mostly formed by foraminifera fossils and skeletal radiolarians as well as rock fragments and mineral particles, which are associated This unit is part of the Synorogenic Mélange and is primarily composed of exotic blocks within a scaly clay matrix.The matrix lithology is similar to that of mudstone of the Wailuli Formation.Exotic blocks of Permian to Cretaceous ages are common and widely distributed, although absent in several areas.(g)Oligocene-Miocene Cablac Limestone (Cablac Formation)This formation is largely composed of oolitic and peloidal limestones and pelagic carbonates with small amounts of intraformational conglomerates, calcilutites, calcarenites, agglomerates, and tuffaceous rocks.The most common clasts of conglomerates are volcanic rocks, calcilutites containing foraminifera, radiolarian chert, biomicarenites, and detrital minerals, such as quartz and magnetite.Several rocks have been affected by alteration processes such as dolomitization and a few have undergone partial silicification, desilicification, and dedolomitization.This formation belongs to the Banda Terrane unit.(h) Oligocene Barique Formation This formation belongs to the Banda Terrane unit and is largely composed of mafic to acidic lavas and tuffs, with minor amounts of serpentinites, volcanic conglomerates, and sandstones.Pillow lava is also observed.Significant alterations are observed in most volcanic rocks.Volcanic rocks are considered to have formed at mid-oceanic ridges and volcanic arcs.
with carbonate, metamorphic, volcanic, and other sedimentary rocks found in the Lolotoi Complex, Maubisse Formation, Aitutu Formation, Wailuli Formation, and Bobonaro Complex.This formation is included in the Synorogenic Megasequence.(f) Middle Miocene Bobonaro Complex (Bobonaro Formation or Bobonaro Scaly-Clay or Bobonaro Mélange or Synorogenic Mélange) (i) Middle-Upper Eocene Dartollu Limestone (Dartolu Formation) This limestone is part of the Banda Terrane unit and its primary constituents are algal and alveolina biomicarenites with minor proportions of calcilutites, siliceous shales, and siltstones.Dolomitization or silicification was not observed.(j) Lower-Upper Cretaceous Waibua Formation This formation is part of the Australian-Margin Megasequence and is largely composed of radiolarites, radiolarian cherts, marls, and shales, with several percentages of calcilutites, marls, and calcarenites.Radiolaria and pelagic foraminifera are important components of these rocks and most limestones are completely or partially silicified.Radiolarian shales, marls, and radiolarites often occur in association with Mn nodules and ferromanganiferous rocks.The formation processes of radiolarites and cherts are closely associated with Mn-rich strata.(k) Late Triassic-Middle Jurassic Wailuli Formation This formation predominantly consists of gray shales and blue-gray marls with minor amounts of sandstones, mudstones, quartz-arenites, coarse polymictic conglomerates, calcarenites, and calcilutites.Most shales are composed of fine micaceous

Table 2 .
Pearson correlation coefficient between major elements and SiO 2 /Al 2 O 3 ratio in river sand samples from selected drainage basin within the study area.(A) Soin; (B) Laclo do Sul; (C) Clerec; and (D) Sahe.The positive correlations of SiO 2 , Al 2 O 3 , and CaO with other major elements, as well as the positive correlation observed between TiO 2 , Fe 2 O 3 , MgO, and other major elements, are denoted by boldface type.Emphasized values with shaded cells indicate the negative associations between TiO 2 , MgO, Fe 2 O 3 , and P 2 O 5 with the SiO 2 /Al 2 O 3 ratio.

Table 3 .
[45]istical results of enrichment factor (EF) calculated with respect to average concentrations of the Post-Archean Australian Shale (PAAS)[45]for river sands from the study area.
The concentrations of Al 2 O 3 , K 2 O, TiO 2 , Fe 2 O 3 , MgO, and P 2 O 5 tended to decrease toward downstream, and on the contrary, CaO and Na 2 O seemed to increase downstream.However, the high content value of Na 2 O was recorded near the upstream area.The SiO 2 and MnO concentration values showed an irregular distribution pattern from upstream to downstream.The presence of CaO and Na 2 O might be attributed to the presence of the calcite and halite minerals associated with carbonate components.Near the upstream area, Na 2 O could be partly attributed to the existence of silicate and aluminosilicate minerals such as plagioclases, clays, and micas.Meanwhile, TiO 2 , Fe 2 O 3 , MgO, K 2 O, and P 2 O 5 concentrations were mostly concentrated in clay minerals such as illite and others.In the upstream area, the composition of river sand is mostly controlled by clay contents and lithic fragments associated with clastic sedimentary rocks.On the other hand, carbonate components related to the carbonate sedimentary rocks make major contributions in the downstream region.(2)Turon River catchment The highest content of CaO and the lowest concentrations of SiO 2 , TiO 2 , Al 2 O 3 , Fe 2 O 3 , MnO, K 2 O, and P 2 O 5 were reported in this river catchment (Table 1), which is mostly covered by the Viqueque and Ainaro Formations (Figures 1C and 3).The average value of the highest concentration was measured to be 20.22 wt% (ranging from 16.19 to 24.25 wt%), and the lowest contents were 43.15 wt% (ranging from 39.00 to 47.30 wt%), 0.78 wt% (ranging from 0.67 to 0.88 wt%), 11.35 wt% (ranging from 10.40 to 12.30 wt%), 5.33 wt% (ranging from 4.75 to 5.91 wt%), 0.

Table 1 )
2 and low CaO, P 2 O 5 , and MnO concentrations, with average values of 51.36 wt% (ranging from 47.70 to 55.01 wt%), 10.61 wt% (ranging from 6.99 to 14.22 wt%), 0.16 wt% (ranging from 0.15 to 0.18 wt%), and 0.16 wt% (ranging from 0.14 to 0.18 wt%) (Table1).SiO 2 , Al 2 O 3 , Fe 2 O 3 , MgO, Na 2 O, and K 2 O contents appeared to increase in the upstream, while the concentrations of TiO 2 , MnO, CaO, and P 2 O 5 tended to decrease toward the downstream (Figure3).TiO 2 , MnO, CaO, and P 2 O 5 could be attributed to the contribution of carbonate components along with manganese and apatite minerals, while clay and silicate minerals could be responsible for SiO 2 , Al 2 O 3 , Fe 2 O 3 , MgO, Na 2 O, and K 2 O. (4) Holarua River catchment The Holarua River catchment is covered by the Lolotoi Complex, Wailuli Formation, and Ainaro Formation (Figures 1C and 3).This drainage basin also recorded the highest concentration of CaO and the lowest contents of SiO 2 , Al 2 O 3 , MnO, Na 2 O, and K 2 O, with the average values reported at 18.11 wt% (ranging from 13.43 to 23.21 wt%), .TiO 2 , MnO, Na 2 O, and P 2 O 5 contents showed an increasing trend toward upstream, and SiO 2 , K 2 O, and Al 2 O 3 concentrations appeared to increase toward downstream, but their high content values were observed close to the upstream area.Although the concentrations of CaO, Fe 2 O 3 , and MgO seemed to decrease toward downstream, their low measured values were found close to the upstream area (Figure 3).SiO 2 , K 2 O, and Al 2 O 3 could be derived from the destruction of feldspar and/or muscovite minerals.TiO 2 , MnO, Na 2 O, and P 2 O 5 were mostly accumulated in clay minerals, and these elements can be partly attributed to the presence of silicate and accessory minerals.The abundance of CaO, Fe 2 O 3 , and MgO could be sourced from calcite, dolomite, and pyrite minerals associated with carbonate components.The upstream and downstream areas are mostly influenced by clay and mica contents, as well as lithic fragments affiliated with clastic sedimentary rocks, which contribute to the composition of river sands.However, there were contributions from silicate and accessory minerals associated with metamorphic rocks in the upstream area.The midstream region is characterized by significant contributions from carbonate components, particularly those associated with carbonate sedimentary rocks.The concentration values of CaO almost showed a decreasing trend toward downstream, and Al 2 O 3 contents also nearly appeared to be increasing in the upstream direction.Near the upstream region, mostly high values of SiO 2 , TiO 2 , Fe 2 O 3 , MnO, MgO, Na 2 O, K 2 O, and P 2 O 5 concentrations were recorded.A high content of CaO and low concentrations of TiO 2 , Al 2 O 3 , K 2 O, P 2 O 5 , Fe 2 O 3 , MnO, and Na 2 O were reported in the downstream region.However, low Al 2 O 3 content and high concentration values of SiO 2 , K 2 O, MnO, P 2 O 5 , and MgO were recorded in the sampling location SS09 (Figure 3).The significant contribution of clay, mica, silicate, and accessory minerals, as well as lithic fragments associated with clastic sedimentary and metamorphic rocks in the upstream region.In the downstream area, carbonate components associated with carbonate sedimentary rocks have important contributions.On the other hand, sample location SS09 is characterized by notable inputs from the silicate and accessory minerals, such as quartz, feldspar, muscovite, garnet, apatite, and chlorite, which are associated with metamorphic rocks.(6) Soin River catchment This catchment area is covered by the Lolotoi Complex, Wailuli Formation, Bobonaro Complex, Viqueque Formation, Dilor Conglomerate, and Ainaro Formation (Figures 1C and 3).The highest P 2 O 5 content was recorded in this drainage basin, with the average value reported at 0.22 wt% (ranging from 0.20 to 0.27 wt%) (Table 1).As shown in Table 2A, SiO 2 showed a positive correlation with Al 2 O 3 , MgO, Na 2 O, and K 2 O, as opposed to Al 2 O 3 , which was negatively associated with MnO and CaO and positively correlated with all other major elements.CaO had a strong positive association with MnO and a negative relationship with all other elements.Positive relationships were observed between TiO 2 , Fe 2 O 3 , MgO, Na 2 O, K 2 O, and P 2 O 5 .High amounts of TiO 2 , Al 2 O 3 , MgO, Na 2 O, K 2 O, P 2 O 5 , Fe 2 O 3 , and SiO 2 (5)Karau Ulun River catchment This river catchment is mostly covered by the Lolotoi Complex, Wailuli Formation, Dartollu Limestone, Viqueque Formation, Dilor Conglomerate, and Ainaro Formation (Figures1C and 3).The lowest concentration of Na 2 O was registered in this catchment area, with an average value of 1.82 wt% (range 1.06-2.47wt%)(Table1).

Table 1 )
. As shown in Table 2B, SiO 2 was positively correlated with Na 2 O, K 2 O, P 2 O 5 , and Al 2 O 3 .Al 2 O 3 was positively associated with SiO 2 , TiO 2 , Fe 2 O 3 , Na 2 O, K 2 O, and P 2 O 5 .Moderate-to-very strong correlations were observed between TiO 2 , Fe 2 O 3 , MnO, MgO, and CaO.The Al 2 O 3 concentration appeared to increase in the upstream direction.High SiO 2 , Al 2 O 3 , K 2 O, P 2 O 5 , Na 2 O, TiO 2 , Fe 2 O 3 , and MgO contents and low CaO and MnO concentrations were recorded near the upstream areas.Downstream and midstream regions, specifically sample location SS045, had high concentrations of Fe 2 O 3 , MgO, TiO 2 , and MnO and low Na 2 O and K 2 O contents.In addition, high values of SiO 2 , Na 2 O, and CaO, along with low concentrations of TiO 2 , Fe 2 O 3 , and Al 2 O 3 , were also observed in the downstream area (Figure (10)3 , Na 2 O, and P 2 O 5 and the lowest concentrations of CaO and MnO were recorded.The reported distribution of the elemental concentrations near the downstream regions showed high measured values of TiO 2 , Al 2 O 3 , Fe 2 O 3 , MgO, MnO, and P 2 O 5 .In contrast, elevated values of SiO 2 , K 2 O, and MnO were recorded at sample site SS049.The downstream regions also appeared to have high concentrations of SiO 2 , Na 2 O, and Al 2 O 3 and low TiO 2 and CaO contents.These findings suggest that the midstream areas were largely influenced by quartz, clay, and mica contents as well as lithic fragments, which are associated with clastic sedimentary rocks.Near the downstream regions, clay, muscovite, quartz, plagioclase, chlorite, sphene, garnet and apatite minerals affiliated with clastic sedimentary rocks and their altered rocks owing to metamorphic processes made major contributions.However, sample location SS049 was characterized by notable inputs from quartz, muscovite, garnet, calcium carbonate, and manganese minerals, which were associated with clastic sedimentary interbedded with carbonate sedimentary rocks (such as micaceous shales and sandstones, marls and calcilutites) of the Wailuli Formation and their altered rocks due to metamorphic processes.In addition, there were major contributions from silicate and clay contents related to clastic sedimentary rocks in the downstream areas.(10)SaheRiver catchment The Sahe River catchment is covered by the Lolotoi Complex, Aitutu Formation, Wailuli Formation, Dartollu Limestone, Cablac Limestone, Viqueque Formation, Dilor Conglomerate, and Suai Formation (Figures 1C and 3).The highest concentrations of CaO and K 2 O, as well as the lowest values of Al 2 O 3 , SiO 2 , MgO, Na 2 O, and P 2 O 5 , were also registered in this river catchment (Table 1).CaO had a positive correlation with TiO 2 and MnO, and SiO 2 showed a very strong association with Al 2 O 3 , and they had a moderate-to-strong positive relationship with MgO, Na 2 O, and K 2 O (Table 2D).Positive associations were observed between TiO 2 , Fe 2 O 3 , MnO, MgO, Na 2 O, and P 2 O 5 .The presence of TiO 2 in carbonate sedimentary rocks associated with Mn minerals is not common.TiO 2 could be associated with Mn minerals in carbonate sedimentary rocks, suggesting the presence of secondary Ti-bearing minerals owing to certain alteration processes.In this river catchment, the highest concentrations of SiO 2 and K 2 O appeared near the upstream areas, along with the lowest concentrations of TiO 2 , Fe 2 O 3 , Na 2 O, MnO, and CaO.Although high measured concentrations of SiO 2 , K 2 O, and MnO were reported in the midstream regions, sample location SS015 also appeared to have high measured contents of CaO and MnO and low concentrations of SiO 2 , Al 2 O 3 , K 2 O, and Na 2 O.The downstream regions were also reported to have high SiO 2 , Na 2 O, K 2 O, MgO, and MnO contents.

Table 4 .
Variable results on first three components of PCA.Variables selected are represented by loadings in bold and shaded cells.

Table 5 .
Sample results on first three components of principal compenent analysis (PCA).The samples selected are represented by loadings in bold and shaded cells.