Analysis of recent climate change over the Arctic using ERA-Interim reanalysis data

This study investigates recent climate change over the Arctic and its link to the mid-latitudes using the ERA-Interim global atmospheric reanalysis data from the European Center for Medium-Range Weather Forecast （ECMWF）. Since 1979, sub- stantial surface warming, associated with the increase in anthropogenic greenhouse gases, has occurred over the Arctic. The great- est warming in winter has taken place offshore in the Kara-Barents Sea, and is associated with the increase in turbulent heat fluxes from the marginal ice zone. In contrast to the marked warming over the Arctic Ocean in winter, substantial cooling appears over Siberia and eastern Asia, linked to the reduction of Arctic sea ice during the freezing season （September-March）. However, in summer, very little change is observed in surface air temperature over the Arctic because increased radiative heat melts the sea ice and the amount of turbulent heat gain from the ocean is relatively small. The heat stored in the upper ocean mixed layer in summer with the opening of the Arctic Ocean is released back to the atmosphere as turbulent heat fluxes during the autumn and through to the following spring. This warming of the Arctic and the reduced sea ice amplifies surface cooling over Siberia and eastern Asia in winter.     

Interpretation of sea surface wind interannual vector EOFs over the China seas

Using interpolation and averaging methods, we analyzed the sea surface wind data obtained from December 1992 to November 2008 by the scatterometers ERS-1, ERS-2, and QuikSCAT in the area of 2°N–39 °N, 105°E–130°E, and we reported the monthly mean distributions of the sea surface wind field. A vector empirical orthogonal function (VEOF) method was employed to study the data and three temporal and spatial patterns were obtained. The first interannual VEOF accounts for 26% of the interannual variance and displays the interannual variability of the East Asian monsoon. The second interannual VEOF accounts for 21% of the variance and reflects the response of China sea winds to El Niño events. The temporal mode of VEOF-2 is in good agreement with the curve of the Niño 3.4 index with a four-month lag. The spatial mode of VEOF-2 indicates that four months after an El Niño event, the southwesterly anomalous winds over the northern South China Sea, the East China Sea, the Yellow Sea, and the Bohai Sea can weaken the prevailing winds in winter, and can strengthen the prevailing winds in summer. The third interannual VEOF accounts for 10% of the variance and also reflects the influence of the ENSO events to China Sea winds. The temporal mode of VEOF-3 is similar to the curve of the Southern Oscillation Index. The spatial mode of VEOF-3 shows that the northeasterly anomalous winds over the South China Sea and the southern part of the East China Sea can weaken the prevailing winds, and southwesterly anomalous winds over the northern part of the East China Sea, the Yellow Sea, and the Bohai Sea can strengthen the prevailing winds when El Niño occurs in winter. If El Niño happens in summer, the reverse is true.     

Features of visibility variation at Great Wall Station, Antarctica

The variation of visibility at Great Wall Station （GWS） was analyzed using manual observational data for the period of 1986 to 2012. Results show that the frequencies of occurrence of high （≥ 10 km） and low visibility （0-1 km） are 61.0% and 8.0%, respectively. Visibility at GWS shows an evident seasonal variation： The highest visibility between November and March, and the lowest visibility from June to October. Sea fog and precipitation are the main factors for low visibility during summer, whereas frequent adverse weather, such as falling snow, blowing snow, or blizzards, are responsible for low visibility in winter. The frequency of occurrence of low visibility has decreased significantly from 1986 to 2012. Conversely, the frequency of occurrence of high visibility has shown a significant increasing trend, especially during winter. The decreasing tendencies of fog, blowing snow, and snowfall have contributed to the increasing trend of high visibility during winter. Visibility at GWS exhibits significant synoptic-scale （2.1 to 8.3 d）, annual, and inter-annual periods （2 a, 4.1 a, and 6.9 a to 8.2 a）, among which the most significant period is 4.1 a. The visibility observed during 2012 indicates that instrumental observation can be applied in the continuous monitoring of visibility at GWS.     

Population distribution,structure and growth condition of Antarctic krill （Euphausia superba Dana） during the austral summer in the Southern Ocean

