Modeling complex aquifer systems: a case study in Baton Rouge,Louisiana (USA)

Hydrogeology Journal - This study targets two challenges in groundwater model development: grid generation and model calibration for aquifer systems that are fluvial in origin. Realistic...     

Trace element remobilization following the resuspension of sediments under controlled redox conditions: City Park Lake,Baton Rouge,LA

The sediments of City Park Lake (Baton Rouge, Louisiana, USA) are characterized by Pb concentrations that exceed the corresponding National Oceanic and Atmospheric Administration probable effect level by more than a factor of two and concentrations of As, Cd, Cu, Ni and Zn that exceed the corresponding threshold effect levels. The sediments of the lake are likely to be dredged because of concern over shoaling of the lake. There is concern that exposing these sediments to an oxidizing environment could mobilize toxic metal(loid)s that are effectively sequestered under reducing conditions. To address this concern, mixtures of sediments collected from seven locations in City Park Lake were incubated in microcosms under controlled oxidizing and reducing conditions. The slurries were sampled at time intervals over a period of nearly 750 h, and the dissolved Fe, Mn and trace metal concentrations determined. Eh and pH were negatively correlated, with Eh declining by 600 mV per unit increase in pH. Of the toxic metals, only the concentration of Cu was positively correlated with Eh, apparently because of the sequestration of Cu as a very insoluble sulfide precipitate under low Eh conditions. Precipitation of the other metals as sulfides was evidently insignificant because of the low S content of the sediments. Concentrations of the remaining potentially toxic metals were negatively correlated with Eh, apparently because of their sequestration as oxides and/or adsorption to Fe oxides or oxy-hydroxides. Concentrations of Ca, Mg, and Sr were positively correlated with Eh, presumably because of their sequestration in carbonate precipitates under low Eh/high pH conditions.     

Theoretical and measured aeolian sand transport on a barrier island, Louisiana, USA

Over the past 100 years, the Isles Dernieres, a low lying barrier island chain along the coast of central Louisiana, Usa , has undergone more than 1 km of northward beach face retreat with the loss of 70% of its surface area. The erosion results from a long term relative sea level rise coupled with day to day wind and wave action that ultimately favours erosion over deposition. At a site in the central Isles Dernieres, 8 days of wind and beach profile measurements during the passage of one winter cold front documented aeolian erosion and deposition patterns under both onshore and offshore winds. For offshore winds, the theoretical erosion rate, based on wind shear velocity, closely matched the measured erosion rate; for onshore winds, the theoretical rate matched the measured rate only after being corrected by a factor that accounted for beach face morphology. In late February 1989, a strong cold front moved into coastal Louisiana. That cold front stalled over the Gulf of Mexico, resulting in 4 days of strong northerly winds at a study site on the Isles Dernieres. During those 4 days, the wind moved sand from the backshore to the upper beach face. When the cold front finally moved out of the area, the wind shifted to the south and decreased in strength. The onshore wind then restored some of the upper beach face sand to the backshore while increased wave activity moved the rest into the nearshore. The theoretical estimate of 1·28 m³ m^?1 for the rate of sand transport by the northerly wind compares well with the measured backshore erosion rate of 1·26 m³ m^?1, which was determined by comparing beach profiles from the start and end of the period of northerly winds. The theoretical estimate of 0·04 m³ m^?1 for the rate of sand transport by the southerly wind, however, is notably less than the measured rate of 0·45 m³ m^?1. The large discrepancy between the two rates can be explained by a difference in the shear velocity of the wind between the beach face, where the erosion occurred, and the backshore, where the wind stress was measured. Using an empirical relationship for the wind shear drag coefficient as a function of coastal environment, the theoretical estimate for the rate of sand transport by the southerly wind becomes 0·44 m³ m^?1     

Sources of elevated salinity in the Mississippi River Alluvial Aquifer, south-central Louisiana, USA

