Now showing 1 - 10 of 14
  • Publication
    From Richardson to early numerical weather prediction
    (Cambridge University Press, 2010-12)
    The development of computer models for numerical simulation of the atmosphere and oceans is one of the great scientific triumphs of the past fifty years. These models have added enormously to our understanding of the complex processes in the atmosphere and oceans. The consequences for humankind of ongoing climate change will be far-reaching. Earth system models are the best means we have of predicting the future of our climate. The basic ideas of numerical forecasting and climate modeling were developed about a century ago, long before the first electronic computer was constructed. However, advances on several fronts were necessary before numerical prediction could be put into practice. A fuller understanding of atmospheric dynamics allowed the development of simplified systems of equations; regular observations of the free atmosphere provided the initial conditions; stable finite difference schemes were developed; and powerful electronic computers provided a practical means of carrying out the calculations required to predict the changes in the weather. In this chapter, we trace the history of computer forecasting from Richardson’s prodigious manual computation, through the ENIAC (Electronic Numerical Integrator and Computer) integrations to the early days of operational numerical weather prediction and climate modeling. The useful range of deterministic prediction is increasing by about one day each decade. We set the scene for the story of the remarkable progress in weather forecasting and in climate modeling over the past fifty years, which will be treated in subsequent chapters.
      371
  • Publication
    The two-child paradox : dichotomy and ambiguity
    (Irish Mathematical Society, 2011-07)
    Given that one of the children in a two-child family is a boy, what are the chances that the other is also a boy. The intuitive answer is 50 : 50. More careful investigation leads us to a 1-in-3 chance. We investigate circumstances under which these answers are correct. The imposition of further conditions yields some very surprising results.
      404
  • Publication
    The Fractal Boundary of the Power Tower Function
    (Associacao Ludus, 2017-08-23)
    We consider the function called the power tower function, defined by iterated exponentiation (or tetration) of the complex variable z. For real values x, it converges on the interval exp(−e)
      333
  • Publication
    Stokes's Fundamental Contributions to Fluid Dynamics
    (Oxford University Press, 2019-06-27)
    George Gabriel Stokes was one of the giants of hydrodynamics in the nineteenth century. He made fundamental mathematical contributions to fluid dynamics that had profound practical consequences. The basic equations formulated by him, the Navier-Stokes equations, are capable of describing fluid flows over a vast range of magnitudes. They play a central role in numerical weather prediction, in the simulation of blood flow in the body and in countless other important applications. In this chapter we put the primary focus on the two most important areas of Stokes’s work on fluid dynamics, the derivation of the Navier-Stokes equations and the theory of finite amplitude oscillatory water waves. Stokes became an undergraduate at Cambridge in 1837. He was coached by the ‘Senior Wrangler-maker’, William Hopkins and, in 1841, Stokes was Senior Wrangler and first Smith’s Prizeman. It was following a suggestion of Hopkins that Stokes took up the study of hydrodynamics, which was at that time a neglected area of study in Cambridge. Stokes was to make profound contributions to hydrodynamics, his most important being the rigorous establishment of the mathematical equations for fluid motions, and the theoretical explanation of a wide range of phenomena relating to wave motions in water.
      344
  • Publication
    Laplace transform integration of the shallow-water equations. Part 1: Eulerian formulation and Kelvin waves
    (Wiley, 2011-04) ;
    A filtering integration scheme is developed, using a modification of the contour used to invert the Laplace transform (LT). It is shown to eliminate components with frequencies higher than a specified cut-off value. Thus it is valuable for integrations of the equations governing atmospheric flow. The scheme is implemented in a shallow water model with an Eulerian treatment of advection. It is compared to a reference model using the semi-implicit (SI) scheme. The LT scheme is shown to treat dynamically important Kelvin waves more accurately than the SI scheme.
      323Scopus© Citations 8
  • Publication
    Adaptive post-processing of short-term wind forecasts for energy applications
    (Wiley, 2011-04) ;
    We present a new method of reducing the error in predicted wind speed, thus enabling better management of wind energy facilities. A numerical weather prediction model, COSMO, was used to produce 48 h forecast data every day in 2008 at horizontal resolutions of 10 and 3 km. A new adaptive statistical method was applied to the model output to improve the forecast skill. The method applied corrective weights to a set of forecasts generated using several post-processing methods. The weights were calculated based on the recent skill of the different forecasts. The resulting forecast data were compared with observed data, and skill scores were calculated to allow comparison between different post-processing methods. The total root mean square error performance of the composite forecast is superior to that of any of the individual methods.
      777Scopus© Citations 13
  • Publication
    Initialization
    (Springer, 2010-08) ;
    The spectrum of atmospheric motions is vast, encompassing phenomena having periods ranging from seconds to millennia. The motions of interest to the forecaster typically have time-scales of a day or longer, but the mathematical models used for numerical prediction describe a broader span of dynamical features than those of direct concern. For many purposes these higher frequency components can be regarded as noise contaminating the motions of meteorological interest. The elimination of this noise is achieved by adjustment of the initial fields, a process called initialization.
      177Scopus© Citations 16
  • Publication
    Precession and recession of the rock'n'roller
    (IOP Publishing, 2009-09-30) ;
    We study the dynamics of a spherical rigid body that rocks and rolls on a plane under the effect of gravity. The distribution of mass is non-uniform and the centre of mass does not coincide with the geometric centre. The symmetric case, with moments of inertia I1 = I2 < I3, is integrable and themotion is completely regular. Three known conservation laws are the total energy E, Jellett’s quantity QJ and Routh’s quantity QR. When the inertial symmetry I1 = I2 is broken, even slightly, the character of the solutions is profoundly changed and new types of motion become possible. We derive the equations governing the general motion and present analytical and numerical evidence of the recession, or reversal of precession, that has been observed in physical experiments. We present an analysis of recession in terms of critical lines dividing the (QR,QJ ) plane into four dynamically disjoint zones. We prove that recession implies the lack of conservation of Jellett’s and Routh’s quantities, by identifying individual reversals as crossings of the orbit (QR(t ),QJ (t)) through the critical lines. Consequently, a method is found to produce a large number of initial conditions so that the system will exhibit recession.
      433Scopus© Citations 11
  • Publication
    On resonant Rossby-Haurwitz triads
    (Wiley-Blackwell, 2009-05)
    The dynamics of non-divergent flow on a rotating sphere are described by the conservation of absolute vorticity. The analytical study of the non-linear barotropic vorticity equation is greatly facilitated by the expansion of the solution in spherical harmonics and truncation at low order. The normal modes are the well-known Rossby–Haurwitz (RH) waves, which represent the natural oscillations of the system. Triads of RH waves, which satisfy conditions for resonance, are of critical importance for the distribution of energy in the atmosphere. We show how non-linear interactions of resonant RH triads may result in dynamic instability of large-scale components. We also demonstrate a mathematical equivalence between the equations for an orographically forced triad and a simple mechanical system, the forced-damped swinging spring. This equivalence yields insight concerning the bounded response to a constant forcing in the absence of damping. An examination of triad interactions in atmospheric reanalysis data would be of great interest.
      404Scopus© Citations 22
  • Publication
    Reducing errors of wind speed forecasts by an optimal combination of post-processing methods
    (Wiley-Blackwell, 2011-09-13) ; ;
    Seven adaptive approaches to post-processing wind speed forecasts are discussed and compared. 48-hour forecasts are run at horizontal resolutions of 7 km and 3 km for a domain centred over Ireland. Forecast wind speeds over a two year period are compared to observed wind speeds at seven synoptic stations around Ireland and skill scores calculated. Two automatic methods for combining forecast streams are applied. The forecasts produced by the combined methods give bias and root mean squared errors that are better than the numerical weather prediction forecasts at all station locations. One of the combined forecast methods results in skill scores that are equal to or better than all of its component forecast streams. This method is straightforward to apply and should prove beneficial in operational wind forecasting.
      981Scopus© Citations 58