Partial Melting and Crustal Deformation during the Early Paleozoic Wuyi–Yunkai Orogeny: Insights from Zircon U-Pb Geochronology and Structural Analysis of the Fuhuling Migmatites in the Yunkai Region, South China

: Migmatites record crucial information about the rheology and tectonothermal evolutionof the deep crust during orogenesis. In the Wuyi–Yunkai orogen in South China, migmatites at Fuhuling record Early Paleozoic high temperatures and associated partial melting. However, the absolute timing and implications for the rheology of the deep crust during orogenesis are poorly constrained. In this contribution, we used spatial analysis of migmatitic leucosomes, structural analysis, and U-Pb geochronology of zircon to elucidate the absolute timing of crustal partial melting, the degree of partial melting, and the role of partial melting on the rheology of the crust during the Wuyi–Yunkai orogeny. Partial melting of the Fuhuling migmatites occurred at c. 440 Ma during Early Paleozoic Wuyi–Yunkai orogenesis. Subsequent lower temperature metamorphism associated with Indosinian movement that caused minor zircon recrystallization was temporally associated with the crystallization of nearby biotite monzogranites, but it did not inﬂuence the morphology of the Fuhuling migmatites. The migmatites preserve a morphological transition from metatexite to diatexite with an increasing proportion of leucosome. This transition preserves di ﬀ erent structural characteristics that represent the response of the solid framework and melt network to variable melt fractions during partial melting. The large proportion of in situ or in source leucosome in the Fuhuling migmatites suggests that it was a melt-rich crustal horizon during orogenesis, and that a substantial proportion of anatectic melt was retained in the deep crust. The rheological transition documented in the Fuhuling migmatites was caused by changes in the melt fraction, and it is an analogue for the rheological transition characteristics of melt-rich crustal horizons in the Yunkai region during Early Paleozoic Wuyi–Yunkai orogenesis and subsequent orogenic collapse.

Abstract: Migmatites record crucial information about the rheology and tectonothermal evolutionof the deep crust during orogenesis. In the Wuyi-Yunkai orogen in South China, migmatites at Fuhuling record Early Paleozoic high temperatures and associated partial melting. However, the absolute timing and implications for the rheology of the deep crust during orogenesis are poorly constrained. In this contribution, we used spatial analysis of migmatitic leucosomes, structural analysis, and U-Pb geochronology of zircon to elucidate the absolute timing of crustal partial melting, the degree of partial melting, and the role of partial melting on the rheology of the crust during the Wuyi-Yunkai orogeny. Partial melting of the Fuhuling migmatites occurred at c. 440 Ma during Early Paleozoic Wuyi-Yunkai orogenesis. Subsequent lower temperature metamorphism associated with Indosinian movement that caused minor zircon recrystallization was temporally associated with the crystallization of nearby biotite monzogranites, but it did not influence the morphology of the Fuhuling migmatites. The migmatites preserve a morphological transition from metatexite to diatexite with an increasing proportion of leucosome. This transition preserves different structural characteristics that represent the response of the solid framework and melt network to variable melt fractions during partial melting. The large proportion of in situ or in source leucosome in the Fuhuling migmatites suggests that it was a melt-rich crustal horizon during orogenesis, and that a substantial proportion of anatectic melt was retained in the deep crust. The rheological transition documented in the Fuhuling migmatites was caused by changes in the melt fraction, and it is an analogue for the rheological transition characteristics of melt-rich crustal horizons in the Yunkai region during Early Paleozoic Wuyi-Yunkai orogenesis and subsequent orogenic collapse.

Introduction
Geological and geophysical studies in modern and ancient orogens demonstrate that the middle to lower continental crust can be partially molten during orogenesis [1][2][3][4][5]. Partial melting drastically weakens the crust and plays an essential role in the rheological evolution of orogenic systems [5][6][7]. Migmatites are rocks that have partially melted at high temperatures in the crust [5] and have two broad morphologies, metatexite and diatexite, based on increasing melt fraction and the continuity of the solid framework [8][9][10][11]. Metatexite generally preserves relict (i.e., pre-anatectic) structures and form as the melt is extracted from intergranular space and the solid framework fuses together [5,12,13]. Diatexite rarely preserves pre-partial melting structures, and it forms when the melt proportion increases to the extent that the solid framework is no longer connected, resulting in strengths close to that of magmas [14][15][16]. These two types of migmatites are structurally and compositionally heterogeneous, and they preserve deformation features that arise from rheological transitions related to melt-dominated to melt-absent behavior [5,11,17,18]. Therefore, migmatites act as an analogue for the rheological behavior of partially molten crust during orogenesis [8,10,19,20]. High-grade gneiss terranes typically undergo multiple deformation and metamorphic episodes, and, consequently, migmatites at the outcrop scale generally have complex geometries and structures that overprint the topology of the melt network that existed during crustal parting melting [10,19,20]. Migmatites can also be overprinted by brittle-ductile shear zones, which further complicates the interpretation of migmatitic structures [19,21]. Considering the complex geological history of most migmatites, detailed field mapping and structural analyses are needed to understand the anatectic melt network of migmatites [11,12,18,[22][23][24][25][26] and are crucial for providing insights into the mechanics of melt production, melt segregation, melt migration, and melt redistribution during crustal partial melting accompanying orogenesis [8,10,11,27,28]. However, there is still a poor understanding of the relative roles of deformation versus melt fraction during the transition from metatexite to diatexite in migmatite terranes, but this information is crucial for linking structures in migmatites to the broader evolution of orogenic belts.
The Wuyi-Yunkai orogen is an intracontinental orogen (traditionally called the "Caledonian orogen" or "Kwangsian orogen" in Chinese literature [29][30][31]). The orogen covers the southeastern half of the South China block (Figure 1b), extends for~2000 km in a northeasterly direction, and contains widespread migmatites formed by the partial melting of continental crust during the Early Paleozoic [29][30][31]. The Yunkai region is a crucial exposure of the deep crust that records Early Paleozoic crustal partial melting in the Wuyi-Yunkai orogen [29][30][31][32]. However, whether partial melting in the Yunkai region occurred in a single event or in multiple pulses is still debated [32][33][34][35][36]. Previous studies on migmatites in the Yunkai region focused on the timing and nature of partial melting on the evolution of the Wuyi-Yunkai orogen [30][31][32]37]. However, little work has been done on the relationships between crustal partial melting, rheological transitions, and structural transitions from metatexites to diatexites, which are all crucial for understanding the nature of partial melting that operated in the Wuyi-Yunkai orogeny and the broader significance of metatexite to diatexite transitions in migmatite terranes.
The Fuhuling migmatites in the Yunkai region are situated on the southwestern margin of the Wuyi-Yunkai orogen (Figure 2a) [38]. A complete metatexite-diatexite profile with distinct lithological changes and variations in migmatite morphology is exposed on the southeast coast. Therefore, the Fuhuling migmatites provide a rare opportunity to investigate the relationship between crustal partial melting, rheological transitions, and migmatite structures for the Yunkai region. Here, we present a detailed structural analysis of the migmatites as well as new zircon U-Pb geochronology and trace elements analysis from metatexite, diatexite, and biotite monzogranite in the Fuhuling region to determine the timing of crustal partial melting and the implications for the tectonic evolution of the Yunkai region. In addition, detailed spatial and structural analysis of the Fuhuling migmatites is used to constrain the relationship of deformation versus melt fraction during the metatexite-diatexite transitions in migmatite terranes. Finally, using the Fuhuling migmatites as an analogue for the deep crust coupled with the results of U-Pb zircon geochronology, we discuss the rheological transition that occurred during crustal partial melting in the Yunkai region during Early Paleozoic Wuyi-Yunkai orogeny.

