National Marine Environmental Forecasting Center

Over the past decade it has become apparent that the atmosphere is a significant pathway for the transport of many natural and pollutant materials from the continents to the ocean. The atmospheric input of many of these species can have an impact (either positive or negative) on biological processes in the sea and on marine chemical cycling. For example, there is now evidence that the atmosphere may be an important transport path for such essential nutrients as iron and nitrogen in some regions. In this report we assess current data in this area, develop global scale estimates of the atmospheric fluxes of trace elements, mineral aerosol, nitrogen species, and synthetic organic compounds to the ocean; and compare the atmospheric input rates of these substances to their input via rivers. Trace elements considered were Pb, Cd, Zn, Cu, Ni, As, Hg, Sn, Al, Fe, Si, and P. Oxidized and reduced forms of nitrogen were considered, including nitrate and ammonium ions and the gaseous species NO, NO 2 , HNO 3 , and NH 3 . Synthetic organic compounds considered included polychlorinated biphenyls (PCBs), hexachlorocyclohexanes (HCHs), DDTs, chlordane, dieldrin, and hexachlorobenzenes (HCBs). Making this assessment was difficult because there are very few actual measurements of deposition rates of these substances to the ocean. However, there are considerably more data on the atmospheric concentrations of these species in aerosol and gaseous form. Mean concentration data for 10° × 10° ocean areas were determined from the available concentration data or from extrapolation of these data into other regions. These concentration distributions were then combined with appropriate exchange coefficients and precipitation fields to obtain the global wet and dry deposition fluxes. Careful consideration was given to atmospheric transport processes as well as to removal mechanisms and the physical and physicochemical properties of aerosols and gases. Only annual values were calculated. On a global scale atmospheric inputs are generally equal to or greater than riverine inputs, and for most species atmospheric input to the ocean is significantly greater in the northern hemisphere than in the southern hemisphere. For dissolved trace metals in seawater, global atmospheric input dominates riverine input for Pb, Cd, and Zn, and the two transport paths are roughly equal for Cu, Ni, As, and Fe. Fluxes and basin‐wide deposition of trace metals are generally a factor of 5‐10 higher in the North Atlantic and North Pacific regions than in the South Atlantic and South Pacific. Global input of oxidized and reduced nitrogen species are roughly equal to each other, although the major fraction of oxidized nitrogen enters the ocean in the northern hemisphere, primarily as a result of pollution sources. Reduced nitrogen species are much more uniformly distributed, suggesting that the ocean itself may be a significant source. The global atmospheric input of such synthetic organic species as HCH,PCBs, DDT, and HCB completely dominates their input via rivers.

Numerical model studies have suggested that the ongoing global climate change will likely affect tropical cyclone activity. Since the global warming has been underway, it is meaningful to ask: Are there evidences of observed changes in tropical cyclone activity? Using best‐track data from 1965 to 2003, we show for the first time that over the past four decades the two prevailing typhoon tracks in the western North Pacific (WNP) have shifted westward significantly; thus the subtropical East Asia has experienced increasing typhoon influence; but the typhoon influence over the South China Sea has considerably decreased. Our trajectory model simulation indicates that the long‐term shifts in the typhoon tracks result primarily from the changes in the mean translation velocity of typhoons or the large‐scale steering flow, which is associated with the westward expansion and strengthening of the WNP subtropical high.

Abstract The polar regions have been attracting more and more attention in recent years, fueled by the perceptible impacts of anthropogenic climate change. Polar climate change provides new opportunities, such as shorter shipping routes between Europe and East Asia, but also new risks such as the potential for industrial accidents or emergencies in ice-covered seas. Here, it is argued that environmental prediction systems for the polar regions are less developed than elsewhere. There are many reasons for this situation, including the polar regions being (historically) lower priority, with fewer in situ observations, and with numerous local physical processes that are less well represented by models. By contrasting the relative importance of different physical processes in polar and lower latitudes, the need for a dedicated polar prediction effort is illustrated. Research priorities are identified that will help to advance environmental polar prediction capabilities. Examples include an improvement of the polar observing system; the use of coupled atmosphere–sea ice–ocean models, even for short-term prediction; and insight into polar–lower-latitude linkages and their role for forecasting. Given the enormity of some of the challenges ahead, in a harsh and remote environment such as the polar regions, it is argued that rapid progress will only be possible with a coordinated international effort. More specifically, it is proposed to hold a Year of Polar Prediction (YOPP) from mid-2017 to mid-2019 in which the international research and operational forecasting communites will work together with stakeholders in a period of intensive observing, modeling, prediction, verification, user engagement, and educational activities.

