Coastal China Research Articles

Unforeseen green tide of floating tubular Ulva meridionalis, a lethal threat to oyster farming, along the west coast of Taiwan: tracking its origin and ecophysiological insights

Abstract Ulva green tides have adversely affected coastal ecosystems. In June 2023 in Changhua County, Taiwan, an unprecedented floating tubular Ulva bloom spanning about 30 km of coastline caused about 30 % reduction in oyster yield as reported by local farmers. Understanding its taxonomic and ecological basis is crucial for preemptive and remedial measures. Based on molecular (rbcL and tufA) and morphological evidence, U lva meridionalis was the cause of this green tide. Haplotype network analysis (based on ITS) suggests that this green tide originated from northern China via the China Coastal Current. Historical survey data indicate that U. meridionalis arrived in Taoyuan Algal Reef (about 150 km north of the bloom area) as far back as 2018. Our ecophysiological experiments revealed that U. meridionalis exhibited a rapid daily growth rate with biomass increment up to 13–21 % when subjected to local nutrient-rich waters under lower salinity and spring-like conditions. Although historical ecological and poultry/livestock data analyses showed no noticeable change in sea surface temperature and rainfall over the past decade, a gradual rise in agricultural nitrogen and phosphorus output was observed. This taxonomic and ecological background lays the groundwork for long-term ecological monitoring. Moreover, this study exemplifies the detrimental impact of an unforeseen Ulva bloom on oyster farming.

Effects of the comprehensive elimination of Spartina alterniflora along China's coast on blue carbon and scenario prediction after ecological restoration

Salt marshes cover the largest area among the three types of traditional blue carbon ecosystems in China's coastal zone, with the introduced smooth cordgrass (Spartina alterniflora Loisel.) being dominant in these marshes. The effects of eradicating S. alterniflora nationwide and the subsequent ecological restoration on blue carbon are unclear. This paper evaluates the variation in blue carbon during the national S. alterniflora eradication campaign, which involves mechanical tillage from 2022 to 2025, and proposes three scenarios for blue carbon changes after native vegetation is reestablished by 2050. The results show that, in 2025, plant carbon stock and soil carbon stock will decrease by 1.38 Tg C and 1.21 Tg C, respectively, in the areas where S. alterniflora has been removed and managed. Although blue carbon is reduced in coastal wetlands in 2025, carbon stock is expected to increase in restored native vegetated wetlands by 2050. S. alterniflora is resilient and competitive, posing a high risk in secondary invasion. Scenario Ⅰ suggests that S. alterniflora marshes could almost recover to their original state from 2022, with 7.70 Tg C stored in plant and soil carbon stocks. Scenario Ⅱ involves native vegetated wetlands coexisting with invasive S. alterniflora marshlands, with a total carbon stock estimated at 7.15 Tg C, reflecting a decrease of 0.39 Tg C in soil carbon stock and by 0.16 Tg C in plant carbon stock. In Scenario Ⅲ, mudflats dominant and native vegetated habitats are reestablished only in suitable sites, with the total carbon stock estimated at 5.63 Tg C, a 26.9% decrease compared to 2022 levels. While eradicating invasive S. alterniflora and restoring native vegetation in China's coast enhance the ecosystem services, it reduces blue carbon stocks. Therefore, developing additional strategies to increase carbon storage in coastal wetlands is necessary.

China coasts facing more tropical cyclone risks during the second decaying summer of double-year La Niña events

Long-lasting La Niña events (including double-year and triple-year La Niña events) have become more frequent in recent years. How the multi-year La Niña events affect tropical cyclone (TC) activities in the western North Pacific (WNP) and whether they differ from single-year La Niña events are unknown. Here we show that TCs are more active over the far-WNP (FWNP, 110°–150°E), leading to marked high risks at China coasts during the second decaying summer of double-year La Niña events. The anomalous TC activities are directly related to the enhanced cyclonic anomaly over the FWNP, possibly a result of large-scale remote forcing initiated by the tropical North Atlantic (TNA) cooling. The persistent TNA cooling from the decaying winter to summer of double-year La Niña events drives westerlies over the Indo-western Pacific through Kelvin waves, which induce the cooling over the north Indian Ocean via the wind-evaporation-sea surface temperature effect, favoring the asymmetric heat distribution pattern and stimulating an anomalous vertical circulation over the eastern Indian Ocean to FWNP. The cooling over the north Indian Ocean also excites Gill responses, magnifying the TNA-induced westerlies and boosting the anomalous vertical circulation, and thus gives rise to the strong cyclonic circulation anomaly over the FWNP in summer. We suggest that the key point of the process is the strong TNA cooling related to the persistent negative Pacific-North American pattern (PNA) and positive North Atlantic Oscillation (NAO) while double-year La Niña events decay, distinct from the rapid decline of PNA and NAO during single-year La Niña events. The work provides a unique perspective on understanding TC activities over the WNP related to the El Niño-Southern Oscillation.