Regional Geology
The South China block (SCB) consists of the Yangtze Block in the northwest and the Cathaysia Block in the southeast. The present boundary between these two blocks in eastern South China is the approximately NE-trending Jiang-Shao Fault, but the extension of this boundary in the southwestern SCB is uncertain because of its poor exposure and younger tectonic overprinting [31,39]. The Wuyi-Yunkai orogen covers the southeastern half of the SCB, stretches for~2000 km in a northeasterly direction, and is unconformable with Devonian strata [29,31]. Metamorphic rocks in the orogen are Paleozoic and older and are intruded by granites with crystallization ages mostly between 464 and 400 Ma ( [40] and reference therein). The Yunkai region is a part of the Wuyi-Yunkai orogen and is located in the western Cathaysia block (Figure 1a). The metamorphosed basement in the Yunkai region is mainly composed of the Gaozhou Complex and Yunkai Group, which are unconformably overlain by weakly metamorphosed to unmetamorphosed Cambrian to Devonian strata ( Figure 2a) [31,37,41,42]. The Gaozhou Complex is composed of amphibolite facies (locally granulite facies) paragneiss, schist, quartzite, marble, and migmatites, whereas the Yunkai Group is composed of lower-grade schist, slate, and phyllite as well as less common paragneiss, amphibolite, and marble [34,37]. Late Paleozoic to Cenozoic sedimentary rocks and Paleozoic to Early Mesozoic granites also occur in the Yunkai region ( Figure 2a) [41][42][43].
Fuhuling is located in Yangxi County of Yangjiang City in Guangdong Province on the edge of the South China Sea (Figure 2a). It is tectonically situated in the southeastern margin of the Yunkai region ( Figure 2a) [38]. The topography of Fuhuling is dominated by a peninsula-like hill surrounded by seas and covered with vegetation; it has a maximum elevation of 138 m above sea level. The dominant rock type in Fuhuling is migmatite, which includes metatexite and diatexite varieties [38]. Here, metatexite defines a migmatite where gneissic banding is evident, and diatexite contains no continuous gneissic banding (e.g., [14]). Neoproterozoic-Early Paleozoic metasedimentary rocks mainly consist of metagreywacke and biotite-quartz schist [38,44,45] that outcrop on the top and the hillside of Fuhuling (Figure 2b). A small area of biotite monzogranite is exposed on the southwestern beach of Fuhuling ( Figure 2b) [38].  [30,31,41,47]).

Outcrop-Scale Structure
Metatexite and diatexite are only exposed on cliffs and beaches along the coast of Fuhuling [38]. There is a complete vertical metatexite-diatexite profile (<15 m in height) exposed in the southeastern portion of the research area (Figures 2 and 3). Based on compositional and structural variations, we divided this profile into four layers from top to bottom ( Figure 3). Below, we use the terms leucosome to refer to the light-colored portion of a migmatite and melanosome for the dark-colored portion [14].
Metatexite and diatexite are only exposed on cliffs and beaches along the coast of Fuhuling [38]. There is a complete vertical metatexite-diatexite profile (<15 m in height) exposed in the southeastern portion of the research area (Figures 2 and 3). Based on compositional and structural variations, we divided this profile into four layers from top to bottom ( Figure 3). Below, we use the terms leucosome to refer to the light-colored portion of a migmatite and melanosome for the dark-colored portion [14].
Layer II is also a stromatic metatexite but with a larger proportion (~10 area %) of leucosome compared to that of Layer I. Leucosomes in Layer II generally occur as stroma <2 cm thick, and some are slightly bent (Figure 4b). The orientations of 20 stromatic leucosomes have poles between 139° and 225° with dips of 20-79° ( Figure 3).  Layer I is a metatexite with~5 area % (inferred to be equivalent to vol. %) of leucosome, which occurs as elongate lenses (<1 cm in diameter) and intermittent stroma (<2 mm in thickness; Figure 4a). These are typical stromatic structures of metatexite [25]. The orientations of eight stromatic leucosomes have poles between 150 • to 252 • with dips of 29 • to 54 • (Figure 3).
Layer II is also a stromatic metatexite but with a larger proportion (~10 area %) of leucosome compared to that of Layer I. Leucosomes in Layer II generally occur as stroma <2 cm thick, and some are slightly bent (Figure 4b). The orientations of 20 stromatic leucosomes have poles between 139 • and 225 • with dips of 20-79 • (Figure 3).
Layer III is a metatexite with folded leucosomes that change in wavelength and amplitude over short distances. Folds vary in morphology and include asymmetric, tight, ptygmatic, parasitic, rootless, and superimposed folds (Figure 4c,d). Leucosomes and melanosomes are thicker in fold hinges than in fold limbs. In the lower part of layer III, metatexite is in transitional contact with diatexite and has fewer stromatic leucosomes (Figure 4g). Stromatic leucosomes have variable orientations; eighty-eight orientation measurements yield variable poles that are mostly 104 • to 241 • (Figure 3).   Layer IV is a diatexite, and leucosomes and melanosomes in this layer mainly form schollen structures as well as schlieren and nebulitic structures (Figure 4d,e). In the diatexite, there are some lenticular, hook-like, and trailing paleosomes (i.e., unmelted rock) or restites (i.e., melt-depleted residues) wrapped by leucosome (Figure 4f). The diatexite exhibits schollen structures as well as schlieren and nebulitic structures (Figure 4d,e) [18,48]. The orientation of fifty-two leucosomes have poles between 61 • to 354 • with dips of 24 • to 85 • , and there is no significant directionality ( Figure 3).

