Corrosion of Longquan celadons in the marine environment: study on the celadons from the Dalian Island shipwreck of the Yuan Dynasty

The Dalian () Island shipwreck of the Yuan Dynasty (1271–1368 CE) was located in the sea area of Pingtan (), Fujian () Province, China. A total of 603 Longquan () celadon wares were excavated, some of which have been severely corroded. In this study, two celadon specimens with severe corrosion were selected to investigate the corrosion mechanism in the marine environment. Optical microscopy (OM), scanning electron microscopy equipped with energy dispersive X-ray spectrometry (SEM–EDS), X-ray diffraction analysis (XRD), Micro-Raman spectroscopy, and thermal expansion analysis were applied to analyze the microstructure of the glaze, corrosion morphology, and composition of corrosion products, and the corrosion mechanism was discussed in depth. The results showed that these two celadon wares are opaque matt glaze, in which the anorthite crystallization-phase separation structure with low chemical stability was found, which was more likely to be corroded. There are three possible stages in the corrosion process of the glaze. The first stage is mainly the long-term cation exchange reaction. Then, because of the corrosion of the CaO-rich droplet phase and the crystal-glass phase interface layer, the porous structure and many microcracks may be formed, promoting the corrosion process. Finally, many corrosion craters, caused by the shedding of anorthite crystals, may appear on the glaze surface. The quintinite group minerals in the glaze were probably formed during the corrosion process, which may be related to the marine buried environment and the corrosion of the glaze. This study discussed the influencing factors of the corrosion process of celadon wares in the marine environment from the perspective of microstructure for the first time, which is of great significance to the study of the corrosion mechanism of ancient ceramics.


Introduction
Chinese porcelain was exported overseas as an important commodity as early as the Tang Dynasty (618-907 CE) and was loved by people all over the world in history [1]. However, in the history of the East-West porcelain trade for thousands of years, many merchant ships loaded with Chinese porcelain wares, unfortunately, sank into the sea due to natural and man-made disasters and were buried by marine sediment [2,3]. These ceramics have experienced hundreds or even thousands of years of alteration process caused by various physical, chemical, and biological factors in the marine environment. Various corrosion phenomena have seriously affected the beauty of these precious cultural relics, and have brought severe challenges to the treatment and preservation after the Open Access *Correspondence: liwd@mail.sic.ac.cn 1 Shanghai Institute of Ceramics, Chinese Academy of Sciences, Shanghai 201899, China Full list of author information is available at the end of the article excavation. At present, more than 90% of the cultural relics excavated by underwater archaeology in China were ceramic wares [4], but there is still a lack of in-depth scientific cognition of the damage to the ceramics in the marine environment, and a systematic protection scheme has not yet been formed. This had become an important issue that urgently needs to be resolved in the protection of cultural relics excavated underwater.
At present, ion exchange (H + /alkali and alkaline earth metal ions), Si-O-Si network hydrolysis, and the resulting microcrack and the loss of protonation of the network were considered to be the important causes of glass corrosion in the seawater [5][6][7][8]. In the protection field of cultural relics, some researchers studied the effect of soluble salt on ceramics and the method of desalination [9][10][11][12][13]. There have also been some preliminary studies on the formation and removal of different types of concretion and contaminant [14][15][16], and some researchers have discussed the secondary damage of ceramics in the process of conservation and treatment after excavation [17,18]. Although the scientific cognition of ceramics excavated underwater has made considerable progress, the research on the corrosion of ancient ceramics in the marine environment is still not enough, especially the porcelains that are generally considered to have good corrosion resistance. The corrosion mechanism of the celadon glaze in the natural environment has not yet reached a consensus, and there is a lack of a comprehensive understanding of the various factors affecting the corrosion process, which is obviously not conducive to the subsequent treatment and protection work.
In this study, the ancient celadon wares from the Dalian Island shipwreck of the Yuan Dynasty were selected as a typical case, aiming to explore the corrosion process of celadons in the marine environment and its influencing factors. The shipwreck is located in the coastal waters of Pingtan, Fujian Province, China ( Fig. 1), with a water depth of 15-18 m. The underwater archaeological excavations of the shipwreck were carried out in 2007. A total of 603 celadon wares were excavated from the shipwreck site, all of which were exported celadons from Longquan kiln [19]. In this paper, the corrosion of the ancient celadons was discussed from the perspectives of chemical composition and microstructure, which is very important for the protection of cultural relics excavated from the marine environment.

