Efficiency Comparison of Reverse Circulation and Diamond Drilling for Phosphate Exploration in Nigeria’s Sokoto Basin

Abstract

This study compares reverse circulation (RC) and diamond core (CD) drilling for exploring phosphate nodules in Nigeria’s Sokoto Basin. A programme of 15 RC wells (1000 m) and 2 CD wells (151 m) assessed performance metrics. RC drilling was significantly faster, progressing three times quicker than CD, and successfully identified phosphate layers 93% of the time. However, its bulk samples caused nodule breakage, reducing grade accuracy. Conversely, CD provided intact cores for precise stratigraphy with 78% recovery but was slower and 40% more expensive. Mineralisation was found primarily in the Dange formation’s gypsiferous shales, with optimal depths varying regionally from under 15 to over 100 metres. The analysis concludes RC is best for rapid lateral resource scoping, while CD is essential for detailed vertical delineation. A cost-time trade-off exists, favouring a hybrid strategy: initial RC mapping followed by targeted CD for reserve classification. The findings offer an optimisation framework for sedimentary phosphate exploration, though future work should integrate advanced coring and 3D modelling to better account for the basin’s geological trends.

Share and Cite:

Salati, L.K. and Adeyemo, J.T. (2026) Efficiency Comparison of Reverse Circulation and Diamond Drilling for Phosphate Exploration in Nigeria’s Sokoto Basin. International Journal of Geosciences, 17, 288-312. doi: 10.4236/ijg.2026.174014.

1. Introduction

Exploring Nigeria’s phosphate deposits in the Sokoto Basin is vital for economic diversification and agricultural development, but faces unique geological challenges requiring optimal appraisal methods [1]-[3]. Phosphate mineralisation there, found in Paleocene formations, shows variable geometry and grade distribution [4]-[6].

Effective exploration requires accurate, cost-effective drilling methods, as methodology selection is a primary determinant of data quality, cost, and viability [7]-[11]. Methods must provide representative samples and precise geological logging [12]-[18]. The choice between reverse circulation (RC) and diamond core (CD) drilling is critical; RC offers speed and lower cost for reconnaissance [19] [20], while CD provides intact core for detailed geology and metallurgy despite higher cost [21].

Evaluating techniques for phosphate involves balancing key factors: penetration rates, cost per metre, sample quality and representativeness (critical for nodular deposits), data fidelity, depth capabilities, logistics, and environmental impact [22]-[26]. While speed and cost are important [27]-[31], sample quality is paramount to avoid distorting resource estimates in variable deposits [32]-[38]. Sample outcomes directly link to geological data fidelity [26] [39]-[41], and practical considerations like terrain and climate are also key [42] [43]. The study synthesises field data to establish guidelines for selecting the optimal drilling method for different Sokoto Basin exploration stages, aiming for effective exploration with mitigated financial risk.

2. Materials and Methods

2.1. Location and Description of the Study Area

The Sokoto Basin in north-western Nigeria, covering approximately 64,000 km2 [44]-[49], contains a prospective phosphate area of about 9000 km2 [2] [50] [51]. Phosphate mineralisation is hosted in the west-dipping shale horizons of the Dange Formation [1]-[3] [6] [52]. Due to the basin’s gentle westward dip, the mineralisation depth increases from east (15 m) to west (>100 m), extending into Kebbi State [52]-[54]. Figure 1 shows the location map of Sokoto Basin with its drainage pattern.

Figure 1. Location map of the Sokoto Basin with drainage pattern (Inset: Map of West Africa showing the location of Nigeria with reference to the position of Sokoto Basin, [44]).

2.2. Field Studies

A drilling programme was planned in the Sokoto Basin to determine the continuity and estimate the size of phosphate resources. Since the phosphate is not in a continuous layer, the plan combined two drilling methods: Reverse Circulation (RC) drilling to collect bulk samples for analysis, and Diamond Core (CD) drilling for detailed geological logging and future research.

The basin was divided into three zones (eastern, central, western), with wells spaced progressively farther apart (10 km, 20 km, and 40 km, respectively) based on the known characteristics of the mineralisation in each area.

Three rigs were used: a Comacchio Geo 602 Explorer with a riffle splitter (Figure 2) for all RC drilling, and two HYDX-5A rigs (Figure 3) for all CD drilling. The RC programme used a pneumatic hammer, drilling 15 wells to 1000 m total. The CD programme used HQ/NQ cores with water-based drilling, completing 2 wells to 151 m total.

Figure 2. Multipurpose Comacchio Geo 602 Explorer complete with riffle splitter.

Figure 3. HYDX-5A Multipurpose core drilling rig.

The disparity in well count—RC (15) versus CD (2)—was deliberate. The RC programme provided rapid, cost-effective reconnaissance across the vast 9000 km2 basin, while the two CD wells delivered high-fidelity calibration and grade data at key locations. Drilling followed a chronological sequence across the Gwadabawa and Binji areas, with all three rigs deployed to accelerate the work. Vertical drilling was used throughout for CD due to the mineralisation’s nature.

2.3. Sample Collection, Labeling and Logging

The study employed a random sampling technique across all locations. Two methods were used: RC (Reverse Circulation) drilling, which collected cyclone-separated cuttings logged via “Sailform” software, and CD (Core Drilling) for continuous cores. For RC drilling, samples were collected from the cyclone base for every meter of drilling depth, as illustrated in Figure 4.

Figure 4. Sample collection at the base of cyclone.

These samples were placed in systematically labelled bags, which were marked with key identifiers including the study name, drilling type, well ID, and sample depth, an example of which is shown in Figure 5. Sample logging was conducted in two phases, with initial real-time logging of weight and lithological data into “Sailform” using a comprehensive template that captured well ID, depth, location, and other meta data.

Figure 5. Sample packaging and labelling.

Samples were collected per meter and placed in labelled trays for initial review. Suspected phosphate horizons were then re-examined after being split with a ripple splitter (Figure 6) to ensure representative sampling and avoid bias. For these horizons, 2 - 3 kg samples were systematically collected, tagged, and sent to the lab, with duplicate samples stored and their tag numbers recorded in the “Sailform” software.

Figure 6. Ripple splitter, splitting and tagging of suspected phosphate horizons.

2.4. Sample Movement and Storage

After drilling, samples were moved to a secure store, documented with check forms, and bagged separately for assaying with accompanying registers. Duplicates and non-tagged samples were systematically arranged in a well-ventilated, leak-proof store (Figure 7).