Sample Selection and Analytical Methods
Six samples from Fuhuling were selected for zircon U-Pb geochronology and trace element analysis, including two metatexite samples (sample No.: SPM and FHM), two diatexite samples (sample No.: SPMG and FHMG), and two biotite monzogranite samples (sample No.: DG2 and DG3). SPM and SPMG were collected from layer III and layer IV of the metatexite-diatexite profile, respectively. FHM, FHMG, DG2, and DG3 were all collected from the southwestern beach of Fuhuling ( Figure 2).
Zircon separation and cathodoluminescence (CL) imaging were both completed by the Honesty Geological Service Company in Langfang, Hebei Province. The CL images were acquired using a JSM6510 (Japan) scanning electron microscope fitted with a GatanMini CL detection system. LA-ICPMS zircon U-Pb isotopic dating and trace element analyses were synchronously carried out using an Agilent 7500a ICP-MS coupled with a GeoLas2005 laser ablation system at the State Key Laboratory of Geological Processes and Mineral Resources, China University of Geosciences (Wuhan). Laser ablation was conducted with a beam diameter of 32 µm for all analyses. Each analysis incorporated an approximately 20 s background acquisition (i.e., gas blank) followed by 50 s of ablation. The software ICPMSDataCal was used to perform off-line selection and integration of background and analyte signals, time-drift corrections, quantitative calibration for trace element analyses, and U-Pb geochronology [50,51]. Specific analytical conditions are detailed in   [51]. Common Pb corrections were undertaken using ComPbCorr#3_17 of Andersen (2002) [52]. Count rates for 29 Si, 204 Pb, 206 Pb, 207 Pb, 208 Pb, 232 Th, and 238 U were collected for age determination. Trace element compositions of zircon were calibrated against NIST SRM 610 and using Si as an internal standard; analytical accuracies are generally better than 5-10% for trace elements [50]. Zircon weighted mean ages and U-Pb concordia diagrams were calculated with the software ISOPLOT (version 3.0) [53]. Uncertainties on individual LA-ICPMS analyses were reported at the 1σ level, and uncertainties on weighted mean ages were quoted at a level of confidence of 95% (2σ). Given the low abundance of 207 Pb in young zircon, we used 206 Pb/ 238 U ages to calculate weighted mean ages for zircon <1000 Ma, and 207 Pb/ 206 Pb ages were preferentially adopted for zircons for >1000 Ma. We excluded zircon U-Pb analyses with discordance values >10% for inherited zircon cores and >5% for metamorphic zircon and magmatic zircon from weighted mean ages, histograms, density probability plots, and U-Pb concordia diagrams. The complete data set, including the excluded discordant analyses, is reported in Supplementary Table S1.

Zircon U-Pb Geochronology
Cathodoluminescence images of representative zircons together with spot ages are shown in Figure 6. U-Pb data for these samples are summarized on Wetherill concordia diagrams in Figure 7.
Given the low abundance of 207 Pb in young zircon, we used 206 Pb/ 238 U ages to calculate weighted mean ages for zircon <1000 Ma, and 207 Pb/ 206 Pb ages were preferentially adopted for zircons for >1000 Ma. We excluded zircon U-Pb analyses with discordance values >10% for inherited zircon cores and >5% for metamorphic zircon and magmatic zircon from weighted mean ages, histograms, density probability plots, and U-Pb concordia diagrams. The complete data set, including the excluded discordant analyses, is reported in supplementary Table S1.

Zircon U-Pb Geochronology
Cathodoluminescence images of representative zircons together with spot ages are shown in Figure 6. U-Pb data for these samples are summarized on Wetherill concordia diagrams in Figure 7.  oscillatory-zoned cores and rims with weak oscillatory zoning ( Figure 6). Nineteen spots were analyzed from 19 zircon grains. Of three analyses of cores, one yielded a 207 Pb/ 206 Pb age of 1128 ± 69 Ma, and the other two analyses gave 206 Pb/ 238 U ages of 437 ± 4 and 443 ± 6 Ma (Table S1). The other 16 analyses on oscillatory zones yielded 206 Pb/ 238 U ages ranging from 234 ± 3 to 255 ± 4 Ma. After excluding ten analyses with discordance values >5%, the remaining six concordant analyses yielded a weighted mean 206 Pb/ 238 U age of 241.0 ± 4.3 (MSWD = 1.3) (Figure 7f), which is interpreted as the crystallization age of the biotite monzogranite.

Diatexite
Zircon from diatexite SPMG was mainly subhedral to euhedral, long prismatic, and 70 to 200 µm long. CL images revealed that there were two types of zircon. One type had core-rim structures; cores exhibited oscillatory, planar, or nebulous zoning, whereas rims showed weak oscillatory zoning or were unzoned ( Figure 6). The second type of zircon was dark under CL and had homogeneous internal structures with weak or no oscillatory zoning. Thirty-eight spots were analyzed from 38 zircon grains in SPMG. Three analyses of inherited cores yielded 206 Pb/ 238 U ages between 573 ± 9 to 639 ± 4 Ma (Table S1). Of the thirty-five analyses on rims of zircons with core-rim structures or on homogeneous zircon, one analysis had a 206 Pb/ 238 U age of 240 ± 3 Ma, and the other 34 analyses yielded 206 Pb/ 238 U ages between 424 ± 5 to 454 ± 5 Ma. After eliminating 11 analyses with discordance values >5%, the other 23 yielded a weighted mean 206 Pb/ 238 U age of 435.4 ± 1.8 Ma (MSWD = 0.77) (Figure 7c).