Abstract Both the Arctic and Antarctic sea ice extents (SIEs) from 44 coupled models in the Coupled Model Intercomparison Project Phase 6 (CMIP6) are evaluated by comparing them with observations and CMIP5 results. The CMIP6 multimodel mean can adequately reproduce the seasonal cycles of both the Arctic and Antarctic SIE. The observed Arctic September SIE declining trend (−0.82 ± 0.18 million km 2 per decade) between 1979 and 2014 is slightly underestimated in CMIP6 models (−0.70 ± 0.06 million km 2 per decade). The observed weak but significant upward trend of the Antarctic SIE is not captured, which was an issue already in the CMIP5 phase. Compared with CMIP5 models, CMIP6 models have lower intermodel spreads in SIE mean values and trends, although their SIE biases are relatively larger. The CMIP6 models did not reproduce the new summer tendencies after 2000, including the faster decline of Arctic SIE and the larger interannual variability in Antarctic SIE.

Pattern and origins of Kuroshio branches in the bottom water of southern East China Sea (ECS), were carefully examined by numerical simulations based on the Regional Ocean Modeling System (ROMS) together with observations. Model results show that in the bottom water of ECS, the intrusion pattern of Kuroshio is mainly composed of an Offshore Kuroshio Branch Current (OKBC) which, bifurcated from the Kuroshio northeast of Taiwan, flows nearly along the isobath of ∼100 m, and a Nearshore Kuroshio Branch Current (NKBC) which, originated from the Kuroshio northeast of Taiwan, upwells northwestward gradually from ∼250 m to ∼60 m, then turns to northeast around 27.5°N, 122°E, thereafter flows northeastward along the isobath of ∼60 m, and finally reaches at 30.5°N where it turns to east. Furthermore, we found that the NKBC mostly originated in the Kuroshio subsurface water (120–250 m) east of Taiwan, whereas the OKBC mainly stemmed from the Kuroshio water (60–120 m) east of Taiwan. This pattern and origins of OKBC and NKBC well addressed the observational phenomena that off the coast of Zhejiang province, China, there were colder, less saline, and more phosphate‐rich bottom water near the isobath of ∼60 m rather than near the isobath of ∼100 m in August 2009. Finally, it is proposed that on southern ECS continental shelf, Kuroshio exhibits its intrusion branches by an anticyclonical stair structure: bottom stair NKBC, middle stair OKBC, and top stair Kuroshio surface branch (KBC).

BACKGROUND: The fall armyworm (FAW), an invasive pest from the Americas, is rapidly spreading through the Old World, and has recently invaded the Indochinese Peninsula and southern China. In the Americas, FAW migrates from winter-breeding areas in the south into summer-breeding areas throughout North America where it is a major pest of corn. Asian populations are also likely to evolve migrations into the corn-producing regions of eastern China, where they will pose a serious threat to food security. RESULTS: To evaluate the invasion risk in eastern China, the rate of expansion and future migratory range was modelled by a trajectory simulation approach, combined with flight behavior and meteorological data. Our results predict that FAW will migrate from its new year-round breeding regions into the two main corn-producing regions of eastern China (Huang-Huai-Hai Summer Corn and Northeast Spring Corn Regions), via two pathways. The western pathway originates in Myanmar and Yunnan, and FAW will take four migration steps (i.e. four generations) to reach the Huang-Huai-Hai Region by July. Migration along the eastern pathway from Indochina and southern China progresses faster, with FAW reaching the Huang-Huai-Hai Region in three steps by June and reaching the Northeast Spring Region in July. CONCLUSION: Our results indicate that there is a high risk that FAW will invade the major corn-producing areas of eastern China via two migration pathways, and cause significant impacts to agricultural productivity. Information on migration pathways and timings can be used to inform integrated pest management strategies for this emerging pest. © 2019 Society of Chemical Industry.

