Origin and impact of initialisation shocks in coupled atmosphere-ocean forecasts

Current methods for initialising coupled atmosphere-ocean forecasts often rely on the use of separate atmosphere and ocean analyses, the combination of which can leave the coupled system imbalanced at the beginning of the forecast, potentially accelerating the development of errors. Using a series of experiments with the European Centre for Medium-range Weather Forecasts coupled system, the magnitude and extent of these so-called initialisation shocks is quantified, and their impact on forecast skill measured. It is found that forecasts initialised by separate ocean and atmospheric analyses do exhibit initialisation shocks in lower atmospheric temperature, when compared to forecasts initialised using a coupled data assimilation method. These shocks result in as much as a doubling of root-mean-square error on the first day of the forecast in some regions, and in increases that are sustained for the duration of the 10-day forecasts performed here. However, the impacts of this choice of initialisation on forecast skill, assessed using independent datasets, were found to be negligible, at least over the limited period studied. Larger initialisation shocks are found to follow a change in either the atmospheric or ocean model component between the analysis and forecast phases: changes in the ocean component can lead to sea surface temperature shocks of more than 0.5K in some equatorial regions during the first day of the forecast. Implications for the development of coupled forecast systems, particularly with respect to coupled data assimilation methods, are discussed.

Tracer Leakage from Modeled Agulhas Rings

In a numerical, isopycnal, ocean model the mixing is investigated with the environment of two idealized Agulhas rings, one that splits, and one that remains coherent. The evolution of a passive tracer , initially contained within the rings, shows that tracer leakage is associated with the for mation of filaments in the early stage of ring evolution. These filaments reach down to the ther mocline. In the deepest layers leakage occurs on a larger scale. Self-advection of the rings is ver y irregular , and it is not possible to compute a Lagrangian boundar y i n order to estimate the transport of leakage from the rings. T o describe the processes that gover n tracer leakage, in a coordinate frame moving with the ring a kinematic separatrix is defined in the streamfunction field for the nondivergent flow . Initially , filaments arise because of the elongation of the ring, which is mainly gover ned by an m 5 2 instability that is collaborating with differential rotation. Because of beta, the symmetr y i s destroyed related to the separatrix associated with a stagnation point in the flow . The filament upstream of the stagnation point grows much faster and is associated with the bulk of tracer leakage. Mixing is enhanced by time dependence of the separatrix. As a result, there are no large differences between the leakage from a coherent ring, where the m 5 2 instability equilibrates, and from a splitting ring, where the m 5 2 instability keeps growing, which confir ms that the amount of leakage is mainly gover ned by the ring’ s initial defor mation combined with unsteady self-advection of the ring and not by the splitting of the ring. The decay of tracer content in the ther mocline shows that in the first months up to 40% of the ring water can be mixed with the environment. In deeper layers the decay of tracer content may reach up to 90%.

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.

Projected changes in medicanes in the HadGEM3 N512 high-resolution global climate model

Medicanes or “Mediterranean hurricanes” represent a rare and physically unique type of Mediterranean mesoscale cyclone. There are similarities with tropical cyclones with regard to their development (based on the thermodynamical disequilibrium between the warm sea and the overlying troposphere) and their kinematic and thermodynamical properties (medicanes are intense vortices with a warm core and even a cloud-free eye). Although medicanes are smaller and their wind speeds are lower than in tropical cyclones, the severity of their winds can cause substantial damage to islands and coastal areas. Concern about how human-induced climate change will affect extreme events is increasing. This includes the future impacts on medicanes due to the warming of the Mediterranean waters and the projected changes in regional atmospheric circulation. However, most global climate models do not have high enough spatial resolution to adequately represent small features such as medicanes. In this study, a cyclone tracking algorithm is applied to high resolution global climate model data with a horizontal grid resolution of approximately 25 km over the Mediterranean region. After a validation of the climatology of general Mediterranean mesoscale cyclones, changes in medicanes are determined using climate model experiments with present and future forcing. The magnitude of the changes in the winds, frequency and location of medicanes is assessed. While no significant changes in the total number of Mediterranean mesoscale cyclones are found, medicanes tend to decrease in number but increase in intensity. The model simulation suggests that medicanes tend to form more frequently in the Gulf of Lion–Genoa and South of Sicily.

Evaluation of a polynya flux model by means of thermal infrared satellite estimates

We test the ability of a two-dimensional flux model to simulate polynya events with narrow open-water zones by comparing model results to ice-thickness and ice-production estimates derived from thermal infrared Moderate Resolution Imaging Spectroradiometer (MODIS) observations in conjunction with an atmospheric dataset. Given a polynya boundary and an atmospheric dataset, the model correctly reproduces the shape of an 11 day long event, using only a few simple conservation laws. Ice production is slightly overestimated by the model, owing to an underestimated ice thickness. We achieved best model results with the consolidation thickness parameterization developed by Biggs and others (2000). Observed regional discrepancies between model and satellite estimates might be a consequence of the missing representation of the dynamic of the thin-ice thickening (e.g. rafting). We conclude that this simplified polynya model is a valuable tool for studying polynya dynamics and estimating associated fluxes of single polynya events.