Biotite Monzogranite
Zircon from biotite monzogranite DG2 was mainly long or short prismatic, subhedral to euhedral, with lengths of 50 to 200 µm. In CL images, most zircon had no inherited cores, was dark under CL, and exhibited weak oscillatory zoning. Some zircon showed core-rim structures with oscillatory-zoned cores, whereas rims exhibited weak oscillatory zoning or no zoning ( Figure 6). Twenty-nine spots were analyzed from 29 zircon grains. Nine of 12 analyses of cores had 206 Pb/ 238 U ages of 413 ± 8 to 449 ± 9 Ma with a weighted average of 435.1 ± 7.0 Ma (MSWD = 2.1), and the other 3 analyses yielded 206 Pb/ 238 U ages of 670 ± 28, 718 ± 14, and 741 ± 7 Ma. The other 17 of 29 analyses were on oscillatory zones and gave 206 Pb/ 238 U ages between 233 ± 3 to 252 ± 5 Ma. After excluding eight analyses with discordance values of >5%, the remaining nine analyses yielded a weighted mean 206 Pb/ 238 U age of 235.3 ± 2.3 Ma (MSWD = 0.13) (Figure 7e), which is interpreted as the timing of magmatic crystallization.
Zircon from biotite monzogranite DG3 was mainly long or short prismatic, subhedral to euhedral, and had lengths of 50 to 220 µm. Most zircon had no inherited cores, was dark under CL, and had strong to weak oscillatory zoning; some zircons showed core-rim structures with oscillatory-zoned cores and rims with weak oscillatory zoning ( Figure 6). Nineteen spots were analyzed from 19 zircon grains. Of three analyses of cores, one yielded a 207 Pb/ 206 Pb age of 1128 ± 69 Ma, and the other two analyses gave 206 Pb/ 238 U ages of 437 ± 4 and 443 ± 6 Ma (Table S1). The other 16 analyses on oscillatory zones yielded 206 Pb/ 238 U ages ranging from 234 ± 3 to 255 ± 4 Ma. After excluding ten analyses with discordance values >5%, the remaining six concordant analyses yielded a weighted mean 206 Pb/ 238 U age of 241.0 ± 4.3 (MSWD = 1.3) (Figure 7f), which is interpreted as the crystallization age of the biotite monzogranite.

Summary of U-Pb Zircon Ages from Migmatites
Analyses of inherited cores in zircon from migmatites had a wide range of ages (c. 2480-553 Ma) with the largest population at c. 1000 Ma. This is distinctly different from the age of inherited zircons from the biotite monzogranites, which ranged from c. 413 to 1128 Ma with the largest population at c. 440 Ma. U-Pb ages of metamorphic zircon from migmatites were mainly Caledonian in age, but some Indosinian ages were similar to the timing of biotite monzogranite crystallization. CL images showed that the metamorphic zircon had core-rim or homogeneous zoning, which indicates that they had undergone different degrees of resorption and recrystallization [54][55][56]. Therefore, zircon from the Fuhuling migmatites recorded metamorphic zircon growth during both the Caledonian and Indosinian.

Metatexite
Trace element data of 39 analyses of 39 zircon grains from sample SPM are shown in Table S2

Nature of the Protoliths
All of the metatexite, diatexite, and biotite monzogranite samples collected from Fuhuling contained inherited zircon. The ages of inherited zircon in granitoid rocks and migmatites is often used to elucidate the age and nature of their sources [59][60][61]. Concentrations of the REE in zircon provide additional information on their geological significance [56,62]. Inherited zircon from metatexite and diatexite has a large range of Th/U ratios (0.023 to 1.64), and most of them have negative Eu anomalies and steep chondrite-normalized HREE slopes, whereas a few of them have smaller negative Eu anomalies and flatter HREE slopes; these variations may represent different source compositions [63]. This inherited zircon from metatexites and diatexites gives a wide range of ages (c. 553-2480 Ma) with a dominant population at c. 1000 Ma, which is broadly similar to the age distribution of detrital zircon from basement metasedimentary rocks in the Yunkai region ( Figure 9). This suggests that the protoliths of the migmatites are probably metasedimentary rocks. By contrast, inherited zircons from biotite monzogranite show a wide range of ages (c. 413-1050 Ma) with a dominant population at c. 440 Ma and Th/U ratios greater than 0.1 (0.12-1.00). This zircon has chondrite-normalized REE patterns characterized by steep HREE slopes, significantly negative Eu anomalies, and prominent, positive Ce anomalies, suggesting that they are of magmatic origin [56,57]. This indicates that the migmatites and biotite monzogranites have different sources.
age distribution of detrital zircon from basement metasedimentary rocks in the Yunkai region ( Figure 9). This suggests that the protoliths of the migmatites are probably metasedimentary rocks. By contrast, inherited zircons from biotite monzogranite show a wide range of ages (c. 413-1050 Ma) with a dominant population at c. 440 Ma and Th/U ratios greater than 0.1 (0.12-1.00). This zircon has chondrite-normalized REE patterns characterized by steep HREE slopes, significantly negative Eu anomalies, and prominent, positive Ce anomalies, suggesting that they are of magmatic origin [56,57]. This indicates that the migmatites and biotite monzogranites have different sources.  [36]; gneiss from Luoding [64]; paragneiss from Xiejie [29]; biotite plagioclase gneiss from Xinyi [65]; and metamorphic feldspar quartz-sandstone from Luola Village, southwest of Luoding, and biotite gneiss from the east of Xinyi [42].