Modeling the spatiotemporal dynamics of electric power consumption in China from 2000 to 2020 based on multisource remote sensing data and machine learning

China's rising electricity power consumption (EPC) is strongly correlated with carbon emissions. Timely and accurate analysis of the spatiotemporal dynamics of EPC is important to the realization of carbon peaking and carbon neutrality goals in China. However, due to data quality problems and limitations of modeling methods, the estimation accuracy warrants further improvement. In this study, the EPC in China from 2000 to 2020 was estimated based on multisource remote sensing data using the Random Forest (RF) model. Compared with previous studies, the accuracy of this study was improved by 39%–47%. The reasons were combining multisource remote sensing data can mitigate the quality issues of nighttime light (NTL) data, and the RF can capture the nonlinear relations between remote sensing data and EPC. In addition, the spatial pattern of the average EPC in China was dominated by the low-level EPC, as well as showing an obvious increasing trend. We also found that in the middle reaches of the Yellow River and northern coastal China, the low-speed increase in EPC led to high carbon emissions and emission intensity. We suggest optimizing the fuel fix of energy and adjusting the industrial structure, combining them with scientific and rational spatial planning.

Characterization and toxicological effects of microplastics in PM2.5 in background sites of the southern coast of China

Microplastics (MPs) are widely present in atmosphere without understanding about their influences on human health. Spatial variations of MPs in fine particulate matter (PM2.5), along with source contribution and health effects analyses, were investigated at two distinct background air quality locations in Hong Kong: TMS, a coastal site, and HT, a high-altitude site, both remote from anthropogenic activities. The higher concentration of total MPs concentration (∑MPs) was observed at HT (31.67 ± 10.25 pg/m3) compared to TMS (11.00 ± 4.69 pg/m3). Polyethylene (PE) was the dominant MP in both locations, and other MPs showed inconsistent concentration patterns in the analysis. Source identification analysis indicated that the MPs at HT was predominantly affected by marine sources, whereas the MPs at TMS was potentially due to continental/local sources. The HT was not an ideal location for monitoring atmospheric background concentration of MP due to high contributions from marine source. Cellular exposure to MPs experiments (A549 cells) demonstrated stronger correlations (p < 0.05) between exposure to MPs and oxidative stress (determined by lactate dehydrogenase (LDH) and malondialdehyde (MDA)) than inflammation factors. The polyvinyl chloride, nylon 6,6 and polystyrene were the chemical compounds that can induce oxidative stress, and interactions among MPs. The results of this study will aid in the development of a new MPs control policy.

University-related Factors in Technology Integration: A Quantitative Study in East Coast of China

University-related Factors in Technology Integration: A Quantitative Study in East Coast of China

Plutonium in the coastal seawater around Chinese nuclear power plants: Sources, distribution, and environmental implications

Coastal surface seawater samples within 30 km around ten Chinese nuclear power plants (NPPs) were systematically investigated. The 239+240Pu activity concentration in the samples varied from 0.226 mBq/m3 to 3.098 mBq/m3, meanwhile the 240Pu/239Pu atom ratios ranged from 0.151 to 0.353. Besides, the Pacific Proving Ground (PPG) close-in fallout and the global fallout were recognized as two primary sources of Pu in these samples. The 239+240Pu activity concentration as well as the PPG contribution showed similar trends as the Kuroshio intrusion path and the coastal currents in the China Seas, illustrating long-range transport and consuming of PPG derived Pu in the coastal China Seas. Moreover, accumulation of PPG sourced Pu in the Beibu Gulf were observed and was attributed to the continuous invasion of the high isotopic Pu that remobilized from the South China Sea (SCS).

Separation and spatial variations of typhoon and non‐typhoon rainfall at different timescales in typical region of southeast China