Re-examining the U.K.'s greatest tornado outbreak: forecasting the limited extent of tornadoes along a cold front

On 23 November 1981, a strong cold front swept across the U.K., producing tornadoes from the west to the east coasts. An extensive campaign to collect tornado reports by the Tornado and Storm Research Organisation (TORRO) resulted in 104 reports, the largest U.K. outbreak. The front was simulated with a convection-permitting numerical model down to 200-m horizontal grid spacing to better understand its evolution and meteorological environment. The event was typical of tornadoes in the U.K., with convective available potential energy (CAPE) less than 150 J kg-1, 0-1-km wind shear of 10-20 m s-1, and a narrow cold-frontal rainband forming precipitation cores and gaps. A line of cyclonic absolute vorticity existed along the front, with maxima as large as 0.04 s-1. Some hook-shaped misovortices bore kinematic similarity to supercells. The narrow swath along which the line was tornadic was bounded on the equatorward side by weak vorticity along the line and on the poleward side by zero CAPE, enclosing a region where the environment was otherwise favorable for tornadogenesis. To determine if the 104 tornado reports were plausible, first possible duplicate reports were eliminated, resulting in as few as 58 tornadoes to as many as 90. Second, the number of possible parent misovortices that may have spawned tornadoes is estimated from model output. The number of plausible tornado reports in the 200-m grid-spacing domain was 22 and as many as 44, whereas the model simulation was used to estimate 30 possible parent misovortices within this domain. These results suggest that 90 reports was plausible.

The teleconnection of the tropical Atlantic to Indo-Pacific sea surface temperatures on inter-annual to centennial time scales: a review of recent findings

In this paper, the teleconnections from the tropical Atlantic to the Indo-Pacific region from inter-annual to centennial time scales will be reviewed. Identified teleconnections and hypotheses on mechanisms at work are reviewed and further explored in a century-long pacemaker coupled ocean-atmosphere simulation ensemble. There is a substantial impact of the tropical Atlantic on the Pacific region at inter-annual time scales. An Atlantic Niño (Niña) event leads to rising (sinking) motion in the Atlantic region, which is compensated by sinking (rising) motion in the central-western Pacific. The sinking (rising) motion in the central-western Pacific induces easterly (westerly) surface wind anomalies just to the west, which alter the thermocline. These perturbations propagate eastward as upwelling (downwelling) Kelvin-waves, where they increase the probability for a La Niña (El Niño) event. Moreover, tropical North Atlantic sea surface temperature anomalies are also able to lead La Niña/El Niño development. At multidecadal time scales, a positive (negative) Atlantic Multidecadal Oscillation leads to a cooling (warming) of the eastern Pacific and a warming (cooling) of the western Pacific and Indian Ocean regions. The physical mechanism for this impact is similar to that at inter-annual time scales. At centennial time scales, the Atlantic warming induces a substantial reduction of the eastern Pacific warming even under CO2 increase and to a strong subsurface cooling.

Annular modes and apparent eddy feedbacks in the Southern Hemisphere

Lagged correlation analysis is often used to infer intraseasonal dynamical effects but is known to be affected by non-stationarity. We highlight a pronounced quasi-two-year peak in the anomalous zonal wind and eddy momentum flux convergence power spectra in the Southern Hemisphere, which is prima facie evidence for non-stationarity. We then investigate the consequences of this non-stationarity for the Southern Annular Mode and for eddy momentum flux convergence. We argue that positive lagged correlations previously attributed to the existence of an eddy feedback are more plausibly attributed to non-stationary interannual variability external to any potential feedback process in the mid-latitude troposphere. The findings have implications for the diagnosis of feedbacks in both models and re-analysis data as well as for understanding the mechanisms underlying variations in the zonal wind.

Cloud streets and land-water interactions in the Amazon

Cloud streets are common feature in the Amazon Basin. They form from the combination of the vertical trade wind stress and moist convection. Here, satellite imagery, data collected during the COBRA-PARA (Caxiuan Observations in the Biosphere, River and Atmosphere of Para) field campaign, and high resolution modeling are used to understand the streets` formation and behavior. The observations show that the streets have an aspect ratio of about 3.5 and they reach their maximum activity around 15:00 UTC when the wind shear is weaker, and the convective boundary layer reaches its maximum height. The simulations reveal that the cloud streets onset is caused by the local circulations and convection produced at the interfaces between forest and rivers of the Amazon. The satellite data and modeling show that the large rivers anchor the cloud streets producing a quasi-stationary horizontal pattern. The streets are associated with horizontal roll vortices parallel to the mean flow that organizes the turbulence causing advection of latent heat flux towards the upward branches. The streets have multiple warm plumes that promote a connection between the rolls. These spatial patterns allow fundamental insights on the interpretation of the Amazon exchanges between surface and atmosphere with important consequences for the climate change understanding.

Density profiles in the raise and peel model with and without a wall; physics and combinatorics

We consider the raise and peel model of a one-dimensional fluctuating interface in the presence of an attractive wall. The model can also describe a pair annihilation process in disordered unquenched media with a source at one end of the system. For the stationary states, several density profiles are studied using Monte Carlo simulations. We point out a deep connection between some profiles seen in the presence of the wall and in its absence. Our results are discussed in the context of conformal invariance ( c = 0 theory). We discover some unexpected values for the critical exponents, which are obtained using combinatorial methods. We have solved known ( Pascal`s hexagon) and new (split-hexagon) bilinear recurrence relations. The solutions of these equations are interesting in their own right since they give information on certain classes of alternating sign matrices.

