Abstract

Using ASCAT, QuikSCAT, and MODIS data, we analyzed the sea surface wind field, temperature, salinity, and chlorophyll concentrations in the mixed zone between the Bohai Sea and Yellow Sea in the winter (the period of winter 2013 included December 2013 and January-February 2014) from 2002 to 2013. We found that the intrusion of the Yellow Sea Warm Current into the Bohai Sea occurred three times in the winters of 2007 (strongest), 2004, and 2013 (weakest) during this 12-year period. We present detailed validation of the intrusion in 2013. This study shows that the intrusion of the Yellow Sea Warm Current into the Bohai Sea occurred when the wind speed, sea surface temperature, and salinity were above (or close to) the multiyear average and the chlorophyll concentration was less than the multiyear average.

1. Introduction

The Yellow Sea Warm Current (YSWC) is a branch of the Tsushima Current, which both branch from the Kuroshio Current [13]. The YSWC originates in the area northeast of the East China Sea and south of the Jeju Island of Korea. It is the only source of sea water from the open sea entering the Yellow Sea and is characterized by high temperature, high salinity (SAL), and low dissolved oxygen content [4].

There are many important issues surrounding the YSWC. For example, existing high-resolution numerical simulations suggest that the path and intensity of the YSWC are altered by the wind [59]. Yuan et al. [8] showed that winter northerly winds could drive the YSWC north to flow across the Yellow Sea trough. Tang et al. [7] showed that the path of the YSWC had clear seasonal and interannual variations. Numerical calculations showed that the variations were closely related to the strength of the northerly wind. Model simulation by Song et al. [4] revealed clear 5-year variability in the YSWC and found a south-north shifting that correlated with the zonal and meridional wind variability. In a strong wind situation, the YSWC intrudes further into the northern Yellow Sea. Using model simulation, Miu and Liu [10] found a northwestward flow along the northern Yellow Sea trough and assumed that it was the YSWC. The flow moves along the Yellow Sea trough against the wind, moving up from the deep area at the north of the Bohai Strait and extending to the central Bohai Sea. Shen and Mao [11] found that sea surface SAL of the eastern Bohai Sea reached a maximum in February and March every year.

From the review above, we found that it is not conclusive how and when the YSWC path enters the Bohai Strait. Furthermore, it is not clear whether the YSWC forms a persistent current into the Bohai Strait.

The Bohai Sea is an inland sea of China, with an average depth of about 18 m. The sea is rich in fishery resources, with busy industrial activities for ocean-shipping, petrochemical, and salt products [12]. The Bohai Strait is the only channel connecting the Bohai Sea to the open ocean, that is, the Yellow Sea. The Yellow Sea is a semiclosed shelf shallow sea (depth 70–100 m and average 44 m). There is a group of sand ridges located to the north of Bohai Strait [13, 14]. Surface observation data in Bohai Strait are not easily obtained, as waves whipped up by the wind roll turbulently in this area.

With the development of remote sensing technology, modern satellite data are frequently collected over the mixed zone between the Yellow Sea and Bohai Sea (Y-B) and are advantageous in monitoring the evolution of the YSWC path. Based on daily sea surface temperature (SST) of NOAA/AVHRR [15, 16], sea surface height (SSH) from AVISO, and the wind data recorded in the Chengshantou Station of China, Zhao et al. [17] identified two tongue-shaped branches of warm waters associated with the YSWC intrusion. Furthermore, statistical results revealed that these two branches were always present from 1981 to 2010. They found that the SSH was lower in the southern Yellow Sea than in the YSWC in the middle of September and from the end of February to early March. Starting in late March, the SSH was higher in the YSWC than in the southern Yellow Sea. Correlation analysis of the variation of the local SSH of the Yellow Sea indicated that the northwestern and northern branches of the YSWC were driven by relatively local SSH changes, caused by regional winds. They analyzed the characteristics of the YSWC using remote sensing data in combination with observational data; however, there was a lack of conclusive data on the SAL to analyze.

