Assimilation of Geosat Altimeter Data for the Agulhas Current Using the Ensemble Kalman Filter with a Quasigeostrophic Model

The ring-shedding process in the Agulhas Current is studied using the ensemble Kalman filter to assimilate geosat altimeter data into a two-layer quasigeostrophic ocean model. The properties of the ensemble Kalman filter are further explored with focus on the analysis scheme and the use of gridded data. The Geosat data consist of 10 fields of gridded sea-surface height anomalies separated 10 days apart that are added to a climatic mean field. This corresponds to a huge number of data values, and a data reduction scheme must be applied to increase the efficiency of the analysis procedure. Further, it is illustrated how one can resolve the rank problem occurring when a too large dataset or a small ensemble is used.

Comparison of terrain following and cut cell grids using a non-hydrostatic model

Terrain following coordinates are widely used in operational models but the cut cell method has been proposed as an alternative that can more accurately represent atmospheric dynamics over steep orography. Because the type of grid is usually chosen during model implementation, it becomes necessary to use different models to compare the accuracy of different grids. In contrast, here a C-grid finite volume model enables a like-for-like comparison of terrain following and cut cell grids. A series of standard two-dimensional tests using idealised terrain are performed: tracer advection in a prescribed horizontal velocity field, a test starting from resting initial conditions, and orographically induced gravity waves described by nonhydrostatic dynamics. In addition, three new tests are formulated: a more challenging resting atmosphere case, and two new advection tests having a velocity field that is everywhere tangential to the terrain following coordinate surfaces. These new tests present a challenge on cut cell grids. The results of the advection tests demonstrate that accuracy depends primarily upon alignment of the flow with the grid rather than grid orthogonality. A resting atmosphere is well-maintained on all grids. In the gravity waves test, results on all grids are in good agreement with existing results from the literature, although terrain following velocity fields lead to errors on cut cell grids. Due to semi-implicit timestepping and an upwind-biased, explicit advection scheme, there are no timestep restrictions associated with small cut cells. We do not find the significant advantages of cut cells or smoothed coordinates that other authors find.

The medieval climate anomaly and Byzantium: a review of the evidence on climatic fluctuations, economic performance and societal change

At the beginning of the Medieval Climate Anomaly, in the ninth and tenth century, the medieval eastern Roman empire, more usually known as Byzantium, was recovering from its early medieval crisis and experiencing favourable climatic conditions for the agricultural and demographic growth. Although in the Balkans and Anatolia such favourable climate conditions were prevalent during the eleventh century, parts of the imperial territories were facing significant challenges as a result of external political/military pressure. The apogee of medieval Byzantine socio-economic development, around AD 1150, coincides with a period of adverse climatic conditions for its economy, so it becomes obvious that the winter dryness and high climate variability at this time did not hinder Byzantine society and economy from achieving that level of expansion. Soon after this peak, towards the end of the twelfth century, the populations of the Byzantine world were experiencing unusual climatic conditions with marked dryness and cooler phases. The weakened Byzantine socio-political system must have contributed to the events leading to the fall of Constantinople in AD 1204 and the sack of the city. The final collapse of the Byzantine political control over western Anatolia took place half century later, thus contemporaneous with the strong cooling effect after a tropical volcanic eruption in AD 1257. We suggest that, regardless of a range of other influential factors, climate change was also an important contributing factor to the socio-economic changes that took place in Byzantium during the Medieval Climate Anomaly. Crucially, therefore, while the relatively sophisticated and complex Byzantine society was certainly influenced by climatic conditions, and while it nevertheless displayed a significant degree of resilience, external pressures as well as tensions within the Byzantine society more broadly contributed to an increasing vulnerability in respect of climate impacts. Our interdisciplinary analysis is based on all available sources of information on the climate and society of Byzantium, that is textual (documentary), archaeological, environmental, climate and climate model-based evidence about the nature and extent of climate variability in the eastern Mediterranean. The key challenge was, therefore, to assess the relative influence to be ascribed to climate variability and change on the one hand, and on the other to the anthropogenic factors in the evolution of Byzantine state and society (such as invasions, changes in international or regional market demand and patterns of production and consumption, etc.). The focus of this interdisciplinary

An idealized study of coupled atmosphere-ocean 4D-Var in the presence of model error

Atmosphere only and ocean only variational data assimilation (DA) schemes are able to use window lengths that are optimal for the error growth rate, non-linearity and observation density of the respective systems. Typical window lengths are 6-12 hours for the atmosphere and 2-10 days for the ocean. However, in the implementation of coupled DA schemes it has been necessary to match the window length of the ocean to that of the atmosphere, which may potentially sacrifice the accuracy of the ocean analysis in order to provide a more balanced coupled state. This paper investigates how extending the window length in the presence of model error affects both the analysis of the coupled state and the initialized forecast when using coupled DA with differing degrees of coupling. Results are illustrated using an idealized single column model of the coupled atmosphere-ocean system. It is found that the analysis error from an uncoupled DA scheme can be smaller than that from a coupled analysis at the initial time, due to faster error growth in the coupled system. However, this does not necessarily lead to a more accurate forecast due to imbalances in the coupled state. Instead coupled DA is more able to update the initial state to reduce the impact of the model error on the accuracy of the forecast. The effect of model error is potentially most detrimental in the weakly coupled formulation due to the inconsistency between the coupled model used in the outer loop and uncoupled models used in the inner loop.