Salinization is a process impacting groundwater quality and availability across much of southern Louisiana, USA. However, a broad divergence of opinion exists regarding the causes of this elevated salinity: updip-migrating marine waters from the Gulf of Mexico, saline fluids migrating up fault planes, movement of water from salt domes, and/or remnant seawater from the last major marine transgression. The Mississippi River Alluvial Aquifer (MRAA) in south-central Louisiana is recharged by the Mississippi River, and there are discharge zones to the west and east. Recharge waters from the Mississippi River are fresh, but Cl⁻ levels in the western portions of the aquifer are as high as 1000 mg/L. The aquifer is an important source of water for several municipalities and industries, but prior to this study the source(s) of the elevated salinity or whether the salinization can be remediated had not been determined.The low Br/Cl ratios in the groundwaters are consistent with a saline endmember produced by subsurface dissolution of salt domes, not a marine source. The H and O isotopic systematics of the aquifer waters indicate meteoric sources for the H₂O, not marine waters or diagenetically-altered deep brines. The westward salinization of aquifer water represents a broad regional process, instead of contamination by point sources. Mapping of spatial variations in salinity has permitted the identification of specific salt domes whose subsurface dissolution is producing waters of elevated salinity in the aquifer. These include the Bayou Choctaw and St. Gabriel domes, and possibly the Bayou Blue dome. Salinization is a natural, on-going process, and the potential for remediation or control is slight, if not non-existent.     

The relative utility of foraminifera and diatoms for reconstructing late Holocene sea-level change in North Carolina, USA

Foraminifera and diatoms preserved in salt-marsh sediments have been used to produce high-resolution records of Holocene relative sea-level (RSL) change. To determine which of these microfossil groups is most appropriate for this purpose we investigated their relative utility from salt marshes in North Carolina, USA. Regional-scale transfer functions were developed using foraminifera, diatoms and a combination of both (multi-proxy) from three salt marshes (Oregon Inlet, Currituck Barrier Island and Pea Island). We evaluated each approach on the basis of transfer-function performance. Foraminifera, diatoms and multi-proxy-based transfer functions all demonstrated a strong relationship between observed and predicted elevations (r²_jack > 0.74 and RMSEP < 0.05 m), suggesting that they have equal utility. Application of the transfer functions to a fossil core from Salvo to reconstruct former sea levels enabled us to consider relative utility in light of ‘paleo-performance’. Fossil foraminifera had strong modern analogues, whilst diatoms had poor modern analogues making them unreliable. This result reflects the high diversity and site-specific distribution of modern diatoms. Consequently, we used foraminifera to reconstruct RSL change for the period since ∼ AD 1800 using a ²¹⁰Pb- and ¹⁴C-based chronology, and we were able to reconcile this with tide-gauge records.     

Three‐dimensional architecture and hydrostratigraphy of cross‐cutting buried valleys using airborne electromagnetics,glaciated Central Lowlands,Nebraska, USA

Buried valleys are characteristic features of glaciated landscapes, and their deposits host important aquifers worldwide. Understanding the stratigraphic architecture of these deposits is essential for protecting groundwater and interpreting sedimentary processes in subglacial and ice‐marginal environments. The relationships between depositional architecture, topography and hydrostratigraphy in dissected, pre‐Illinoian till sheets is poorly understood. Boreholes alone are inadequate to characterize the complex geology of buried valleys, but airborne electromagnetic surveys have proven useful for this purpose. A key question is whether the sedimentary architecture of buried valleys can be interpreted from airborne electromagnetic profiles. This study employs airborne electromagnetic resistivity profiles to interpret the three‐dimensional sedimentary architecture of cross‐cutting buried valleys in a ca 400 km² area along the western margin of Laurentide glaciation in North America. A progenitor bedrock valley is succeeded by at least five generations of tunnel valleys that become progressively younger northward. Tunnel‐valley infills are highly variable, reflecting under‐filled and over‐filled conditions. Under‐filled tunnel valleys are expressed on the modern landscape and contain fine sediments that act as hydraulic barriers. Over‐filled tunnel valleys are not recognized in the modern landscape, but where they are present they form hydraulic windows between deep aquifer units and the land surface. The interpretation of tunnel‐valley genesis herein provides evidence of the relationships between depositional processes and glacial landforms in a dissected, pre‐Illinoian till sheet, and contributes to the understanding of the complex physical hydrology of glacial aquifers in general.     

Application of stable carbon isotopes for reconstructing salt‐marsh floral zones and relative sea level,New Jersey,USA