Figure 7. Sample storing in Sokoto town.

The drilling process involved a strict verification protocol before and after operations. Before drilling, a pre-sign-off form was completed by key personnel to confirm site safety, sample integrity, and capture technical details like the well ID and planned depth. After drilling, a post-sign-off form ensured the site was properly cleaned, rehabilitated, and the borehole was sealed with an engraved identification tag. Finally, the prepared and logged samples were dispatched to laboratory for phosphate analysis.

2.5. Laboratory Analyses and Sample Preparations

A total of 242 phosphate-bearing samples were collected and analyzed. To determine their elemental composition, the samples were prepared as homogeneous pellets and analyzed using X-Ray Fluorescence (XRF). To identify the specific phosphate minerals and their relative abundances, the samples were analyzed using X-Ray Diffraction (XRD). These analyses are essential for understanding the ore’s properties and industrial potential.

2.6. Statistical Analysis

1) Description of phosphate’s thickness, depth, and recovery rates

The wells in the western zone are the deepest, where phosphate is expected to be penetrated at over 100 m, such as in Sakwae, Matankari, and Kalgo areas. In the eastern zone around Chimola, Kagogo, Salame, and Kware, the wells are quite shallow, with phosphate forming in horizons below 15 m. Most of the wells in the eastern zone were bored using the RC method, with samples collected at every 1 m.

2) RC and CD metrics

Drilling reached depths of 40 - 93 m. The rock formation is varied, but shale from the Dange formation dominates and is the primary phosphate-bearing rock, often containing phosphatic nodules and gypsum, as confirmed by spot tests. Other rocks like limestone and marl are generally non-phosphatic. While phosphate mineralization is widespread from shallow depths down to 83 m, the optimal, richest horizon for mineralization typically begins below 15 m.

3. Results and Discussions

The discussions are premised on borehole data (Table 1 and Table 2), results of laboratory analyses (Figures 8-11), summary of RC wells attributes (Table 3) and data from the lithologic logs of wells of the bedding units within the different formations of phosphate mineralisation in the Sokoto Basin (Figures 12(a)-(n)).

3.1. Drilling Efficiency

Drilling across 15 wells indicated widespread phosphate layers that thin southward, consistent with prior research [55] [56]. The local geology is complex, featuring a phosphate-bearing Dange Formation shale layer between sandy and limestone units. The Gamba Formation is often absent due to erosion. Phosphate layer variability is attributed more to local geological features like erosion and unconformities than to elevation, a conclusion supported by multiple studies [57]-[61].

Table 1. Chronological log of SPDD001 (Dange Formation, 71.8 m depth; Long. 5.28268, Lat. 13.42994, and elevation of 261 m).

S/No

Depth (m)

Unit/formation

Mineral

Grain size

Alteration

Texture

Colour

Over form

Lithology

Thickness

Solid core recovery (m)

Barrel core recovery (m)

Total core recovery (m)

Loss (-)

Gain (+)

% Recovery

1

00 - 1.2

Gwandu

Quartz

Gravel

Siliceous

Coarse-grained

Brown

Coarse gravel

Sand and laterite

1.2

0

0.73

0.73

−0.47

60

2

1.2 - 3.77

Gwandu

Quartz

Pebble

Ferruginous sandstone

Coarse-grained to medium-grained

Reddish brown

Ferruginous sandstone, shale, and clay

2.57

0.37

0.33

0.7

−1.87

27

3

3.77 - 6.3

Gwandu

Quartz, feldspar, oolitic iron

Pebble

Carbonate

Fine-grained to medium-grained

Brown, yellowish-brown-white

Clay, ferruginous sandstone, and shale

2.53

1.34

1.19

2

−0.53

79

4

6.3 - 73

Gwandu

Quartz, mica

Siliceous

Medium-grained

Grey

-

Shale

1

0.9

-

0.9

−0.1

90

5

7.3 - 7.8

Gwandu

Calcite and quartz

Carbonate

Fine-grained

Reddish and grey

Clay and laminated limestone

0.5

0.38

0.72

1.1

0.6

220

6

7.8 - 12.69

Gwandu

Quartz and mica

Poorly sorted

Siliceous

Fine-grained and coarse

Brown and dark red

Clay and ferruginous sandstone

4.89

2.62

0.33

2.95

−1.94

60

7

12.69 - 20.52

Gamba

Quartz

Well sorted

Siliceous and carbonate

Fine-grained to coarse-grained

White and brown

Limestone and clay

7.83

0.7

0.5

1.2

−6.63

15

8

20.52 - 20.94

Gamba

Quartz

Siliceous

Medium-grained

Brown

Mudstone

0.42

0.33

-

0.3

−0.12

71

9

20.94 - 23.02

Gamba

Calcite

Carbonate

Fine-grained to medium-grained

White to brown

Limestone and clay

2.08

0.74

-

0.74

−1.34

36

10

23.02 - 24.64

Gamba

Quartz

Pebble

Siliceous

Coarse-grained

Reddish

Coarse gravel

Ferruginous sandstone

1.62

0.05

-

0.05

−1.57

3

11

24.64 - 25.24

Gamba

Calcite

Carbonate

Fine-grained to medium-grained

White

Limestone

0.6

0.23

0.26

0.49

−0.11

81

12

25.24 - 25.66

Gamba

Quartz

Carbonate

Fine-grained to medium-grained

Grey and yellow

-

Shale and mudstone

0.42

0.45

0.19

0.64

0.22

152

13

25.66 - 27.12

Kalambaina

Quartz

Siliceous

Medium-grained

White

-

Mudstone

1.46

0.07

-

0.07

−1.39

5

14

27.12 - 28.05

Kalambaina

Quartz, mica

Siliceous

Coarse-grained

White

Mudstone

0.93

0.15

0.05

0.2

−0.73

21

15

28.05 - 28.85

Kalambaina

Quartz, mica

Carbonate

Fine-grained

Grey

Limestone and shale

0.8

0.34

0.54

0.88

0.08

110

16

28.85 - 30.01

Kalambaina

Calcite and quartz

Carbonate

Fine-grained to medium-grained

White

Limestone

1.16

0.13

0.53

0.66

−0.5

56

17

30.01 - 30.74

Kalambaina

Calcite, quartz

Carbonate

Fine-grained to medium-grained

White

Limestone with traces of fossils

0.73

0.9

0.03