In a different area, Yuan et al. [18] studied the surface Kuroshio path in the Luzon Strait area using satellite ocean color, sea surface temperature, and altimeter data. The results suggested that the dominant path of the surface Kuroshio intrusion in the winter was a direct route from northeast of Luzon to southwest of Taiwan and then westward along the continental slope of the northern South China Sea.

All these issues are related to the YSWC intrusion into the Bohai Strait, which can be investigated using remote sensing ocean color and active microwave data. Many observations have indicated that the Kuroshio loop current is primarily a winter phenomenon and to date its occurrence in summer has not been frequently observed [1921]. Winter is therefore the best season to observe the YSWC [17]. Further, the development of remote sensing technology provides a great possibility of studying the regional characteristics of the YSWC. Research on the YSWC invasion into the Bohai Sea will significantly affect research on its ecological environment and the survival of fisheries spawning grounds. Its invasion clearly related to climate change.

In this study, we investigated the winter YSWC path in the mixed zone Y-B by remote sensing data based on wind field (atmospheric pressure changes cause sea surface wind field) and SST variations, with factor analyses of SAL, chlorophyll concentration (CHL), northerly wind days, and high wind days. We study that if remote sensing data show the intrusion of YSWC into the Bohai Sea, how to form decision criteria by remote sensing data and how to use the in situ data validate the criterion.

2. Materials and Methods

We used remote sensing data from the winter of 2013 as a case study to analyze and validate the in situ and cross section observational data. The sea surface wind field, temperature, SAL, CHL, wind scale, and high wind days were used to analyze characteristics of the mixed zone between the Bohai Sea and Yellow Sea in winters using remote sensing data.

2.1. Data Collection