Samples
In order to systematically study the corrosion of the celadon wares from the Dalian Island shipwreck, 11 pieces of celadon shards with different corrosion degrees, provided by Fujian Museum, were selected as the research object. These shards have been immersed in deionized water for desalination in the museum before they were provided to us. The desalination time is 5 weeks.
The samples were numbered from 01 to 11 and the detailed information can be seen in Table S1 and S2 in the Additional file 1. After observation and analysis, two severely corroded fragments (Nos. 08 and 09) among them were found to have significant representativeness, so they were chosen as the experimental samples in this paper for the key analysis. The photos and appearance descriptions are listed in Table 1.
The characteristic parts of the samples were cut from the shards and then cleaned by deionized water and ethanol in an ultrasonic cleaner (2 times for deionized water, 1 time for ethanol, 15 min each time). After that, all samples were dried in a drying oven at 110 ℃ for 3 h.

Methods
X-ray diffraction (Bruker D8 DISCOVER) using Cu Kα radiation was applied to investigate the phase components of the glaze surface. The voltage is 40 kV, the current is 40 mA, the scan rate is 5(°)/min, and the angle range is 10°-80°.
Micro-Raman spectroscopy (Horiba XploRA one 532 nm) was used to examine the crystals and glass phase in the glaze. The wavenumber calibration of the peak position was at a × 50 objective with a piece of crystal silicon. A confocal microscope (Olympus BX41) was employed to focus the laser beam on the test areas with ×10 and ×50 objectives. The baseline was subtracted using the LABSPEC software. Then, the Raman spectra of glassy phase were de-convoluted and curvefitted with Gaussian functions using the Origin software [20,21].
The macroscale morphology of the glaze surfaces and polished cross-sections were observed with an optical microscope (Keyence VHX-2000). The micromorphology and micro-area chemical compositions of glaze surface and polished cross-section of samples were characterized using field emission scanning electron microscope (FEI Magellan400) equipped with EDS.
The firing temperature of the bodies was examined by dilatometer (NETZSCH DIL 402 C). To prepare the test sample for dilatometer, the glaze layer was ground off first, and the body was made into a strip of 5 mm × 5 mm × 25 mm. The heating rate was 5℃/min and the approximate value of the firing temperature of the sample was obtained according to the inflection point in the thermal expansion curve [22,23].

Results
The microstructure of the original glaze Figure 2a is the OM image of the glaze surface of No. 08. The glaze is almost opaque and the color is uneven. It can be found that many long rod-like crystals are spread on the glaze surface with a length of up to 100 μm and a width of about 5 μm. A cluster of radially interwoven crystals is shown in Fig. 2b. The EDS analysis results show that the long rod-like crystals have a high content of MgO and CaO (Table 2), which are Table 1 The detailed information of the celadon samples

