We will make available to users of the Marine Heat Flow Capability a Matlab program called SlugHeat, developed at UCSC based on the theory and programs described by Villinger and Davis . Slug heat processes individual thermistor penetration records to derive in-situ temperatures, thermal conductivity (when data are collected for this purpose using the multipenetration system), and heat flow. The program can be run interactively or in batch mode, and includes options for discharding individual thermistor records, in case of collection of spurious data, and Monte Carlo analysis for both variations in thermal conductivity and the depths of boundaries between heterogeneous layers. Source code, example input and output files, and instructions will be distributed to researchers using the Marine Heat Flow Capability.
SlugSed is a one-dimensional numerical model of fluid and heat transport coded in Matlab. In its most basic form, the model calculates the temperature distribution within a medium of finite thickness and known thermal conductivity and heat capacity, given appropriate initial and boundary conditions. The initial conditions consist of depth, temperature, rate(s) of radiogenic heat production or loss, and thermal conductivity for discrete points within the medium. The boundary conditions must be specified as temperature at the top of the medium and either temperature or heat flux at the bottom, and both can be time-variant. SlugSed can accommodate multiple distinct basement layers, each having independent property sets. In the case of sedimentation, the depth-variant physical and thermal properties within the accumulating sediment layer are derived from a user-specified porosity versus depth function. The thermal properties of basement layers (when in use) can be specified in bulk (i.e., applied to all points within the layer) or individually. Seepage of fluids within the sediment layer can be driven by compaction and consolidation, fluid overpressures in basement, or imposed explicitly.