The cooling histories of individual meteorites can be empirically reconstructed by using ages from different radioisotopic chronometers with distinct closure temperatures. For a group of meteorites derived from a single parent body such data permit the reconstruction of the cooling history and properties of that body. Particularly suited are H chondrites because precise radiometric ages over a wide range of closure temperatures are available. A thermal evolution model for the H chondrite parent body is constructed by using all H chondrites for which at least three different radiometric ages are available. Several key parameters determining the thermal evolution of the H chondrite parent body and the unknown burial depths of the H chondrites are varied until an optimal fit is obtained. The fit is performed by an 'evolution algorithm'. Empirical data for eight samples are used for which radiometric ages are available for at least three different closure temperatures. A set of parameters for the H chondrite parent body is found that yields excellent agreement (within error bounds) between the thermal evolution model and empirical data of six of the examined eight chondrites. The new thermal model constrains the radius and formation time of the H chondrite parent body (possibly (6) Hebe), the initial burial depths of the individual H chondrites, the average surface temperature of the body, the average initial porosity of the material the body accreted from, and the initial 60Fe content of the H chondrite parent body.

The major aim of this study is to assess the effects of sintering of initially porous material on the thermal evolution of planetesimals, and to constrain the values of basic parameters that determined the structure and evolution of the H chondrite parent body. A new code is presented for modeling the thermal evolution of ordinary chondrite parent bodies that initially are highly porous and undergo sintering by hot pressing as they are heated by decay of radioactive nuclei. The pressure and temperature stratification in the interior of the bodies is calculated by solving the equations of hydrostatic equilibrium and energy transport. The decrease of porosity of the granular material by hot pressing due to self-gravity is followed by solving a set of equations for the sintering of powder materials. For the heat conductivity of granular material we combine recently measured data for highly porous powder materials, relevant for the surface layers of planetesimals, with data for heat conductivity of chondrite material, relevant for the strongly sintered material in deeper layers. To demonstrate the capability of our new model, the thermal evolution of the H chondrite parent body was reconstructed. The model starts with a porous body that is later compacted first by 'cold pressing' at low temperatures and then by 'hot pressing' for temperatures above \approx 700 K, i.e., the threshold temperature for sintering of silicates. The thermal model was fitted to the well constrained cooling histories of the two H chondrites Kernouve (H6) and Richardton (H5).