Air-Sea Interactions During the Passage of a Winter Storm Over the Gulf Stream: A Three-Dimensional Coupled Atmosphere-Ocean Model Study

A three-dimensional, regional coupled atmosphere-ocean model with full physics is developed to study air-sea interactions during winter storms off the U. S. east coast. Because of the scarcity of open ocean observations, models such as this offer valuable opportunities to investigate how oceanic forcing drives atmospheric circulation and vice versa. The study presented here considers conditions of strong atmospheric forcing (high wind speeds) and strong oceanic forcing (significant sea surface temperature (SST) gradients). A simulated atmospheric cyclone evolves in a manner consistent with Eta reanalysis, and the simulated air-sea heat and momentum exchanges strongly affect the circulations in both the atmosphere and the ocean. For the simulated cyclone of 19-20 January 1998, maximum ocean-to-atmosphere heat fluxes first appear over the Gulf Stream in the South Atlantic Bight, and this results in rapid deepening of the cyclone off the Carolina coast. As the cyclone moves eastward, the heat flux maximum shifts into the region near Cape Hatteras and later northeast of Hatteras, where it enhances the wind locally. The oceanic response to the atmospheric forcing is closely related to the wind direction. Southerly and southwesterly winds tend to strengthen surface currents in the Gulf Stream, whereas northeasterly winds weaken the surface currents in the Gulf Stream and generate southwestward flows on the shelf. The oceanic feedback to the atmosphere moderates the cyclone strength. Compared with a simulation in which the oceanic model always passes the initial SST to the atmospheric model, the coupled simulation in which the oceanic model passes the evolving SST to the atmospheric model produces higher ocean-to-atmosphere heat flux near Gulf Stream meander troughs. This is due to wind-driven lateral shifts of the stream, which in turn enhance the local northeasterly winds. Away from the Gulf Stream the coupled simulation produces surface winds that are 5 similar to 10% weaker. Differences in the surface ocean currents between these two experiments are significant on the shelf and in the open ocean.

A Model Study of the Seasonal Circulation in the Gulf of Maine

The Princeton Ocean Model is used to study the circulation in the Gulf of Maine and its seasonal transition in response to wind, surface heat flux, river discharge, and the M-2 tide. The model has an orthogonal-curvature linear grid in the horizontal with variable spacing from 3 km nearshore to 7 km offshore and 19 levels in the vertical. It is initialized and forced at the open boundary with model results from the East Coast Forecast System. The first experiment is forced by monthly climatological wind and heat flux from the Comprehensive Ocean Atmosphere Data Set; discharges from the Saint John, Penobscot, Kennebec, and Merrimack Rivers are added in the second experiment; the semidiurnal lunar tide (M-2) is included as part of the open boundary forcing in the third experiment. It is found that the surface heat flux plays an important role in regulating the annual cycle of the circulation in the Gulf of Maine. The spinup of the cyclonic circulation between April and June is likely caused by the differential heating between the interior gulf and the exterior shelf/slope region. From June to December the cyclonic circulation continues to strengthen, but gradually shrinks in size. When winter cooling erodes the stratification, the cyclonic circulation penetrates deeper into the water column. The circulation quickly spins down from December to February as most of the energy is consumed by bottom friction. While inclusion of river discharge changes details of the circulation pattern, the annual evolution of the circulation is largely unaffected. On the other hand, inclusion of the tide results in not only the anticyclonic circulation on Georges Bank but also modifications to the seasonal circulation.

A Model Study of the Circulation in the Pearl River Estuary (PRE) and Its Adjacent Coastal Waters: 1. Simulations and Comparison with Observations

The Princeton Ocean Model is used to study the circulation features in the Pearl River Estuary and their responses to tide, river discharge, wind, and heat flux in the winter dry and summer wet seasons. The model has an orthogonal curvilinear grid in the horizontal plane with variable spacing from 0.5 km in the estuary to 1 km on the shelf and 15 sigma levels in the vertical direction. The initial conditions and the subtidal open boundary forcing are obtained from an associated larger-scale model of the northern South China Sea. Buoyancy forcing uses the climatological monthly heat fluxes and river discharges, and both the climatological monthly wind and the realistic wind are used in the sensitivity experiments. The tidal forcing is represented by sinusoidal functions with the observed amplitudes and phases. In this paper, the simulated tide is first examined. The simulated seasonal distributions of the salinity, as well as the temporal variations of the salinity and velocity over a tidal cycle are described and then compared with the in situ survey data from July 1999 and January 2000. The model successfully reproduces the main hydrodynamic processes, such as the stratification, mixing, frontal dynamics, summer upwelling, two-layer gravitational circulation, etc., and the distributions of hydrodynamic parameters in the Pearl River Estuary and coastal waters for both the winter and the summer season.

A Model Study of the Circulation in the Pearl River Estuary (PRE) and Its Adjacent Coastal Waters: 2. Sensitivity Experiments

