Scyphozoan jellyfish outbreak events are drawing increasing attentions during the past decade. Elemental compositions of statoliths are helpful to understand jellyfish life history and blooming mechanisms, but very rare endeavor has been focused on the Scyphozoan class. In this work, we explored the feasibility of element analysis of Aurelia aurita (a representative Scyphozoan jellyfish outbreak species in China) which may be used as proxies of environment parameters during jellyfish living and moving. Statolith crystals of Aurelia aurita were found to be a gathering of hexahedron type trigonal needle with size of 10−50 μm long, and 5−10 μm in diameter. By using laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS) downhole profiling method, elements such as Ca, Sr, Mg, Na and P were found to be above the limit of detection and limit of quantification. The epidermis of statocyst could significantly impact the element analyses, so the real statolith element signal section needs to be selected based on elements and Ca profiles together with care. By laser ablated a signal spot repeatedly, the analytical uncertainty was about 3%−4% for Sr/Ca content ratio and Mg/Ca content ratio, but above 10% for other element/Ca content ratios (n=3). Based on the analysis of statolith from temperature-control cultured jellyfish, Sr/Ca content ratios among different statoliths of the same jellyfish were about 6% (n=14), demonstrating biological processes/vital effects causing small variations compared with analytical uncertainties. Therefore, Sr/Ca content ratios may be used as a potential proxy to reveal the living environment variations the Scyphozoan jellyfish has experienced, such as temperature history, which is helpful to understand jellyfish bloom mechanisms.

The small yellow croaker (Larimichthys polyactis) belongs to the family Sciaenidae, which is an offshore warm fish species and widely distributed in the western Pacific. In this study, the variation of genetic diversity and genetic differentiation among L. polyactis populations was analyzed by mitochondrial DNA control region. A total of 110 polymorphic sites were checked, which defined 134 haplotypes. High level of haplotype diversity (h=0.993±0.002) was detected in the examined range. Population genetic structure analyse (analysis of molecular variance, F_st) showed there were high gene flow among L. polyactis populations. The result showed that there were relatively high genetic diversity and low genetic differentiation among the Yellow Sea and the East China Sea populations, which can be attributed to diverse habitats, wide distribution range and high mutation rate of control region. Using phylogenetic methods, coalescent analyses (neutrality tests, mismatch distribution analysis, Bayesian skyline analyses) and molecular dating interpreted in conjunction with paleoclimatic and physiographic evidence, we inferred that the genetic make-up of extant populations of L. polyactis was shaped by Pleistocene environmental impacts on the historical demography of this species. Besides, relatively constant genetic diversity and larger effective population size were detected in recent L. polyactis population. The result showed that the fishing policy certainly, such as the summer closed fishing, played a role in protecting resources of L. polyactis. This study can offer a wealth of biological novelties which indicates genetic structure of L. polyactis population and provides the foundation for resources protection and policy setting.

An ecological optimization project (semi-closed reclamation project) was implemented to control the invasion of Spartina alterniflora, and optimize the habitat of the Chongming Dongtan wetland, in the Changjiang (Yangtze) River Estuary. After project implementation, a macrobenthic ecological survey was conducted in a natural tidal flat and a semi-closed reclamation restoration area within the Chongming Dongtan wetland from 2019 to 2020. Compared with historical data before reclamation, findings showed that the groups, numbers, and species diversity of the macrobenthos increased significantly, and the ecological optimization project resulted in good ecological benefits. In addition, compared to the natural tidal flat, the number of collected macrobenthic phyla, and the macrobenthic density and biomass were significantly lower in the restoration area. Furthermore, the biodiversity index and functional redundancy of natural tidal flats were generally higher, indicating that the community composition and function of natural tidal flats were relatively more stable. Even though the species composition differed between a number of restoration areas and natural tidal flats, there was no difference in functional diversity, indicating that the effect of restoring ecological functions in restoration areas was optimal. Among them, the biodiversity and functional redundancy of Site S2 were significantly reduced, and the ecosystem function was extremely unstable. Habitat heterogeneity, vegetation community and decreasing salinity were the main factors that affected the ecological functions of macrobenthos. The ecological quality was also evaluated; the Transects N3 and N4 showed good quality. The overall ecological quality of the restoration area was generally high, but that of Site S2 was poor and that of Site S2E was merely good, which was mainly due to modifications of the ecological function of macrobenthos. It is suggested that reeds mowing and freshwater species release should be adopted in restoration areas to improve the community function and the environmental disturbance resistance of the macrobenthos.