Pedal : identification of models for duodenal administration of levodopa

Backgound and aims: The main purpose of the PEDAL study is to identify and estimate sample individual pharmacokinetic- pharmacodynamic (PK/PD) models for duodenal infusion of levodopa/carbidopa (Duodopa®) that can be used for in numero simulation of treatment strategies. Other objectives are to study the absorption of Duodopa® and to form a basis for power calculation for a future larger study. PK/PD based on oral levodopa is problematic because of irregular gastric emptying. Preliminary work with data from [Gundert-Remy U et al. Eur J Clin Pharmacol 1983;25:69-72] suggested that levodopa infusion pharmacokinetics can be described by a two-compartment model. Background research led to a hypothesis for an effect model incorporating concentration-unrelated fluctuations, more complex than standard E-max models. Methods: PEDAL involved a few patients already on Duodopa®. A bolus dose (normal morning dose plus 50%) was given after a washout during night. Data collection continued until the clinical effect was back at baseline. The procedure was repeated on two non-consecutive days per patient. The following data were collected in 5 to 15 minutes intervals: i) Accelerometer data. ii) Three e-diary questions about ability to walk, feelings of “off” and “dyskinesia”. iii) Clinical assessment of motor function by a physician. iv) Plasma concentrations of levodopa, carbidopa and the metabolite 3-O-methyldopa. The main effect variable will be the clinical assessment. Results: At date of abstract submission, lab analyses were currently being performed. Modelling results, simulation experiments and conclusions will be presented in our poster.

QoS-aware model driven architecture through the UML and CIM

The specification of Quality of Service (QoS) constraints over software design requires measures that ensure such requirements are met by the delivered product. Achieving this goal is non-trivial, as it involves, at least, identifying how QoS constraint specifications should be checked at the runtime. In this paper we present an implementation of a Model Driven Architecture (MDA) based framework for the runtime monitoring of QoS properties. We incorporate the UML2 superstructure and the UML profile for Quality of Service to provide abstract descriptions of component-and-connector systems. We then define transformations that refine the UML2 models to conform with the Distributed Management Taskforce (DMTF) Common Information Model (CIM) (Distributed Management Task Force Inc. 2006), a schema standard for management and instrumentation of hardware and software. Finally, we provide a mapping the CIM metamodel to a .NET-based metamodel for implementation of the monitoring infrastructure utilising various .NET features including the Windows Management Instrumentation (WMI) interface.

Income inequality in a job-search model with heterogeneous time preferences

This paper explores the use of an intertemporal job-search model in the investigation of within-cohort and between-cohort income inequality, the latter being generated by the heterogeneity of time preferences among cohorts of homogenous workers and the former by the cross-sectional turnover in the job market. It also offers an alternative explanation for the empirically-documented negative correlation between time preference and labor income. Under some speciÖc distributions regarding wage offers and time preferences, we show how the within-cohort and between-cohort Gini coe¢ cients of income distribution can be calculated, and how they vary as a function of the parameters of the model.

Foreign funding to an emerging market: the Monetary Premium Theory and the Brazilian Case, 1991 - 1998

We develop a framework to explain the private capital flows between the rest of the world and an emerging economy. The model, based on the monetary premium theory, relates an endogenous supply of foreign capitals to an endogenous differential of interest rates; its estimation uses the econometric techniques initiated by Heckman. Four questions regarding the capital flows phenomenon are explored, including the statistical process that governs the events of default and the impact of the probability of default on the interest rate differential. Using the methodology, we analyse the dynamics of foreign capital movements in Brazil during the 1991- 1998 period.

On the growth effects of barriers to trade

We study the macroeconomic effects of international trade policy by integrating a Hecksher-Ohlin trade model into an optimal-growth framework. The model predicts that an open economy will have higher factor productivity and faster growth. Also, under protectionist policies there may be “development traps,” or additional steady states with low income. In the last case, higher tariffs imply lower incomes, so that the large cross-country differences in barriers to trade may explain part of the huge dispersion of per capita income observed across countries. The model simulation shows that the link between trade and macroeconomic performance may be quantitatively important.