Timing of Crustal Partial Melting in Yunkai Region
In the Fuhuling migmatite complex, metamorphic zircon from metatexite (SPM and FHM) and diatexite (SPMG and FHMG) has both core-rim and homogeneous structures consistent with growth during the crystallization of anatectic melt [54][55][56][66][67][68]. Zircon with core-rim structures has cores with oscillatory, planar, or nebulous zoning, which is surrounded by overgrowths with weak oscillatory zoning or no zoning. Homogeneous-structured neoblastic zircon has no cores and weak or no oscillatory zoning. These different internal structures suggest that the metamorphic zircons experienced variable degrees of resorption and recrystallization [54,57,69]. Metamorphic zirconformed by dissolution-reprecipitation also has specific trace element characteristics [54,57,70]. U-Pb ages of homogeneous-structured neoblastic zircons or zircon rims with core-rim structures in samples SPM, FHM, SPMG, and FHMG can be divided into two groups.
One group of 78 analyses yields 206 Pb/ 238 U ages that range from 431 ± 5 to 454 ± 5 Ma. Compared with inherited zircon in these samples, these analyses also have chondrite-normalized REE patterns with steep HREE slopes, negative Eu anomalies, and prominent, positive Ce anomalies, which are characteristics of zircon crystallized from melt [54,56,70]. Therefore, the weighted mean 206 Pb/ 238 U ages of 443.3 ±2.7, 44 ±13, 435.4 ± 1.8, and 437.5 ± 3.0 Ma from zircon rims with core-rim structures or homogeneous-structured neoblastic zircon of samples SPM, FHM, SPMG, and FHMG, respectively, are interpreted as the general timings of partial melting. More specifically, because zircon is expected to dissolve during heating and partial melting and grow during cooling and melt crystallization [71][72][73][74][75], these ages most likely represent the timing of neoblastic zircon growth during melt crystallization. The second group of 13 analyses yields 206 Pb/ 238 U ages of 235 ± 3 to 244 ± 4 Ma and highly variable Th/U ratios (0.02-1.92). Compared with inherited zircon of these samples, chondrite-normalized REE patterns of these analyses are characterized by smaller, negative Eu anomalies, positive Ce anomalies, and flatter HREE slopes. These arefeatures similar tometamorphic zircon that grows or is recrystallized in the presence of a metamorphic fluid [54,56,70].
Recently,   [31] obtained zircon U-Pb ages from gneissic migmatites from the Longxiu-Yunlu section of the Gaozhou area, and Ke et al. (2018) [42] obtained zircon U-Pb ages from leucosome in migmatites from Luoding, Guangdong, and Xinyi, Guangdong. Both studies documented metamorphic zircon growth mostly in the Early Paleozoic(~440 Ma) and minor growth in the Early Mesozoic(~230 Ma).   [31] suggested that the Early Paleozoicages represent the timing of crustal partial melting, and the Early Mesozoicages may be related to low-temperature metamorphism during Indosinian movement. Ke et al. (2018) [42] interpreted both the Early Paleozoic ages and the Early Mesozoic ages as the timing of melt crystallization and suggestedthat the Yunkai region wasaffected by crustal partial melting in the Early Paleozoic and the Early Mesozoic.
Compared with the zircon analyzed in these studies, the morphology and trace element compositions of zircon from metatexite and diatexite at Fuhuling are similar to that of zircon from gneissic migmatites described byWang et al. weak or no oscillatory zoning. These different internal structures suggest that the metamorphic zircons experienced variable degrees of resorption and recrystallization [54,57,69]. Metamorphic zirconformed by dissolution-reprecipitation also has specific trace element characteristics [54,57,70]. U-Pb ages of homogeneous-structured neoblastic zircons or zircon rims with core-rim structures in samples SPM, FHM, SPMG, and FHMG can be divided into two groups.
One group of 78 analyses yields 206 Pb/ 238 U ages that range from 431 ± 5 to 454 ± 5 Ma. Compared with inherited zircon in these samples, these analyses also have chondrite-normalized REE patterns with steep HREE slopes, negative Eu anomalies, and prominent, positive Ce anomalies, which are characteristics of zircon crystallized from melt [54,56,70]. Therefore, the weighted mean 206 Pb/ 238 U ages of 443.3 ±2.7, 44 ±13, 435.4 ± 1.8, and 437.5 ± 3.0 Ma from zircon rims with core-rim structures or homogeneous-structured neoblastic zircon of samples SPM, FHM, SPMG, and FHMG, respectively, are interpreted as the general timings of partial melting. More specifically, because zircon is expected to dissolve during heating and partial melting and grow during cooling and melt crystallization [71][72][73][74][75], these ages most likely represent the timing of neoblastic zircon growth during melt crystallization. The second group of 13 analyses yields 206 Pb/ 238 U ages of 235 ± 3 to 244 ± 4 Ma and highly variable Th/U ratios (0.02-1.92). Compared with inherited zircon of these samples, chondrite-normalized REE patterns of these analyses are characterized by smaller, negative Eu anomalies, positive Ce anomalies, and flatter HREE slopes. These arefeatures similar tometamorphic zircon that grows or is recrystallized in the presence of a metamorphic fluid [54,56,70].
Recently,   [31] obtained zircon U-Pb ages from gneissic migmatites from the Longxiu-Yunlu section of the Gaozhou area, and Ke et al. (2018) [42] obtained zircon U-Pb ages from leucosome in migmatites from Luoding, Guangdong, and Xinyi, Guangdong. Both studies documented metamorphic zircon growth mostly in the Early Paleozoic(~440 Ma) and minor growth in the Early Mesozoic(~230 Ma).   [31] suggested that the Early Paleozoicages represent the timing of crustal partial melting, and the Early Mesozoicages may be related to low-temperature metamorphism during Indosinian movement. Ke et al. (2018) [42] interpreted both the Early Paleozoic ages and the Early Mesozoic ages as the timing of melt crystallization and suggestedthat the Yunkai region wasaffected by crustal partial melting in the Early Paleozoic and the Early Mesozoic. Indosinian zircon from this study and from [30]. Chondrite values are from [58]. Figure 10. Comparison of chondrite-normalized REE patterns of (a) Caledonian zircon and (b) Indosinian zircon from this study and from [30]. Chondrite values are from [58].