The CWRF is developed as a climate extension of the Weather Research and Forecasting model (WRF) by incorporating numerous improvements in the representation of physical processes and integration of external (top, surface, lateral) forcings that are crucial to climate scales, including interactions between land, atmosphere, and ocean; convection and microphysics; and cloud, aerosol, and radiation; and system consistency throughout all process modules. This extension inherits all WRF functionalities for numerical weather prediction while enhancing the capability for climate modeling. As such, CWRF can be applied seamlessly to weather forecast and climate prediction. The CWRF is built with a comprehensive ensemble of alternative parameterization schemes for each of the key physical processes, including surface (land, ocean), planetary boundary layer, cumulus (deep, shallow), microphysics, cloud, aerosol, and radiation, and their interactions. This facilitates the use of an optimized physics ensemble approach to improve weather or climate prediction along with a reliable uncertainty estimate. The CWRF also emphasizes the societal service capability to provide impactrelevant information by coupling with detailed models of terrestrial hydrology, coastal ocean, crop growth, air quality, and a recently expanded interactive water quality and ecosystem model. This study provides a general CWRF description and basic skill evaluation based on a continuous integration for the period 1979– 2009 as compared with that of WRF, using a 30-km grid spacing over a domain that includes the contiguous United States plus southern Canada and northern Mexico. In addition to advantages of greater application capability, CWRF improves performance in radiation and terrestrial hydrology over WRF and other regional models. Precipitation simulation, however, remains a challenge for all of the tested models.

Abstract The global physical and biogeochemical environment has been substantially altered in response to increased atmospheric greenhouse gases from human activities. In 2023, the sea surface temperature (SST) and upper 2000 m ocean heat content (OHC) reached record highs. The 0–2000 m OHC in 2023 exceeded that of 2022 by 15 ± 10 ZJ (1 Zetta Joules = 10 21 Joules) (updated IAP/CAS data); 9 ± 5 ZJ (NCEI/NOAA data). The Tropical Atlantic Ocean, the Mediterranean Sea, and southern oceans recorded their highest OHC observed since the 1950s. Associated with the onset of a strong El Niño, the global SST reached its record high in 2023 with an annual mean of ∼0.23°C higher than 2022 and an astounding > 0.3°C above 2022 values for the second half of 2023. The density stratification and spatial temperature inhomogeneity indexes reached their highest values in 2023.

Abstract Typhoons Hato (2017) and Typhoon Mangkhut (2018) struck the Pearl River Delta, China, in two consecutive years, causing intense flooding and severe damage along the coast, especially in Macau. These two typhoons are extreme events because of their record‐breaking storm tides, and they are typical typhoons in terms of their track forward directions, which affect the coast of Guangdong Province, China. This paper analyzes and compares the physical processes of these two major typhoons. We first report the results of our posttyphoon surveys in Macau and the constructed inundation maps for both typhoon events. The inundation maps in conjunction with measured wind, atmospheric pressure, and wave data are used to validate a package of numerical models that simulate fully coupled tide‐surge‐wave fields during a typhoon. We quantify storm surges and wave heights in the Pearl River Delta. The surge around Macau was mainly affected by wind, while atmosphere pressure and ocean waves played a secondary role. In additional to inundation depth, we also investigate inundation durations and flow velocities. Typhoon Mangkhut made its landfall at much lower tidal level than that of Typhoon Hato. However, Typhoon Mangkhut generated much larger storm surge with larger affected areas and significantly longer inundation durations. We attribute the causes mainly to Typhoon Mangkhut's storm size, which was almost three times larger than that of Typhoon Hato.

Atmospheric concentrations of aluminium, an indicator of dust substances, have been determined in a set of high-volume aerosol particle samples collected at different locations over continental China and over the China Sea. High concentrations of dust were observed in northern continental China, and at certain locations such as Beijing dust may include an anthropogenic fraction. The mass particle-size distributions of dust varied depending on its distance from source regions, with the mass median diameter for Al of ~1.6-5.9 μm at Beijing in northern China and ~ 1.9 μm over off-shore areas of the East China Sea. Model-predicted mean dry deposition velocities of dust particles are from 1.4 to 4.8 cm s-1 over northern continental China and from 1.4 to 2.1 cm s-1 over the China Sea. Atmospheric deposition models have been applied to estimate the atmospheric fluxes and deposition of dust at different locations. The estimated atmospheric flux of dust at Xi’an of the Loess Plateau is 25 (4.9 to 44) g m-2 mo-1 which is the highest among the regions we studied. The estimated present-day dust flux is comparable to the late quaternary records of eolian dust accumulation at this site. The total atmospheric deposition of dust to the China Sea is 67 Tg yr-1, accounting for 14% of the total atmospheric deposition of dust to the entire North Pacific. With such a high deposition rate, Asian dust may play an important rôle in biogeochemical cycles of trace substances in the Asia/North Pacific region.