Since 2015, green tides with Ulva prolifera as the dominant species in the Qinhuangdao coastal waters have continued to occur. In this study, the relationship between green tides in Qinhuangdao and the Yellow Sea (setting sites in Rudong and Qingdao) was evaluated by genetic analyses of U. prolifera. Single nucleotide polymorphism (SNP) markers were used to analyze genetic diversity and genetic relationships among groups. Genetic differentiation was lower among floating U. prolifera populations in Rudong and Qingdao than in Qinhuangdao. The floating U. prolifera population had higher genetic diversity and polymorphism levels in Qingdao and Rudong than in Qinhuangdao. Physiological experiments showed that the growth rate and net buoyancy of floating U. prolifera were highest in Qinhuangdao and Qingdao, respectively, under the same environmental conditions (temperature and light). Overall, these findings showed that U. prolifera populations in the Qinhuangdao and Yellow Sea green tides (Rudong and Qingdao) differ significantly at the molecular and physiological levels. Therefore, the Qinhuangdao green tide is not correlated with the Yellow Sea green tide and has a different origin and development mode. This study provides insight into the mechanism underlying green tide blooms in coastal waters of China.

The productivity and health of seagrass depend on the combined inputs of nutrients from the water and sediments in which they grow and the microbiota with which they live intimately. However, little is known about the composition and diversity pattern of single-celled benthic eukaryotes in seagrass meadows. Here, we investigated how the structure and diversity of the benthic microeukaryotic community vary with respect to season, location, and seagrass colonization, by applying 18S rRNA gene amplicon sequencing for 96 surface sediment samples that were collected from three different seagrass habitats through four seasons. We found that benthic microeukaryotic communities associated with seagrass Zostera japonica exhibited remarkable spatial and seasonal variations, as well as differences between vegetated and unvegetated sediments. Diatoms and dinoflagellates predominated in the benthic microeukaryotic communities, but they were inversely correlated and displaced each other as the dominant microbial group in different seasons or habitats. Mucoromycota was more prevalent in vegetated sediments, whereas Lobulomycetales and Chytridiales had higher proportions in unvegetated sites. Total organic carbon and total organic nitrogen were the most important environmental factors in driving the microeukaryotic assemblages and diversity. Our study expands the available knowledge on the biogeographic distribution patterns and niche preferences for benthic microeukaryotes in seagrass systems.

The impacts of the Hong Kong-Zhuhai-Macao Bridge (HKZMB) on suspended sediment content (SSC) were analysed in the Zhujiang River Estuary based on data from HY-1C, which was launched in September 2018 in China, carrying Coastal Zone Imager (CZI) and Chinese Ocean Color and Temperature Scanner on it. A new SSC inversion model was established based on the relationship between in-situ SSC and the remote sensing reflectance in red and near-infrared bands of CZI image. HY-1C satellite data obtained from October to December 2019 were applied to retrieve SSC in the Zhujiang River Estuary. The results show that SSC around the HKZMB is ranging from 20 mg/L to 95 mg/L. SSC change obviously on two sides of the bridge. During flooding and ebbing period, SSC increases obviously downstream of the bridge. SSC difference between upstream and downstream is ranging from 5 mg/L to 20 mg/L. Currents flowing across the HKZMB, the change trend of SSC in most places upstream and downstream is almost the same that SSC downstream of the bridge is higher than SSC upstream. The tidal currents interact with bridge piers, inducing vortexes downstream, leading the sediment to re-suspend downstream of the bridge piers. Other factors, including seafloor topography and wind, can also contribute to the distribution of SSC in the Zhujiang River Estuary.