The Princeton Ocean Model is used to study the circulation in the Pear River Estuary (PRE) and the adjacent coastal waters in the winter and summer seasons. Wong et al. [2003] compares the simulation results with the in situ measurements collected during the Pearl River Estuary Pollution Project (PREPP). In this paper, sensitivity experiments are carried out to examine the plume and the associated frontal dynamics in response to seasonal discharges and monsoon winds. During the winter, convergence between the seaward spreading plume water and the saline coastal water sets up a salinity front that aligns from the northeast to the southwest inside the PRE. During the summer the plume water fills the PRE at the surface and spreads eastward in the coastal waters in response to the prevailing southwesterly monsoon. The overall alignment of the plume is from the northwest to the southeast. The subsurface front is similar to that in the winter and summer except that the summer front is closer to the mouth and the winter front closer to the head of the estuary. Inside the PRE, bottom flows are always toward the head of the estuary, attributed to the density gradient associated with the plume front. In contrast, bottom flows in the shelf change from offshore in winter to onshore in summer, reflecting respectively the wintertime downwelling and summertime upwelling. Wind also plays an essential role in controlling the plume at the surface. An easterly wind drives the plume westward regardless winter or summer. The eastward spreading of the plume during the summer can be attributed to the southerly component of the wind. On the other hand, the surface area of the plume is positively proportional to the amount of discharge.

Three time-slice simulations of the Eemian (eem), the mid-Holocene (hol), and a pre-industrial (pri) control run with the coupled atmosphere-ocean-model ECHAM5/MPI-OM

Orbital forcing does not only exert direct insolation effects, but also alters climate indirectly through feedback mechanisms that modify atmosphere and ocean dynamics and meridional heat and moisture transfers. We investigate the regional effects of these changes by detailed analysis of atmosphere and ocean circulation and heat transports in a coupled atmosphere-ocean-sea ice-biosphere general circulation model (ECHAM5/JSBACH/MPI-OM). We perform long term quasi equilibrium simulations under pre-industrial, mid-Holocene (6000 years before present - yBP), and Eemian (125 000 yBP) orbital boundary conditions. Compared to pre-industrial climate, Eemian and Holocene temperatures show generally warmer conditions at higher and cooler conditions at lower latitudes. Changes in sea-ice cover, ocean heat transports, and atmospheric circulation patterns lead to pronounced regional heterogeneity. Over Europe, the warming is most pronounced over the north-eastern part in accordance with recent reconstructions for the Holocene. We attribute this warming to enhanced ocean circulation in the Nordic Seas and enhanced ocean-atmosphere heat flux over the Barents Shelf in conduction with retreat of sea ice and intensified winter storm tracks over northern Europe.

Results of a Finite-Element Sea-Ice Ocean Model (FESOM) set-up to study the interannual to decadal variability in the deep-water formation rates

The characteristics of a global set-up of the Finite-Element Sea-Ice Ocean Model under forcing of the period 1958-2004 are presented. The model set-up is designed to study the variability in the deep-water mass formation areas and was therefore regionally better resolved in the deep-water formation areas in the Labrador Sea, Greenland Sea, Weddell Sea and Ross Sea. The sea-ice model reproduces realistic sea-ice distributions and variabilities in the sea-ice extent of both hemispheres as well as sea-ice transport that compares well with observational data. Based on a comparison between model and ocean weather ship data in the North Atlantic, we observe that the vertical structure is well captured in areas with a high resolution. In our model set-up, we are able to simulate decadal ocean variability including several salinity anomaly events and corresponding fingerprint in the vertical hydrography. The ocean state of the model set-up features pronounced variability in the Atlantic Meridional Overturning Circulation as well as the associated mixed layer depth pattern in the North Atlantic deep-water formation areas.

On the effects of circulation, sediment resuspension and biological incorporation by diatoms in an ocean model of aluminium, link to model data

The distribution of dissolved aluminium in the West Atlantic Ocean shows a mirror image with that of dissolved silicic acid, hinting at intricate interactions between the ocean cycling of Al and Si. The marine biogeochemistry of Al is of interest because of its potential impact on diatom opal remineralisation, hence Si availability. Furthermore, the dissolved Al concentration at the surface ocean has been used as a tracer for dust input, dust being the most important source of the bio-essential trace element iron to the ocean. Previously, the dissolved concentration of Al was simulated reasonably well with only a dust source, and scavenging by adsorption on settling biogenic debris as the only removal process. Here we explore the impacts of (i) a sediment source of Al in the Northern Hemisphere (especially north of ~ 40° N), (ii) the imposed velocity field, and (iii) biological incorporation of Al on the modelled Al distribution in the ocean. The sediment source clearly improves the model results, and using a different velocity field shows the importance of advection on the simulated Al distribution. Biological incorporation appears to be a potentially important removal process. However, conclusive independent data to constrain the Al / Si incorporation ratio by growing diatoms are missing. Therefore, this study does not provide a definitive answer to the question of the relative importance of Al removal by incorporation compared to removal by adsorptive scavenging.

Regional ocean model for the Mediterranean, links to model files