A teleconnection between the North Atlantic Ocean and the Eurasian continent is suggested by statistical and dynamical analysis of the northern summer 500 hPa geopotential height field. This teleconnection, termed the Atlantic–Eurasian (AEA) teleconnection, has five centers of action, in the subtropical North Atlantic Ocean, northeastern North Atlantic Ocean, Eastern Europe, the Kara Sea, and north China. The AEA index (AEAI) shows that the AEA undergoes a high degree of variability from year to year, and the AEAI has an increasing trend over the last 30 years. Our results suggest that this phenomenon is a large-scale Rossby wave train that originates in the subtropical North Atlantic Ocean. We support this conclusion by the methods of stationary wave ray tracing in non-uniform horizontal basic flow, wave activity flux calculations, and numerical models. The AEA and midlatitude circumglobal teleconnection pattern manifest distinct features at the hemispheric scale, despite the anomalies associated with them bear some similarities in the northeastern North Atlantic and Eastern Europe. Regional climate variations are strongly linked to this AEA along its path through northern Eurasia.

During the last 15 years, operational oceanography systems have emerged in several countries around the world. This emergence has been largely fostered by the GODAE experiment, during which each nation engaged in this activity have organised partnership and constructive competition. This trans-national coordination was very beneficial for the development of operational oceanography, leading to economies of scales and more targeted actions. Today, several systems provide routine real-time ocean analysis and forecast and/or reanalysis products. They are all based on (i) state-of-the-art primitive equation baroclinic Ocean General Circulation Model (OGCM) configurations, either global or regional (basin-scale), with resolutions that range from coarse to eddy resolving and (ii) data assimilation techniques whose complexity ranges from simple analysis correction to advanced 4D variational schemes. They assimilate altimeter sea level anomalies, remotely sensed SST such as GHRSST products and in situ profiles of T and S, including ARGO. Some systems have implemented downscaling capacities in specific regions of interest including shelf/coastal seas. Some also have implemented coupling

Operational evolution of global and regional ocean forecasting systems has been extremely significant in recent years. Global Ocean Data Assimilation Experiment (GODAE) Oceanview supports the national research groups providing them with coordination and sharing expertise among the partners. Several systems have been set up and developed pre-operationally, and the majority of these are now fully operational; at the present time, they provide medium- and long-term forecasts of the most relevant ocean physical variables. These systems are based on ocean general circulation models and data-assimilation techniques that are able to correct the model with the information inferred from different types of observations. A few systems also incorporate a biogeochemical component coupled with the physical system, while others are based on coupled ocean–wave–ice–atmosphere models. The products are routinely validated with observations in order to assess their quality. Data and product implementation and organization, as well as service, for users have been well tried and tested, and most of the products are now available to users. The interaction with different users is an important factor in the development process. This paper provides a synthetic overview of the GODAE OceanView prediction systems.

Abstract The impact of assimilating sea ice thickness data derived from ESA's Soil Moisture and Ocean Salinity (SMOS) satellite together with Special Sensor Microwave Imager/Sounder (SSMIS) sea ice concentration data of the National Snow and Ice Data Center (NSIDC) in a coupled sea ice‐ocean model is examined. A period of 3 months from 1 November 2011 to 31 January 2012 is selected to assess the forecast skill of the assimilation system. The 24 h forecasts and longer forecasts are based on the Massachusetts Institute of Technology general circulation model (MITgcm), and the assimilation is performed by a localized Singular Evolutive Interpolated Kalman (LSEIK) filter. For comparison, the assimilation is repeated only with the SSMIS sea ice concentrations. By running two different assimilation experiments, and comparing with the unassimilated model, independent satellite‐derived data, and in situ observation, it is shown that the SMOS ice thickness assimilation leads to improved thickness forecasts. With SMOS thickness data, the sea ice concentration forecasts also agree better with observations, although this improvement is smaller.

The deep structure of the continental collision between India and Asia and whether India's lower crust is underplated beneath Tibet or subducted into the mantle remain controversial. It is also unknown whether the active normal faults that facilitate orogen-parallel extension of Tibetan upper crust continue into the lower crust and upper mantle. Our receiver-function images collected parallel to the India-Tibet collision zone show the 20-km-thick Indian lower crust that underplates Tibet at 88.5-92°E beneath the Yarlung-Zangbo suture is essentially absent in the vicinity of the Cona-Sangri and Pumqu-Xainza grabens, demonstrating a clear link between upper-crustal and lower-crustal thinning. Satellite gravity data that covary with the thickness of Indian lower crust are consistent with the lower crust being only ∼30% eclogitized so gravitationally stable. Deep earthquakes coincide with Moho offsets and with lateral thinning of the Indian lower crust near the bottom of the partially eclogitized Indian lower crust, suggesting the Indian lower crust is locally foundering or stoping into the mantle. Loss of Indian lower crust by these means implies gravitational instability that can result from localized rapid eclogitization enabled by dehydration reactions in weakly hydrous mafic granulites or by volatile-rich asthenospheric upwelling directly beneath the two grabens. We propose that two competing processes, plateau formation by underplating and continental loss by foundering or stoping, are simultaneously operating beneath the collision zone.