0.93

0.2

127

18

30.74 - 32.77

Kalambaina

Calcite and quartz

Carbonate

Fine-grained to medium-grained

Dull white

Limestone

2.03

0.10

-

0.10

−1.93

5

19

32.77 - 33.69

Kalambaina

Quartz and feldspar

Siliceous

Fine-grained

Brown

-

Clay

0.92

0.62

-

0.62

−0.3

67

20

33.69 - 34.71

Dange

Quartz, gypsum, mica, and glauconite

Siliceous

Fine-grained to medium-grained

Yellowish to greenish

Shale

1.02

0.80

0.20

1

0.02

98

21

34.71 - 35.51

Dange

Quartz, mica

Carbonate

Medium-grained

Greenish to yellow and grey

Shale

0.80

0.60

0.30

0.90

0.1

113

22

35.51 - 36.67

Dange

Quartz, mica

Siliceous

Fine-grained to medium-grained

Yellowish to green

Shale

1.16

0.86

-

0.86

−0.3

74

23

36.67 - 38.07

Dange

Quartz and olivine

Carbonate

Fine-grained

Greenish to yellow

Shale or clay

1.4

1.1

-

1.1

−0.3

78

24

38.07 - 39.51

Dange

Quartz, feldspar

Siliceous

Very fine-grained

Brown-pale green

Clay (bentonite)

1.44

0.23

-

0.23

−1.21

16

25

39.51 - 40.60

Dange

Quartz, gypsum

Carbonate

Fine-grained to medium-grained

Yellowish to green

Shale

1.09

2.02

0.2

2.22

1.13

203

26

40.60 - 42.04

Dange

Quartz and gypsum

Siliceous

Fine-grained to medium-grained

Dark

Shale

1.44

0.62

0.3

0.92

−0.52

156

27

42.04 - 43.25

Wurno

Quartz

Carbonate

Very fine-grained

Grey

Shale

1.21

0.24

-

0.24

−0.97

20

28

43.25 - 45.01

Wurno

Quartz and gypsum

Carbonate

Fine-grained to medium-grained

Dark

Shale and intercession of gypsum

1.79

0.22

0.72

0.96

−0.8

183

29

45.01 - 46.74

Wurno

Quartz

Carbonate

Fine-grained to medium-grained

Dark

Shale

1.73

1.9

-

1.9

0.17

109

30

46.74 - 48.11

Wurno

Quartz, mica

Carbonate

Fine-grained to medium-grained

Dark

Shale with traces of fossils

1.34

1.0

0.2

1.20

−0.14

90

31

48.11 - 49.91

Wurno

Quartz, mica

Carbonate

Fine-grained to medium-grained

Dark

Shale and intercalation of limestone nodules

1.8

1.06

0.44

1.50

−0.3

83

32

49.91 - 51.14

Wurno

Quartz, mica

Carbonate

Medium-grained

Dark

Shale

1.23

1.0

0.21

1.21

−0.02

98

33

51.14 - 53.66

Wurno

Quartz, mica

Siliceous

Fine-grained to medium-grained

Dark

Shale

2.46

1.58

0.28

1.86

−0.6

75

34

53.66 - 54.60

Wurno

Quartz, mica

Quartz, mica

Fine-grained to medium-grained

Grey to greenish

Limestone with traces of fossils and phosphate nodules

0.94

1.09

0.62

1.69

0.74

178

35

54.60 - 56.15

Wurno

Quartz, mica, pyrite

Siliceous

Fine-grained to medium-grained and coarse

Grey, yellowish, and green

Shale, pyrite, and limestone intercalation

1.55

1.10

0.20

1.30

−0.25

84

36

56.15 - 57.46

Wurno

Quartz, mica

Carbonate

Coarse-grained

Dark

Sandstone with voids

1.31

0.06

-

0.06

−1.25

5

37

57.46 - 58.62

Wurno

Quartz, mica

Siliceous

Fine-grained to medium-grained

Grey

Shale

1.16

0.16

0.16

−1

14

38

58.62 - 59.31

Wurno

Quartz, mica

Siliceous

Fine-grained

Grey

Shale

0.69

0.16

-

0.16

−0.53

23

39

59.31 - 60.51

Wurno

Quartz, mica

Siliceous

Fine-grained

Grey

Shale

1.2

1.0

0.1.

1.10

−0.1

92

40

60.51 - 61.57

Wurno

Quartz, mica

Siliceous

Fine-grained

Grey

Shale

1.06

0.36

-

0.36

−0.7

35

41

61.57 - 62.84

Wurno

Quartz, mica

Carbonate

Fine-grained to medium-grained

Dark-grey

Shale with intercalations of limestone nodules

1.27

0.71

0.36

0.89

−0.38

70

42

62.84 - 64.04

Wurno

Quartz, mica, and calcite

Fissile

Carbonate

Medium-grained

Dark, light to brown

Shale, limestone, and clay

1.20

0.6

0.52

1.12

0.08

93

43

64.04 - 65.46

Wurno

Quartz, calcite, and mica

Carbonate

Fine-grained to medium-grained

Grey

Shale

1.42

0.83

0.14

0.97

−0.45

68.3

44

65.46 - 66.86

Wurno

Quartz, mica

Carbonate

Medium-grained

Grey

Shale

1.40

0.98

0.12

1.1

−0.3

78

45

66.86 - 69.52

Wurno

Quartz, mica, calcite

Carbonate

Fine-grained

Grey

Shale with intercalations of limestone and gypsum

2.66

0.17

0.92

1.09

−1.57

40

46

69.52 - 71.83

Wurno

Quartz, mica

Carbonate

Medium-grained

Dark

Shale

2.31

1.5

0.3

1.8

−0.51

78

Table 2. Chronological sequence of SPDD00 w16 borehole (Long. 5.27935, Lat. 13.58550, and elevation of 259 m).

S/No

Depth (m)

Unit or formation

Mineral

Grain size

Alteration

Texture

Colour

Contact

Lithology

Thickness

Solid core recovery (m)

Barrel core recovery (m)

Total core recovery (m)

Loss (-)

Gain (+)

% Recovery

1

00 - 1.16

Gamba

Quartz

Medium-grained

Siliceous

Semi, consolidated

Brown-whitish

Erosional

Sand and laterite

1.16

0.5

0.5

0.66

43.1

2

1.16 - 2.57

Gamba

Quartz

Medium-grained to coarse-grained

