A three-dimensional coupled physical-biological model study in the spring of 1993 in the Bohai Sea of China

A three-dimensional (3-D) coupled physical and biological model was used to investigate the physical processes and their influence on the ecosystem dynamics of the Bohai Sea of China. The physical processes include M-2 tide, time - varying wind forcing and river discharge. Wind records from I to 31 May in 1993 were selected to force the model. The biological model is based on a simple, nitrate and phosphate limited, lower trophic food web system. The simulated results showed that variation of residual currents forced by M, tide, river discharge and time-varying wind had great impact on the distribution of phytoplankton biomass in the Laizhou Bay. High phytoplankton biomass appeared in the upwelling region. Numerical experiments based on the barotropic model and baroclinic model with no wind and water discharge were also conducted. Differences in the results by the baroclinic model and the barotropic model were significant: more patches appeared in the baroclinic model comparing with the barotropic model. And in the baroclinic model, the subsurface maximum phytoplankton biomass patches formed in the stratified water.

Have climate induced changes in zooplankton communities led to poor conditions for larval fish growth? A test on cod larvae on the UK shelf

Climate effects have been shown to be at least partly responsible for the reorganisation in the plankton ecosystem on the shelf seas of NW Europe over the last 50 years. Most fish larvae feed primarily on zooplankton, so changes in zooplankton quantity, quality and seasonal timing have been hypothesized to be a key factor affecting their survival. To investigate this we have implemented a 1-dimensional trophodynamic growth model of cod larvae for the waters around the UK covering the period 1960 to 2003. Larval growth is modelled as the difference between the amount of food absorbed by the larva and its various metabolic costs. Prey availability is based upon the biomass and size of available preys (i.e. adults and nauplii copepods and cladocerans) taken from the Continuous Plankton Recorder dataset. Temperature and wind forcing are also taken into account. Results suggest that observed changes in plankton community structure may have had less impact than previously suggested. This is because changes in prey availability may be compensated for by increased temperatures resulting in little overall impact on potential larval growth. Stock recovery, at least in the short term is likely to be more dependent upon conserving the year classes recruited to allow spawning stock biomass to rebuild. If as our model suggests, the larvae are still able to survive in the changing environment, reduction in fishing on the adults is needed to allow the stock to recover.

Sea-level changes in Iceland and the influence of the North Atlantic Oscillation during the last half millennium

We present a new, diatom-based sea-level reconstruction for Iceland spanning the last -500 years, and investigate the possible mechanisms driving the sea-level changes. A sea-level reconstruction from near the Icelandic low pressure system is important as it can improve understanding of ocean-atmosphere forcing on North Atlantic sea-level variability over multi-decadal to centennial timescales. Our reconstruction is from Viarhólmi salt marsh in Snæfellsnes in western Iceland, a site from where we previously obtained a 2000-yr record based upon less precise sea-level indicators (salt-marsh foraminifera). The 20th century part of our record is corroborated by tide-gauge data from Reykjavik. Overall, the new reconstruction shows ca0.6m rise of relative sea level during the last four centuries, of which ca0.2m occurred during the 20th century. Low-amplitude and high-frequency sea-level variability is super-imposed on the pre-industrial long-term rising trend of 0.65m per 1000 years. Most of the relative sea-level rise occurred in three distinct periods: AD 1620-1650, AD 1780-1850 and AD 1950-2000, with maximum rates of ~3±2mm/yr during the latter two of these periods. Maximum rates were achieved at the end of large shifts (from negative to positive) of the winter North Atlantic Oscillation (NAO) Index as reconstructed from proxy data. Instrumental data demonstrate that a strong and sustained positive NAO (a deep Icelandic Low) generates setup on the west coast of Iceland resulting in rising sea levels. There is no strong evidence that the periods of rapid sea-level rise were caused by ocean mass changes, glacial isostatic adjustment or regional steric change. We suggest that wind forcing plays an important role in causing regional-scale coastal sea-level variability in the North Atlantic, not only on (multi-)annual timescales, but also on multi-decadal to centennial timescales.

Modelling the Minho river plume intrusion into the rias baixas