Migmatite Morphology and Development
The structures of migmatites formed by partial melting are usually heterogeneous. If migmatites are affected by later high-grade metamorphism, crustal partial melting (anatexis), and even crustal magmatism, then the pre-existing structures formed by partial melting may homogenize or disappear depending on the extent of subsequent melting and the amount of deformation [18][19][20]25]. Therefore, the macro-and microstructures preserved by migmatites provide important information to determine if they were substantially affected by multiple high-temperature tectonothermal events. Considering that U-Pb zircon ages indicate that the Fuhuling migmatite did not partially melt during the Early Mesozoic, the leucosome-melanosome framework therefore records Early Paleozoic (c. 435-444 Ma) melt organization and can be directly linked to deformation during Wuyi-Yunkai orogenesis.
The metatexite-diatexite profile at Fuhuling has a complex geometry at the outcrop scale and exhibitsstructural stratification. Leucosomes in the metatexite are preferentially oriented, whereas leucosome in the diatexite show no preferred orientation (Figure 3). This indicates that the configuration of leucosomes in migmatite changes, and their directionality is weakened during the transition from metatexite to diatexite. These phenomena may be caused by changes in melt fraction under static conditions [16,25]. Along the profile, the metatexite of layer I and layer II has a typical stromatic structure; the metatexite of layer III contains a large number of flowing folds with different shapes, and the diatexite of layer IV exhibits schlieren and nebulitic structures. All of these are typical mesostructures of migmatite associated with partial melting [20,25,76]. The mesostructures of metatexite and diatexite are distinctly different, which is manifested by the geometry and abundance of leucosome. Leucosome geometry in migmatites is inferred to record evidence of active melt flow, rather than stagnant melt [18,48]. However, leucosomes represent the final vestiges of melt before solidification of the anatectic system. The high abundance of leucosome in some migmatites may reflect melt accumulation rather than a large degree of partial melting of the host. The morphological differences between the metatexite and diatexite layers are a product of the proportion of melt in the system and the amount of pervasive ductile deformation. Metatexite is characterized by the presence of leucosomes parallel to the dominant foliation in the host or superimposed on a transposed inherited compositional layering of sedimentary or magmatic origin [10,18,48]. There is a higher proportion of leucosome in fold hinges than fold limbs, which suggests that these were areas of melt accumulation and probably represent preferred flow pathways for the anatectic melt [10,18,19,48]. The transition from metatexite to diatexite in the Fuhuling migmatites is associated with an increase in the proportion and structural variety of leucosome [10,48]. Diatexite with schlieren or nebulitic structures contains partially resorbed and dismembered paleosome or residuum within a coarser-grained, diffuse leucosome; the preservation of these features suggests that there was no strong pervasive deformation during partial melting [10,19,48]. Diatexite with a schollen structure is characterized by disruption of the more competent layers (typically paleosome or residuum) into "schollen" or "raft" fragments with different orientations, which indicates a more substantial component of ductile deformation compared with the schlieren-bearing diatexite [10,19,48]. The microstructure in migmatites continually readjusts tochanges of partial melting, which reflects information about the crystallization of anatectic melt rather than the melt-producing reactions [5,8,13]. The microstructures in metatexite include fine-grained interstitial quartz, feldspar, and drop-shaped quartz, which probably represent crystallization of the last vestiges of anatectic melt [17,[22][23][24]. Elongated or preferred orientation of feldspar in diatexite suggests a component of magmatic flow while the diatexite was still molten [17,[22][23][24].

Mechanisms and Causes of Partial Melting at Fuhuling
The Wuyi-Yunkai orogen is an intraplate orogen, and the geodynamic models for its evolution during the Early Paleozoic are still debated [29,31,77]. However, intraplate orogeny probably caused crustal shortening and thickening and resulted in the closure of the pre-existing Nanhua Rift [29,31,77]. The heat needed to initiate and sustain partial melting in the region was probably the result of crustal thickening associated with the Wuyi-Yunkai orogen and concomitant radioactive heat generation from the sedimentary rocks [29,31,77,78] as well as possible residual heat generated by Neoproterozoic plume activity in the region [29,77]. Together, this geodynamic environment of an intraplate orogeny overriding an old plume system is inferred to have led to partial melting of the metasedimentary rocks and some Cathaysian basement rocks [29,31,79].
Partial melting of crustal rocks can proceed through hydrate-breakdown melting, where hydrous minerals are consumed (e.g., biotite, muscovite, amphibole) and anhydrous peritectic minerals (e.g., garnet, orthopyroxene) are generated [80], or through the influx of an externally derived hydrous fluid (e.g., [81,82]). Hydrate-breakdown melting of fertile protoliths can generate up to~15 vol.% melt through muscovite-breakdown reactions and up to 40 vol.% by biotite-breakdown melting at ultra-high-temperature (>900 • C) conditions [83]. At the same pressure-temperature conditions, fluid-present melting can generate more melt than hydrate-breakdown melting as long as there is a large enough influx of H 2 O into the system. In addition, fluid-present partial melting results in a net volume reduction and promotes the retention of melt, as opposed to fluid-absent melting that requires an increase in volume, which is commonly accommodated by melt loss from the system. Hence, the absence of anhydrous peritectic minerals in the Fuhuling migmatites is inconsistent with fluid-absent, hydrate-breakdown melting.
In general, migmatites in the Wuyi-Yunkai orogen generally contain anhydrous peritectic minerals [31,50], which is consistent with hydrate-breakdown melting. However, the absence of anhydrous peritectic minerals in the Fuhuling migmatites in the Yunkai region is inconsistent with hydrate-breakdown melting. Furthermore, the amount of leucosome in the Fuhuling migmatites is >40 vol.% in many cases ( Figure 11); most of this leucosome is inferred to be the product of in situ or in source partial melting. These large volumes of leucosome (i.e., former anatectic melt and entrained minerals) are inconsistent with fluid-absent melting unless ultra-high-temperature conditions were reached (c.f. [84]), these temperatures are not compatible with the abundance of biotite and hornblende in the Fuhuling migmatites. Therefore, migmatites at Fuhuling probably underwent fluid-present melting in the Early Paleozoic, and this mechanism was partly responsible for the morphology of the migmatites.
Minerals 2019, 9, x FOR PEER REVIEW Fuhuling migmatites and proposed that, in addition to MCT and SLT, there is a third rheological threshold termed the framework-melting transition (FMT). The FMT is characterized by a sharp decrease of aggregate strength at melt fractions of ~21 vol.% and can be interpreted as a transition of the solid framework from compact to fused. An increased proportion of melt results in a progressively weaker crust [93] and has significant effects not only on the nature of partially molten rocks(such as the Fuhuling migmatites) but also on the crustal-scale tectonic processes [10,19,20] and the evolution of the Wuyi-Yunkai orogen [29][30][31][32]39,40].
The structural transition from metatexite to diatexite, such as that observed in the Fuhuling migmatites, is related to the evolution of the rheology of partially molten rocks, which is the result of a complex and dynamic interplay between melt production, melt segregation, melt migration, and melt redistribution [10,19,20]. The amount of melt in migmatitic crust and its role in controlling its rheology can be investigated through macro-and micro-structural analyses based on the distribution and proportion of the melt at the outcrop scale [10,[18][19][20]. To assess the influence of the aggregate strength and melt fraction of the partially melted protolith on the structures of Fuhuling migmatite in the Yunkai region, we quantified leucosome proportions along the metatexite-diatexite transect ( Figure 11). The methodology of our approach is the same as described in detail by Zhang et al. (2017) [38] and Chen et al. (2017) [92]. We measured leucosome proportions in each layer of the metatexite-diatexite profile and linked this with the possible rheological transitions controlled by melt fraction (Table S3). Based on the relationship between aggregate strength and melt fraction of the partially melted protolith at Fuhuling [92] and our structural analysis of leucosome networks, we now propose a conceptual model for the variable morphologies observed in the Fuhuling migmatites (Figure 12a-i). This is essential for establishing a link between rheological transitions and crustal partial melting in the Yunkai region during Early Paleozoic Wuyi-Yunkai orogeny.
At melt fractions < 7 vol.%, only a small amount of melt formed at the interstices of mineral grains (Figure 12f). The partially melted protolith remained coherent and banded, and the solid framework and structure of protolith were not compromised (Figure 12b). At the onset of partial

Rheology of the Yunkai Region during the Wuyi-Yunkai Orogeny
Large phanerozoic orogens commonly contain migmatites formed by crustal partial melting, and this process has important implications for the rheology of the crust during orogenesis [3,11,78]. The rheology and strength of partially molten rocks change significantly at two thresholds that are relevant for most migmatite terranes [11,25,85], including the Yunkai region. The first is known as the melt connectivity transition (MCT) [86][87][88][89], which is generally thought to occur at 7 vol.% melt in a static system [88]. A second threshold is known as the solid-to-liquid transition (SLT) [88], which is similar to the rheologically critical melt percentage (RCMP), and is characterized by the less-pronounced aggregate strength drop at a melt fraction around 40 vol.% [90][91][92].
Fuhuling migmatites have variable proportions of leucosome, which is inferred to approximate the variable amount of anatectic melt that was present during orogenesis. Chen et al. (2017) [92] investigated the relationship between aggregate strength versus leucosome fraction in the Fuhuling migmatites and proposed that, in addition to MCT and SLT, there is a third rheological threshold termed the framework-melting transition (FMT). The FMT is characterized by a sharp decrease of aggregate strength at melt fractions of~21 vol.% and can be interpreted as a transition of the solid framework from compact to fused. An increased proportion of melt results in a progressively weaker crust [93] and has significant effects not only on the nature of partially molten rocks (such as the Fuhuling migmatites) but also on the crustal-scale tectonic processes [10,19,20] and the evolution of the Wuyi-Yunkai orogen [29][30][31][32]39,40].
The structural transition from metatexite to diatexite, such as that observed in the Fuhuling migmatites, is related to the evolution of the rheology of partially molten rocks, which is the result of a complex and dynamic interplay between melt production, melt segregation, melt migration, and melt redistribution [10,19,20]. The amount of melt in migmatitic crust and its role in controlling its rheology can be investigated through macro-and micro-structural analyses based on the distribution and proportion of the melt at the outcrop scale [10,[18][19][20].
To assess the influence of the aggregate strength and melt fraction of the partially melted protolith on the structures of Fuhuling migmatite in the Yunkai region, we quantified leucosome proportions along the metatexite-diatexite transect ( Figure 11). The methodology of our approach is the same as described in detail by Zhang et al. (2017) [38] and Chen et al. (2017) [92]. We measured leucosome proportions in each layer of the metatexite-diatexite profile and linked this with the possible rheological transitions controlled by melt fraction (Table S3). Based on the relationship between aggregate strength and melt fraction of the partially melted protolith at Fuhuling [92] and our structural analysis of leucosome networks, we now propose a conceptual model for the variable morphologies observed in the Fuhuling migmatites (Figure 12a-i). This is essential for establishing a link between rheological transitions and crustal partial melting in the Yunkai region during Early Paleozoic Wuyi-Yunkai orogeny.
Minerals 2019, 9, x FOR PEER REVIEW buoyancy-driven magma ascent [97,98], and the heat advected with migrating melt may promote amplification of thermal anti form, extending the zone of partial melting to shallower crust that weakens the crust during orogeny [13,85,99]. The buoyancy contrast between metatexite and relatively melt-rich diatexite may also promote the buoyant rise of partially molten diatexite, which is a common inference from gneiss domes [99][100][101]. The Fuhuling diatexites have complex morphologies with leucosomes that lack preferred orientations, which contrasts with the preferred orientation of leucosomes in the overlying metatexites. This suggests that the loss of cohesion at high melt fractions in the diatexites possibly coupled with buoyancy-driven uplift were important factors in the heterogeneity of the leucosomes in diatexites compared to leucosome in the metatexites. Overall, the Fuhuling migmatites show rheological transitions caused by a change in melt fraction. Given that structural processes are expected to be scale-invariant in migmatite terranes [26,102], the rheological transitions documented in the Fuhuling migmatites are probably reflective of the crustal-scale rheological heterogeneity that existed in the Yunkai region during the Early Paleozoic Wuyi-Yunkai orogen ( Figure 13). Consequently, the high melt fractions inferred in the deep crust may have played a key role in facilitating orogenic collapse in the Yunkai region in the Early Paleozoic [79].  (4), and as the partial melting proceeds, granite will eventually form. The 650 °C isotherm approximates the wet solidus of crustal rocks. (b-e) Macrostructural sketches of the "mixture" composed of partiallymolten rocks and melt, illustrating the rheological transitions due to changesin the melt fraction when going through the processes of melt production, melt segregation, melt migration, and melt redistribution, which correspond to the melt fractions of 0-7 vol.%, 7-21 vol.%, 21-40 vol.%, and > 40 vol.%, respectively (pink = melt, gray or yellow blocks = (meta)sedimentary rocks) (modified after [11,24]).   (4), and as the partial melting proceeds, granite will eventually form. The 650 • C isotherm approximates the wet solidus of crustal rocks. (b-e) Macrostructural sketches of the "mixture" composed of partiallymolten rocks and melt, illustrating the rheological transitions due to changesin the melt fraction when going through the processes of melt production, melt segregation, melt migration, and melt redistribution, which correspond to the melt fractions of 0-7 vol.%, 7-21 vol.%, 21-40 vol.%, and >40 vol.%, respectively (pink = melt, gray or yellow blocks = (meta)sedimentary rocks) (modified after [11,24]). At melt fractions <7 vol.%, only a small amount of melt formed at the interstices of mineral grains (Figure 12f). The partially melted protolith remained coherent and banded, and the solid framework and structure of protolith were not compromised (Figure 12b). At the onset of partial melting, the resulting melt formed fine-grained quartzofeldspathic films, which may have a granular macroscopic appearance, along the grain boundaries between the reactant minerals [10,25]. As the fraction of melt in the system increased due to progressive anatexis, the melt formed elongate mottles and intermittent strips (Figure 3a) [25], which are both typical stromatic structures of metatexites [10,25].
At melt fractions between 7 and 21 vol.%, the interstitial melt became connected along grain boundaries (Figure 12g). The partially melted protolith was still coherent and banded (Figure 12c). The melt then organized into thin parallel and laterally continuous leucosomes, which would give the system a layered structure (Figure 3b) [10,25]. This structure is also typical of stromatic metatexites and probably relates to short-range segregation of melt from its residuum in situ or in source [10,19,25].
At melt fractions between 21 and 40 vol.%, the proportion of melt increased due to framework melting but was still restricted within a coherent solid framework (Figure 12h). The partially melted protolith experienced a significant decrease in aggregate strength, coherency compositional banding gradually weakened, and stromatic structures began to deteriorate (Figure 12d). This formed metatexite with folds of a wide range of leucosome orientations as well as melt accumulation in the hinges of folds (Figure 3c,d), suggesting the local migration of melt during folding [19,25].
At melt fractions >40 vol.%, the solid framework was no longer connected ( Figure 12i). The coherency and banding of the partially melted protolith were destroyed, and pre-partial melting structures in the migmatites were destroyed (Figure 12e). More competent layers were progressively disrupted into tabular fragments, and typically paleosome and residuum were suspended in the melt forming "schollen" or "rafts" in diatexite (Figure 3f) [25]. Elongate or preferentially oriented quartz, plagioclase, and biotite are common in diatexite due to magmatic flow (Figure 4f,h) [20,25]. These characteristics indicate that the continuity of the solid framework was lost [10,88].
If partial melting proceeds and melt cannot be effectively drained, the accumulation of melt and settling of partiallymolten residuum will lead to the transition from a "mixture"-composed of the partially molten protolith and melt-to magmas that may ascend to form granite plutons [11] or remain in the anatectic crust. In the Yunkai region, Early Paleozoic migmatites and granites are considered to have a genetic link because of the synchronicity between crustal partial melting and magmatism [30][31][32]. Current geodynamic models of the Wuyi-Yunkai orogeny include crustal partial melting (anatexis) during the Early Paleozoic [29,31,77], but the rheological effects of a melt-rich crust on this intraplate orogen are unclear.During orogenesis, partial melting of the crust is heterogeneous given the variable fertilities of different lithologies [94]. If the crust incorporates large quantities of fertile (meta)sedimentary rocks, it may be able to generate voluminous anatectic melt [15,85]. In contrast to the fluid-present melting documented at Fuhuling, anhydrous peritectic minerals associated with higher-temperature hydrate-breakdown melting are present elsewhere in the Yunkai region [31]. The preservation of anhydrous peritectic minerals in anatectic rocks requires the loss of anatectic melt [95] and leaves behind a more refractory residue that is expected to be stronger than the unmelted protolith [96]. This contrasts with the large amount of leucosome observed at Fuhuling, which suggests substantial melt retention and possible melt accumulation. Consequently, the Yunkai region may represent a section through a heterogeneous crust with both melt-rich (weak) and melt-poor (strong) horizons, and this suggests that partial melting of the crust and the amount of melt in different crustal horizons was heterogeneous during the Wuyi-Yunkai orogeny ( Figure 13).  Spatial analysis of leucosomes at Fuhuling implies that a large proportion of melt remained in the deep crust, which is evidenced by the presence of diatexites and metatexites with high melt fractions at Fuhuling. Melt retention in the middle to deep crust has also been proposed in other migmatite terranes [26,46], but this contrasts with the general view that the continental crust is differentiated due to the ascent of granitic materials in dykes that have no or minimal interaction with the surrounding rock during ascent (c.f. [80]). Such crustal-scale mass transfer is facilitated by buoyancy-driven magma ascent [97,98], and the heat advected with migrating melt may promote amplification of thermal anti form, extending the zone of partial melting to shallower crust that weakens the crust during orogeny [13,85,99]. The buoyancy contrast between metatexite and relatively melt-rich diatexite may also promote the buoyant rise of partially molten diatexite, which is a common inference from gneiss domes [99][100][101]. The Fuhuling diatexites have complex morphologies with leucosomes that lack preferred orientations, which contrasts with the preferred orientation of leucosomes in the overlying metatexites. This suggests that the loss of cohesion at high melt fractions in the diatexites possibly coupled with buoyancy-driven uplift were important factors in the heterogeneity of the leucosomes in diatexites compared to leucosome in the metatexites. Overall, the Fuhuling migmatites show rheological transitions caused by a change in melt fraction. Given that structural processes are expected to be scale-invariant in migmatite terranes [26,102], the rheological transitions documented in the Fuhuling migmatites are probably reflective of the crustal-scale rheological heterogeneity that existed in the Yunkai region during the Early Paleozoic Wuyi-Yunkai orogen ( Figure 13). Consequently, the high melt fractions inferred in the deep crust may have played a key role in facilitating orogenic collapse in the Yunkai region in the Early Paleozoic [79].

Conclusions
1. Zircon U-Pb geochronology and trace element analysis indicate that biotite monzogranites are not derived from the migmatites at Fuhuling. Partial melting at Fuhuling in the Yunkai region occurred during Early Paleozoic Wuyi-Yunkai orogenesis. Minor zircon dissolution-reprecipitation in the Indosinian was probably related to lower-temperature metamorphism.
2. Migmatites at Fuhuling were produced by fluid-present melting, have morphologies ranging from metatexite to diatexites, and were not affected by deformation associated with Indosinian movement.
3. Morphological differences across the metatexite-diatexite transition in the Fuhuling migmatites may have been caused by a change in melt connectivity and continuity of the solid framework because there were variable proportions of solid and melt during partial melting and deformation.
4. The abundance of leucosome in the Fuhuling migmatites suggests that it acted as a melt-rich and buoyant horizon. The rheological transitions documented in the Fuhuling migmatites may be analogous with the rheological transitions that operated the deep melt-rich crust after crustal thickening associated with the Early Paleozoic Wuyi-Yunkai orogenesis.