We investigated use of δ¹³C in bulk organic sediment to define the botanical origin of samples preserved in coastal sediment as a means to reconstruct relative sea level in New Jersey, USA. Modern transects at three sites demonstrated that low and high salt‐marsh floral zones dominated by C₄ species (Spartina alterniflora and Spartina patens) were associated with sediment δ¹³C values between ?18.9‰ and ?15.8‰ and occurred from mean tide level (MTL) to mean higher high water (MHHW). Brackish transitional settings vegetated by Phragmites australis with Iva fructescens and Typha sp. (C₃ species) and freshwater upland samples (C₃ species) were characterized by bulk sediment δ¹³C values of ?27.0‰ to ?22.0‰ and existed above MHHW. Parallel transects at one site suggested that intra‐site variability was not discernible. The utility of δ¹³C values for reconstructing relative sea level in New Jersey is limited by an inability to differentiate between brackish sediments related to sea level and freshwater upland samples. To facilitate this distinction in a 4.4 m core, we used a multi‐proxy approach (δ¹³C values with presence or absence of agglutinated foraminifera) to recognize indicative meanings for four sample types. Sediment with δ¹³C values greater than ?18.9‰ was derived from a vegetated salt‐marsh and formed between MTL and MHHW. Sediment with δ¹³C values less than ?22.0‰ and containing agglutinated foraminifera formed in a brackish transitional zone between MHHW and highest astronomical tide (HAT). This is the narrowest elevational range of the four sample types and most precise sea‐level indicator. Sediment with δ¹³C values less than ?22.0‰ and lacking foraminifera can only constrain the upper bound of former sea level. Samples with intermediate values (?22.0‰ to ?18.9‰) formed between MTL and HAT. Using these indicative meanings and radiocarbon dates, we suggest that a transition from brackish to salt‐marsh δ¹³C values recorded in the core took approximately 350 years (from 1800 to 1450 cal. a BP). Copyright © 2011 John Wiley & Sons, Ltd.     

Seasonal and spatial water quality changes in the outflow plume of the Atchafalaya River,Louisiana, USA

The objective of this study was to examine the interaction between the Atchafalaya River and the Atchafalaya Delta estuarine complex. Measurements of suspended sediments, inorganic nutrients (NO₃ ^?, NH₄ ⁺, PO₄ ^3?), chlorophylla (chla), and-salinity were taken monthly from December 1996 to January 1998. These data were compiled by season, and the Atchafalaya River plume data were also analyzed using the Generalized Additive Model technique. There were significant decreases in NO₃ ^? concentrations during summer, fall, and winter as river water passed through the estuary, that were attributable to chemical and biological processes rather than dilution with ambient water. In some regions there were higher chla concentrations during summer and fall compared to winter and spring, when river discharge and the introduction of inorganic nutrients were highest, suggesting biological processes were active during this study. The presence of NH₄ ⁺, as a percentage of available dissolved inorganic nitrogen, increased with distance from the Atchafalaya River, indicative of remineralization processes and NO₃ ^? reduction. Mean PO₄ ^3? concentrations were often higher in the estuarine regions compared to the Atchafalaya River. During summer total suspended solid (TSS) concentrations increased with distance from the river mouth, suggesting a turbidity maximum. Highest chla concentrations were found in the bayous and shallow water bodies of the Terrebonne marshes, as were the lowest TSS concentrations. The low chla concentrations found in other areas of this study, despite high inorganic nutrient concentrations, suggest light limitation as the major control of phytoplankton growth. Salinity reached near seawater concentrations at the outer edge of the Atchafalaya River plume, but much lower salinities (<10 psu) were observed at all other regions. The Atchafalaya Delta estuarine complex buffers the impact of the Atchafalaya River on the Louisiana coastal shelf zone, with a 41% of 47% decrease in Atchafalaya River NO₃ ^? concentrations before reaching Gulf waters.     

Geophysical data integration,stochastic simulation and significance analysis of groundwater responses using ANOVA in the Chicot Aquifer system,Louisiana, USA

Data integration is challenging where there are different levels of support between primary and secondary data that need to be correlated in various ways. A geostatistical method is described, which integrates the hydraulic conductivity (K) measurements and electrical resistivity data to better estimate the K distribution in the Upper Chicot Aquifer of southwestern Louisiana, USA. The K measurements were obtained from pumping tests and represent the primary (hard) data. Borehole electrical resistivity data from electrical logs were regarded as the secondary (soft) data, and were used to infer K values through Archie’s law and the Kozeny-Carman equation. A pseudo cross-semivariogram was developed to cope with the resistivity data non-collocation. Uncertainties in the auto-semivariograms and pseudo cross-semivariogram were quantified. The groundwater flow model responses by the regionalized and coregionalized models of K were compared using analysis of variance (ANOVA). The results indicate that non-collocated secondary data may improve estimates of K and affect groundwater flow responses of practical interest, including specific capacity and drawdown.