Antarctic krill （Euphausia superba Dana） was collected using a High Speed Collector and an Isaac-Kidd midwater trawl （IKMT） net during the austral summer of 2007/2008 and 2008/2009 in the circumpolar and Prydz Bay regions of the Southern Ocean, respectively. Combined with the simultaneous recording of environmental factors, spatial distribution, population structure and growth condition of E. superba were studied. The abundance of E. superba in the Weddell Sea was higher than in Prydz Bay. However, the abundance of E. superba in both the Weddell Sea and Prydz Bay was lower than figures reported in previous krill surveys for the same time period. With respect to the total study area, E. superba displayed a normal growing state during the two expeditions. E. superba grew relatively poorly in some stations, which may be due to the late retreat of sea ice or lower chlorophyll a concentrations. The number of juvenile E. superba collected using the High Speed Collector was proportionally greater in stations located at the edge of the sea ice, while adults dominated in long-term non-ice regions. This phenomenon reflects the different distribution pattern between juvenile and adult krill. The population structure of E. superba differed between sea regions, which may affect recruitment.     

Comparisons of several evaporation/precipitation datasets for the Bohai Sea based on salinity simulation

Evaporation (E) rate and precipitation (P) rate are two significant meteorological elements required in the ocean baroclinic modeling as external forcings. However, there are some uncertainties in the currently used E/P rates datasets, especially in terms of the data quality. In this study, we collected E/P rates data from ERA-40, NCEP/NCAR Reanalysis, HOAPS for the Bohai Sea and nine routine stations around Laizhou Bay, and made comparisons among them. It was found that the differences in E/P rates between land and sea are remarkable, which was due to the difference in underlying surfaces. Therefore, the traditional way of using E/P rates acquired on land directly at sea is not correct. Since no final conclusion has been reached concerning the net water transport between the Bohai Sea and the Yellow Sea, it is unfeasible to judge the adequacy of the four kinds of data by using the water budget equation. However, the E/P rates at ERA-40 sea points were considered to be the optimal in terms of temporal/spatial coverage and resolution for the hindcast of salinity variation in the Bohai Sea. Besides, using the 3-D hydrodynamic model HAMSOM (HAMburg Shelf Ocean Model), we performed numerical experiments with different E/P datasets and found that the E/P rates at sea points from ERA-40 dataset are better than those from NCEP/NCAR Reanalysis dataset. If NCEP/NCAR Reanalysis E/P rates are to be used, they need to be adjusted and tested prior to simulation so that more close-to-reality salinity values can be reproduced.     

Ice concentration assimilation in a regional ice-ocean cou- pled model and its application in sea ice forecasting

A reasonable initial state of ice concentration is essential for accurate short-term forecasts of sea ice using ice-ocean coupled models. In this study, sea ice concentration data are assimilated into an operational ice forecast system based on a com- bined optimal interpolation and nudging scheme. The scheme produces a modeled sea ice concentration at every time step, based on the difference between observational and forecast data and on the ratio of observational error to modeled error. The impact and the effectiveness of data assimilation are investigated. Significant improvements to predictions of sea ice extent were obtained through the assimilation of ice concentration, and minor improvements through the adjustment of the upper ocean properties. The assimilation of ice thickness data did not significantly improve predictions. Forecast experiments show that the forecast accuracy is higher in summer, and that the errors on five-day forecasts occur mainly around the marginal ice zone.     

Influence of sea ice and sea surface temperature on the abundance of Antarctic krill in Area 48.2

In this paper we examine the relationship between Antarctic krill catch, sea ice concentration, and sea surface tempera- ture （SST）. Data on the Antarctic krill catch from 2003 to 2010 in CCAMLR Area 48.2 were combined with sea ice and SST data. Results showed that krill fishing in Area 48.2 took place from February to August each year but the catch was concentrated from March to July, with production during this period accounting for about 99.3% of the annual catch. Regression analysis showed that the catch per unit effort （CPUE） was clearly related to sea ice concentration and SST intervals. CPUE was negatively correlated with the area of sea ice among years （R2=0.64）, and the correlation was strongest （R2=0.71） when sea ice concentration was greater than 90%. Over the months the CPUE initially increased, then decreased as the area of sea ice increased. The relationship was strongest （R2=0.88） when the concentration of sea ice was 60%--70%. There was no negative correlation among years between CPUE and the ice-free area when S ST was between -2 ℃ and 3 ℃ （R2=0.21）, but there was a significant negative correlation when SST was between 1 ℃ and 2℃ （R2=0.82）. Over the months, CPUE initially increased then decreased with increasing sea ice-free area, and the relationship was strongest （R2=0.94） when SST was between 0℃and 1 ℃. This study shows that sea ice concentra- tion and SST have significant effects on the abundance of krill in Area 48.2, and the findings have practical significance for the use and conservation of Antarctic krill resources.     