Siliceous

Consolidated

Yellow-reddish

Gradational

Sandstone and laterite

1.41

0.90

0.90

0.51

63.8

3

2.57 - 5.01

Gamba

Quartz

Coarse-grained

Siliceous

Consolidated

Dark reddish yellow

Ferruginous sandstone

2.44

0.84

0.84

1.6

34.4

4

5.01 - 6.73

Gamba

Quartz

Coarse-grained to medium-grained

Siliceous

Semi-consolidated

Dark red dish-yellow

Sharp

Iron, sandstone, and siltstone

1.72

0.96

0.96

0.76

55.8

5

6.73 - 7.03

Gamba

Quartz

Fine-grained

Siliceous

Semi-consolidated

Yellow

Gradational

Clay

0.3

0.28

0.28

0.02

93.3

6

7.03 - 7.83

Kalambaina

Quartz

Fine-grained

Fine-grained

Fine-grained

Yellow- light grey

Gradational

Clay and shale

0.8

0.79

0.79

0.01

98.6

7

7.83 - 10.13

Kalambaina

Quartz

Fine-grained

Fine-grained

Fissile

Light grey-reddish brown

Shale, laterite, and clay

2.30

1.67

1.67

0.63

72.6

8

10.13 - 12.35

Kalambaina

Quartz

Fine-grained to medium-grained

Siliceous

Consolidated

Dark reddish-brownish

Sharp

Ferruginous sandstone, clay, and intercalation of limestone

2.22

1.65

1.65

0.57

74.3

9

12.35 - 13.63

Kalambaina

Quartz

Fine-grained

Fissile

Carbonate

Dark grey

Sharp

Shale

1.28

1.26

1.26

0.02

98.4

10

13.63 - 14.17

Kalambaina

Quartz

Fine-grained

Fissile

Carbonate

Dark grey

Shale

0.54

0.48

0.48

0.06

88.8

11

14.17 - 16.13

Kalambaina

Quartz

Fine-grained

Fine-grained

Carbonate

Light grey

Gradational

Shale

1.96

1.85

1.85

0.11

94.4

12

16.13 - 16.84

Kalambaina

Quartz, pyrite, phosphate nodule

Fine-grained

Fissile

Carbonate

Light grey

Shale with phosphate nodules and pyrite

0.71

0.7

0.7

0.01

98.6

13

16.84 - 17.99

Kalambaina

Pyrite and quartz

Fine-grained

Fissile

Carbonaate

Dark grey

Shale with pyrite

1.15

1.14

1.14

0.01

99.1

14

17.99 - 20.02

Kalambaina

Quartz, pyrite

Fine-grained

Fissile

Carbonate

Light grey

Shale

2.03

2.0

2.0

0.03

98.5

15

20.02 - 21.52

Kalambaina

Quartz, gypsum

Fine-grained

Fissile

Carbonate

Light grey

Shale with lenses of gypsum

1.5

1.46

1.46

0.04

97.3

16

21.52 - 22.73

Kalambaina

Quartz, calcite

Fine-grained

Consolidated

Carbonate

White

Sharp

Limestone with traces of fossils

1.21

1.08

1.08

0.13

89.3

17

22.73 - 25.52

Kalambaina

Quartz, calcite

Fine-grained

Consolidated

Carbonate

Grey

Limestone

2.79

2.74

2.79

0.05

98.2

18

25.52 - 28.02

Kalambaina

Quartz, calcite

Fine-grained

Consolidated

Carbonate

White

Limestone

2.5

2.45

2.45

0.05

98.0

19

28.02 - 32.94

Kalambaina

Calcite

Fine-grained

Consolidated

Carbonate

White

Limestone

4.92

3.68

3.68

1.24

74.8

20

32.94 - 36.62

Kalambaina

Calcite, quartz

Fine-grained

Consolidated

Calcareous

Light grey

Gradational

Limestone with shell

3.68

2.10

2.10

1.58

57.1

21

36.62 - 42.12

Kalambaina

Quartz, calcite

Fine-grained

Consolidated

Calcareous

Light grey

Limestone

5.5

2.72

2.73

2.77

50

22

42.12 - 45.85

Kalambaina

Quartz, calcite

Fine-grained

Consolidated

Calcareous

Light grey

Limestone with shale

3.73

1.9

1.9

1.83

51

23

45.85 - 47.85

Dange

Quartz, pyrite, gypsum, and calcite

Fine-grained

Fissile

Carbonate

Light grey

Gradational

Shale

2.0

1.65

1.65

0.35

82.5

24

47.85 - 51.32

Dange

Quartz, pyrite, phosphate, and calcite

Fine-grained

Fissile

Carbonaceous

Dark

Shale with phosphate nodules, pyrite. Start

3.47

2.52

2.52

0.95

72.6

25

51.32 - 55.37

Dange

Quartz, phosphate

Fine-grained

Fissile

Carbonaceous

Dark

Shale with phosphate nodules

4.05

3.77

3.77

0.35

93.1

26

55.37 - 57.51

Dange

Quartz, gypsum, phosphate

Fine-grained

Fissile

Carbonaceous

Dark

Shale with phosphate nodules

2.14

1.8

1.8

0.34

84.1

27

57.51 - 60.79

Dange

Quartz, phosphate

Fine-grained

Fissile

Carbonaceous

Dark

Shale with phosphate nodules. End

3.28

1.83

1.83

1.45

55.8

28

60.79 - 63.01

Wurno

Quartz

Fine-grained

Fissile

Carbonaceous

Dark

Shale

2.22

1.21

1.12

1.19

54.5

29

63.01 - 64.35

Wurno

Quartz

Fine-grained

Fissile

Carbonaceous

Dark

Shale

1.34

1.33

1.33

0.01

99.3

30

64.35 - 68.72

Wurno

Quartz

Fine-grained

Fissile

Carbonaceous

Light grey and dark

Shale

4.37

2.73

2.73

1.66

62.5

31

68.72 - 72.32

Wurno

Quartz

Fine-grained

Fissile

Carbonaceous

Dark

Shale

3.6

2.41

2.41

1.19

67.0

32

72.32 - 75.72

Wurno

Quartz, pyrite

Fine-grained

Fissile

Carbonaceous

Dark

Shale

3.4

3.10

3.10

0.3

91.1

33

75.72 - 80.00

Wurno

4.28

Figure 8. Elevation map of RC drilling boreholes showing wells distribution in relation to topography around Dange area of Sokoto Basin.

