The 1500 m thick sequence of Huronian continental volcanics at Thessalon, Ontario is subdivided into 4 volcanic cycles, each of which includes abundant early mafic end-members, central intermediate flows, and late rhyolite units. Major and trace element concentrations are dominated by extensive gabbroic fractionation trends that ultimately produced two types of felsic flows: (1) rhyolites with high light rare earth element (LREE) and relatively low large-ion lithophile element (LILE) concentrations (high-LREE, low-LILE rhyolites), and (2) following late separation of REE-rich accessory phases, rhyolites depleted in LREE (low-LREE, high-LILE rhyolites). Mafic end-members of individual volcanic cycle are progressively less siliceous and less enriched in LILE and LREE with height in the stratigraphic section. Ti/Zr ratios gradually rise from 35 in early mafic flows to stabilize at about 85 in late units, while average SiO2 contents decrease from 56 to about 50% and Mg# rises from about 48 to 52. ɛ-Nd values are consistently negative, indicating variable degrees of pre-fractionation crustal contamination of the end-member magmas during their uprise through the crust. Mixing models are consistent with up to 50% contamination by crustal material of tonalitic hornblende-gneiss composition. A progressive increase in ɛ-Nd, from about-5.0 to-0.5 upward in the volcanic succession, reflects a decreasing degree of crustal contamination due to development of insulating layers along margins of the feeder system. Detailed stratigraphic variations suggest that successive magmas batches were intercepted by a progressively fractionating, periodically replenished magma source, giving rise to open-system magmatism. Despite the prevalence of crustal assimilation in the Huronian lavas, (La/Sr)N ratios are too low in least contaminated end-members to be explained by contamination of tholeiitic magmas. The late basalts resemble instead modern island are basalts, and it is suggested that the subcontinental mantle source was enriched by subduction-related processes during crustal formation. Within individual volcanic cycles gabbroic fractionation trends systematically deviate from calculated factors toward compositions characteristic of hornblende-gneiss. Such relations suggest that further crustal contamination of the magmas occurred simultaneous with crystal fractionation. probably within undulating sills at upper crustal levels. Quantitative analysis suggests assimilation/fractional crystallization (A/FC) ratios of about 0.45. As a result of extensive two-stage contamination, rhyolites from the initial volcanic cycle incorporate a total of over 60% of crust.