Thickness distribution and extent of sea ice and snow in Prydz Bay and surrounding waters observed in 2012/2013 austral summer

Three ship-based observational campaigns were conducted to survey sea ice and snow in Prydz Bay and the surrounding waters(64.40°S–69.40°S, 76.11°E–81.29°E) from 28 November 2012 to 3 February 2013. In this paper, we present the sea ice extent and its variation, and the ice and snow thickness distributions and their variations with time in the observed zone. In the pack ice zone, the southern edge of the pack ice changed little, whereas the northern edge retreated significantly during the two earlier observation periods. Compared with the pack ice, the fast ice exhibited a significantly slower variation in extent with its northernmost edge retreating southwards by 6.7 km at a rate of 0.37 km?d-1. Generally, ice showed an increment in thickness with increasing latitude from the end of November to the middle of December. Ice and snow thickness followed an approximate normal distribution during the two earlier observations(79.7±28.9 cm, 79.1±19.1 cm for ice thickness, and 11.6±6.1 cm, 9.6±3.4 cm for snow thickness, respectively), and the distribution tended to be more concentrated in mid-December than in late November. The expected value of ice thickness decreased by 0.6 cm, whereas that of snow thickness decreased by 2 cm from 28 November to 18 December 2012. Ice thickness distribution showed no obvious regularity between 31 January and 3 February, 2013.     

Seasonal and inter-annual variations of the primary types of the Arctic sea-ice drifting patterns

Abstract Monthly mean sea ice motion vectors and monthly mean sea level pressure （SLP） for the period of 1979-2006 are investigated to understand the spatial and temporal changes of Arctic sea-ice drift. According to the distinct differences in monthly mean ice velocity field as well as in the distribution of SLP, there are four primary types in the Arctic Ocean： Beaufort Gyre＋Transpolar Drift, Anticyclonic Drift, Cyclonic Drift and Double Gyre Drift. These four types account for 81% of the total, and reveal distinct seasonal variations. The Cyclonic Drift with a large-scale anticlockwise ice motion pattern trends to prevail in summer while the Anticyclonic Drift with an opposite pattern trends to prevail in winter and spring. The prevailing seasons for the Beaufort Gyre＋Transpolar Drift are spring and autumn, while the Double Gyre Drift trends to prevail in winter, especially in Feb- ruary. The annual occurring times of the Anticyclonic Drift and the Cyclonic Drift are closely correlated with the yearly mean Arc- tic Oscillation （AO） index, with a correlation coefficient of -0.54 and 0.54 （both significant with the confident level of 99%）, re- spectively. When the AO index stays in a high positive （negative） condition, the sea-ice motion in the Arctic Ocean demonstrates a more anticlockwise （clockwise） drifting pattern as a whole. When the AO index stays in a neutral condition, the sea-ice motion becomes much more complicated and more transitional types trend to take place.     

Interaction between sea ice/iceberg and ship structures: A review

For ship structural design and good maneuverability in an ice-covered sea, the local and global load of ice cover on ships should be well understood. This paper reviews the extensive work done on ice loads on ships, including： （a） Ice pressure and local load determination based on field and model tests; （b） Global ice loads on ships from full-scale field observations, model tests and numerical models under different ice conditions （level ice and pack ice） and ship operations （maneuvering and mooring）. Spe- cial attention is paid to the discrete element simulation of global ice loads on ships; and （c） Analytical solutions and numerical models of impact loads of icebergs on ships for polar navigation. Finally, research potential in these areas is discussed.     

Composition and distribution of diatom assemblages in the surface sediments of the Bering Sea