---

Résumé  L’intégration de données entre en jeu lorsque plusieurs niveaux intermédiaires d’assistance sont nécessaires pour corréler données primaires et secondaires de diverses manières. Le présent article décrit une méthode géostatistique qui intègre les mesures de conductivité hydraulique (K) et les données de résistivité électrique, afin d’estimer plus efficacement la distribution de K dans l’Aquifère Supérieur de Chicot, au sud-ouest de la Louisiane (Etats-Unis). Les mesures de K sont issues des pompages d’essai et représentent les données primaires (“dures”). Les données des diagraphies de résistivité électrique ont été considérées comme des données secondaires (“molles”), à partir desquelles les valeurs de K ont été déduites, par la loi d’Archie et l’équation de Kozeny-Carman. Un pseudo semi-variogramme croisé a été développé afin de pallier à l’absence de colocalisation des données de résistivité. Les incertitudes sur les semi-variogrammes automatiques et sur les pseudo semi-variogrammes croisés ont été quantifiées. Les réponses du modèle d’écoulement des eaux souterraines aux modèles régionalisés et co-régionalisés de K ont été comparés, par les analyses de variance (ANOVA). Les résultats montrent que les données secondaires non-colocalisées peuvent améliorer les estimations de K, et affecter efficacement les réponses des écoulements souterrains, y compris les débits spécifiques et les rabattements.

---

Resumen  La integración de datos es un gran desafío cuando existen diferentes niveles de apoyo entre datos primarios y secundarios que es necesario correlacionar de varias maneras. Se describe un método geoestadístico el cual integra mediciones de conductividad hidráulica (K) y datos de resistividad eléctrica para tener una mejor estimación de la distribución de K en el Acuífero Chicot Superior del suroeste de Luisiana, Estados Unidos de América. Las mediciones de K se obtuvieron de pruebas de bombeo y representan los datos primarios (duros). Los datos de sondeos de resistividad eléctrica se consideraron como datos secundarios (suaves) y se usaron para inferir valores de K a través de la ley de Archie y la ecuación de Carman-Kozeny. Se desarrolló un pseudo semivariograma cruzado para enfrentar la falta de colocación de datos de resistividad. Se cuantificaron las incertidumbres en los auto-semivariogramas y en los semivariogramas cruzados. Las respuestas del modelo de flujo de agua subterránea por los modelos coregionalizados y regionalizados de K se compararon usando el análisis de varianza (ANOVA). Los resultados indican que los datos secundarios no colocados pueden mejorar los estimados de K y afectar las respuestas de flujo de agua subterránea de interés práctico, incluyendo capacidad específica y descenso.

     

Hierarchical dynamic stratigraphy in various Quaternary gravel deposits,Rhine glacier area (SW Germany): implications for hydrostratigraphy

