Introduction
Bentonite is used for various applications (e.g. Harvey and Lagaly, Reference Harvey, Lagaly, Bergaya, Theng and Lagaly2006) and mined worldwide. The properties of bentonites from different mining districts, however, vary with the geological settings throughout bentonitization (Grim and Güven, Reference Grim and Güven2011; Gilg et al., Reference Gilg, Kaufhold, Ufer and Kaufhold2022). Significant variability of the bentonite materials can sometimes be observed within a single outcrop ranging from appearance to compositional variation. Bentonites from Wyoming (USA), for example, show similar mineralogical compositions but variable amounts of exchangeable Na+ (Williams et al., Reference Williams, Elsley and Weintritt1953).
In Milos (Greece), both the type of dominating exchangeable cation and mineralogical composition show variability at least partly caused by post-formational alteration (Decher, Reference Decher1996; Kaufhold, Reference Kaufhold2001; Kaufhold and Decher, Reference Kaufhold and Decher2003). The type of compositional variability of the different deposits translates directly into the required quality control. Despite the mineralogical compositions being significantly variable, most Milos bentonites have a high smectite content. In Milos, the most important parameter for quality control is the amount of naturally present exchangeable Na+, as it determines the amount of Na2CO3 that has to be added technically during production (for most products). Conversely, in the Bavarian mining district, the exchangeable cation population was found to be quite constant (Kaufhold, Reference Kaufhold, Decher and Dohrmann1998; Köster et al., Reference Köster, Williams, Kudejova and Gilg2019; Supplementary material) while the smectite content was observed to vary significantly (Kaufhold et al., Reference Kaufhold, Dohrmann and Decher2003). In Bavarian deposits, therefore, the smectite content is an important quality-determining parameter while other technologically important compounds, e.g. carbonate content, are usually easier to measure accurately. Typically, it is sufficient to test the smectite content, e.g. by the traditional methylene blue (MB) method or by more straightforward cation exchange capacity (CEC) methods such as Cu-trien (Kaufhold and Dohrmann, Reference Kaufhold and Dohrmann2003). For the calculation of the smectite content (Swt, in wt.%) from CEC values, the molar weight (MFU) of one formula unit (FU), and the smectite layer charge and the portion of variable charge, have to be known (Vogt and Köster, Reference Vogt and Köster1978; Lagaly, Reference Lagaly and Mermut1993; Kaufhold et al., Reference Kaufhold, Dohrmann, Ufer and Meyer2002):

where CECmeas is the experimentally measured CEC value, whereas CECS is a theoretical CEC of pure smectite, assuming that smectite is the only compound contributing to bulk sample CEC:

In this study, layer charge (LC, given per FU) is used here as the charge generated only by substitution of heteroatoms in the 2:1 layer and neutralized by interlayer cations, which can then be measured by a probe specific to the interlayer space (Christidis et al., Reference Christidis, Chryssikos, Derkowski, Dohrmann, Eberl, Joussein and Kaufhold2023). The variable charge (VC) occurring on the smectite crystallite edges, is pH dependent, and contributes to the CECS. The accuracy of CEC-based methods for the determination of the smectite content (Eqns 1 and 2) for bentonite mining quality control, therefore, depends heavily on the variability of the layer charge (Kaufhold, Reference Kaufhold2005). In smectites, LC varies from 0.2 to 0.6 eq/FU. Köster (Reference Köster1977) suggested use of the value 0.34 eq/FU for smectites if the LC was not determined analytically. With only a few exceptions (e.g. Christidis and Dunham, Reference Christidis and Dunham1993; Christidis and Dunham, Reference Christidis and Dunham1997; Christidis, Reference Christidis2006), a statistically sound study of the smectite LC distribution in a bentonite mine or even across a mining district, however, has not been conducted so far because both classical methods, the structural formula method (SFM) and the alkylammonium method (AAM), are laborious and/or erroneous (e.g. Laird, Reference Laird and Mermut1994; Christidis et al., Reference Christidis, Chryssikos, Derkowski, Dohrmann, Eberl, Joussein and Kaufhold2023). One of few exceptions was provided by Christidis (Reference Christidis2006) who used microprobe analysis to investigate Cyprus bentonites. Based on the chemical data, structural formulae were derived from which the LC was deduced. Those authors found significant variation of the LC even within different deposits (overall ranging from mean 0.38 to mean 0.59 eq/FU). Those authors, however, used a method that probes single smectite particles, in contrast to most other methods which provide a bulk-sample parameter. If considering the individual particle measurements, the LC variability is even broader (0.26–0.66 eq/FU). In contrast, Chiou et al. (Reference Chiou, Dohrmann, Kaufhold, Plötze, Stucki, Ufer and Kaufhold2025) used AAM to investigate five bentonites from different quarries of the Milos deposit and found quite similar values (0.31–0.33 eq/FU). In the same study, five different Wyoming bentonites were additionally investigated, which are all low charged (as widely accepted for Wyoming bentonites) but still showed a large range of LC (0.18, 0.21, 0.25, 0.25, and 0.27 eq/FU). These results indicate that the LC can be more or less constant in a bentonite mining district or show comparatively large variation. The first comprehensive study on LC of Bavarian bentonites was published by Vogt and Köster (Reference Vogt and Köster1978) who, based on 14 Bavarian samples, found low LC variability ranging from 0.36 to 0.39 eq/FU (determined by SFM) with one sample showing a lower LC accounting for 0.33 eq/FU.
Kuligiewicz et al. (Reference Kuligiewicz, Derkowski, Emmerich, Christidis, Tsiantos, Gionis and Chryssikos2015) introduced an efficient method for measuring the LC of smectite which is based on the determination of the position of a mid-infrared (MIR) band of OD stretching after a sample is saturated with D2O. The D2O saturation enables the use of water as an intrinsic probe of the wettable interlayer. The band position depends directly on the LC and hence the method can either be calibrated against SFM or AAM (LC (AAM) or (LC (SFM)). The systematic difference between values derived by both methods was previously discussed (see Kaufhold et al., Reference Kaufhold, Dohrmann, Stucki and Anastácio2011). However, in contrast to SFM, neither the O-D- nor the AAM-derived LC values are affected by the variable charge (Kuligiewicz et al., Reference Kuligiewicz, Derkowski, Emmerich, Christidis, Tsiantos, Gionis and Chryssikos2015). Determining the LC of one sample using the O-D method takes only ~40 min. Using this method, therefore, enables a significant set of bentonite samples to be measured quickly and precisely from one or several bentonite mines within a bentonite mining district. Hence the aim of the present study was based on the possibility of measuring the LC of large sample sets resulting from the recently established O-D method. None of the other laborious LC methods could facilitate the study of this quantity of samples. The aim of the present study was as follows: (1) use of the O-D method in a comprehensive investigation of the LC distribution over the Bavarian bentonite mining district; and (2) assess the error of smectite content calculations based on CEC data resulting from the actual range of LC values.
Geological setting
The Bavarian bentonite mining district is located within the eastern part of the Molasse basin also called the North Alpine Foreland Basin (NAFB; Schlunegger and Kissling, Reference Schlunegger, Kissling, Rosenberg and Bellahsen2024). The NAFB is an asymmetric basin across the northern margin of the Alps and encompasses an area 1000 km long and 10–200 km wide (Rocholl et al., Reference Rocholl, Schaltegger, Gilg, Wijbrans and Böhme2018). It was formed as a result of the Alpine orogenesis in the mid-Cenozoic and ranges from the Rhone valley to the alpine foothills of southern Germany and Austria (Homewood et al., Reference Homewood, Allen, Williams, Allen and Homewood1986; Schlunegger et al., Reference Schlunegger, Matter, Burbank, Leu, Mange and Matyas1997). Significant uplift during the Alpine orogeny resulted in a lowering of the Molasse basin and the sedimentation regime changed from marine to terrigenous in the Upper Oligocene (Probst, Reference Probst1986), then back to marine, and back again to terrigenous in the early Pannonian (Vogt, Reference Vogt1980). The Obere Süßwassermolasse (OSM) represents the end of the second marine to terrigenous clastic sedimentation cycle. The complete retreat of the Paratethys from the west side of the basin by the Middle Miocene (~16.3 Ma) (Reichenbacher et al., Reference Reichenbacher, Krijgsman, Lataster, Pipperr, Van Baak, Chang and Bachtadse2013) caused a transition to fluvial sedimentation at the axis of the basin and alluvial fan sedimentation in the southern rim (Kuhlemann and Kempf, Reference Kuhlemann and Kempf2002). Multiple methods which include litho-, magneto-, and bio-stratigraphic and isotopic dating (e.g. Dehm, Reference Dehm1951; Gubler et al., Reference Gubler, Meier and Oberli1992; Bolliger, Reference Bolliger1994; Heissig, Reference Heissig1997; Kälin and Kempf, Reference Kälin and Kempf2009; Abdul Aziz et al., Reference Abdul Aziz, Böhme, Rocholl, Zwing, Prieto, Wijbrans and Bachtadse2008; Gubler, Reference Gubler2009; Abdul Aziz et al., Reference Abdul Aziz, Böhme, Rocholl, Prieto, Wijbrans, Bachtadse and Ulbig2010) have been used to obtain stratigraphic subdivisions of the OSM sediments in order to overcome the lateral variability in sedimentation dynamics caused by the complex geometry of the basin (Reichenbacher et al., Reference Reichenbacher, Krijgsman, Lataster, Pipperr, Van Baak, Chang and Bachtadse2013). For a recent and detailed analysis of the OSM stratigraphic subdivisions the reader is referred to Abdul Aziz et al. (Reference Abdul Aziz, Böhme, Rocholl, Zwing, Prieto, Wijbrans and Bachtadse2008), Köster and Gilg (Reference Köster and Gilg2015), Rocholl et al. (Reference Rocholl, Schaltegger, Gilg, Wijbrans and Böhme2018), and references within, as only the main subdivision features relevant for the context of this study are discussed below.
Briefly, the oldest sediments of the OSM are termed the Süßwasserschichten and consist of fine gravels, sand, and silt (Fig. 1). They are overlain by the thick (10–50 m) coarse gravel of the Nördlicher Vollschotter (NV) which in turn is overlain by the Sand-Mergel-Decke (SMD), an informal local lithostratigraphic unit ~15–40 m thick, consisting of fine fluvial to lacustrine gravel, sand, and carbonate-rich silts (sandy marls). The Brock horizon is an ejecta layer situated roughly between the NV and the SMD, containing up to 20 cm angular Jurassic limestone blocks that originated from the ‘Nördlinger Ries’ impact crater area ~100 km from the bentonite mining district (Unger et al., Reference Unger, Fiest and Niemeyer1990; Ulbig, Reference Ulbig1999; Abdul Aziz et al., Reference Abdul Aziz, Böhme, Rocholl, Zwing, Prieto, Wijbrans and Bachtadse2008). The bentonite deposits are situated over the irregular erosional surface of the SMD. The bentonite deposits are covered by Miocene Hangend-Serie or Jüngere Grobkies coarse clastic sediments and are unconformably covered by Pleistocene loess.

Figure 1. Schematic representation of the informal lithostratigraphy of the ‘Obere Süßwassermolasse’ (Upper Freshwater Molasse) near Mainburg-Landshut, modified after Ulbig (Reference Ulbig1994).
In Bavaria, bentonites have been mined for more than 100 years and over four different regions: Malgersdorf, Landshut, and East and West of Augsburg (Gilg and Ulbig, Reference Gilg and Ulbig2017). Samples investigated in the present study were all from the Landshut region.
The formation of the Bavarian bentonites was investigated by Unger and Niemeyer (Reference Unger and Niemeyer1985), Unger et al. (Reference Unger, Fiest and Niemeyer1990), and Ulbig (Reference Ulbig1994). Bavarian bentonite mines are distributed over an area of about 400 km2 (Köster, Reference Köster2018). They consist of flat, lenticular clay bodies varying in thickness from 0.5 to 2 m (may reach 8 m). Beds with a thickness of up to 4 m are rare but of particular interest for mining. Occasionally, even thicker beds are observed in which incompletely altered tuffite, sandy clay, or silicified bentonite can be found (‘Harte Platte’ = ‘hard plate’) (Vogt and Köster, Reference Vogt and Köster1978; Köster and Gilg, Reference Köster and Gilg2015). Horizontal areas of the clay bodies range from 100 m2 to 2 km2. The tuffs/ash are considered a precursor of the bentonites (Ulbig, Reference Ulbig1994, Abdul Aziz et al., Reference Abdul Aziz, Böhme, Rocholl, Zwing, Prieto, Wijbrans and Bachtadse2008, Rocholl et al., Reference Rocholl, Schaltegger, Gilg, Wijbrans and Böhme2018) and Köster et al. (Reference Köster, Hölzl and Gilg2017) found the age of bentonitzation to be ~14.7±4.1 Ma. The thick bentonite beds formed within the braided river system of the upper fluviatile NV sand and gravel, and above the SMD floodplain sediments (Ulbig, Reference Ulbig1994; Schmid, Reference Schmid2003; Maurer and Buchner, Reference Maurer and Buchner2007) resulting in the formation of oxbow lakes into which ash accumulated and resedimented (Rocholl et al., Reference Rocholl, Schaltegger, Gilg, Wijbrans and Böhme2018; Köster et al., Reference Köster, Hölzl and Gilg2017). The rhyolitic volcanic ash then altered in situ to smectitic clay minerals (Ulbig, Reference Ulbig1994; Abdul Aziz et al., Reference Abdul Aziz, Böhme, Rocholl, Zwing, Prieto, Wijbrans and Bachtadse2008; Abdul Aziz et al., Reference Abdul Aziz, Böhme, Rocholl, Prieto, Wijbrans, Bachtadse and Ulbig2010; Köster and Gilg Reference Köster and Gilg2015; Bauer et al., Reference Bauer, Vennemann and Gilg2016). Vogt (Reference Vogt1980) analyzed the chemical composition of the smectites from different locations and found that despite the spatial variability, the smectite samples had a typical montmorillonitic character with similar formulae. From a geochemical point of view and based on the current knowledge about the formation of Bavarian bentonites, this low variability can be explained by similar conditions throughout bentonitization and the absence of post-formational alteration. The protolith material of the Landshut bentonites is dominated by distal rhyolitic volcanic glassy ash from a single Plinian eruption, and a Ca-Mg-rich reduced meteoric water at low temperatures (Ulbig, Reference Ulbig1994; Gilg and Ulbig, Reference Gilg and Ulbig2017; Bauer et al., Reference Bauer, Vennemann and Gilg2016). Nevertheless, different carbonate contents and colors of the bentonite point towards locally different geochemical environments (including Eh) which in turn affects Fe mobility and charge in the octahedral sheet, which may also lead to different LC. The composition and, thus, LC of smectite from different locations are therefore expected to be similar. As such, any difference in the CEC between bentonite deposits would originate exclusively from the content of detrital material in the bentonite. The average temperature in this region is currently about 10°C with 1022 mm of annual precipitation. The soils of this region are used for livestock farming and growing hops.
Materials and methods
Materials
The Landshut region of the Bavarian bentonite mining district was selected for the investigation of the LC distribution as several samples from eight different mines were collected over several years, which is the reason for the different sample assignment in Fig. 2. These samples were probed to establish a new in situ quality control method that is based on measuring the water content and the specific electrical resistivity, which were correlated with methylene blue (MB) (Kaufhold et al., Reference Kaufhold, Decher and Dohrmann1998) and CEC (Kaufhold and Penner, Reference Kaufhold and Penner2006; Kaufhold et al., Reference Kaufhold, Grissemann, Dohrmann, Klinkenberg and Decher2014) values in different clay raw materials. An example for sampling in the Osterwaal deposit is given in Fig. 3.

Figure 2. List of samples and locations of the mines (green circles) in the Landshut bentonite mining district. Larger towns are shown for orientation (red circles).

Figure 3. Example for sampling a small profile in the Osterwaal deposit. (A) Sketch of the different materials, sample numbers, and smectite contents based on MB (blue numbers), and samples included in the present study are marked with a red circle. (B) Example of a sampling spot which was used for geoelectrical investigation and sample collection (about 50 g wet state).
Mineralogical and geochemical data were not available for all samples but CEC values were available. The LC values measured with AAM were available for five samples and were published in Kaufhold (Reference Kaufhold2001). The latter were used for the comparison of the AAM and O-D method.
Layer-charge determination
The O-D method (Kuligiewicz et al., Reference Kuligiewicz, Derkowski, Emmerich, Christidis, Tsiantos, Gionis and Chryssikos2015) is based on a linear regression between the position of the sharp O-D stretching band (νOD) at ~2700 cm–1 and the layer charge. The O-D method is insensitive to interlayer chemistry as long as it is filled with Na+, Li+, Ca2+, and/or Mg2+ cations which allow full interlayer expansion to two-water layers at high humidity (not suitable for use with K+ or Cs+ cations) and is independent of tetrahedral or octahedral charge locations. The best precision, however, is obtained when comparing samples with the same interlayer cation composition. As the O-D method measures only the wettable smectitic surfaces, it is not affected by bulk rock impurities or variable charge at crystallite edges. Therefore, no influence is observed from particle edges or the illitic interlayer charge as only the basal surface charges of fundamental particles are probed by the O-D method (Kuligiewicz et al., Reference Kuligiewicz, Derkowski, Środoń, Gionis and Chryssikos2018). This means that the LC component is not measured for any 2:1 layers that are not fully expandable (only smectitic layers). In the case of partially collapsed interlayer and interstratified illite-smectite phase, the LC determined should be considered as the wettable surface charge.
Approximately 2 mg of sample material was added to a 1.5 mL vial, saturated with ~1 mL 99.99% pure D2O water (Sigma-Aldrich), and ultrasonicated for 4 s. A pipette was used to place a drop of the suspension onto the ATR crystal of a Nicolet 6700 by ThermoScientific (Waltham, MA, USA) FT-IR instrument (Bukas et al., Reference Bukas, Tsampodimou, Gionis and Chryssikos2013). A custom-made cap was used to cover the suspension during drying, to a thin film, by N2 gas flowing through the cap. An N2 gas and D2O vapor mix, controlled at between 60 and 80% (±2%) relative humidity (RH), was delivered to the sample by a custom gas blending system. FT-IR spectra were collected in the range of 4000 to 580 cm–1 at RH values of 60, 65, 70, 75, and 80%. The LC for each of the RH values was calculated from the IR measurements with LC values reported as calibrated against the AAM (Eqn 3; = LC (AAM)), although calibration for the SFM is also available (see Kuligiewicz et al., Reference Kuligiewicz, Derkowski, Emmerich, Christidis, Tsiantos, Gionis and Chryssikos2015). In the present study, however, all O-D derived LC values were calibrated against the alkylammonium method. Calculation of the final LC value and the standard deviation of the measurement for each sample was done by averaging the LC values collected at each step differing by 5% RH:

where LC is expressed in electrons per formula unit (FU) and νOD is the measured νO-Dw in cm–1. The σ of 0.02 and R 2 of 0.96 were calculated from this study, as that reported by Kuligiewicz et al. (Reference Kuligiewicz, Derkowski, Emmerich, Christidis, Tsiantos, Gionis and Chryssikos2015) are σ = 0.01 and R 2 = 0.92.
Spectra processing was done using OMNIC software following Kuligiewicz et al. (Reference Kuligiewicz, Derkowski, Emmerich, Christidis, Tsiantos, Gionis and Chryssikos2015).
Results and Discussion
Layer-charge range and distribution
Examples of spectra for the samples with the highest and lowest LC values are shown in Fig. 4. The spectra at RH of 80 and 60% represent, respectively, the highest and lowest RHs of the five spectra collected for calculation of one LC value. The band intensity maximum between 2650 and 2700 cm–1 shows the position of the O-D stretching band correlated with LC and the inset depicts its second derivative (Fig. 4).

Figure 4. Spectra of the sample with the highest and lowest LC values of the 40 samples analyzed. Spectra taken at two of the five relative humidities (RHs) are shown for each sample. The inset shows the spectra of the sample measured with the O-D method and the highest LC at 80% RH along with its second derivative.
For five of the samples used in the present study, LC values determined previously by AAM were available (Kaufhold, Reference Kaufhold2001; Kaufhold et al., Reference Kaufhold, Dohrmann, Ufer and Meyer2002). These values were compared with the results obtained by the O-D method in the present study (Table 1). The AAM values were slightly but systematically higher: 0.01 eq/FU for most of the samples. Although the sample from the Rehbach mine showed a slightly higher deviation of ~0.02 eq/FU, for other samples of the Rehbach deposit (RAB, RAN; Kaufhold et al., Reference Kaufhold, Dohrmann, Ufer and Meyer2002) a value of 0.32 eq/FU was published, which would fit perfectly the trend shown by the others. Given that the analysis was done with different methods and by different operators, and also included a two-decade time difference between the analyses, the LC values measured for these samples are quite consistent.
Table 1. Comparison of LC values from five samples previously published by Kaufhold et al. (Reference Kaufhold, Dohrmann, Ufer and Meyer2002) with the LC values measured by the O-D (AAM) method (O-D calibrated against the alkylammonium method) two decades later (current study). Values for LC O-D (AAM) are shown here with 3 decimal places for ease of individual comparison but standard reporting is recognized to 2 decimal places, as the O-D methodological work of Kuligiewicz et al. (Reference Kuligiewicz, Derkowski, Emmerich, Christidis, Tsiantos, Gionis and Chryssikos2015) reports a sigma value of 0.01 (AMM). As seen below, the precision of the O-D method for the bentonite samples here, was high.

All the measured LC values are summarized in Fig. 5. The lowest measured value was 0.29 eq/FU; however, most of the samples showed values between 0.29 and 0.30 eq/FU and only about 10% of the samples showed LC values above 0.31 eq/FU. The range of the measured values was surprisingly narrow but is in good agreement with the values published by Kaufhold et al. (Reference Kaufhold, Dohrmann, Ufer and Meyer2002) for eight Bavarian bentonites (0.29–0.33 eq/FU) and Wolters et al. (Reference Wolters, Lagaly, Kahr, Nueesch and Emmerich2009) for two Bavarian bentonites (AAM, 0.30 and 0.32 eq/FU). Based on AAM measurements, Vogt and Köster (Reference Vogt and Köster1978) reported higher LC values for Bavarian montmorillonites ranging from 0.33 to 0.39 eq/FU; however, the reason for this deviation from other reported values remains unclear. Müller-Vonmoos and Kahr (Reference Müller-Vonmoos and Kahr1983) reported a LC for Montigel (a technical grade Bavarian bentonite from former Südchemie AG) of 0.28 eq/FU, which – based on the present study – is not representative of the average Bavarian bentonite but instead represents a low LC end member of Bavarian bentonite. The bentonite sample set described in Chiou et al. (Reference Chiou, Dohrmann, Kaufhold, Plötze, Stucki, Ufer and Kaufhold2025) contained two Bavarian bentonites. Sample B25 has a LC of 0.30 eq/FU, which based on the present study can be considered typical for Bavarian bentonites. Sample B16, however, showed a LC value of 0.27 eq/FU which is comparatively low. This suggests that materials with LC values outside the range measured in this study may exist but should not be considered representative of the entire bentonite mining district.

Figure 5. The distribution of smectite LC (calibrated against AAM) measured in the present study (A); compared with the CEC distribution in the same samples, from Kaufhold et al. (Reference Kaufhold, Dohrmann, Ufer and Meyer2002) (B); and the comparison of both the values on axes covering a typical range of LC and CEC values in bentonite (C).
In Fig. 6, all LC results are averaged according to the different bentonite mines (numbers represent the number of samples/measurements per mine). The hypothesis tested was that different mines could show different typical LC ranges. Most of the mines showed values of ~0.30 eq/FU (Bergham, Berg, Obergolzaberg, Osterwaal 1, Oberpriel, Rehbach). Two mines showed lower values, closer to 0.29 eq/FU (Osterwaal 2 and Priel) and only Martinszell showed a value closer to 0.32 eq/FU. The standard deviation over the samples analyzed ranged from 0.001 to 0.007 eq/FU. For the location which consisted of only a single sample analysis, the average standard deviation of 0.005 was assumed (Fig. 6). Overall, the differences among locations are minimal, and most of them overlap within the measurement range. Moreover, the highest value that was measured in the Martinszell deposit is based on only one sample and hence should not be interpreted further. The series Osterwaal 1 and Osterwaal 2, however, were both obtained in the same mine while the Osterwaal 2 sample series was collected several years later than the other. Although the differences between average values may represent different extents of bentonitization and possible post-formational overprints (Christidis, Reference Christidis2008), the overlapping measurement errors of LC for the two Osterwaal series indicate that LC differences within a single mine can be as significant as between different mines. In conclusion, no LC values or ranges can be given for a specific bentonite mine in the Bavarian mining district. Fortunately, the greater quantity of measurements possible with the use of the O-D method now enables the typical LC range of Bavarian bentonites to be determined based on solid statistics.

Figure 6. Comparison of LC values (O-D calibrated against AAM) from the different bentonite mines. The numbers above the bars denote the number of samples/measurements (n). Standard deviation of mines with n>4 was 0.005 eq/FU on average.
Smectite content calculation
In the Bavarian bentonite mining district, quality control is based mainly on estimation of the smectite content, which in turn is based on MB or CEC values (Kaufhold and Dohrmann, Reference Kaufhold and Dohrmann2003) and reference samples of known composition. CEC values can, however, also be used for the calculation of absolute smectite content values without the need to use reference materials, in contrast to the procedure used to obtain MB values. The distribution of CEC values in Bavarian bentonites (Fig. 5B) compared with the variability of LC (Fig. 5C) implies that smectite crystallochemical composition varies far less than CEC; hence, it is not the smectitic LC but the smectite content which controls bulk rock CEC in the bentonites studied.
In order to calculate the smectite content from CEC data, information or assumptions about the LC, the molar mass of a formula unit, and the portion of variable charge are required (Eqns (1) and (2); Lagaly, Reference Lagaly and Mermut1993). In a test shown in Table 2, these parameters were varied systematically in order to assess the error of smectite content calculations based on the CEC data from Bavarian bentonites. Based on the results of the present study, one can assume smectites of most Bavarian bentonites to be represented by a LC range from 0.29 to 0.31 eq/FU. The CEC of Bavarian bentonites typically ranges from 60 to 80 meq 100 g–1 (Fig. 5B,C), which corresponds to the reported average of 70 meq 100 g–1 in the first rows of Table 2. The molar mass of a formula unit (MFU) depends largely on the Fe content (Środoń and McCarty, Reference Środoń and McCarty2008) but also some single remaining water molecules after 105°C drying in an oven for 4 days (which is the reference state of the CEC measurements) must also be considered. Based on the moderate Fe content in the Bavarian bentonites and, thus, low Fe in the montmorillonite, structure MFU results in about 370 g mol–1 and is considered to be a suitable reference value for Bavarian bentonites (Kaufhold, Reference Kaufhold2005). To be able to compare the CEC with LC values that represent permanent charge sites only, the CEC has to be reduced by the amount of CEC index cation adsorbed to variable charge sites (Eqn 2). This amount depends on the pH throughout the CEC measurements, which is discussed in detail by Kaufhold and Dohrmann (Reference Kaufhold and Dohrmann2013). Based on the present study in Table 2, a 10% variable charge was used as a reference value. Based on the reference values and using Eqns 1 and 3, a variation of the resulting smectite content of ±3 mass% was calculated. In the case of higher CEC values (90 meq 100 g–1, rows 4–6 in Table 2), this value may increase to 3.2–3.4%. In rows 7–9 of Table 2, the MFU was varied and the calculations based on the average LC of 0.30 eq/FU. Varying Fe content and/or amount of relict water may cause some variation of the MFU. The effect of lower and higher MFU on the smectite content resulted in an error of ±2 mass%. The MFU range of 360–380 g mol–1, however, is large compared with the actual variability. More detailed discussions on the variability of MFU can be found elsewhere (Kaufhold et al., Reference Kaufhold, Dohrmann, Ufer and Meyer2002; Kaufhold et al., Reference Kaufhold, Plötze, Klinkenberg and Dohrmann2013). The effect of variable MFU on the calculated smectite content is considered to be low. On the other hand, a significant effect of the variable charge on the calculated smectite content was observed (±4 mass%) which shows the importance of considering the pH throughout CEC measurements. Overall, the surprisingly narrow range of LC values can be considered beneficial with respect to the accuracy of smectite content calculations for the Bavarian bentonites based on CEC values. The error in smectite content obtained based on the narrow LC range measured and realistically assumed from the other components of Eqn 2 is no greater than the error observed in the best analytical approaches for smectite quantification in clay-bearing rocks (Omotoso et al., Reference Omotoso, McCarty, Hiller and Kleeberg2006).
Table 2. Effect of relevant parameter variation on the smectite content as calculated based on CEC values

The values do not correspond to specific samples but were selected in order to represent the realistic range.
The assumed 10% of variable charge and 370 MFU value were used to calculate the smectite content in all the samples analyzed, taking their individually measured LC (this study) and CEC from Kaufhold et al. (Reference Kaufhold, Dohrmann, Ufer and Meyer2002). With one exception (Oberpriel), all other samples showed standard deviation of smectite content in 5–8% range (Fig. 7). The deposits with only one sample analyzed were assumed to have a standard deviation of 6%.

Figure 7. Comparison of calculated per cent smectite content for samples from the different locations. The numbers above the bars denote the number of samples (n).
Conclusions
The new O-D method can be used for comprehensive characterizations of LC distribution in bentonite mining districts and was demonstrated in this study based on the Bavarian bentonites mined in eight different open pits or parts of mines. The O-D-LC values showed low variation (0.29–0.32 eq/FU) and the values are in agreement with published data of single Bavarian bentonite samples. While the previously reported range was slightly larger, it should be noted that these values were analyzed by AAM in different laboratories and the measurement itself is very laborious (Janek and Smrčok, Reference Janek and Smrčok1999). The narrow LC range has positive implications for the accuracy of determining smectite content calculated from CEC data, as routinely done for quality control of Bavarian bentonites. The remaining error of the CEC based smectite contents resulting from LC variations accounts for, on average, ±3 wt.%, which is as good or better than smectite quantification using structural analysis methods. The similarity of LC among the deposits implies a similar origin, i.e. tephra composition and formation conditions of smectite, regardless of the smectite content (Vogt, Reference Vogt1980).
Supplementary material
The supplementary material for this article can be found at http://doi.org/10.1017/cmn.2025.10013.
Data availability statement
All data are available in the online Supplementary Materials and by request.
Acknowledgements
None.
Author contribution
N.J. Kanik: writing, data acquisition, analysis, critical review; S. Kaufhold: concept, writing, data acquisition, analysis, critical review, resources; A. Derkowski: concept, critical review, resources; R. Dohrmann: concept, critical review, resources.
Financial support
Project funding from BGR funding No. 204-10121999
Competing interests
The authors declare none.