Sample No
Digital photo Appearance and state of conservation  likely to be diopside crystals based on the stoichiometric ratio, and some Mg 2+ and Si 4+ ions in the crystal are partially replaced by Al 3+ ions and a few Fe 3+ ions. The Raman spectrum of the long rod-like crystal is in good agreement with the reference spectrum of the diopside crystal (Fig. 3a), which further confirms its existence. Besides, the broad peaks around 500 cm −1 , 767 cm −1 , and 1000 cm −1 are the characteristics of the silicate glass. It can be seen in the OM image of the polished crosssection of No. 08 (Fig. 2c) that many cracks penetrate the glaze layer. The SEM image shows (Fig. 2d) that many unmelted quartz grains (the Raman spectrum is shown in Fig. 3b) and bubbles with different sizes are densely distributed in the glaze. A large number of column-like anorthite and crystals with high P 2 O 5 content are distributed in the glaze ( Fig. 2e and f ), the Raman spectrum of the column-like crystals is shown in Fig. 3c, and the EDS analysis results are listed in Table 2. The presence of phosphorus-containing crystals may be related to the use of botanic ash, a glaze flux widely used in ancient China [24,25].
The glaze layer of sample No. 09 has been corroded severely, and the remaining opaque glaze was loose and porous. Figure 4a is the OM image of the glaze surface of No. 09. There are plenty of yellow substances on the surface, which are mixed with the remaining glaze. The SEM image of the glaze surface (Fig. 4b) shows that the remaining glaze on the left is already full of interconnected holes. Under further magnification (Fig. 4c), it is found that there are many long rod-like diopside crystals and some quartz grains in the glaze (see Table 3 for the EDS analysis results), which has a certain similarity with sample No. 08.

The micromorphology of the corroded glaze of No.08
There are many corrosion craters on the glaze surface, and the edges of the craters are opaque white. The cracks and corrosion craters on the glaze surface of No. 08 can be observed more clearly from the SEM image (Fig. 5a). There are a large number of craters with different depths near the crystals. The craters near the cracks on the glaze surface are particularly prominent (Fig. 5b), and they are connected to each other in a large area. The magnified SEM images of the craters ( Fig. 5c and d) show that the column-like anorthite crystals (Table 4) in the glaze are exposed outside, and the phase-separated droplets in the nearby glass phase were corroded into a porous structure. At the edge of the crater near glaze surface, it can be found that a large number of column-like crystals with lower contrast and rhomboid crystals with higher contrast (Fig. 5e), which may be anorthite crystals and diopside crystals according to the morphology and previous EDS analysis results. There are obvious interstices between the column-like anorthite crystals and the nearby glass phase. The magnified SEM images ( Fig. 5f ) showed that there is a sub-micron phase-separation structure in the glass phase near the anorthite crystals, and the original phase-separated droplets with higher contrast have been corroded, forming a porous structure. The phase-separated droplets at the interface between the anorthite crystals and glass phase were also corroded by seawater, forming interstices.
To characterize the alteration of the glass phase, the Raman spectra of the glass phase were recorded at the surface craters (corroded region) and cross-section (original glaze), and the representative results were shown in the Fig. 6a and b.Two main broad peaks at ∼480 cm −1 and ∼1000 cm −1 originated from the SiO 4

Fig. 3
Raman spectra recorded in the crystalline zone of the glaze of sample No. 08, a long rod-like diopside crystals and glass phase on the glaze surface; b quartz grains; c column-like crystals and glass phase, and for comparison, the reference spectra of diopside (RRUFF ID: R040009) [56], quartz (RRUFF ID: R040031) [57], and anorthite (RRUFF ID: R040059) [58] are cited tetrahedron bending and stretching modes, respectively [21]. Four typical components (Q 0 , Q 1 , Q 2 , and Q 3−4 ) could be extracted from the stretching envelopes, according to the description with a set of more or less connected SiO 4 tetrahedron [21], and their positions were listed in Table 5. The area ratios of bending peak over stretching peak (A 500 /A 1000 ) were calculated to analyze the polymerization degree of the glaze [20]. The results (Table 5) showed that A 500 / A 1000 ratio of the surface crater is lower than that of the cross-section, indicating a different polymerization degree. And the relatively lower A 500 /A 1000 ratio (~ 1.5) may also correspond to the lower firing temperature [20,26]. The wavenumber downshift of the Q 3−4 component of the surface crater (Table 5) may be related to a lengthening of the Si-O bond because of the interaction of the oxygen atoms with newly inserted protonic species [6]. The area ratios of three components (Q 1 , Q 2 , and Q 3−4 ) over stretching peak subtracting Q 0 ( Table 5) were also calculated to analyze the alteration of each component [27]. Figure 6c shows that A 3-4 /(A 1000 -A 0 ) ratios were larger than 0.45, indicating the sheet-like region and tectosilicates (Q 3−4 ) mainly preserve the main connected framework [20]. The A 1 /(A 1000 -A 0 ) ratio of surface crater increased a lot while A 2 /(A 1000 -A 0 ) and A 3-4 /(A 1000 -A 0 ) decreased compared to that of the cross-section,  indicating the breaking of Si-O bonds associated with corrosion [6,27]. Therefore, the integrity of the glass network seemed to be altered after seawater corrosion.