AbstractRainfall in East Asia is affected by two rain‐bearing systems: tropical cyclones and monsoon‐related frontal systems. Distinguishing typhoon rainfall (TR) and non‐typhoon rainfall (non‐TR) helps to understand the evolution process of regional rainfall at different timescales. Taking Fujian Province in the southeast coast of China as an example, based on the fixed box approach of separating TR, the method of determining the size of fixed box is explored. TR and non‐TR are separated in Fujian Province, and the spatial variations of TR and non‐TR at different timescales (annual, monthly, day of annual‐maximum‐rain) are analysed. The results showed that (1) according to the relationship between the sizes of fixed box and the rate of change of TR, the size of fixed box could be reasonably determined; thus, the expended size of fixed box suitable for separating TR in Fujian was 3.5°, namely the range of 20°–31.8°N and 112.3°–124.2°E. (2) The spatial variations of TR at different timescales in Fujian were similar: TR decreased from the coast to the inland, and the northeast coast of Fujian was the high‐value region. Due to the difference of water vapour sources, non‐TR in March–June increased from the coast to the inland, but the high value of non‐TR in July–September was distributed in the eastern coast. (3) The contribution rates of average TR to the total rainfall of the year, July–September and 1 day were 12.8%, 34.6% and 35.7%, respectively. In eastern coast, TR in July–September accounted for 1/3–1/2 of total rainfall, and the rainfall in a day was mainly affected by TR; while in western inland, rainfall was mainly non‐TR and the influence of TR was less than 1/4.

Whole-Genome Resequencing Reveals Signatures of Adaptive Evolution in Acanthopagrus latus and Rhabdosargus sarba.

Acanthopagrus latus and Rhabdosargus sarba are economically important marine species along the coast of China, with similar external morphological characteristics and living habits, with wide distribution and strong adaptability. To investigate the molecular mechanisms underlying the adaptive evolution of these two species, we conducted whole-genome resequencing of 10 individuals of both species from the coastal waters of Wuyu Island, Fujian, China, using high-throughput sequencing technology. We obtained SNP, InDel, CNV, and SV variation information and annotated these variations, constructing a genomic variation database for both species. By comparing the resequencing data with reference genomes, we identified 9,829,511 SNP loci in the population of A. latus and 34,051,056 SNP loci in the population of R. sarba. Using whole-genome SNP data, we employed Fst and ROD methods to identify candidate genomic regions under selection. Functional annotation and enrichment analysis using GO and KEGG databases revealed potential adaptive evolution in R. sarba associated with immune response, feeding, growth and development, and locomotion, while A. latus showed potential adaptive evolution associated with immune response, nervous system, growth and development, and metabolism.

Optimizing entropy weight model to accurately predict variation of pharmaceuticals and personal care products in the estuaries and coasts of the South China

Pharmaceuticals and personal care products (PPCPs) have raised increasing concern worldwide due to their continuous release and potential hazards to the ecosystem and human health. This study optimized the entropy weight model (EW-WRSR) that combines entropy weight with multi-criteria decision analysis to investigate pollution patterns of PPCPs in the coasts and estuaries. The results revealed that occurrences of PPCPs from the 1940s to the present were consistent with using PPCPs, different types of human activities, and local urban development. This helped better understand the history of PPCP contamination and evaluate the uncertainty of EW-WRSR. The model predicted hotspots of PPCPs that were consistent with the actual situation, indicating that PPCPs mainly enter the nearshore ecosystem by the form of sewage discharge and residual aquaculture. This study can provide method that identifying highly contaminated regions on a global scale.

Decoding the RNA viromes in shrew lungs along the eastern coast of China

Shrews being insectivores, serve as natural reservoirs for a wide array of zoonotic viruses, including the recently discovered Langya henipavirus (LayV) in China in 2018. It is crucial to understand the shrew-associated virome, viral diversity, and new viruses. In the current study, we conducted high-throughput sequencing on lung samples obtained from 398 shrews captured along the eastern coast of China, and characterized the high-depth virome of 6 common shrew species (Anourosorex squamipes, Crocidura lasiura, Crocidura shantungensis, Crocidura tanakae, Sorex caecutiens, and Suncus murinus). Our analysis revealed numerous shrew-associated viruses comprising 54 known viruses and 72 new viruses that significantly enhance our understanding of mammalian viruses. Notably, 34 identified viruses possess spillover-risk potential and six were human pathogenic viruses: LayV, influenza A virus (H5N6), rotavirus A, rabies virus, avian paramyxovirus 1, and rat hepatitis E virus. Moreover, ten previously unreported viruses in China were discovered, six among them have spillover-risk potential. Additionally, all 54 known viruses and 12 new viruses had the ability to cross species boundaries. Our data underscore the diversity of shrew-associated viruses and provide a foundation for further studies into tracing and predicting emerging infectious diseases originated from shrews.