Diatoms collected from the top 1 cm of the surface sediment layer at seven Bering Sea stations during the Fourth Chi- nese National Arctic Research Expedition （4th CH1NARE-Arctic） in 2010 were studied. In total, 101 taxa belonging to 38 genera were found. The species were divided into four groups： eurythermal species, boreal-temperate species, polar species, and sea-ice species. The diatom assemblages at the stations in the southwestern basin area were primarily composed of boreal-temperate species, such as Neodenticula seminae, Thalassiosira trifulta, Rhizosolenia hebetata f. hiemalis, and Actinocyclus curvatulus. The northeastern shelf stations were dominated by polar species, including Fragilariopsis oceanica, Thalassiosira antarctica spora, Thalassiosira nordensldoeldii, and Thalassiosira hyalina. The overall abundance was highest at the basin stations with 3.7 ×10^6 cells·g^-1 of wet sediment, whereas the eastern shelf stations had the lowest abundance of 0.7×10^6 cells·g^-1 of wet sediment （excluding the resting spores of Chaetoceros spp.）. The relationship between the distribution of the surface sediment diatom assemblages and the environment is discussed.     

Variability of thermohaline structure at 4°S, 156°E during TOGA COARE IOP

The thermohaline structure at 4°S, 156°E was analyzed based on CTD data acquired during the TOGA COARE Intensive Observing Period (IOP) from November, 1992 to February, 1993. The ocean responses during two Madden-Julian Oscillation (MJO) events were preliminarily studied based on meteorological field observation. The main water masses at the observation point were Tropical Surface Water, Southern Subtropical Lower Water and Southern Intermediate Water from surface downward. There was good correlation of sea surface temperature with the wind field, and of the surface salinity with wind speed and rainfalls. Both of the two surface variables were also modulated by upwelling caused by westerly winds at the observation point. The isohaline layer was not always shallower than the isothemal layer in this observation and could be considered as the lower limit of the diurnal variation of the isothernal layers in most cases. The existence of large variations of the maximum salinity core is suggested to be related to the meridional motion in that depth. Contribution No. 2264 from the Institute of Oceanology, Chinese Academy of Sciences. This project was supported by NSFC (No. 49176255).     

A physics-based dual-frequency approach for altimeter wind speed retrieval

The altimeter normalized radar cross section(NRCS) has been used to retrieve the sea surface wind speed for decades, and more than a dozen of wind speed retrieval algorithms have been proposed. Despite the continuing efforts to improve the wind speed measurements, a bias dependence on wave state persists in all wind algorithms. On the basis of recent evidence that short waves are essentially modulated by local winds and much less affected by wave state, we proposed a physics-based approach to retrieve the wind speed from the dual-frequency difference in terms of the mean square slope of short waves. A collocated dataset of coincident altimeter/buoy measurements were used to develop and validate the approach. Validation against buoy measurements indicates that the approach is almost unbiased and has an overall root mean square error of 1.24 m s-1, which is 5.3% lower than the single-parameter algorithm in operational use(Witter and Chelton, 1991) and 2.4% lower than another dual-frequency approach(Chen et al., 2002). Furthermore, the results indicate that the new approach significantly improves the wave-dependent bias compared to the single-parameter algorithm. The capacity of altimeter to retrieve sea surface wind speed appears to be limited for the case of winds below 3 m s-1. The validity of the approach at high winds needs to be further examined in the future study.     

Assimilation of GMS-5 satellite winds using nudging method with MM5

With the aid of Meteorological Information Composite and Processing System （MICAPS）, satellite wind vectors derived from the Geostationary Meteorological Statellite-5 （GMS-5） and retrieved by National Satellite Meteorology Center of China （NSMC） can be obtained. Based on the nudging method built in the fifth-generation Mesoscale Model （MM5） of Pennsylvania State University and National Center for Atmospheric Research, a data preprocessor is developed to convert these satellite wind vectors to those with specified format required in MM5. To examine the data preprocessor and evaluate the impact of satellite winds from GMS-5 on MM5 simulations, a series of numerical experimental forecasts consisting of four typhoon cases in 2002 are designed and implemented. The results show that the preprocessor can process satellite winds smoothly and MM5 model runs successfully with a little extra computational load during ingesting these winds, and that assimilation of satellite winds by MM5 nudging method can obviously improve typhoon track forecast but contributes a little to typhoon intensity forecast. The impact of the satellite winds depends heavily upon whether the typhoon bogussing scheme in MM5 was turned on or not. The data preprocessor developed in this paper not only can treat GMS-5 satellite winds but also has capability with little modification to process derived winds from other geostationary satellites.     