A regional ocean general circulation model of the Mediterranean is used to study the climate of the Last Glacial Maximum. The atmospheric forcing for these simulations has been derived from simulations with an atmospheric general circulation model, which in turn was forced with surface conditions from a coarse resolution earth system model. The model is successful in reproducing the general patterns of reconstructed sea surface temperature anomalies with the strongest cooling in summer in the northwestern Mediterranean and weak cooling in the Levantine, although the model underestimates the extent of the summer cooling in the western Mediterranean. However, there is a strong vertical gradient associated with this pattern of summer cooling, which makes the comparison with reconstructions complicated. The exchange with the Atlantic is decreased to roughly one half of its present value, which can be explained by the shallower Strait of Gibraltar as a consequence of lower global sea level. This reduced exchange causes a strong increase of salinity in the Mediterranean in spite of reduced net evaporation.

Efficient parallelization of a regional ocean model for the western Mediterranean Sea

This paper focuses on the parallelization of an ocean model applying current multicore processor-based cluster architectures to an irregular computational mesh. The aim is to maximize the efficiency of the computational resources used. To make the best use of the resources offered by these architectures, this parallelization has been addressed at all the hardware levels of modern supercomputers: firstly, exploiting the internal parallelism of the CPU through vectorization; secondly, taking advantage of the multiple cores of each node using OpenMP; and finally, using the cluster nodes to distribute the computational mesh, using MPI for communication within the nodes. The speedup obtained with each parallelization technique as well as the combined overall speedup have been measured for the western Mediterranean Sea for different cluster configurations, achieving a speedup factor of 73.3 using 256 processors. The results also show the efficiency achieved in the different cluster nodes and the advantages obtained by combining OpenMP and MPI versus using only OpenMP or MPI. Finally, the scalability of the model has been analysed by examining computation and communication times as well as the communication and synchronization overhead due to parallelization.

Holocene experiment with coupled atmosphere-ocean-model ECHAM5/MPI-OM

Changes in the Earth's orbit lead to changes in the seasonal and meridional distribution of insolation. We quantify the influence of orbitally induced changes on the seasonal temperature cycle in a transient simulation of the last 6000 years - from the mid-Holocene to today - using a coupled atmosphere-ocean general circulation model (ECHAM5/MPI-OM) including a land surface model (JSBACH). The seasonal temperature cycle responds directly to the insolation changes almost everywhere. In the Northern Hemisphere, its amplitude decreases according to an increase in winter insolation and a decrease in summer insolation. In the Southern Hemisphere, the opposite is true. Over the Arctic Ocean, decreasing summer insolation leads to an increase in sea-ice cover. The insulating effect of sea ice between the ocean and the atmosphere leads to decreasing heat flux and favors more "continental" conditions over the Arctic Ocean in winter, resulting in strongly decreasing temperatures. Consequently, there are two competing effects: the direct response to insolation changes and a sea-ice insulation effect. The sea-ice insulation effect is stronger, and thus an increase in the amplitude of the seasonal temperature cycle over the Arctic Ocean occurs. This increase is strongest over the Barents Shelf and influences the temperature response over northern Europe. We compare our modeled seasonal temperatures over Europe to paleo reconstructions. We find better agreements in winter temperatures than in summer temperatures and better agreements in northern Europe than in southern Europe, since the model does not reproduce the southern European Holocene summer cooling inferred from the paleo reconstructions. The temperature reconstructions for northern Europe support the notion of the influence of the sea-ice insulation effect on the evolution of the seasonal temperature cycle.

Experimental and diagnostic protocol for the physical component of the CMIP6 Ocean Model Intercomparison Project (OMIP)

The Ocean Model Intercomparison Project (OMIP) aims to provide a framework for evaluating, understanding, and improving the ocean and sea-ice components of global climate and earth system models contributing to the Coupled Model Intercomparison Project Phase 6 (CMIP6). OMIP addresses these aims in two complementary manners: (A) by providing an experimental protocol for global ocean/sea-ice models run with a prescribed atmospheric forcing, (B) by providing a protocol for ocean diagnostics to be saved as part of CMIP6. We focus here on the physical component of OMIP, with a companion paper (Orr et al., 2016) offering details for the inert chemistry and interactive biogeochemistry. The physical portion of the OMIP experimental protocol follows that of the interannual Coordinated Ocean-ice Reference Experiments (CORE-II). Since 2009, CORE-I (Normal Year Forcing) and CORE-II have become the standard method to evaluate global ocean/sea-ice simulations and to examine mechanisms for forced ocean climate variability. The OMIP diagnostic protocol is relevant for any ocean model component of CMIP6, including the DECK (Diagnostic, Evaluation and Characterization of Klima experiments), historical simulations, FAFMIP (Flux Anomaly Forced MIP), C4MIP (Coupled Carbon Cycle Climate MIP), DAMIP (Detection and Attribution MIP), DCPP (Decadal Climate Prediction Project), ScenarioMIP (Scenario MIP), as well as the ocean-sea ice OMIP simulations. The bulk of this paper offers scientific rationale for saving these diagnostics.