The characteristic of HLA-A, HLA-B, HLA-C, HLA-DRB1, HLA-DRB3/4/5, HLA-DQA1, HLA-DQB1, HLA-DPA1, and HLA-DPB1 alleles in Zhejiang Han population.

The Zhejiang Han population, a subgroup of the Southern Han ethnic group, resides in Zhejiang Province, situated on the southeast coast of China. In this study, we conducted HLA genotyping for 813 voluntary umbilical cord blood donors from the Zhejiang Han population, targeting 11 HLA loci, namely HLA-A, HLA-B, HLA-C, HLA-DRB1, HLA-DRB3/4/5, HLA-DQA1, HLA-DQB1, HLA-DPA1, and HLA-DPB1, using the next-generation sequencing method. Our analysis of the alleles and haplotypes revealed a high degree of polymorphism within these loci. A total of 289 unique HLA alleles were identified, with the HLA-B locus exhibiting the most significant diversity, while HLA-DRB4 displayed the lowest variation. Due to the inherent limitations of the sequencing method, some unresolvable alleles in the specific loci, such as HLA-DRB1, HLA-DPA1, and HLA-DPB1, were assigned as G group designation. In our comprehensive analysis across all 11 HLA loci, a total of 1204 haplotypes were estimated. The distribution of these alleles was similar to those of the Chinese Southern Han population while highly different from the Caucasian population. These findings contribute to a deeper understanding of the genetic characteristics of HLA loci within the Chinese Southern Han population.

Role of Topographic Forcing in the Maritime Continent Region in Interaction Between the East Asian Summer Monsoon and the Australian Winter Monsoon

AbstractThe topography in the Maritime Continent (MC) has significant impact on climate anomalies in the Asian‐Australian monsoons region. In the present study, the Regional Climate Model Version 4.6 (RegCM4.6) is applied for the simulation of climate over the Asian‐Australian monsoon region during the boreal summer. Results demonstrate that the RegcM4.6 is able to well reproduce precipitation, temperature and low‐level and upper‐level circulation patterns over the Asian‐Australian monsoons region. A sensitivity experiment with zero topographic height in the MC region shows that the intensity of western Pacific subtropical high (WPSH) and Australian cold high both weaken simultaneously, while the cross‐equatorial flows also decline and the East Asian‐Australian monsoons become weaker. Meanwhile, the anomalous cyclonic circulation with significant convergence prevails in lower levels over the western Pacific, leading to more precipitation and higher temperature. In the MC region, there are more precipitation and high temperature in the north while there are less precipitation and low temperature in the south. Temperature increases over a large area from the Yunnan‐Guizhou Plateau to the Loess Plateau but decreases in the southeastern coast of China and eastern India. These results have important implications for better understanding the topographic impact of the MC region on the interaction between the east Asian‐Australian monsoons.

Clinicopathologic features and new factors affecting constitution of eyelid neoplasms: a 10-year multicentre experience in southern coastal China

ObjectiveThis study aims to delineate the key characteristics of eyelid neoplasms in central and southern coastal China and to explore potential factors affecting tumour development. MethodsThis was a retrospective, multicentre observational study. In total, 1302 cases of eyelid neoplasms diagnosed in 2 tertiary hospitals were reviewed from January 2013 to May 2023. Histogenesis, pathologic diagnosis, and potential risk factors were investigated, and the findings were compared with data from various regions and countries. Age- and sex-specific incidence rates, time trends, distributions across age groups, and the relationship with the Sociodemographic Index, latitude, and altitude were calculated. ResultsThere were 1162 benign (89.2%) and 140 malignant (10.8%) cases, representing 63 pathologic types. The median age of benign cases (50 years) was significantly lower than that of malignant cases (72 years) (P < 0.05). The most common benign and malignant lesions were intradermal nevus (28.83%) and basal cell carcinoma (BCC) (74.29%). A relatively younger median age (68 years) for sebaceous gland carcinoma (SGC) was observed, with a male preponderance (sex ratio 2:1). Analyses of the proportions of BCC, SGC, and squamous cell carcinoma in eyelid malignancies revealed significant regression equations (P < 0.05). ConclusionsMost eyelid neoplasms are of epithelial and melanocytic origin, with benign tumours being predominant. The regions with greater Sociodemographic Index exhibit a greater incidence of BCC and a lower incidence of squamous cell carcinoma and SGC. BCC appears to have a greater incidence in regions located at higher latitudes.

