AbstractThis study integrates geochemical modeling, spatial analysis and several statistical methods including principal component analysis, multivariate regression and cluster analysis to investigate hydrogeologic controls of arsenic and uranium contamination within groundwater of the Arikaree aquifer on the Pine Ridge Reservation (PRR). Located in southwestern South Dakota, the PRR is largely rural and many people rely on domestic supply wells completed in the Arikaree aquifer as their primary drinking water source. Locally, the White River Group, which unconformably underlies the Arikaree Group, is enriched in arsenic and uranium related to volcanic ash deposits and acts as a geogenic metal source. Geochemical data from over 250 groundwater samples were obtained through collaboration with the Oglala Sioux Tribe. Cluster spatial statistics analyses delineated four regions of statistically significant variations in groundwater chemistry that represent upgradient, intermediate, and downgradient portions of the Arikaree aquifer. Groundwater evolves as it flows through the Arikaree aquifer with increasing alkalinity, sodium, and pH along flow paths. These chemical changes are likely due to dissolution of carbonate minerals and volcanic ash. Thermodynamic calculations suggest increasing supersaturation of the groundwater with respect to calcite; thus, volcanic ash dissolution may be an important secondary source of alkalinity. Elevated alkalinity and pH levels were found to be the driving factors of arsenic and uranium mobility, and downgradient sections of the aquifer in the northern portions of the PRR are most likely to be impacted by metal(loid) contamination, with 73% of the wells in this grouping failing a USEPA maximum contaminant level for arsenic, uranium, and/or gross alpha.