On the recursive sequence Xn = axn-1+bxn-2/c+dxn-1xn-2

In this Paper, we study the invariant intervals, the global attractivity of the equilibrium points, and the asymptotic behavior of the solutions of the difference equation x(n) = ax(n-1) + bx(n-2) / c + dx(n-1)x(n-2), n =1, 2, ..., where a greater than or equal to 0, b, c, d > 0. (C) 2004 Elsevier Inc. All rights reserved.

A weak formulation of Roe's approximate Riemann Solver applied to 'Barotropic' flows

A weak formulation of Roe's approximate Riemann solver is applied to the equations of ‘barotropic’ flow, including the shallow water equations, and it is shown that this leads to an approximate Riemann solver recently presented for such flows.

Simplified recursive identifier for ARMA processes

Symmetrical behaviour of the covariance matrix and the positive-definite criterion are used to simplify identification of single-input/single-output systems using recursive least squares. Simulation results are obtained and these are compared with ordinary recursive least squares. The adaptive nature of the identifier is verified by varying the system parameters on convergence.

Evaluating a new LISFLOOD-FP formulation with data from the summer 2007 floods in Tewkesbury, UK

Two-dimensional flood inundation modelling is a widely used tool to aid flood risk management. In urban areas, the model spatial resolution required to represent flows through a typical street network often results in an impractical computational cost at the city scale. This paper presents the calibration and evaluation of a recently developed formulation of the LISFLOOD-FP model, which is more computationally efficient at these resolutions. Aerial photography was available for model evaluation on 3 days from the 24 to the 31 of July. The new formulation was benchmarked against the original version of the model at 20 and 40 m resolutions, demonstrating equally accurate simulation, given the evaluation data but at a 67 times faster computation time. The July event was then simulated at the 2 m resolution of the available airborne LiDAR DEM. This resulted in more accurate simulation of the floodplain drying dynamics compared with the coarse resolution models, although maximum inundation levels were simulated equally well at all resolutions tested.

Adaptive recursive control of a soil testing load frame

Recursive Learning Control (RLC) has the potential to significantly reduce the tracking error in many repetitive trajectory applications. This paper presents an application of RLC to a soil testing load frame where non-adaptive techniques struggle with the highly nonlinear nature of soil. The main purpose of the controller is to apply a sinusoidal force reference trajectory on a soil sample with a high degree of accuracy and repeatability. The controller uses a feedforward control structure, recursive least squares adaptation algorithm and RLC to compensate for periodic errors. Tracking error is reduced and stability is maintained across various soil sample responses.

Prediction of visibility and aerosol within the operational Met Office unified model. I: Model formulation and variational assimilation

The formulation and performance of the Met Office visibility analysis and prediction system are described. The visibility diagnostic within the limited-area Unified Model is a function of humidity and a prognostic aerosol content. The aerosol model includes advection, industrial and general urban sources, plus boundary-layer mixing and removal by rain. The assimilation is a 3-dimensional variational scheme in which the visibility observation operator is a very nonlinear function of humidity, aerosol and temperature. A quality control scheme for visibility data is included. Visibility observations can give rise to humidity increments of significant magnitude compared with the direct impact of humidity observations. We present the results of sensitivity studies which show the contribution of different components of the system to improved skill in visibility forecasts. Visibility assimilation is most important within the first 6-12 hours of the forecast and for visibilities below 1 km, while modelling of aerosol sources and advection is important for slightly higher visibilities (1-5 km) and is still significant at longer forecast times

Review of the formulation of present‐generation stratospheric chemistry‐climate models and associated external forcings

The goal of the Chemistry‐Climate Model Validation (CCMVal) activity is to improve understanding of chemistry‐climate models (CCMs) through process‐oriented evaluation and to provide reliable projections of stratospheric ozone and its impact on climate. An appreciation of the details of model formulations is essential for understanding how models respond to the changing external forcings of greenhouse gases and ozonedepleting substances, and hence for understanding the ozone and climate forecasts produced by the models participating in this activity. Here we introduce and review the models used for the second round (CCMVal‐2) of this intercomparison, regarding the implementation of chemical, transport, radiative, and dynamical processes in these models. In particular, we review the advantages and problems associated with approaches used to model processes of relevance to stratospheric dynamics and chemistry. Furthermore, we state the definitions of the reference simulations performed, and describe the forcing data used in these simulations. We identify some developments in chemistry‐climate modeling that make models more physically based or more comprehensive, including the introduction of an interactive ocean, online photolysis, troposphere‐stratosphere chemistry, and non‐orographic gravity‐wave deposition as linked to tropospheric convection. The relatively new developments indicate that stratospheric CCM modeling is becoming more consistent with our physically based understanding of the atmosphere.

Origin of regional climate differences: role of boundary conditions and model formulation in two GCMs

Model differences in projections of extratropical regional climate change due to increasing greenhouse gases are investigated using two atmospheric general circulation models (AGCMs): ECHAM4 (Max Planck Institute, version 4) and CCM3 (National Center for Atmospheric Research Community Climate Model version 3). Sea-surface temperature (SST) fields calculated from observations and coupled versions of the two models are used to force each AGCM in experiments based on time-slice methodology. Results from the forced AGCMs are then compared to coupled model results from the Coupled Model Intercomparison Project 2 (CMIP2) database. The time-slice methodology is verified by showing that the response of each model to doubled CO2 and SST forcing from the CMIP2 experiments is consistent with the results of the coupled GCMs. The differences in the responses of the models are attributed to (1) the different tropical SST warmings in the coupled simulations and (2) the different atmospheric model responses to the same tropical SST warmings. Both are found to have important contributions to differences in implied Northern Hemisphere (NH) winter extratropical regional 500 mb height and tropical precipitation climate changes. Forced teleconnection patterns from tropical SST differences are primarily responsible for sensitivity differences in the extratropical North Pacific, but have relatively little impact on the North Atlantic. There are also significant differences in the extratropical response of the models to the same tropical SST anomalies due to differences in numerical and physical parameterizations. Differences due to parameterizations dominate in the North Atlantic. Differences in the control climates of the two coupled models from the current climate, in particular for the coupled model containing CCM3, are also demonstrated to be important in leading to differences in extratropical regional sensitivity.