Figure 9. Isopach map showing thicknesses of phosphate bearing horizons recovered from 15 RC drilled wells within the Sokoto basin.

Figure 10. Elevation map showing drilled RC well locations.

Figure 11. Inferred thicknesses of phosphate bearing horizons around Dange area where RC drilling was conducted.

Figure 12. (a): Lithologic log of well 35 (N 13.22873, E 005.27907) kware, Kware LGA, Sokoto; (b): Lithologic log of well 53 (N 12. 71278. E 005.13532) Gidan Dandala, Shagari LGA, Sokoto State; (c): Lithologic log of well 37 (N 13.06292, E 005.18015) Kalambaina, Wamako LGA, Sokoto State; (d): Lithologic log of well 38 (N 12.97233, E 005.14647), Bodinga LGA, Sokoto State; (e): Lithologic log of well 39 (N 12.91047, E 005.13461) Gadan Dadinkai, Bodinga LGA, Sokoto State; (f): Lithologic log of well 40 (N 12.83771, E 005.06572) Birnin Ruwa, Bodinga LGA, Sokoto State; (g): Lithologic log of well 41 (N 12.75044, E 005.02237) Dabagi, Yabo LGA, Sokoto State; (h): Lithologic log of well 43 (N 12.59028, E 004.97746) Jazomo, Shagari LGA, Sokoto State; (i): Lithologic log of well 44 (N 12.51159, E 004.89166); (j): Lithologic log of well 45 (N 12.45996, E 004.77916); (k): Lithologic log of well 51 (N 13.20801, E 005.45325) Wurno, Wurno LGA, Sokoto State; (l): Lithologic log of well 52 (N 12.92562, E 005.30902) Amana, Dange Shuni LGA, Sokoto State; (m): Lithologic log of well 58 (N 13.14804, E 005.34524) Wurno LGA, Sokoto State; (n): Lithologic log of well 56 (N 12.39171, E 005.01839) Murdowu, Tambuwal LGA, Sokoto State.

Table 3. Summary of RC wells attributes in the study area.

S/N

Well id

Easting

Northing

Elv(m)

EOH(m)

Suspected phosphate horizon

No. of samples for assay

1

SP-RC W030 (kebbi)

4.202632

12395300

205

66

No phosphate bearing formation

-

2

SP-RC W035

5.27907

13.22873

259

50

9 - 21 m

15

3

SP-RC W051

5.45325

13.20801

334

72

30 - 41 m

14

4

SP-RC W058

5.34524

13.14804

310

74

30 - 57 m

29

5

SP-RC W037

5.18015

13.06292

261

49

15 - 26 m

14

6

SP-RC W038

5.14647

12.97233

272

58

27 - 43 m

19

7

SP-RC W052

5.30902

12.92562

321

68

8 - 26 m

21

8

SP-RC W039

5.13461

12.91047

277

80

23 - 41 m

21

9

SP-RC W040

5.06572

12.83771

312

91

63 - 83 m

23

10

SP-RC W053

5.13532

12.71278

300

41

11 - 21 m

13

11

SP-RC W041

5.02237

12.74044

299

93

54 - 72 m

21

12

SP-RC W043

4.97746

12.59028

265

60

6 - 15 m

12

13

SP-RC W056

5.01839

12.39171

323

40

0 - 7 m

9

14

SP-RC W044

4.89166

12.51159

290

81

26 - 39 m

16

15

SP-RC W045

4.77916

12.45996

265

77

7 - 19 m

15

Total = 1000

Total = 242

A comparison of drilling methods found RC faster and cheaper but with lower sample recovery, while DC was slower and more expensive but provided superior sample quality and recovery. The study acknowledges a potential limitation, as its DC dataset was much smaller than its RC dataset, possibly not capturing the basin’s full complexity, a limitation also noted by few studies [62] [63].

3.2. Mineralisation Characteristics

The phosphate occurrences in the Dange Formation are found as nodules and pellets within shale layers 6.9 - 27 meters thick, associated with limestone, marl, black shale, and gypsum in the Chadawa area, consistent with prior Sokoto Basin studies [3] [64]. However, RC drilling is unsuitable for sampling these nodules as it fragments them, producing unreliable mineral content data and preventing accurate measurements, resulting in low-confidence data inadequate for detailed mine planning or financial analysis.

3.3. Challenges

RC drilling was primarily limited by sample contamination, where unconsolidated materials from overlying formations (e.g., the Gwandu Formation) infiltrated samples, obscuring the target phosphate shale and complicating the definition of mineralized zones. Conversely, CD was constrained by high costs and slow progress, exacerbated by a crew capacity gap, the inherently slow coring method, and difficult ground conditions that reduced progress to 2 - 2.7 meters per day and necessitated drilling mud for stability. These operational challenges find support in the literature [65] [66].

3.4. Spatial Distribution

The southward thinning of phosphate horizons shown by isopach mapping indicates that mineralization is not topographically controlled, particularly near the Wurno and Bodinga wells. This supports the views of [53] [67]-[69] and is a critical pathfinder for phosphate exploration in the Sokoto Basin.

4. Conclusion

RC drilling is the fastest and most cost-effective method for a quick reconnaissance. However, to build a detailed resource model, CD is essential, though it requires more budget and planning. Because the phosphate mineralisation is discontinuous, the best strategy is to use a hybrid of both drilling approaches. This optimises costs and ensures the collection of data needed for a reliable and bankable resource. Further research in the study area is capable of adopting a two-phase drilling approach to locate targets with broad RC drilling and define the resource using focused CD. Improved efficiency and data quality can also be achieved with expanded drilling westward, enhanced CD training, and integration of new technologies like geophysical surveys.

Acknowledgements

The Nigerian Geological Survey Agency (NGSA) is gratefully acknowledged by the authors for funding this research, while the collaborative support rendered by Geocardinal Engineering Services Limited, AG Vision Company Limited and Sokoto State Government is greatly appreciated.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

References