The Minho River, situated 30 km south of the Rias Baixas is the most important freshwater source flowing into the Western Galician Coast (NW of the Iberian Peninsula). This discharge is important to determine the hydrological patterns adjacent to its mouth, particularly close to the Galician coastal region. The buoyancy generated by the Minho plume can flood the Rias Baixas for long periods, reversing the normal estuarine density gradients. Thus, it becomes important to analyse its dynamics as well as the thermohaline patterns of the areas affected by the freshwater spreading. Thus, the main aim of this work was to study the propagation of the Minho estuarine plume to the Rias Baixas, establishing the conditions in which this plume affects the circulation and hydrographic features of these coastal systems, through the development and application of the numerical model MOHID. For this purpose, the hydrographic features of the Rias Baixas mouths were studied. It was observed that at the northern mouths, due to their shallowness, the heat fluxes between the atmosphere and ocean are the major forcing, influencing the water temperature, while at the southern mouths the influence of the upwelling events and the Minho River discharge were more frequent. The salinity increases from south to north, revealing that the observed low values may be caused by the Minho River freshwater discharge. An assessment of wind data along the Galician coast was carried out, in order to evaluate the applicability of the study to the dispersal of the Minho estuarine plume. Firstly, a comparative analysis between winds obtained from land meteorological stations and offshore QuikSCAT satellite were performed. This comparison revealed that satellite data constitute a good approach to study wind induced coastal phenomena. However, since the numerical model MOHID requires wind data with high spatial and temporal resolution close to the coast, results of the forecasted model WRF were added to the previous study. The analyses revealed that the WRF model data is a consistent tool to obtain representative wind data near the coast, showing good results when comparing with in situ wind observations from oceanographic buoys. To study the influence of the Minho buoyant discharge influence on the Rias Baixas, a set of three one-way nested models was developed and implemented, using the numerical model MOHID. The first model domain is a barotropic model and includes the whole Iberian Peninsula coast. The second and third domains are baroclinic models, where the second domain is a coarse representation of the Rias Baixas and adjacent coastal area, while the third includes the same area with a higher resolution. A bi-dimensional model was also implemented in the Minho estuary, in order to quantify the flow (and its properties) that the estuary injects into the ocean. The chosen period for the Minho estuarine plume propagation validation was the spring of 1998, since a high Minho River discharge was reported, as well as favourable wind patterns to advect the estuarine plume towards the Rias Baixas, and there was field data available to compare with the model predictions. The obtained results show that the adopted nesting methodology was successful implemented. Model predictions reproduce accurately the hydrodynamics and thermohaline patterns on the Minho estuary and Rias Baixas. The importance of the Minho river discharge and the wind forcing in the event of May 1998 was also studied. The model results showed that a continuous moderate Minho River discharge combined with southerly winds is enough to reverse the Rias Baixas circulation pattern, reducing the importance of the occurrence of specific events of high runoff values. The conditions in which the estuarine plume Minho affects circulation and hydrography of the Rias Baixas were evaluated. The numerical results revealed that the Minho estuarine plume responds rapidly to wind variations and is also influenced by the bathymetry and morphology of the coastline. Without wind forcing, the plume expands offshore, creating a bulge in front of the river mouth. When the wind blows southwards, the main feature is the offshore extension of the plume. Otherwise, northward wind spreads the river plume towards the Rias Baixas. The plume is confined close to the coast, reaching the Rias Baixas after 1.5 days. However, for Minho River discharges higher than 800 m3 s-1, the Minho estuarine plume reverses the circulation patterns in the Rias Baixas. It was also observed that the wind stress and Minho River discharge are the most important factors influencing the size and shape of the Minho estuarine plume. Under the same conditions, the water exchange between Rias Baixas was analysed following the trajectories particles released close to the Minho River mouth. Over 5 days, under Minho River discharges higher than 2100 m3 s-1 combined with southerly winds of 6 m s-1, an intense water exchange between Rias was observed. However, only 20% of the particles found in Ria de Pontevedra come directly from the Minho River. In summary, the model application developed in this study contributed to the characterization and understanding of the influence of the Minho River on the Rias Baixas circulation and hydrography, highlighting that this methodology can be replicated to other coastal systems.

Ocean Boundary Layer Dynamics and Air-Sea Interaction

Thesis (Ph.D.)--University of Washington, 2005-12