Sustainable urban planning to control flood exposure in the coastal zones of China

ContextSustainable development in coastal zones faces escalating flood risk in the context of climate change and urbanization, and the rapid urban growth in flood zones has been one of the key drivers. Therefore, understanding the Urban Exposure to Flooding (UEF) and its future scenarios is important in coastal zones.ObjectivesThe objectives of this study were: (1) to assess the future dynamics of UEFs in China's coastal zones, and (2) to identify a sustainable way of urban planning in controlling the growth of UEFs.MethodsFuture UEFs in coastal China were assessed during 2020–2050 by combining urban expansion model, scenario analysis, and flood exposure assessment. Alternative scenarios were considered of shared socioeconomic pathways (SSPs), representative concentration pathways (RCPs), strategies of urban planning.ResultsThe results show that the 1000-year flood UEFs along coastal China was expected to grow under SSP2-RCP4.5 from 9,879 km2 in 2020 to 13,424 (12,997–13,981) km2 in 2050, representing an increment of 35.88% (31.56%–41.52%). Alternatively, the strategy of sustainable development planning could reduce the newly added UEF by 16.98% (15.63%–18.67%) in a 1000-year flood scenario.ConclusionsThe findings proved that the ways of urban growth matters in terms of affecting food exposure and risk and flood risk should be incorporated into urban planning for a sustainable landscape. The study could offer methodology and support for sustainable development strategies in reducing future urban flood risk.

Spatial and temporal characteristics of warm season convective storms over the southeastern coast of China

Utilizing observational data from the Xiamen radar station during the warm season of 2015–2018 (May to September), this study employs a convective storm identification algorithm to statistically analyze the spatiotemporal distribution, diurnal propagation, and seasonal variability of convective storms over the southeastern coast of China. Key findings include: (1) Significant monthly variations in convective storm frequency, with peaks in August over land areas of Zhangzhou and Xiamen's northwest, and offshore southeast of Xiamen. Seasonal circulation patterns, particularly the subtropical high and northern high-level troughs, drive these variations. (2) Large-scale convective storms are most frequent, while small-scale ones are less common. Mid-deep convective storms dominate, particularly in southern Zhangzhou and southwestern mountainous Quanzhou, whereas shallow convective storms are rare and scattered. High-frequency areas correlate with higher terrain, underscoring the influence of topography on storm occurrence and development. (3) Hovmöller plots reveal a bimodal diurnal pattern in propagation of storms for July and August, with peaks in the daytime and late night. Daytime storms propagate from coastal to higher terrain areas, while nighttime storms maybe driven by enhanced vertical wind shear. These findings enhance the understanding of convective storms in the region and highlight the crucial role of the circulation background, terrain, and prevailing wind directions in the spatiotemporal characteristics of warm season convective storms in southeastern China.

Development and validation of a 40 K liquid SNP array for the mud crab (Scylla paramamosain)

Genomic selection (GS) and genome-wide association studies (GWAS) have been widely implemented to accelerate the genetic improvement of important economic traits. Developing an efficient and reliable high-throughput genotyping tool is the premise and basis of implementing GS and GWAS on the mud crab (Scylla paramamosain). In this study, leveraging the whole-genome resequencing data of 146 mud crabs from the southeast coast of China, a total of 16,714,943 high-quality SNPs were discovered. Based on polymorphic information content and evenly distribution on 49 chromosomes, a total of 40,352 SNPs was selected to develop a 40 K liquid SNP array (named “Xiexin No. 1”) using the genotyping by target sequencing (GBTS) technology. The genotyping performance test showed that this SNP array had a high site polymorphism (over 92.70%) and individual calling rate (97.02%). Furthermore, the mean of SNP concordance, genotype concordance, non-reference sensitivity, and non-reference discrepancy between “Xiexin No. 1” SNP array and resequencing were identified to be 0.986, 0.954, 0.942, and 0.097, respectively. Using this SNP array data, the genomic heritabilities were evaluated as 0.887 ± 0.065, 0.643 ± 0.108, 0.599 ± 0.111, and 0.624 ± 0.111 for body weight (BW), carapace length (CL), carapace length (CW), and body height (BH) in a cultured population (233 crabs). Furthermore, five significant SNPs associated with growth traits were identified by GWAS. After annotation, one candidate gene (XP_027221997.1) associated with BW and one candidate gene (XP_027210446.1) associated with CW and BH were identified. Most importantly, the prediction accuracy of GBLUP for BW, BH, CL, and CW were 0.535 ± 0.064, 0.544 ± 0.101, 0.574 ± 0.087, and 0.585 ± 0.083, with the unbiasedness values of 0.988 ± 0.22, 1.015 ± 0.274, 1.011 ± 0.196 and 1.037 ± 0.235. In conclusion, the “Xiexin No. 1” SNP array provides an efficient and reliable platform for high-quality genotyping, which would greatly promote the genetic studies and thus accelerate the genetic improvement of mud crabs.