The wave-induced setup and circulation in a two dimensional horizontal (2DH) reef-lagoon-channel system is investigated by a non-hydrostatic model. The simulated results agree well with observations from the laboratory experiments, revealing that the model is valid in simulating wave transformation and currents over reefs. The effects of incident wave height, period, and reef flat water depth on the mean sea level and wave-driven currents are examined. Results show that the distributions of mean sea level and current velocities on the reef flat adjacent to the channel vary significantly from those in the area close to the side walls. From the wave averaged current field, an obvious alongshore flux flowing from the reef flat to the channel is captured. The flux from the reef flat composes the second source of the offshore rip current, while the first source is from the lagoon. A detailed momentum balance analysis shows that the alongshore current is mainly induced by the pressure gradient between the reef flat and the channel. In the lagoon, the momentum balances are between the pressure and radiation stress gradient, which drives flow towards the channel. Along the channel, the offshore current is mainly driven by the pressure gradient.

We report field measurements of vertical profiles of the turbulent diffusivity and temperature at different stations in the South China Sea (SCS). Our study shows that the measured turbulent diffusivity follows a power-law distribution with a varying exponent in water layers. Similar multiple-layer scaling regimes were also observed from the temperature fluctuations. Combining turbulent diffusivity and temperature fluctuations, the vertical structure of temperature was revealed. Furthermore, we discussed the temperature profiles in each layer. A constant function of a dimensionless temperature profile was found in water layers that have identical turbulence conditions. Our results reveal the multiple-layer structure of temperature in the SCS. This study contributes to the understanding of the vertical structure of multiple layers in the SCS and provides clues for exploring the physical mechanism for maintaining the temperature structure.

We introduced the Coupled Model Intercomparison Project Phase 6 (CMIP6) Ocean Model Intercomparison Project CORE2-forced (OMIP-1) experiment by using the First Institute of Oceanography Earth System Model version 2.0 (FIO-ESM v2.0), and comprehensively evaluated the simulation results. Unlike other OMIP models, FIO-ESM v2.0 includes a coupled ocean surface wave component model that takes into account non-breaking surface wave-induced vertical mixing in the ocean and effect of surface wave Stokes drift on air-sea momentum and heat fluxes in the climate system. A sub-layer sea surface temperature (SST) diurnal cycle parameterization was also employed to take into account effect of SST diurnal cycle on air-sea heat fluxes to improve simulations of air-sea interactions. Evaluations show that mean values and long-term trends of significant wave height were adequately reproduced in the FIO-ESM v2.0 OMIP-1 simulations, and there is a reasonable fit between the SST diurnal cycle obtained from in situ observations and that parameterized by FIO-ESM v2.0. Evaluations of model drift, temperature, salinity, mixed layer depth, and the Atlantic Meridional Overturning Circulation show that the model performs well in the FIO-ESM v2.0 OMIP-1 simulation. However, the summer sea ice extent of the Arctic and Antarctic is underestimated.

In the Xiangshan Bay at the east coast of China, coastal marine pollution is conspicuous and severe in recent years. As transport of the pollutants is closely related to the coastal circulation, there is a great practical significance to investigate the circulation in this area. In this work, the surface pattern and vertical profiles of Lagrangian residual velocity (LRV) were studied based on field observation data from the inner Xiangshan Bay. By tracking GPS-GPRS drifters’ trajectories, the surface LRV pattern is going out in the central deep trough and flowing inwards near the shoreside. Combined with data from two mooring stations, vertical profiles of LRV is flowing out at surface and flowing in at the bottom, consistent with the gravitational circulation induced by baroclinic effects at the estuary. However, according to the diagnostic analysis, the main mechanism driving the residual current is barotropic rather than baroclinic. The LRV equation is controlled by the tidally-averaged barotropic pressure gradient force, tidal body force and tidally-averaged turbulent stress, while the tidally-averaged baroclinic pressure gradient force is one order of magnitude less than other forces. Additionally, the tidally mean eddy viscosity coefficient which is used in the expression of tidally-averaged turbulent stress might be not adequate and requires further studies.