Abstract The 168 year trends of summer (July–September) sea ice area (SIA) variations in six Arctic regions during 1850–2017 are analyzed. SIA has been significantly decreasing in most Arctic regions since 1850. The rate of retreat for the period of 1948–2017 accelerated multi-fold. For the nearly four decades since 1979, most Arctic regions are experiencing the highest reduction rate. Besides the increasing surface air temperature, the key drivers to the accelerated summer Arctic sea ice decline are found to be the combined global warming and the regional Arctic warming exerted simultaneously by the Arctic Oscillation, North Atlantic Oscillation, Atlantic Multidecadal Oscillation and Pacific Decadal Oscillation during the last several decades. The dynamical and thermodynamical warming, driven by the internal variability of the teleconnection patterns, occurred in the last several decades, in particular on the multidecadal timescales. This leads to Arctic amplification that accelerates the positive ice/ocean albedo feedback loop, resulting in accelerating summer sea ice decline.

Abstract Changes in ocean heat content (OHC), salinity, and stratification provide critical indicators for changes in Earth’s energy and water cycles. These cycles have been profoundly altered due to the emission of greenhouse gasses and other anthropogenic substances by human activities, driving pervasive changes in Earth’s climate system. In 2022, the world’s oceans, as given by OHC, were again the hottest in the historical record and exceeded the previous 2021 record maximum. According to IAP/CAS data, the 0–2000 m OHC in 2022 exceeded that of 2021 by 10.9 ± 8.3 ZJ (1 Zetta Joules = 10 21 Joules); and according to NCEI/NOAA data, by 9.1 ± 8.7 ZJ. Among seven regions, four basins (the North Pacific, North Atlantic, the Mediterranean Sea, and southern oceans) recorded their highest OHC since the 1950s. The salinity-contrast index, a quantification of the “salty gets saltier—fresh gets fresher” pattern, also reached its highest level on record in 2022, implying continued amplification of the global hydrological cycle. Regional OHC and salinity changes in 2022 were dominated by a strong La Niña event. Global upper-ocean stratification continued its increasing trend and was among the top seven in 2022.

The frequent occurrences of the second-year surface cooling condition in the eastern equatorial Pacific, as observed in late 2021, are attributed to decadal changes in the thermocline depth, which determine the relative dominances of local cooling effect in the east and subsurface warming effect remotely from the west. Coupled models need to adequately represent these processes in a balanced way, thus being able to successfully predict the observed sea surface temperature evolution in late 2021.

An extended Stochastic Impacts by Regression on Population, Affluence and Technology (STIRPAT) model, incorporating factors that drive carbon emissions, is built from the regional perspective. A spatial Durbin model is applied to investigate the factors, including population, urbanization level, economic development, energy intensity, industrial structure, energy consumption structure, energy price, and openness, that impact both the scale and intensity of carbon emissions. After performing the model, we find that the revealed negative and significant impact of spatial-lagged variables suggests that the carbon emissions among regions are highly correlated. Therefore, the empirical results suggest that the provinces are doing an exemplary job of lowering carbon emissions. The driving factors, with the exception of energy prices, significantly impact carbon emissions both directly and indirectly. We, thus, argue that spatial correlation, endogeneity and externality should be taken into account in formulating polices that seek to reduce carbon emissions in China. Carbon emissions will not be met by controlling economic development, but by energy consumption and low-carbon path.

The world's largest macroalgal blooms caused by Ulva prolifera have occurred in the Yellow Sea for 11 consecutive years. The area covered by blooms has been approximately 500 km2 in previous years, while in 2017, the maximum area decreased significantly to 312 km2. In this study, we concluded that species competition between Ulva and Sargassum (fast rise of the golden tides), extreme high sea surface temperature and harvest for floating Ulva macroalgae were the three critical factors influencing the sharp reduction in covered area for blooms in 2017. In addition, analysis of annual variations of Pyropia aquaculture area in the Southern Yellow Sea over the past two decades revealed that a great expansion in "Sansha" regions was mainly responsible for the initial blooms in 2007, and that this expansion supported the great biomass of the blooms in following years. Based on these findings, we suggest comprehensive utilization of the macroalgal blooms is a feasible way to control them.