Exploring the feasibility of multi-elements coupled with chemometrics for discriminating the geographical origins of oysters (Crassostrea ariakensis)

This study explored the efficacy of multi-elements combined with chemometrics to discriminate the geographical origins of oysters (Crassostrea ariakensi). We determined 52 elements in 166 samples from four regions along the southeast coast of China. Significant regional variations of 51 elements were revealed (P < 0.05), while the principal component analysis (PCA) provided no clear regional delineations. The training models (n = 117) established on linear discriminant analysis (LDA), partial least square discriminant analysis (PLS-DA), and random forest (RF) uniformly achieved 100% predictive accuracy. The cross-validation accuracies of the final models (n = 166) derived from LDA, PLS-DA, and RF were 100%, 100%, and 99.4%, respectively. Even with the models simplified to 8 elements (Zn, Al, K, Cd, Cu, Rb, B, and Ag), high predictive and cross-validation accuracies were maintained, underscoring the robustness and algorithm flexibility of elemental profiling for accurately identifying the geographical origins of oysters.

A comprehensive assessment of the importance of the “Maritime Silk Road” route based on multi-dimensional data-driven approaches

Worldwide straits are crucial pathways for global trade and energy transportation, occupying key strategic positions. Therefore, this study utilizes automatic identification system (AIS) data and the structure of the maritime transportation network formed by ship navigation, combined with geospatial big data analysis, to fully consider factors including the natural navigational conditions of routes, the support capabilities of surrounding ports, the geographic advantages of neighboring countries, and the distribution of the shipping network importance in the CCAS evaluation model. To assess the navigability of maritime routes, a data-driven model was constructed using three layers of this model. Empirical analyses were conducted on 22 major routes along China's coast and the “Maritime Silk Road”. This method has significantly enhanced the accuracy and reliability of route assessments, and the improved assessment model demonstrated better navigational conditions in the South China Sea routes, identifying internal factors hindering the development of these critical maritime routes. The findings of this study provide a comprehensive assessment and analysis framework for global maritime transportation, offering crucial reference information for key stakeholders, including policy makers, port and berth infrastructure investors, maritime economists, and large shipping companies. Additionally, the study's results can be practically applied in the actual shipping trade, emphasizing the role of our framework in facilitating effective decision-making and strategic planning within the maritime industry.

Impact of Marine Shipping Emissions on Ozone Pollution During the Warm Seasons in China

AbstractAs China's land‐based anthropogenic emissions are decreasing, the impact of marine shipping emissions (MSEs) on the atmosphere, especially in coastal areas, deserves further attention. This study investigates the impact of MSEs on MDA8 ozone (O3) levels during the warm seasons of 2017 in China, considering different seasons and synoptic patterns. The results indicate that the average impact of MSEs on O3 decreases from offshore to inland, peaking at over 29.0 ppb at sea and 13.8 ppb along the coast of mainland China. Influenced by precursor emissions, meteorology and other factors, MSEs contribute differently to O3 in Beijing‐Tianjin‐Hebei (BTH), Yangtze River Delta (YRD) and Pearl River Delta (PRD), with contributions of 3.0, 5.2, and 4.9 ppb, respectively, and ranging from 2.7 to 7.3 ppb in 13 coastal port cities. The O3 impacts of MSEs are higher on polluted days than on clean days, especially during onshore winds. In the BTH, MSEs increase O3 by 5.5 ppb on polluted days and 3.0 ppb on clean days with northeast winds from the Bohai Sea. In the YRD, MSEs increase O3 by 9.4 ppb on polluted days and 7.3 ppb on clean days with southeast winds. MSEs significantly increase O3 levels in the PRD by 11.0 ppb on polluted days and 5.0 ppb on clean days with southeast winds. Although the emission inventories, initial and boundary conditions, etc. may introduce uncertainties, our results still provide useful information for O3 pollution management in coastal cities as a reasonable way to track mass contributions.