Multiple processes affecting surface seawater N20 satura- tion anomalies in tropical oceans and Prydz Bay, Antarc- tica

Abstract We analyzed the N20 content of surface seawater sampled from Prydz Bay, Antarctica, on a cruise track between 30°S and 30°N during the twenty-second Chinese National Antarctic Research Expedition during austral summer, 2006. The surface water showed an average pN20 value of 311.9a：7.6 nL.L-1 （14.1：t：0.4 nmol-Ll）, which was slightly undersaturated. The air-sea N2O flux in the region was -0.3±0.8 p.mol.mE,dl; howeveh N20 in the surface water was oversaturated in most stations along the cruise track. Saturation anomalies were greater than 10%, with a maximum of 54.7% being observed at the Eqaator, followed by 31% at 10°N in the Sulu Sea. The air-sea fluxes at these locations were 12.4 and 4 μmol-mLd-1, respectively. Overall, the results indicated that surface water in Prydz Bay was near equilibrium with atmospheric N20, and that ocean waters in lower latitudes acted as a N20 source. Physical processes such as stratification, ice-melt water dilution, and solar radiation dominate the factors leading to N20 saturation of surface water of Prydz Bay, while biological production and upwelling are primarily responsible for the N20 oversaturation of surface water observed in subtropical and tropical regions along the cruise track.     

Evaluation of the biological toxicity of lfuorine in Antarctic krill

Antarctic krill is a potentially nutritious food source for humans, but lfuorine （F） toxicity is a matter of concern. To evaluate the toxicity of F in Antarctic krill, 30 Wistar rats were divided into three groups with different dietary regimens：a control group, a krill treatment group （150 mg＆#183;kg-1 F）, and a sodium lfuoride （NaF） treatment group （150 mg＆#183;kg-1 F）. After three months, F concentrations in feces, plasma, and bone were determined, and the degree of dental and skeletal lfuorosis was assessed. The F concentrations in plasma and bone from the krill treatment group were 0.167 0＆#177;0.020 4 mg.L-1 and 2 709.8＆#177;301.9 mg＆#183;kg-1, respectively, compared with 0.043 8＆#177;0.005 5 mg＆#183;L-1 and 442.4＆#177;60.7 mg＆#183;kg-1, respectively, in samples from the control group. Concentrations of F in plasma and bone in the krill treatment group were higher than in the control group, but lower than in the NaF treatment group. The degree of dental lfuorosis in the krill treatment group was moderate, compared with severe in the NaF treatment group and normal in the control group. The degree of skeletal lfuorosis did not change signiifcantly in any group. These results showed that the toxicity of F in Antarctic krill was lower than for an equivalent concentration of F in NaF, but it was toxic for rats consuming krill in large quantities. To conclude, we discuss possible reasons for the reduced toxicity of F in Antarctic krill. The present study provides a direct toxicological reference for the consideration of Antarctic krill for human consumption.     

Sea Salt Sodium Record in a Shallow Ice Core from East Antarctica as a Potential Proxy of the Antarctic Sea Ice Extent in Southern Indian Ocean

Antarctic sea ice has experienced an increasing trend in recent decades, especially in the Ross Sea and Indian Ocean sectors. Sea ice variability affects greatly the maritime airmass transport from high latitude to Antarctic continent. Here we present a new ice core record of sea salt sodium(ssNa+) concentration at annual-resolution in the Princess Elizabeth Land spanning from 1990 to 2016, showing that this marker could be used as a potential proxy for reconstructing the sea ice extent(SIE) in the Southern Indian Ocean(SIO) given their significant correlation(R =-0.6, P 0.01) over the past 27 years. The correlation and composite analyses results show that the ssNa~+ at the 202 km inland from Zhongshan Station and the SIE changes in SIO are closely related to the Indian Ocean Dipole(IOD) and Southern Annular Mode(SAM). The northward wind in central SIO occurs during positive IOD and the strengthened westerlies occurs during positive SAM, both of which favor increased sea ice in SIO and lead to the decreased ssNa~+ concentration at the coastal site.     