Inertial oscillation forced by the September 1997 cyclone in the Bay of Bengal

Time-series measurements from a moored data buoy located in the Bay of Bengal captured signals of inertial oscillation forced by the September 1997 cyclone. The progressive vector diagram showed mean northeastward current with well-defined clockwise circulation. Spectral analysis exhibited inertial peak at 0.67 cpd with blue shift and high rotary coefficient of –0.99, which signifies strong circular inertial oscillation. The wind and SST also exhibited spectral peak at inertial band (0.69 cpd) with higher blue shift. The inertial amplitude of 148.8 cm/s corresponding to a wind stress of 0.99 N/m2 and spectral peak near the local inertial frequency (0.653 cpd) indicate that the transfer of momentum was high.

Coastal Upwelling of the South Eastern Arabian Sea - An Integrated Approach

Upwelling regions occupies only a small portion of the global ocean surface. However it accounts for a large fraction of the oceanic primary production as well as fishery. Therefore understanding and quantifying the upwelling is of great importance for the marine resources management. Most of the coastal upwelling zones in the Arabian Sea are wind driven uniform systems. Mesoscale studies along the southwest coast of India have shown high spatial and temporal variability in the forcing mechanism and intensity of upwelling. There exists an equatorward component of wind stress as similar to the most upwelling zones along the eastern oceanic boundaries. Therefore an offshore component of surface Ekman transport is expected throughout the year. But several studies supported with in situ evidences have revealed that the process is purely recurring on seasonal basis. The explanation merely based on local wind forcing alone is not sufficient to support the observations. So, it is assumed that upwelling along the South Eastern Arabian Sea is an effect of basin wide wind forcing rather than local wind forcing. In the present study an integrated approach has been made to understand the process of upwelling of the South Eastern Arabian Sea. The latitudinal and seasonal variations (based on Sea Surface Temperature, wind forcing, Chlorophyll a and primary production), forcing mechanisms (local wind and remote forcing) and the factors influencing the system (Arabian Sea High Saline Water, Bay of Bengal water, runoff, coastal geomorphology) are addressed herewith.

Basin-scale hydrodynamics in a Mediterranean reservoir. Implications for the phytoplankton dynamics

Procesos hidrodinámicos determinan, en un alto grado la calidad del agua en embalse, sin embargo dichos procesos han sido tradicionalmente olvidados en la gestión de embalse. En esta tesis se presentan evidencias de los principales procesos hidrodinámicos que ocurren en un embalse Mediterráneo a escala de cuenca a través de campañas experimentales y modelización numérica; y su influencia en la dinámica de poblaciones de fitoplancton. Dichos procesos son principalmente la generación de ondas internas o secas y la intrusión del río. La presencia de viento periódico genera secas forzadas, amplificando los modos cercanos al periodo del viento, de manera que modos verticales altos, considerados como raros en la naturaleza, tienden a dominar en el sistema.

Comparison between observations and models of the Mozambique Channel transport: Seasonal cycle and eddy frequencies

A time series of the observed transport through an array of moorings across the Mozambique Channel is compared with that of six model runs with ocean general circulation models. In the observations, the seasonal cycle cannot be distinguished from red noise, while this cycle is dominant in the transport of the numerical models. It is found, however, that the seasonal cycles of the observations and numerical models are similar in strength and phase. These cycles have an amplitude of 5 Sv and a maximum in September, and can be explained by the yearly variation of the wind forcing. The seasonal cycle in the models is dominant because the spectral density at other frequencies is underrepresented. Main deviations from the observations are found at depths shallower than 1500 m and in the 5/y–6/y frequency range. Nevertheless, the structure of eddies in the models is close to the observed eddy structure. The discrepancy is found to be related to the formation mechanism and the formation position of the eddies. In the observations, eddies are frequently formed from an overshooting current near the mooring section, as proposed by Ridderinkhof and de Ruijter (2003) and Harlander et al. (2009). This causes an alternation of events at the mooring section, varying between a strong southward current, and the formation and passing of an eddy. This results in a large variation of transport in the frequency range of 5/y–6/y. In the models, the eddies are formed further north and propagate through the section. No alternation similar to the observations is observed, resulting in a more constant transport.

