A Preliminary Heat Flow Model for Cooling a Batholith near Ica, Peru

This research models the cooling of a suite in the Linga Super-unit located at the north end of the Arequipa segment in the Cretaceous Peruvian Coastal Batholith. The monzogabbro to granite Linga suite is approximately 50 km long and 15 km wide, with an estimated vertical extent of about 5 km originally intruded to a depth of 3 km. The emplacement was in andesitic volcanics on the west and the Pampahuasi diorite Super-unit on the east and has incorporated earlier gabbroic bodies. The Linga suite is thought to be the result of a sequence of three pulses: an elongate unit to the west then two elliptical units to the northeast and southeast. The data for modeling comes from field observations on internal and external contacts, some well-defined magma chamber walls and roof, pendant and stoped blocks, magma chamber zoning, the nature and distribution of enclaves and xenoliths, magmatic fabric, evidences of magma mingling, rock porosity, mineralogical associations in metamorphic aureoles, extensive mineralization and brecciated conduits, and the types of hydrothermal alteration varying with distance from contacts. More than forty hand samples, thin sections, and geochemical analyses were used to estimate water content, magma and country rock temperature, liquid density, and viscosity. Further data will come from: zircon U-Pb ages for country rock and magma batch timeframes, fluid inclusions for magma pressure and temperature, and δ18O data for source of hydrothermal fluids. Simple heat conduction calculations using MATLAB and HEAT 3D for a single tabular intrusion estimated a cooling time to solidus of about 300 k.y. More complex modeling includes magma convection and multiple intrusions. Extensive veining and pervasive alteration suggested the use of HYDROTHERM to model possible additional heat flow effects from hydrothermal fluids. Extensive propylitic and significant potassic alteration were observed and, with TerraSpec infrared spectroscopy to identify mineralogical content of altered outcrops, fluid flux and temperatures were estimated. Rough calculations for cooling times are made using the modeling of Hayba and Ingebritsen (1997).