An Observational Study of Typhoon Imbudo in 2003

Typhoon Imbudo was a super-typhoon over the northwestern Pacific in 2003. It caused tremendous damage when it made landfalls in the Philippines and China. This paper documents observational analyses of Typhoon Imbudo during its landfall in China. All available observations are used to study its motion, intensity changes, convection, structure and precipitation. Best-track data indicate that Imbudo moved west-northwestward until 1800 UTC 23 July and then turned northwestward. FNL （final） analysis data show that the motion of Imbudo is dominated by changes of the subtropical high. At Imbudo＇s mature stage, the minimum sea level pressure dropped to 910 hPa and the maximum sustained winds were as high as 67 m s^-1, which is the intensity of a super-typhoon. The surface wind field exhibited asymmetric characteristics. Polar-orbiting satellite imagery also manifested convective asymmetry before Imbudo made landfall in China. Analyzed the vertical wind shear, it is shown that the convection has a downshear-left pattern. All kinds of precipitation data were used to identify the asymmetric characteristic of the rainfall associated with the Imbudo. The maximum rainfalls were located in the southern boundary area between Guangxi and Guangdong. However, the lack of in situ observations limited further analyses of this typhoon.     

Major mid-late Holocene cooling in the East China Sea revealed by an alkenone sea surface temperature record

Although the mid-late Holocene cold and dry event about 4000 years ago(the 4 ka event) has been observed almost globally, it was most prominent in terrestrial climate proxies from the lower latitudes. Here we evaluate the oceanic response to this event in terms of a Holocene sea surface temperature(SST) record reconstructed using the K'37U index for Core B3 on the continental shelf of the East China Sea. The record reveals a large temperature drop of about 5℃ from the mid-Holocene(24.7℃ at 5.6 ka) to the 4 ka event(19.2℃ at 3.8 ka). This mid-late Holocene cooling period in Core B3 correlated with(i) decreases in the East Asia summer monsoon intensity and(ii) the transition period with increased El Nino/Southern Oscillation activities in the Equatorial Pacific. Our SST record provides oceanic evidence for a more global nature of the mid-late Holocene climate change, which was most likely caused by a southward migration of the Intertropical Converge Zone in response to the decreasing summer solar insolation in the Northern Hemisphere. However, the large SST drop around Core B3 indicates that the mid-late Holocene cooling was regionally amplified by the initiation/strengthening of eddy circulation/cold front which caused upwelling and resulted in additional SST decrease. Upwelling during the mid-late Holocene also enhanced with surface productivity in the East China Sea as reflected by higher alkenone content around Core B3.     

Holocene Temperature Records from the East China Sea Mud Area Southwest of the Cheju Island Reconstructed by the U^K_37 and TEX86 Paleothermometers

As an important marginal sea under the influences of both the Changjiang River and the Kuroshio, the East China Sea （ECS） environment is sensitive to both continental and oceanic forcing. Paleoenvironmental records are essential for understanding the long-term environmental evolution of the ECS and adjacent areas. However, paleo-temperature records from the ECS shelf are currently very limited. In this study, the U^K_37 and TEX86 paleothermometers were used to reconstruct surface and subsurface temperature changes of the mud area southwest of the Cheju Island （Site F10B） in the ECS during the Holocene. The results indicate that temperature changes of F 10B during the early Holocene （11.6-6.2 kyr） are associated with global climate change. During the period of 6.2-2.5 kyr, the similar variability trends of smoothing average of AT （the difference between surface and subsurface temperature） of Site F10B and the strength of the Kuroshio suggest that the Kuroshio influence on the site started around 6.2kyr when the Kuroshio entered the Yellow Sea and continued to 2.5 kyr. During the late Holocene （2.5-1.45 kyr）, apparent decreases of U^K_37 sea surface temperature （SST） and AT imply that the direct influence of the Kuroshio was reduced while cold eddy induced by the Kuroshio gradually controlled hydrological conditions of this region around 2.5 kyr.