The origin of the Messinian Hoya de la Sima (HS) gypsum (Betic foreland) is constrained using Sr-87/Sr-86, delta S-34, Sr concentration, and petrographic data. The Lower and Middle HS units consist of subaqueous vertically-aligned and stromatolitic selenites, the latter containing unusual microbial depositional textures. The Upper Unit consists of very-shallow-water bioturbated lenticular gypsum with Paracamelus ichnites. Sr-87/Sr-86 and delta S-34 indicate precipitation from predominantly marine waters, with upward increasing continental influence. Mixing models between Messinian seawater and continental water that dissolved Triassic evaporites show that the percentages of seawater required to explain the measured Sr-87/Sr-86 are analogous to the percentages obtained using delta S-34, supporting precipitation from such mixtures. Sr-87/Sr-86 and delta S-34 of Lower HS selenites resemble those of the Primary Lower Gypsum (PLG) of the Messinian Salinity Crisis (MSC), in both cases indicating precipitation from seawater-continental water mixtures in which most Sr and SO4 were supplied by Messinian seawater. In the Lower HS selenites, Sr concentrations indicate contributing continental waters with Sr/Ca ratios similar to seawater. However, Sr concentrations of PLC selenites from other Betic basins (Bajo Segura, Sorbas indicate parent waters with Sr/Ca ratios lower than seawater. If the Sr contents of the betic PLG selenites are representative, it is unlikely that the Lower HS selenites represent the PLC. However, we cannot completely discard that option since different LPG subbasins could have had variable Sr/Ca. The HS gypsums formed coevally to diapirism of Triassic evaporites, in a restricted lagoonal basin developed during or slightly after a phase of strike-slip faulting in the Betic Cordillera. More general implications of this work are that Sr concentrations, combined with Sr-87/Sr-86 and delta S-34 data, provide key constraints on the origin of parent brines, and using Sr concentrations as paleosalinity proxy in gypsum precipitated from seawater-continental water mixtures is inconclusive because they largely depend on the initial Sr/Ca ratio of the brine. This study also illustrates that the electron microprobe is ideal to obtain reliable, spatially-resolved, Sr concentration data in selenitic gypsum. These data, in combination with fluorescence petrography, provide valuable criteria to discriminate primary depositional from diagenetic gypsum, a critical step in the interpretation of geochemical data. (C) 2015 Elsevier B.V. All rights reserved.