The micromorphology of the corroded glaze of No.09
To further study the corrosion process of sample No. 09, an area where the glaze remained more was selected for in-depth characterization. The OM image of the polished cross-section of this area is shown in Fig. 7a. The body of the area has a concave face, the glaze layer is obviously thicker than other areas, and the alteration of the glaze layer in the concave area is relatively slight. In the SEM image (Fig. 7b), a crack penetrating from the glaze surface to the glaze-body interface can be observed. There  is a corrosion region near the crack, where phase separation structure and many column-like anorthite crystals are found (P1 and P2 in Fig. 7d, EDS analysis results are listed in Table 6). The corrosion morphology is similar to sample No. 08. The morphology of porous structure in the craters near the crack can be clearly seen in the secondary electron SEM image (Fig. 7f ). In addition, a few crystals with higher P 2 O 5 content were also found in the nearby glaze (P3 in Fig. 7d). Figure 7i is a SEM image of a severely corroded area in the glaze, which shows that some yellow substances with lower contrast compared with the glaze are distributed in the pores and cracks of the glaze except for the glaze surface, showing a trend of extending from outside to inside. Figure 7g and h are the OM and SEM images of the corrosion area near a crack. It can be found that the corroded area presenting as yellow (area A) extends along the crack to the inside of the glaze layer. A small number of substances with lower contrast can be found in nearby pores (P4 in Fig. 7h). EDS analysis results show that the MgO and Al 2 O 3 contents are very high, reaching 38.91 wt% and 38.31 wt%, respectively, and the Fe 2 O 3 content is also high ( Table 6). The morphology and composition of the yellow substances will be described in detail in the next chapter.

The morphology and composition of the corrosion products
To identify the phase components of the yellow area on the surface of No. 09, XRD analysis was applied to examine two different areas (area 09-1 with yellow substances and gray area 09-2) of No. 09. The results show that the main crystal phases in the gray area are anorthite    Table 6 (See figure on next page.) and quartz, while diopside, quartz, and quintinite-3 T (Mg 4 Al 2 (CO 3 )(OH) 12 ·3H 2 O) are detected in the yellow area (in Fig. 8). Figure 9a is the OM image of the yellow corrosion area on the glaze surface of No. 09 and Fig. 9b and c are the magnified SEM images. It can be found that the yellow corrosion area is mainly composed of a large number of nubbly anorthite crystals and the MgO-rich phase with low contrast (EDS analysis results are listed in Table 7). The MgO-rich phase filled the interspaces between the crystals, while the glass phase in the original glaze cannot be observed. In another position (Fig. 9d), many interwoven lamellar substances and a large number of densely spaced long rod-like diopside crystals can be found ( Fig. 9e and f ). The SiO 2 content of the lamellar substances is lower than the glaze, while the MgO and Al 2 O 3 contents are very high, reaching 60.98 wt% and 38.13 wt%, respectively (Table 7). Similar morphology can also be observed in the SEM images of the polished cross-section of the yellow area ( Fig. 9g and i). Figure 9i is the SEM image of the interwoven substances in a pore, and its composition has higher Fe 2 O 3 content, reaching 35.23 wt% (Table 7), which is different from the corrosion product on the glaze surface. The MgO/Al 2 O 3 (molar ratio) of the lamellar substances in sample No. 09, calculated using the EDS results, is approximately 2:1, which is consistent with the composition of the quintinite phase, supporting the XRD results.