The approach of dynamic stratigraphy aims to understand genetic processes that form stratigraphic units in a hierarchy of spatial and temporal scales. This approach was used to investigate Quaternary gravel deposits in terms of their sedimentology and in order to characterize the various sedimentary units in terms of their hydrogeological properties. Facies analysis within 62 gravel pits, laboratory permeability measurements of field samples and geophysical surveys (3-D georadar, 2-D seismic reflection) led to the detection and classification of sedimentary heterogeneity according to the following six scales whereby each scale can be translated into defined hydrostratigraphic units. (1) Particles and pores (micro scale) that reflect depositional and diagenetic fluid dynamics as well as source material behaviour (e.g. grain-size, roundness, lithological composition). This was found to be important for the hydrogeochemistry of groundwater in gravel aquifers (e.g. higher sorption capacity of carbon-rich limestone particles for organic pollutants). (2) Strata (meso scale) contain the recognition of sorting, fabric, texture and stratinomic features, which can give an indication of transport and depositional dynamics. Five major lithofacies groups, for example, were distinguished within fluvial gravel-bed deposits. Their variable hydraulic properties led to their subdivision into 12 hydrofacies types. They form the smallest mappable hydrostratigraphic units, which may result in either preferred pathways for fluid flow or flow barriers. (3) Depositional elements (macro scale) enable reconstruction of sedimentary/geomorphic elements and their dynamics within a depositional system (e.g. gravel-bed braided river systems are dominated by gravel sheet, gravel dunes and scour pool depositional elements). Hydrostratigraphically, the architecture of depositional elements influences the hydraulic connectivity and local permeability structure/distribution within an aquifer body. Five types of depositional elements in fluvial gravel-bed deposits were distinguished and their geometries/dimensions quantified. (4) Facies bodies (mega scale) composed of a stack of depositional elements and strata recording distinct environmental systems and their dynamics (e.g. a coarse-grained prograding delta system). Hydrostratigraphically, facies bodies represent major compartments of an aquifer. Six major types of meltwater-controlled facies bodies were identified in the study area. (5) Genetic sequences (mega scale) reflect the shifts of depositional environments caused by allocyclic changes (e.g. glacial advance recorded by a coarsening upward sequence) or autocyclic changes of landscape shaping events. These sequences may form separate hydrostratigraphic units or aquifer storeys. (6) Basin fill (giga scale) comprising the lateral and vertical stacking of facies bodies and genetic sequences controlled by either long-term glacier dynamics or short term flood events. The regional distribution of permeable gravel units and, for example, less permeable diamicts builds the larger scale hydrostratigraphy.     

Tropical Storm and Hurricane wind effects on water level, salinity, and sediment transport in the river-influenced Atchafalaya-Vermilion Bay System, Louisiana, USA

Changes in circulation, water level, salinity, suspended sediments, and sediment flux resulted from Tropical Storm Frances and Hurricane Georges in the Vermilion-Atchafalaya Bay region during September 1998. Tropical Storm Frances made landfall near Port Aransas, Texas, 400 km west of the study area, and yet the strong and long-lived southeasterly winds resulted in the highest water levels and salinity values of the year at one station in West Cote Blanche Bay. Water levels were abnormally high across this coastal bay system, although salinity impacts varied spatially. Over 24 h, salinity increased from 5 to 20 psu at Site 1 on the east side of West Cote Blanche Bay. Abnormally high salinities were recorded in Atchafalaya Bay but not at stations in Vermilion Bay. On September 28, 1998, Hurricane Georges made landfall near Biloxi, Mississippi, 240 km east of the study area. On the west side of the storm, wind stress was from the north and maximum winds locally reached 14 m s⁻¹. The wind forcing and physical responses of the bay system were analogous to those experienced during a winter cold-front passage. During the strong, north wind stress period, coastal water levels fell, salinity decreased, and sediment-laden bay water was transported onto the inner shelf. As the north wind stress subsided, a pulse of relatively saline water entered Vermilion Bay through Southwest Pass increasing salinity from 5 to 20 psu over a 24-h period. National Oceanic and Atmospheric Administration (NOAA)-14 reflectance imagery revealed the regional impacts of wind-wave resuspension and the bay-shelf exchange of waters. During both storm events, suspended solid concentrations increased by an order of magnitude from 75 to over 750 mg l⁻¹. The measurements demonstrated that even remote storm systems can have marked impacts on the physical processes that affect ecological processes in shallow coastal bay systems.     

Variation in the composition and partitioning of adsorbed cations at a brine-contaminated crude oil production facility in southeastern Louisiana,USA

The first geological materials impacted by oil field wastes released into near-surface environments in southern Louisiana, USA, are typically clays and silts. Clay minerals within these siliciclastic sediments have the potential for altering the composition of produced water wastes through cation exchange. The general relations between the composition of adsorbed cations and interstitial water salinity in brine-contaminated samples from a site in southeastern Louisiana are consistent with previous studies of multicomponent exchange in groundwater systems of varying salinity. The divalent cations Ca and Mg dominate as adsorbed cations at low salinities (<1200 mg/L), but Na is dominant at moderate to high salinities (up to 53,000 mg/L). The change in the proportions of adsorbed cations is a non-linear function of salinity, and the transition from Ca-dominated adsorption to Na-dominated adsorption occurs over a narrow range of salinities.     

Integrated geostatistics and GIS techniques for assessing groundwater contamination in Al Arish area, Sinai, Egypt

