Provenance of the Southeastern South China Block in the Late Triassic and Initiation of Paleo-Pacific Subduction:Evidence from Detrital Zircon U-Pb Geochronology

During the Late Paleozoic-Early Mesozoic Era,the sediment transport system and tectonic regime in the southeastern margin of the South China Block(SESCB)all changed,significantly affected by the Paleo-Pacific subduction.However,controversy exists about the Paleo-Pacific subduction's initiation time.This study uses detrital zircon U-Pb ages to discuss the Late Triassic source-tosink system in the SESCB.It provides some references for the Paleo-Pacific subduction process based on crucial age information and zircons'trace elements.The paleogeography and similarity of detrital zircon age distribution indicate that three sinks were found in the SESCB during the Late Triassic:1.the Yangchun-Kaiping-Gaoming area,comprising major age ranges of 260-220,460-400,and 1200-800 Ma,which might be sourced from the Yunkai terrane;2.the Jiexi-Kanshi-Nanjing area,characterized by the significant age component of 2000-1800 Ma,which corresponded to the Wuyi terrane;3.the Xinan area,consisting of significant age groups of 290-250 and 380-320 Ma,which might be sourced from the magmatic rocks formed by the Huinan Movement and Paleo-Pacific subduction.Note that 290-250 Ma zircons were widely distributed in the Upper Triassic strata,and their trace elements suggested the existence of a magmatic arc near the SESCB during the 290-250 Ma.Thus,we propose that the Paleo-Pacific subduction might have begun in the Early Permian.

Quantifying the age peaks,age ranges and weights of detrital ages based on the EM algorithm

Detrital geochronology fundamentally involves the quantification of major age ranges and their weights winthin an age distribution.This study presents a streamlined approach,modeling the age distribution of detrital zircons using a normal mixture model,and employs the Expectation-Maximization(EM)algorithm for precise estimations.A method is introduced to automatically select appropriate initial mean values for EM algorithm,enhancing its efficacy in detrital geochronology.This process entails multiple trials with varying numbers of age components leading to diverse k-component models.The model with the lowest Bayesian Information Criterion(BIC)is identified as the most suitable.For accurate component number and weight determination,a substantial sample size(n>200)is advisable. Our findings based on both synthetic and empirical datasets confirm that the normal mixture model,refined by the EM algorithm,reliably identifies key age parameters with minimal error.As a kind of probability density estimator,the normal mixture model offers a novel visualization tool for detrital data and an alternative foundation for KDE in calculating existing similarity metrics.Another focus of this study is the critical examination of quantitative metrics for comparing detrital zircon age patterns.Through a case study,this study demonstrates that metrics based on empirical cumulative probability distribution(such as K-S and Kuiper statistics)may lead to erroneous conclusions.The employment of the Kullback-Leibler(KL)divergence,a metric grounded in probability density estimation,is proposed.Reference critical values,simulated via the Monte Carlo method,provide more objective benchmarks for these quantitative metrics. All methodologies discussed are encapsulated in a series of MATLAB scripts,available as open-source code and a standalone application,facilitating wider adoption and application in the field.