Corrosion process of celadon
The glaze is a thin glass layer covering the ceramic body [28]. Generally, the glaze is considered to be chemically stable and able to resist most environmental erosion except strong acids and alkalis. The chemical composition of traditional Chinese ancient high-temperature glaze mainly includes SiO 2 , Al 2 O 3 , CaO, MgO, K 2 O, Na 2 O, Fe 2 O 3 , etc., and it can be approximatively represented by SiO 2 -Al 2 O 3 -CaO(MgO)-K 2 O(Na 2 O) glass system. Affected by the chemical composition and the external environment, the corrosion mechanism of the glass presents different characteristics. There are usually two reaction processes, including interdiffusion (or ion exchange) and hydrolysis, but the specific kinetics depends on the glass composition, temperature, pH, composition and renewal rate of the solution, and other factors [29,30]. In the acidic or neutral conditions, ionexchange (or interdiffusion) between hydrogen-bearing species (H + and H 3 O + ) in aqueous solutions and modified ions in glass (such as alkali and alkaline earth metal ions) (Eq. 1 and Eq. 2) leads to selective leaching of glass elements, and a dealkalinised hydrated layer is usually formed on the glass surface [6,[31][32][33][34][35]. As the pH value increases (pH > 9), the hydrated glass network gradually dissolves by hydrolysis of ionic covalent bonds (Si-O-Si) (Eq. 3), and local condensation reactions of insoluble species can occur to form a gel [29,36]. The sea area where the Dalian Island shipwreck is located is close to the mainland, and the average salinity of the seawater is 30.46-31.89‰ [19]. The seawater in Fujian, China is a weakly alkaline environment, with a pH of about 8.2, and the pH value of seawater at different depths has little difference [37]. Therefore, seawater can be approximately regarded as a weakly alkaline salt solution containing Cl − , Na + , Mg 2+ , SO 4 2− , Ca 2+ , K + , etc. In the corrosion process of celadon glaze, the ion exchange seems to be insignificant, and there is no dealkalinised layer formed on the glaze surface, which may be related to the higher pH and salinity of seawater. The dissolution process of the glass network may be the main reaction in the corrosion process, and although the ion exchange reaction also occurs at the same time, there is no obvious feature that can be distinguished.
From the perspective of microstructure, however, in addition to the glass phase, there are a certain number of unmelted quartz grains and pores in the ancient glaze. During the firing process of the high-temperature calcium glaze, crystallization, phase-separation, or phase separation-crystallization structure can be formed, which (1)  [38,39]. The study on ancient Jun glazes buried in the soil environment has shown that the wollastonite crystalline phase in the glaze will be preferentially corroded [40].
The chemical stability of the glass is changed to a certain extent because of the phase-separation structure [41,42], in which the composition of each phase and the shape and distribution of the phase-separation structure are important influencing factors [42]. The phase-separation structure in the typical crystalline phase-separated glaze in ancient China is usually accompanied by the crystallization of anorthite, which is different from the relatively uniform phase-separation structure in glass. Table 6 Chemical compositions of the marked positions in Fig. 7d and i   Three different morphologies of the phase-separation structure are presented due to the inhomogeneity of the glaze composition. When the glass composition is located in the low CaO content region within the liquid-liquid immiscibility boundary above 950 ℃ in the CaO-Al 2 O 3 -SiO 2 equilibrium phase diagram, the discrete droplet phase should be CaO-rich phase, and the glass matrix is the SiO 2 -rich phase [43]. At the same time, other network modifiers such as iron ions and magnesium ions also tend to be enriched in the CaO-rich phase with a large number of non-bridging oxygen atoms (NBO).
The tests revealed that the firing temperature of Nos. 08 and 09 are approximately 1215 ℃ and 1190 ℃, respectively (Fig. 10), which are lower than Longquan celadon wares with proper firing temperature (1220-1280 ℃ [44,45]). The typical representative is sample No. 06 (1265 ℃) excavated from the same shipwreck. Previously the influence of firing temperature and chemical composition on the microstructure of celadon glaze have been discussed, these celadon glazes with lower firing temperatures and high MgO content have a special microstructure different from transparent glazes. In addition to many crystals on the glaze surface, there are also a large number of crystals in the glaze layer as a whole, resulting in opaque and matt glaze. In the anorthite crystallization-phase separation structure in the glaze, the CaO-rich droplets with higher contrast are discretely distributed in the SiO 2 -rich matrix with lower contrast (Figs. 5d and f). Due to the existence of a large number of NBO, the CaO-rich phase has lower corrosion resistance than the SiO 2 -rich phase. On the other hand, the glass phase near the crystal interface is also more likely to be corroded (Fig. 5f), which may be explained by the less tight arrangement of atoms and more vacancies at the crystal-glass phase interface. At the same time, the interstices created by corrosion at the phase interface provide a channel for external solution to enter the glaze layer, thereby promoting the corrosion process. With the corrosion of the interface layer, the combination of the anorthite crystals and the glass phase is gradually destroyed, which causes the shedding of the local crystals. Combining the previous research on the celadon wares with transparent glaze that is well preserved in the same ship, it can be found that the crystallization-phase separation structure with lower chemical stability does significantly affect the corrosion process in the seawater.
For the glaze with the characteristics of crystallizationphase separation structure and the specific marine buried environment, we propose a possible corrosion process in the form of a schematic diagram, as shown in Fig. 11. First of all, ion exchange occurred between the glaze and seawater, giving rise to the minor alteration of the Si-O-Si network. The long-term reaction may lead to Table 7 Chemical compositions of the marked positions in Fig. 9c, e, f, and i (in wt%) microcracks because of the network contraction, which increased the surface available for the reaction [6,8].
The CaO-rich phase-separated droplets are more likely to react with seawater and the phase-separated droplets were gradually corroded, forming a local porous structure. Secondly, the porous structure, increasing the local specific surface area, which also accelerates the corrosion process of the local SiO 2 -rich matrix to a certain extent. At the same time, the glass-crystal interface becomes loose and porous due to the corrosion of the phase-separated droplets. Finally, as the hydrolysis of SiO 2 -rich matrix, the interior CaO-rich phase-separated droplets are gradually exposed. In the process of alternate dissolution of the two separated glass phases, anorthite crystals are also lost and a hole with a size of tens of microns is formed. The above steps are not absolute but change dynamically under different actual conditions. The corrosion process will continue to change dynamically as the local corrosion environment changes [5]. Finally, different types of the secondary phase precipitated due to the reaction between the released cations from glass and the exogenous anions, or the result of intense hydrolysis-reprecipitation [30,46]. The large number of corrosion products appearing in the celadon from the Dalian Island shipwreck also reflects the complexity of the marine burial environment and corrosion process. The formation process of corrosion products will be discussed further in the next chapter.

The formation of corrosion products
The buried environment of celadon wares can be indirectly reflected by the corrosion products in the body and glaze, which is of great significance to deeply understand the influence of the marine environment on the corrosion process. The original appearance of sample No. 09 is significantly changed because of a lot of yellow substances (mainly quintinite phase) on the surface. It is worth noting that the yellow substances with a similar composition was also found in the cracks of sample No. 6 [47][48][49]. Among them, quintinite belongs to the quintinite group minerals, and the divalent ions and trivalent ions are Mg and Al, respectively. Because of the extensive solid-solution between the Mg and Fe end-members in the quintinite group minerals [47], the various compositions of the lamellar crystals in different pores ( Fig. 9e and i) are also related to their formation process.
Since the quintinite phase usually starts to undergo irreversible transformation above 400 °C, decomposing into the periclase and the magnesia-aluminum spinel phase [50][51][52], which means that the substances are impossibly formed during the firing process of the ceramics. Therefore, the more possibility is that the substance was formed during the long-term alteration of ceramics in seawater. As far as we know, no quintinite phase has been found as a corrosion product in ceramics or glass before. It is worth noting that the quintinite phase is formed during the alkali-activation process of stone wools with high Mg content (Si-rich amorphous phase) [53], but the pH value of the solution (NaOH solution, 6.25 mol/L) is very high. A recent study found that the quintinite phase was detected in the ceramic floor tiles (Akko Tower Wreck, Akko Port, Israel, nineteenth century) in a shipwreck under the sea [54], although its source was not specified and it was not directly observed, Fig. 11 Schematic diagram of the corrosion process of crystalline phase-separated structure in the glaze in the marine environment it is probably also due to the long-term alteration. Therefore, under alkaline conditions, materials containing Mg and Al may be transformed into quintinite phases under certain circumstances, although the specific reaction conditions and reaction process still need to be studied.
In our research, during the corrosion process of the crystalline phase-separated structure of No. 09, the corrosion of the phase-separated droplets leads to the appearance of the sub-micron porous structure. These small holes may gradually connect after further corrosion to form porous channels or cracks in the glaze. As seawater enters, a relatively stable micro-environment (relative to the glaze surface contacting with the flowing seawater directly) is formed inside these pores [7,8]. The internal corrosive micro-environment has a higher S/V value (S = exposed surface of glass and V = volume of solution), and the material exchange with the external seawater is relatively slow, which is conducive to the saturation and precipitation of insoluble salts [55]. During the corrosion process, the pH value of the local corrosion environment probably increases, providing alkaline conditions, while Mg 2+ and HCO 3 − ions are provided by the seawater, Al, Si, or Fe are from the dissolution of the glaze, Eq. 4 may occur under certain conditions to form quintinite crystals.
The formation process of corrosion products varies in different regions, which may be related to the space for crystal growth. The space is limited in the small holes of the porous structure left by the corrosion of phase-separated droplets, causing the crystals to not develop into good morphology. These corrosion products fill in the porous channels and gradually wrap the nearby anorthite crystals to form a composite corrosion structure. However, the loose and porous structure cannot completely prevent the ingress of seawater, and seawater could also directly enter the interior of the glaze through the cracks, which further intensifies the corrosion process. While the larger pores in the glaze provide a good growth space. After seawater enters and fills the pores, the ions in the solution are gradually saturated with the corrosion process, leading to the formation of a well-developed lamellar structure. In Fig. 7i we have found that defects such as pores and cracks are usually fast corrosion paths in the glaze layer, which provide spaces with different volumes for the formation of corrosion products. (4) 4Mg 2+ + 2Al 3+ + HCO − 3 + 13OH − + 2H 2 O → Mg 4 Al 2 (OH) 12 (CO 3 ) · 3H 2 O Generally speaking, the corrosion products in the glaze are formed under specific conditions, which are closely related to the marine buried environment and the original glaze composition. The porous structure formed by the corrosion of crystalline phase-separated structure, as well as cracks, pores, and other specific structures in the glaze also play an important role in the formation of corrosion products.

Conclusions
The two Longquan celadon wares from the Dalian Island shipwreck of the Yuan Dynasty have been buried under the seabed for more than 600 years and corroded severely. This paper analyzed the microstructure and corrosion products of the samples, and the following conclusions are obtained: 1. The severely corroded opaque matte glazes belong to the products with lower firing temperatures in Longquan celadon. They have a special anorthite crystallization-phase separation structure with lower chemical stability, which is an important reason affecting the corrosion process of the glaze in the marine environment. 2. The corrosion process of the glaze with anorthite crystallization-phase separation structure can be roughly divided into three stages, including ion exchange, hydrolysis of silicate glass network, and the shedding of the crystals. The CaO-rich phase-separated droplets and the crystal-glass interface layer are preferentially corroded, forming the porous structure, which accelerates the hydrolysis of the SiO 2 -rich matrix. Eventually, the anorthite crystals are also lost and a large-sized hole is formed. 3. The corrosion products in the glaze are mainly quintinite group minerals. The formation process is closely related to the corrosion of the celadon and the external marine environment. The porous structure and cracks in the glaze promote the process, and pores provide enough space for crystal growth.