[1] Awalla, C.O., Ezeh, C.C. and Nwabufo-Ene, K.E. (2002) Review of Sokoto Phosphate Deposit Sokoto, Nigeria. Global Journal of Pure and Applied Sciences, 8, 89-93. [CrossRef]
[2] Mustapha, S.O., Olatunji, A.S., Ajayi, F.F. and Isibor, R.A. (2022) Mineralogy and Geochemistry of Some Phosphate Deposits for Possible Rare Earth Elements Mineralization Potentials within Sokoto Basin, Northwestern Nigeria. Indian Journal of Science and Technology, 15, 2570-2578. [CrossRef]
[3] Okosun, E.A. and Alkali, Y.B. (2013) The Geochemistry, Origin and Reserve Evaluation of Sokoto Phosphate Deposit, North Western, Nigeria. Earth Science Research, 2, 111-121. [CrossRef]
[4] Obaje, N.G., Faruq, U.Z., Bomai, A., Dabai Moses, S., Ali, M., Adamu, S., et al. (2020) A Short Note on the Petroleum Potential of the Sokoto Basin in North-Western Nigeria. Petroleum Science and Engineering, 4, 34-38. [CrossRef]
[5] Obaje, S.O., Ojutalayo, A., Ogedengbe, O. and Okosun, E.A. (2014) Nigeria’s Phosphate and Uranium Mineral Occurrences: Implication for Mineral Investment. Journal of Environmental & Earth Sciences, 4, 83-89.
[6] Onuigbo, E.N., Okoro, A.U. and Nwokeabia, C.N. (2017) Clay Minerals as Indicator of Phosphatogenesis: A Case Study of Sokoto Basin, Northwestern Nigeria. Journal of Natural Sciences Research, 7, 21-27.
[7] Amer, M.Y., Salem, S.K., Farahat, M.S. and Salem, A.M. (2025) Reducing Drilling Cost of Geothermal Wells by Optimizing Drilling Operations: Cost Effective Study. Unconventional Resources, 7, Article ID: 100196. [CrossRef]
[8] Babayeju, O.A., Jambol, D.D. and Esiri, A.E. (2025) Reducing Drilling Risks through Enhanced Reservoir Characterization for Safer Oil and Gas Operations. GSC Advanced Research and Reviews, 19, 86-101. [CrossRef]
[9] Ganat, T.A.O. (2020) Understanding Exploration, Appraisal and Production Stages. In: Al-Arbi Omar Ganat, T., Ed., Technical Guidance for Petroleum Exploration and Production Plans, Springer International Publishing, 21-28. [CrossRef]
[10] Iwe, K.A., Isong, D.E., Daramola, G.O. and Agho, M.O. (2021) Sustainable Pore Pressure Prediction and Its Impact on Geo-Mechanical Modelling for Enhanced Drilling Operations. World Journal of Advanced Research and Reviews, 12, 540-557. [CrossRef]
[11] Kaiser, M.J. (2021) A Review of Exploration, Development, and Production Cost Offshore Newfoundland. Natural Resources Research, 30, 1253-1290. [CrossRef]
[12] Bruno, P.P.G. (2023) Seismic Exploration Methods for Structural Studies and for Active Fault Characterization: A Review. Applied Sciences, 13, Article No. 9473. [CrossRef]
[13] Groves, D.I., Santosh, M. and Zhang, L. (2020) A Scale-Integrated Exploration Model for Orogenic Gold Deposits Based on a Mineral System Approach. Geoscience Frontiers, 11, 719-738. [CrossRef]
[14] Hill, E.J., Fabris, A., Uvarova, Y. and Tiddy, C. (2021) Improving Geological Logging of Drill Holes Using Geochemical Data and Data Analytics for Mineral Exploration in the Gawler Ranges, South Australia. Australian Journal of Earth Sciences, 70, 1067-1093. [CrossRef]
[15] Liu, L., Li, T. and Ma, C. (2024) Research on 3D Geological Modeling Method Based on Deep Neural Networks for Drilling Data. Applied Sciences, 14, Article No. 423. [CrossRef]
[16] Tshanga M, M., Ncube, L. and van Niekerk, E. (2024) Remote Sensing Insights into Subsurface-Surface Relationships: Land Cover Analysis and Copper Deposits Exploration. Earth Science Informatics, 17, 3979-4000. [CrossRef]
[17] Zhao, J., He, L., Gong, J., He, Z., Feng, Z., Pang, J., et al. (2025) Predicting Manganese Mineralization Using Multi-Source Remote Sensing and Machine Learning: A Case Study from the Malkansu Manganese Belt, Western Kunlun. Minerals, 15, Article No. 113. [CrossRef]
[18] Zhao, T., Wang, S., Ouyang, C., Chen, M., Liu, C., Zhang, J., et al. (2024) Artificial Intelligence for Geoscience: Progress, Challenges, and Perspectives. The Innovation, 5, Article ID: 100691. [CrossRef] [PubMed]
[19] Sarath, H. and Manuel, J.R. (2010) Techniques Used and Problems Encountered in Normal and Reverse Circulation Drilling at the Mt Roseby Project, Cloncurry, Northwest Queensland. XXV International Mineral Processing Congress (IMPC) 2010, Brisbane, 6-10 September 2010.
[20] Li, L., Ru, D., Li, K. and Yishan, W. (2006) Research and Application of Reverse Circulation Drilling Technology. International Oil & Gas Conference and Exhibition in China, Beijing, December 2006, SPE-104412-MS. [CrossRef]
[21] Orpen, J. and Orpen, D. (2020) Error-Proofing Diamond Drilling and Drill Core Measurements. SEG Discovery, 123, 23-34. [CrossRef]
[22] Aghlan, M.H., Asghari, O. and Emery, X. (2025) Optimizing Infill Drill Hole Decisions While Capturing the Spatial Continuity of Geochemical and Geometallurgical Properties: Application to Gol Gohar Iron Ore Mine, Iran. Minerals, 15, Article No. 478. [CrossRef]
[23] Nicholls, L.V. (2025) Multi-Scale Geometallurgical Characterisation of Legacy Tail-ings for Sustainable Critical Metal Valorisation: A Case Study Exploring the Secondary Prospectivity of Cobalt. PhD Thesis, University of Queensland, 576 p.
[24] Kalpande, V. and Wahhab, A. (2019) Data Mining Analytics of Exploration Drilling Well Cost for Petroleum Economics in Western Offshore India. Petro Tech Conference 2019, New Delhi, 10-12 February 2019, 2233.
[25] Riazi, K., Azhari, H., Charette, J.H., Underwood, F.E., King, J.A., Afshar, E.E., et al. (2022) The Prevalence and Incidence of NAFLD Worldwide: A Systematic Review and Meta-Analysis. The Lancet Gastroenterology & Hepatology, 7, 851-861. [CrossRef] [PubMed]
[26] Rock, O.M., Omonona, V.O. and Abraham, E.M. (2017) Exploring and Reserve Estimation for Industrial Mineral Potential in Parts of Calabar Area (Ewen/Iwuru/Agbangana Axis) Southern Nigeria. Journal of Geology & Geophysics, 6, Article ID: 1000317.
[27] Sterk, R. and Franey, D. (2017) RC Drilling and Sampling: Optimising the Quantity-Quality Balance. Book/Conference Proceedings for the Drilling for Geology Conference, Brisbane, May 2017, 17 p.
[28] Allaire, A., Dorval, A., Beaulieu, C., Girard, M., Paquette, M.-H., Millet, E., Farhangi, S. and Khoury, T. (2024) Pre-Feasibility Study for the Kamistiatusset (Kami) Iron Ore Property Newfoundland and Labrador, Canada. NI 43-101 Technical Report, BBA Document No.: 3813101-000000-40-ERA-0001-R00.
[29] Ghaffari, H., Huang, J.J., Li, J.J. and O’brien, M.F. (2024) Oxide Tailings Project Prefeasibility Study for the Avino Property, Durango, Mexico, NI 43-101 Technical Report. Tetra Tech Canada Inc.
[30] Hatton, M., Lewis, M., Jones, L., Howery, J., Franceschini, M., McIntyre, A. and Kalenchuk, K. (2022) Technical Report on the Feasibility Study for the McIlvenna Bay Project, Saskatchewan Canada. Foran Mining Corporation, NI 43-101, Project Number: 169521520.
[31] Oghenebrume, U.L. and Abiodun, O.A. (2025) Optimization of Drilling Bit Performance: A Cost-per-Foot and Breakeven Analysis of Well-A and Well-B at Pan Ocean Oil Corporation Nigeria. British Journal of Multidisciplinary and Advanced Studies, 6, 9-26. [CrossRef]
[32] Wang, X., Dong, J., Liu, J. and Lyu, J. (2021) Brief Analysis of Cost per Meter Management in Drilling and Well Completion Operations in Bohai Oilfield. Academic Journal of Business & Management, 3, 59-63.
[33] Hapugoda, S. and Manuel, J. (2010) Techniques Used and Problems Encountered in Normal and Reverse Circulation Drilling at the Mt Roseby Project, Cloncurry, Northwest Queensland. XXV International Mineral Processing Congress (IMPC) 2010 Proceedings, Brisbane, 6-10 September 2010.
[34] Stanić, T. (2023) Assessing the Feasibility and Effectiveness of Reverse Circulation in Drilling Operation. Unpublished Master’s Thesis, Montan Universitat Leoben, 136.
https://pure.unileoben.ac.at/ws/portalfiles/portal/24898652/AC17033493.pdf
[35] Engström, K. (2017) A Comprehensive Literature Review Reflecting Fifteen Years of Debate Regarding the Representativity of Reverse Circulation vs Blast Hole Drill Sampling. TOS Forum, 2013, 36-46. [CrossRef]
[36] Annels, A.E. and Dominy, S.C. (2003) Core Recovery and Quality: Important Factors in Mineral Resource Estimation. Applied Earth Science, 112, 305-312. [CrossRef]
[37] Szilas, C. (2022) The Tanzanian Minjingu Phosphate Rock—Possibilities and Limitations for Direct Application. Unpublished PhD Thesis, Royal Veterinary and Agri-cultural University, 178.
https://www.gpsagro.dk/thesis.pdf
[38] Qin, X.Q. (2025) Geochemistry of Trace Metals Associated with Fe-Mn Nodules (FMNs) of Red Soil Profiles in China. Unpublished PhD Thesis, Université Paris Cité, NNT: 2024UNIP7087.
https://theses.hal.science/tel-04989843v1/file/va_Qin_Xiaoquan.pdf
[39] Raabe, D., Sun, B., Da Silva, A.K., Gault, B., Yen, H., Sedighiani, K., et al. (2020) Current Challenges and Opportunities in Microstructure-Related Properties of Advanced High-Strength Steels. Metallurgical and Materials Transactions A, 51, 5517-5586. [CrossRef]
[40] Amon, D.J., Gollner, S., Morato, T., Smith, C.R., Chen, C., Christiansen, S., et al. (2022) Assessment of Scientific Gaps Related to the Effective Environmental Management of Deep-Seabed Mining. Marine Policy, 138, 105006. [CrossRef]
[41] Dominy, S.C. and Glass, H.J. (2025) Geometallurgical Sampling and Testwork for Gold Mineralisation: General Considerations and a Case Study. Minerals, 15, Article No. 370. [CrossRef]
[42] Sketchley, D.A. (1998) Gold Deposits: Establishing Sampling Protocols and Monitoring Quality Control. Exploration and Mining Geology, 7, 129-138.
[43] Heddoun, O., Ait Lemkademe, A. and Benzaazoua, M. (2025) Comprehensive Hydrogeological Analysis of the Beni Amir Deposit in the Oulad Abdoun Phosphate Basin, Morocco. Mine Water and the Environment, 44, 482-501. [CrossRef]
[44] Dehkordi, M.M., Nodeh, Z.P., Dehkordi, K.S., salmanvandi, H., Khorjestan, R.R. and Ghaffarzadeh, M. (2024) Soil, Air, and Water Pollution from Mining and Industrial Activities: Sources of Pollution, Environmental Impacts, and Prevention and Control Methods. Results in Engineering, 23, Article ID: 102729. [CrossRef]
[45] Adelana, S.M.A., Olasehinde, P.I. and Vrbka, P. (2006) A Quantitative Estimation of Groundwater Recharge in Part of the Sokoto Basin, Nigeria. Journal of Environmental Hydrology, 14, 1-14.
[46] Obaje, N.G., Aduku, M. and Yusuf, I. (2013) The Sokoto Basin of Northwestern Nigeria: A Preliminary Assessment of the Hydrocarbon Prospectivity. Petroleum Technology Development Journal, 3, 71-86.
[47] Ezeh, O.C., Egwuonwu, G.N., Amadansun, C.V.O. and Umego, M.N. (2021) Source Parameter Imaging of Aeromagnetic Map of Sokoto Basin, Northwestern Nigeria. International Journal of Mechanical and Civil Engineering, 4, 1-11.
[48] Abubakar, A., Likkason, O.K., Abdulganiyu, Y., Tsafe, H.U., Mahmood, U., Balarabe, B., et al. (2024) Investigation into a Deep-Rooted Crustal Framework Deduced from Potential Field Data in Cretaceous and Tertiary Strata, Sokoto Basin NW, Nigeria. Science World Journal, 18, 526-540. [CrossRef]
[49] Wali, S.U., Gada, M.A., Dankani, I.M. and Hamisu, I. (2022) Relating Groundwater to Geology in Sokoto Basin, North-Western Nigeria Using Multivariate and Regression Analysis: Implications for Groundwater Availability. International Journal of Hydrology, 6, 57-65. [CrossRef]
[50] Bonde, D.S., Udensi, E.E. and Rai, J.K. (2014) Spectral Depth Analysis of Sokoto Basin. IOSR Journal of Applied Physics, 6, 15-21. [CrossRef]
[51] Christopher, O.A., Alaga, T.A. and Ogbole, O.J. (2016) Geological Site Suitability in Sokoto State, Northwestern Nigeria, for Phosphates Prospecting. Greener Journal of Geology and Earth Sciences, 4, 40-50. [CrossRef]
[52] Obaje, N.G., Faruq, U.Z., Bomai, A., Moses, S.D., Ali, M.S., Adamu, S., et al. (2020) Advances on the Geology and Evaluation of Potential Petroleum Systems of the Sokoto Basin in NW Nigeria. African Journal of Engineering and Environment Research, 1, 78-92.
[53] Ogbesejana, A.B., Radda, F.D., Gobirawa, N.S., Muazu, L.A., Okunola, O.J. and Bello O.M. (2024) Elemental Geochemistry of the Late Paleocene Shales from Dange and Gamba Formations, Sokoto Basin, Northwest Nigeria. Nigerian Research Journal of Chemical Sciences, 12, 23-50.
http://www.unn.edu.ng/nigerian-research-journal-of-chemical-sciences
[54] Ozumba, B.M. (2018) Potential Reservoir Rock Characteristics of the Sokoto Sector of the Iullemmeden Basin: Geological Similarities with the Niger Delta Petroleum Province and Prospects for Hydrocarbon Discovery. Archives of Petroleum & Environmental Biotechnology, 2018, 1-15.
[55] Magawata, U.Z., Olasunkanmi, N.K., Egbeyale, G.B. and Saleh, A. (2025) Structure and Mineralization Potential of the Southwestern Part of Sokoto Basin in Nigeria. Frontiers in Earth Science, 13, Article ID: 1500162. [CrossRef]
[56] Aisabokhae, J.E. (2021) Geophysical Mapping and Mineralisation Characterisation of the Mesothermal Auriferous Basement Complex in Southern Kebbi, NW Nigeria. NRIAG Journal of Astronomy and Geophysics, 10, 443-465. [CrossRef]
[57] Brookfield, M.E., Hemmings, D.P. and Van Straaten, P. (2009) Paleoenvironments and Origin of the Sedimentary Phosphorites of the Napo Formation (Late Cretaceous, Oriente Basin, Ecuador). Journal of South American Earth Sciences, 28, 180-192. [Google Scholar] [CrossRef]
[58] Davies, J.R., Wilson, D. and Williamson, I.T. (2004) Geology of the Country around Flint. Memoir of the British Geological Survey, Sheet 108 (England and Wales).
[59] Akinlalu, A.A. and Afolabi, D.O. (2018) Borehole Depth Determination to Freshwater and Well Design Using Geophysical Logs in Coastal Regions of Lagos, Southwestern Nigeria. Applied Water Science, 8, Article No. 152. [CrossRef]
[60] Gómez-Escalonilla, Diancoumba, O., Traoré, D.Y., Montero, E., Martín-Loeches, M. and Martínez-Santos, P. (2022) Multiclass Spatial Predictions of Borehole Yield in Southern Mali by Means of Machine Learning Classifiers. Journal of Hydrology: Regional Studies, 44, Article ID: 101245. [CrossRef]
[61] Kariuki, D.W. (2021) Borehole Geology and Hydrothermal Alteration Mineralogy of OW-807A, Olkaria Geothermal Field, Reports 2021 (10). United Nations Education-al, Scientific and Cultural Organisation.
[62] Kovalevsky, V.S., Kruseman, G.P. and Rushton, K.R. (2004) Groundwater Studies: An International Guide for Hydrogeological Investigations. IHP-VI Series on Groundwater, No. 3. UNESCO.
[63] Raji, B.A. and Alagbe, S.A. (2000) A Topo-Geochemical Sequence Study of Groundwater in Asa Drainage Basin, Kwara State, Nigeria. Environmental Geology, 39, 544-548. [CrossRef]
[64] Zaki, M.M., Chen, S., Zhang, J., Feng, F., Khoreshok, A.A., Mahdy, M.A., et al. (2022) A Novel Approach for Resource Estimation of Highly Skewed Gold Using Machine Learning Algorithms. Minerals, 12, Article No. 900. [CrossRef]
[65] Simoni, M.U., Drielsma, J.A., Ericsson, M., Gunn, A.G., Heiberg, S., Heldal, T.A., et al. (2024) Mass-Balance-Consistent Geological Stock Accounting: A New Approach toward Sustainable Management of Mineral Resources. Environmental Science & Technology, 58, 971-990. [CrossRef] [PubMed]
[66] Adekeye, O.A. and Akande, S.O. (2005) Mineralogical and Geochemical Studies of Phosphorite Nodules in the Dange Formation Sokoto Basin, Northwestern Niveria. Journal of Mining and Geology, 40, 101-106. [CrossRef]
[67] García-Rincón, J., Gatsios, E., Lenhard, R.J., Atekwana, E.A. and Naidu, R. (2024) Advances in the Characterisation and Remediation of Sites Contaminated with Petroleum Hydrocarbons. Environmental Contamination Remediation and Management, Switzerland. Springer.
[68] Hossain, M.E. and Islam, M.R. (2018) Drilling Engineering Problems and Solutions: A Field Guide for Engineers and Students. John Wiley & Sons, 627. [CrossRef]
[69] Kennedy, W.J. and Garrison, R.E. (1975) Morphology and Genesis of Nodular Phosphates in the Cenomanian Glauconitic Marl of Southeast England. Lethaia, 8, 339-360. [CrossRef]

Copyright © 2026 by authors and Scientific Research Publishing Inc.

Creative Commons License

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.