The primary satellite data (Table 1) used to study the YSWC intrusion into the Bohai Sea in the winter was sea surface wind field data. Data were obtained from ASCAT_L3 daily average data on a global grid of about 25 km resolution (ftp://ftp.ifremer.fr/ifremer/cersat/products/gridded/MWF/L3/ASCAT/Daily/) and QuikSCAT_V4 daily average data in the Ku band on a global grid of about 25 km resolution (http://www.remss.com/missions/qscat). To reduce the cloudiness of the MODIS data, the monthly averaged sea surface SST and CHL data from MODIS_L3/Terra at 9 km resolution were used in this study (http://oceancolor.gsfc.nasa.gov/cgi/l3). The period of winter 2013 included December 2013 and January-February 2014. Scatterometer wind had an accuracy of less than 1 m/s [22, 23].

2.2. Methods

We computed the monthly average for all parameters. Validation data were from in situ and cross section observational data from ships. The SAL-retrieving algorithm for this study is from Qing et al. [24]. The model was recalibrated based on remote sensing reflectance data derived from MERIS and from in situ SAL data, yielding a root mean square error (RMSE) in modeled SAL of 1.31 psu (relative to in situ data). This recalibrated model was applied to MODIS data following a linear correction for SAL in the mixed zone Y-B: With winter remote sensing data from nearly 12 years, from 2002 to 2013, we analyzed the characteristics of the entry of the YSWC into the Bohai Sea, in the study area (Figure 1, spans 37.5–39°N/120–124°E in Y-B to avoid the land).

3. Results and Discussion

3.1. Wind Field Inversion of Y-B Mixed Zone from ASCAT Data

Wind field data were obtained from ASCAT_L3 data and displayed in Table 1. There were 90 wind field data points retrieved from ASCAT data, in the mixed zone, Y-B. In winter 2013, the average wind speed was 7.83 m/s (Figure 2(a)) compared with the average value of 7.71 m/s for 2002–2012 (Figure 2(b)). The wind speeds (2002–2012) were a fusion of ASCAT (2007–2012) and QuikSCAT (2002–2006) data, through Kriging interpolation in Figure 2(b). In winter 2013, the numbers of windy days were defined on the Beaufort scale 1-2, 2-3, 3-4, 4-5, 5-6, and 6-7, with 1, 12, 28, 24, 19, and 6 days, respectively (Table 2). On the Beaufort scale, a wind with average maximum speed above level 6 (≥10.8 m/s) for 10 min is termed a high wind and above level 8 (≥17.2 m/s) is a strong wind [25].

Table 2 shows the number of days with each of the different wind directions: east (E) (occurring in 2 days), south (S), west (W), north (N), NNE, NNW, SSW, SSE, NWW, WSW, ESE, and NEE. We found that northerly winds occurred in 55 days (61%) in winter 2013.

3.2. SST Inversion Revealed by MODIS Data

The retrieved SST from MODIS_L3 monthly average data in winters of 2013 and 2002–2012 showed that temperature varied from 4°C to 8°C in the study region. There were 174 SST data points retrieved from MODIS_L3 data in the mixed zone Y-B. The average temperature was 7.23°C in winter 2013 (Figure 3(a)), which was higher than the average value of 6.71°C (Figure 3(b)) from winter of 2002–2012.

3.3. CHL Inversion from MODIS Data

The CHL distribution retrieved from MODIS_L3 monthly average data in winters of 2013 and 2002–2012 varied from 0 to 4 mg/m3. There were 428 CHL data points retrieved from MODIS_L3 data in the mixed zone Y-B. The average CHL in winter 2013 was 3.02 mg/m3 (Figure 4(a)) and 3.06 mg/m3 for 2002–2012 (Figure 4(b)). There was no more variation in CHL in 2013 than in 2002–2012.

3.4. Salinity Inversion from MODIS

MODIS_L3 observations and inversion SAL from formula (1) showed that the SAL ranged from 31.03 to 32.93 psu in the winter of 2013 (Figure 5(a)). There were 779 SAL data points retrieved from MODIS_L3 data in the mixed zone Y-B. In the study region, the averaged value ranged from 28.97 to 32.82 psu (Figure 5(b)) in 2002–2012.

3.5. The In Situ and Cross Section Observation Data Analysis in Winter 2013

We analyzed the in situ and cross section observational SST and SAL data measured by shipboard CTD collected at the 42 stations shown in Figure 1 (+ symbol display site). These observations were made in Y-B on December 17–24, 2013, at depths of 0 m and 10 m and at the bottom (Figures 6 and 7).

As shown in Figure 6, the SST data came from the sea surface (Figure 6(a)), from 10 m water depth (Figure 6(b)), and at the bottom (Figure 6(c)). The SAL data came from the sea surface (Figure 7(a)), 10 m water depth (Figure 7(b)), and the bottom (Figure 7(c)).

3.6. Results

Through the aforementioned remote sensing inversion average data for years contrast, we found that, in winter 2013, winds were mostly northerly (Table 2). The average wind speeds of the two periods were similar, only 0.12 m/s higher in 2013 than in 2002–2012, with an RMSE of 0.26 m/s. On average high wind days were 2 days shorter in winter 2013 compared with 2002–2012, which had an average of 8.97 days. However, the average SST was 0.52°C in 2013, which was much higher than the average of 2002–2012; the RMSE was 1.85°C. The average CHL in the mixed zone in winter 2013 was similar to the average from 2002 to 2012, only 0.04 mg/m3 lower in 2013 than in 2002–2012; the RMSE was 0.59 mg/m3. Figure 4 shows a low average CHL at 0–2 mg/m3 near the mixing zone in winter 2013, which is not existent in the previous decade to see Figure 4(b). The average SAL in winter 2013 was 32.01 and 31.29 psu in 2002–2012 (see Figure 5); RMSE was 1.84 psu. Thus, we hypothesized that the YSWC entered the Bohai Sea in winter 2013.

From in situ and cross section observational data, we found that a warm water mass (SST: 7–10°C; SAL: 30.1–31.5 psu) invaded the Bohai Sea from the Yellow Sea on December 17–24, 2013, at all depths. We found that, in December 2013, the SAL and SST from observational data did not vary greatly between different depths. Because of the shallowness of the research area and the wind and wave mixing in winter, the region formed a homogeneous mixed layer across different depths of seawater in Y-B, which confirms the winter 2013 remote sensing analysis. The YSWC invaded the Bohai Sea in winter 2013.

3.7. Discussion

To enable an understanding of the movement of the YSWC into the Bohai Sea, it is important to establish decision criteria by remote sensing data.

The YSWC is characterized by high temperature, high SAL, and low dissolved oxygen content. Figures 8(a), 8(b), 8(c), 9(a), 9(b), and 9(c) show that the average , average , average high wind days, average SST, average SAL, and average CHL over 12 winters were 7.72 m/s, 195.4°, 25.17 days, 6.76°C, 31.37 psu, and 3.05 mg/m3, respectively. Table 3 shows that northerly winds dominated in the mixed zone Y-B in the winters of 2002–2013. There are no laws to follow from the average high wind days.

We found that the average V, average SST, and average SAL in winter 2013, retrieved by remote sensing data, were higher than the 12-year average. High wind days and average CHL were lower than the 12-year average (Table 4). The northerly winds dominated in winter 2013. By combining in situ, cross section observation and remote sensing data analysis, we found that the movement of the YSWC into the Bohai Sea is caused by positive anomalies of northerly wind speeds, SST, and SAL and a negative anomaly of CHL over multiyear average in the mixed zone Y-B.

According to the above criteria, the average SST in winter did not conform to the conditions in 2002, 2005, 2008, 2009, 2010, and 2012. The average SAL in winter did not conform to the conditions in 2006, 2009, 2010, 2011, and 2012. The average CHL in winter did not conform to the conditions in 2003, 2005, 2008, 2010, 2011, and 2012. The average wind speed in winter did not conform to conditions in 2002, 2003, 2006, 2011, and 2012 (see Tables 3 and 4).

4. Conclusions

In winters of 2004, 2007, and 2013, when northerly winds dominated, the average wind speeds were 8.18 m/s, 8.01 m/s, and 7.83 m/s, respectively, which were all greater than the 12-year average (7.72 m/s); the SAL and SST were higher than the 12-year average (6.76°C and 31.37 psu), while the CHL was less than the 12-year average (3.05 mg/m3, the average CHL is close to this value in winter 2007). The YSWC intrusion happened in 2007, 2004, and 2013, ranked by the intrusion strength. This conclusion was validated by the in situ and cross section observational data from winter 2013. We conclude that the YSWC intrusion into the Bohai Sea occurred when the wind speed, sea surface temperature, and SAL were above (or close to) the multiyear average and the CHL was less than (or close to) the multiyear average. This finding agrees with the model simulation results published in the literature [4, 7, 8, 17].

The study area is fully mixed in the winter owing to the wind field and ocean currents, and thus it is the ideal season to observe the intrusion of the YSWC into the Bohai Sea, by means of remote sensing; this conclusion is the same as Zhao et al. [17]. Therefore, the conclusions obtained from remote sensing data in this study are reasonable in this area, in the winter.

We also found that the YSWC intrusion into the Bohai Sea did not happen every year and was closely related to the wind field.

Future work shall increase in situ and cross section observational data from the study area and continue to verify the conclusions from remote sensing data. We will focus on the conditions and mechanism of the YSWC intrusion into the Bohai Sea during the autumn and winter and expand the study area to the north Yellow Sea and the Kuroshio area.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This project is financially supported by the “Strategic Priority Research Program” of the Chinese Academy of Sciences (no. XDA01020304) and by the National Natural Science Foundation of China (no. 41176160 and no. 4141101049). In addition, the authors thank ECMWF, EUMETSAT, KNMI, IFREMER/CERSAT, GODIVA, NOAA, Météo-France, NDBC, O&SI SAF, PMEL, and the UK Met Office for the courtesy of providing buoy, numerical, and satellite data used for daily ASCAT data production, evaluation, and validation. MODIS data were provided by NASA. Helpful advice from Dr. Xiaofeng Li NOAA is greatly appreciated.