On the response of the Antarctic Circumpolar Current transport to climate change in coupled climate models

The CMIP3 (IPCC AR4) models show a consistent intensification and poleward shift of the westerly winds over the Southern Ocean during the 21st century. However, the responses of the Antarctic Circumpolar Currents (ACC) show great diversity in these models, with many even showing reductions in transport. To obtain some understanding of diverse responses in the ACC transport, we investigate both external atmospheric and internal oceanic processes that control the ACC transport responses in these models. While the strengthened westerlies act to increase the tilt of isopycnal surfaces and hence the ACC transport through Ekman pumping effects, the associated changes in buoyancy forcing generally tend to reduce the surface meridional density gradient. The steepening of isopycnal surfaces induced by increased wind forcing leads to enhanced (parameterized) eddy-induced transports that act to reduce the isopycnal slopes. There is also considerable narrowing of the ACC that tends to reduce the ACC transport, caused mainly by the poleward shifts of the subtropical gyres and to a lesser extent by the equatorward expansions of the subpolar gyres in some models. If the combined effect of these retarding processes is larger than that of enhanced Ekman pumping, the ACC transport will be reduced. In addition, the effect of Ekman pumping on the ACC is reduced in weakly stratified models. These findings give insight into the reliability of IPCC-class model predictions of the Southern Ocean circulation, and into the observed decadal-scale steady ACC transport.

Spin-up and adjustment of the Antarctic Circumpolar Current and global pycnocline

A theory is presented for the adjustment of the Antarctic Circumpolar Current (ACC) and global pycnocline to a sudden and sustained change in wind forcing. The adjustment timescale is controlled by the mesoscale eddy diffusivity across the ACC, the mean width of the ACC, the surface area of the ocean basins to the north, and deep water formation in the North Atlantic. In particular, northern sinking may have the potential to shorten the timescale and reduce its sensitivity to Southern Ocean eddies, but the relative importance of northern sinking and Southern Ocean eddies cannot be determined precisely, largely due to limitations in the parameterization of northern sinking. Although it is clear that the main processes that control the adjustment timescale are those which counteract the deepening of the global pycnocline, the theory also suggests that the timescale can be subtly modified by wind forcing over the ACC and global diapycnal mixing. Results from calculations with a reduced-gravity model compare well with the theory. The multidecadal-centennial adjustment timescale implies that long observational time series will be required to detect dynamic change in the ACC due to anthropogenic forcing. The potential role of Southern Ocean mesoscale eddy activity in determining both the equilibrium state of the ACC and the timescale over which it adjusts suggests that the response to anthropogenic forcing may be different in coupled ocean-atmosphere climate models that parameterize and resolve mesoscale eddies.

A global perspective on Langmuir turbulence in the ocean surface boundary layer.

The turbulent mixing in thin ocean surface boundary layers (OSBL), which occupy the upper 100 m or so of the ocean, control the exchange of heat and trace gases between the atmosphere and ocean. Here we show that current parameterizations of this turbulent mixing lead to systematic and substantial errors in the depth of the OSBL in global climate models, which then leads to biases in sea surface temperature. One reason, we argue, is that current parameterizations are missing key surface-wave processes that force Langmuir turbulence that deepens the OSBL more rapidly than steady wind forcing. Scaling arguments are presented to identify two dimensionless parameters that measure the importance of wave forcing against wind forcing, and against buoyancy forcing. A global perspective on the occurrence of waveforced turbulence is developed using re-analysis data to compute these parameters globally. The diagnostic study developed here suggests that turbulent energy available for mixing the OSBL is under-estimated without forcing by surface waves. Wave-forcing and hence Langmuir turbulence could be important over wide areas of the ocean and in all seasons in the Southern Ocean. We conclude that surfacewave- forced Langmuir turbulence is an important process in the OSBL that requires parameterization.

Ocean heat uptake processes: a model intercomparison