Sustainable development in El Arish area of North Sinai, Egypt, is retarded by serious environmental problems, where the land-use and land cover of the region is changing over present time. The impact of human activities in the study area is accompanied by the destruction and over-exploitation of the environment. This study applies multivariate statistics (factor and cluster analyses) and GIS techniques to identify both anthropogenic and natural processes affecting the groundwater quality in the Quaternary sands aquifer. The aim of this study was to investigate the impacts on groundwater resources, the potential pollution sources, and to identify the main anthropogenic inputs of both nutrients and trace metal. Since the depth to the water table is shallow especially in the northern part (<4?m), and the aquifer was exposed on the ground surface, it has poor buffering capacity and the pollution risk is very high. Groundwater chemistry in this coastal region has complex contaminant sources, where intensive farming activities and untreated wastes put stress on groundwater quality. Several areal distribution maps were constructed for correlating water quality with possible contributing factors such as location, land-use, and aquifer depth. These maps identified both anthropogenic and natural processes affecting groundwater quality of the studied aquifer. Cluster analysis was used to classify water chemistry and determine the hydrochemical groups, Q-mode dendrogram is interpreted and there are three main clusters. Factor analyses identify the potential contamination sources affecting groundwater hydrochemistry such as: nitrate, sulfate, phosphate and potassium fertilizers, pesticides, sewage pond wastes, and salinization due to circulation of dissolved salts in the irrigation water itself.     

Sources of Nutrients and Fecal Indicator Bacteria to Nearshore Waters on the North Shore of Kaua`i (Hawai`i,USA)

Water quality monitoring in Hanalei Bay, Kaua`i (Hawai`i, USA) has documented intermittent high concentrations of nutrients (nitrate, phosphate, silica, and ammonium) and fecal indicator bacteria (FIB, i.e., enterococci and Escherichia coli) in nearshore waters and spurred concern that contaminated groundwater might be discharging into the bay. The present study sought to identify and track sources of nutrients and FIB to four beaches in Hanalei Bay and one beach outside the bay, together representing a wide range of land uses. ²²³Ra and ²²⁴Ra activity, salinity, nutrient and FIB concentrations were measured in samples from the coastal aquifer, the nearshore ocean, springs, the Hanalei River, and smaller streams. In addition, FIB concentrations in beach sands were measured at each site, and the enterococcal surface protein (esp) gene assay was used to investigate whether the observed FIB originated from a human source. Nutrient concentrations in groundwater were significantly higher than in nearshore water, inversely correlated to salinity, and highly site specific, indicating local controls on groundwater quality. Fluxes of groundwater into Hanalei Bay were calculated using a mass-balance approach and represented at least 2–10% of river discharges. However, submarine groundwater discharge (SGD) may provide 2.7 times as much nitrate + nitrite to Hanalei Bay as does the Hanalei River. It may also provide significant fluxes of phosphate and ammonium, comprising 15% and 20% of Hanalei River inputs, respectively. SGD-derived silica inputs to the bay comprised less than 3% of Hanalei River inputs. FIB concentrations in groundwater were typically lower than those in nearshore water, suggesting that significant FIB inputs from SGD are unlikely. Positive esp gene assays suggested that some enterococci in environmental samples were of human fecal origin. Identifying how nutrients and FIB enter nearshore waters will help environmental managers address pressing water quality issues, including exceedances of the state Enterococcus water quality standard and nutrient loading to coral reefs.     

Ultraphyllonite origin for slate,mid-Atlantic Piedmont,USA

Structural, petrographic, and geochemical criteria are applied in combination to determine the petrologic origins of slaty rocks in polydeformed metasedimentary terranes. Multiple lines of evidence suggest that the Peach Bottom slate (PBs), mined over several centuries, is an ultraphyllonite. It occurs within an Alleghanian dextral retraining bend within the Appalachian Piedmont, as a segment of the Pleasant Grove-Huntington Valley shear zone system (>200 km long). The PBs differs markedly from other slates, in that, it contains abundant, micro-scale remnant porphyroclasts of higher-grade minerals. The PBs has major, trace and rare earth element geochemistry and ¹⁴⁷Sm/¹⁴³Nd ratios that are similar to the surrounding pelitic rocks. When compared to average slates, the geochemistry of the PBs exhibits depletion in alkali earth and alkali elements.