We compare the quasi-equilibrium heat balances, as well as their responses to 4×CO2 perturbation, among three global climate models with the aim to identify and explain inter-model differences in ocean heat uptake (OHU) processes. We find that, in quasi-equilibrium, convective and mixed layer processes, as well as eddy-related processes, cause cooling of the subsurface ocean. The cooling is balanced by warming caused by advective and diapycnally diffusive processes. We also find that in the CO2-perturbed climates the largest contribution to OHU comes from changes in vertical mixing processes and the mean circulation, particularly in the extra-tropics, caused both by changes in wind forcing, and by changes in high-latitude buoyancy forcing. There is a substantial warming in the tropics, a significant part of which occurs because of changes in horizontal advection in extra-tropics. Diapycnal diffusion makes only a weak contribution to the OHU, mainly in the tropics, due to increased stratification. There are important qualitative differences in the contribution of eddy-induced advection and isopycnal diffusion to the OHU among the models. The former is related to the different values of the coefficients used in the corresponding scheme. The latter is related to the different tapering formulations of the isopycnal diffusion scheme. These differences affect the OHU in the deep ocean, which is substantial in two of the models, the dominant region of deep warming being the Southern Ocean. However, most of the OHU takes place above 2000 m, and the three models are quantitatively similar in their global OHU efficiency and its breakdown among processes and as a function of latitude.

Coordinated ground-based, low altitude satellite and Cluster observations on global and local scales during a transient post-noon sector excursion of the magnetospheric cusp

On 14 January 2001, the four Cluster spacecraft passed through the northern magnetospheric mantle in close conjunction to the EISCAT Svalbard Radar (ESR) and approached the post-noon dayside magnetopause over Greenland between 13:00 and 14:00 UT During that interval, a sudden reorganisation of the high-latitude dayside convection pattern accurred after 13:20 UT most likely caused by a direction change of the Solar wind magnetic field. The result was an eastward and poleward directed flow-channel, as monitored by the SuperDARN radar network and also by arrays of ground-based magnetometers in Canada, Greenland and Scandinavia. After an initial eastward and later poleward expansion of the flow-channel between 13:20 and 13:40 UT, the four Cluster spacecraft, and the field line footprints covered by the eastward looking scan cycle of the Sondre Stromfjord incoherent scatter radar were engulfed by cusp-like precipitation with transient magnetic and electric field signatures. In addition, the EISCAT Svalbard Radar detected strong transient effects of the convection reorganisation, a poleward moving precipitation, and a fast ion flow-channel in association with the auroral structures that suddenly formed to the west and north of the radar. From a detailed analysis of the coordinated Cluster and ground-based data, it was found that this extraordinary transient convection pattern, indeed, had moved the cusp precipitation from its former pre-noon position into the late post-noon sector, allowing for the first and quite unexpected encounter of the cusp by the Cluster spacecraft. Our findings illustrate the large amplitude of cusp dynamics even in response to moderate solar wind forcing. The global ground-based data proves to be an invaluable tool to monitor the dynamics and width of the affected magnetospheric regions.

An assessment of Indian monsoon seasonal forecasts and mechanisms underlying monsoon interannual variability in the Met Office GloSea5-GC2 system

We assess Indian summer monsoon seasonal forecasts in GloSea5-GC2, the Met Office fully coupled subseasonal to seasonal ensemble forecasting system. Using several metrics, GloSea5-GC2 shows similar skill to other state-of-the-art forecast systems. The prediction skill of the large-scale South Asian monsoon circulation is higher than that of Indian monsoon rainfall. Using multiple linear regression analysis we evaluate relationships between Indian monsoon rainfall and five possible drivers of monsoon interannual variability. Over the time period studied (1992-2011), the El Nino-Southern Oscillation (ENSO) and the Indian Ocean dipole (IOD) are the most important of these drivers in both observations and GloSea5-GC2. Our analysis indicates that ENSO and its teleconnection with the Indian rainfall are well represented in GloSea5-GC2. However, the relationship between the IOD and Indian rainfall anomalies is too weak in GloSea5-GC2, which may be limiting the prediction skill of the local monsoon circulation and Indian rainfall. We show that this weak relationship likely results from a coupled mean state bias that limits the impact of anomalous wind forcing on SST variability, resulting in erroneous IOD SST anomalies. Known difficulties in representing convective precipitation over India may also play a role. Since Indian rainfall responds weakly to the IOD, it responds more consistently to ENSO than in observations. Our assessment identifies specific coupled biases that are likely limiting GloSea5-GC2 prediction skill, providing targets for model improvement.