On increasing global temperatures: 75 years after Callendar

In 1938, Guy Stewart Callendar was the first to demonstrate that the Earth’s land surface was warming. Callendar also suggested that the production of carbon dioxide by the combustion of fossil fuels was responsible for much of this modern change in climate. This short note marks the 75th anniversary of Callendar’s landmark study and demonstrates that his global land temperature estimates agree remarkably well with more recent analyses.

An exact local conservation theorem for finite-amplitude disturbances to non-parallel shear flows, with remarks on Hamiltonian structure and on Arnol'd's stability theorems

Disturbances of arbitrary amplitude are superposed on a basic flow which is assumed to be steady and either (a) two-dimensional, homogeneous, and incompressible (rotating or non-rotating) or (b) stably stratified and quasi-geostrophic. Flow over shallow topography is allowed in either case. The basic flow, as well as the disturbance, is assumed to be subject neither to external forcing nor to dissipative processes like viscosity. An exact, local ‘wave-activity conservation theorem’ is derived in which the density A and flux F are second-order ‘wave properties’ or ‘disturbance properties’, meaning that they are O(a2) in magnitude as disturbance amplitude a [rightward arrow] 0, and that they are evaluable correct to O(a2) from linear theory, to O(a3) from second-order theory, and so on to higher orders in a. For a disturbance in the form of a single, slowly varying, non-stationary Rossby wavetrain, $\overline{F}/\overline{A}$ reduces approximately to the Rossby-wave group velocity, where (${}^{-}$) is an appropriate averaging operator. F and A have the formal appearance of Eulerian quantities, but generally involve a multivalued function the correct branch of which requires a certain amount of Lagrangian information for its determination. It is shown that, in a certain sense, the construction of conservable, quasi-Eulerian wave properties like A is unique and that the multivaluedness is inescapable in general. The connection with the concepts of pseudoenergy (quasi-energy), pseudomomentum (quasi-momentum), and ‘Eliassen-Palm wave activity’ is noted. The relationship of this and similar conservation theorems to dynamical fundamentals and to Arnol'd's nonlinear stability theorems is discussed in the light of recent advances in Hamiltonian dynamics. These show where such conservation theorems come from and how to construct them in other cases. An elementary proof of the Hamiltonian structure of two-dimensional Eulerian vortex dynamics is put on record, with explicit attention to the boundary conditions. The connection between Arnol'd's second stability theorem and the suppression of shear and self-tuning resonant instabilities by boundary constraints is discussed, and a finite-amplitude counterpart to Rayleigh's inflection-point theorem noted

First order k-th moment finite element analysis of nonlinear operator equations with stochastic data

We develop and analyze a class of efficient Galerkin approximation methods for uncertainty quantification of nonlinear operator equations. The algorithms are based on sparse Galerkin discretizations of tensorized linearizations at nominal parameters. Specifically, we consider abstract, nonlinear, parametric operator equations J(\alpha ,u)=0 for random input \alpha (\omega ) with almost sure realizations in a neighborhood of a nominal input parameter \alpha _0. Under some structural assumptions on the parameter dependence, we prove existence and uniqueness of a random solution, u(\omega ) = S(\alpha (\omega )). We derive a multilinear, tensorized operator equation for the deterministic computation of k-th order statistical moments of the random solution's fluctuations u(\omega ) - S(\alpha _0). We introduce and analyse sparse tensor Galerkin discretization schemes for the efficient, deterministic computation of the k-th statistical moment equation. We prove a shift theorem for the k-point correlation equation in anisotropic smoothness scales and deduce that sparse tensor Galerkin discretizations of this equation converge in accuracy vs. complexity which equals, up to logarithmic terms, that of the Galerkin discretization of a single instance of the mean field problem. We illustrate the abstract theory for nonstationary diffusion problems in random domains.

An improved sum-product estimate for general finite fields

An improved sum-product estimate for subsets of a finite field whose order is not prime is provided. It is shown, under certain conditions, that max{∣∣∣A+A∣∣∣,∣∣∣A⋅A∣∣∣}≫∣∣A∣∣12/11(log2∣∣A∣∣)5/11. This new estimate matches, up to a logarithmic factor, the current best known bound obtained over prime fields by Rudnev

Proteolysis of milk heated at high temperatures by native enzymes analysed by trinitrobenzene sulphonic acid (TNBS) method

Native enzymes play a significant role in proteolysis of milk during storage. This is significant for heat resistant native enzymes. Plasmin is one of the most heat resistant enzymes found in milk. It has been reported to survive several heat treatments, causing spoilage during storage. The aim of this study was to assess susceptibility of high temperature heated milk to proteolysis by native enzymes. The trinitrobenzene sulphonic acid (TNBS) method was used for this purpose. Raw milk was heated at 110, 120, 130,142°C for 2 s and 85°C for 15 s and milk processed at low temperature (85°C /15s) was selected to mimic pasteurisation. TNBS method confirmed that raw milk and milk processed at 85°C /15s were the most proteolysed, whereas treatment of milk at high temperatures (110, 120, 130 and 142°C for 2 s) inactivated the native enzymes. It may thus be concluded that high temperature processing positively affects proteolysis by lowering its susceptibility to spoilage during storage.

Gradient recovery in adaptive finite-element methods for parabolic problems

We derive energy-norm a posteriori error bounds, using gradient recovery (ZZ) estimators to control the spatial error, for fully discrete schemes for the linear heat equation. This appears to be the �rst completely rigorous derivation of ZZ estimators for fully discrete schemes for evolution problems, without any restrictive assumption on the timestep size. An essential tool for the analysis is the elliptic reconstruction technique.Our theoretical results are backed with extensive numerical experimentation aimed at (a) testing the practical sharpness and asymptotic behaviour of the error estimator against the error, and (b) deriving an adaptive method based on our estimators. An extra novelty provided is an implementation of a coarsening error "preindicator", with a complete implementation guide in ALBERTA in the appendix.

A finite element method for nonlinear elliptic problems

We present a Galerkin method with piecewise polynomial continuous elements for fully nonlinear elliptic equations. A key tool is the discretization proposed in Lakkis and Pryer, 2011, allowing us to work directly on the strong form of a linear PDE. An added benefit to making use of this discretization method is that a recovered (finite element) Hessian is a byproduct of the solution process. We build on the linear method and ultimately construct two different methodologies for the solution of second order fully nonlinear PDEs. Benchmark numerical results illustrate the convergence properties of the scheme for some test problems as well as the Monge–Amp`ere equation and the Pucci equation.

A finite element method for second order nonvariational elliptic problems

We propose a numerical method to approximate the solution of second order elliptic problems in nonvariational form. The method is of Galerkin type using conforming finite elements and applied directly to the nonvariational (nondivergence) form of a second order linear elliptic problem. The key tools are an appropriate concept of “finite element Hessian” and a Schur complement approach to solving the resulting linear algebra problem. The method is illustrated with computational experiments on three linear and one quasi-linear PDE, all in nonvariational form.

A finite element-spherical harmonics model for radiative transfer in inhomogeneous clouds. Part II. some applications

EVENT has been used to examine the effects of 3D cloud structure, distribution, and inhomogeneity on the scattering of visible solar radiation and the resulting 3D radiation field. Large eddy simulation and aircraft measurements are used to create realistic cloud fields which are continuous or broken with smooth or uneven tops. The values, patterns and variance in the resulting downwelling and upwelling radiation from incident visible solar radiation at different angles are then examined and compared to measurements. The results from EVENT confirm that 3D cloud structure is important in determining the visible radiation field, and that these results are strongly influenced by the solar zenith angle. The results match those from other models using visible solar radiation, and are supported by aircraft measurements of visible radiation, providing confidence in the new model.

Nature of vegetation and building morphology characteristics across a city: influence on shadow patterns and mean radiant temperatures in London

Vegetation and building morphology characteristics are investigated at 19 sites on a north-south LiDAR transect across the megacity of London. Local maxima of mean building height and building plan area density at the city centre are evident. Surprisingly, the mean vegetation height (zv3) is also found to be highest in the city centre. From the LiDAR data various morphological parameters are derived as well as shadow patterns. Continuous images of the effects of buildings and of buildings plus vegetationon sky view factor (Ψ) are derived. A general reduction of Ψ is found, indicating the importance of including vegetation when deriving Ψ in urban areas. The contribution of vegetation to the shadowing at ground level is higher during summer than in autumn. Using these 3D data the influence on urban climate and mean radiant temperature (T mrt ) is calculated with SOLWEIG. The results from these simulations highlight that vegetation can be most effective at reducing heat stress within dense urban environments in summer. The daytime average T mrt is found to be lowest in the densest urban environments due to shadowing; foremost from buildings but also from trees. It is clearly shown that this method could be used to quantify the influence of vegetation on T mrt within the urban environment. The results presented in this paper highlight a number of possible climate sensitive planning practices for urban areas at the local scale (i.e. 102- 5 × 103 m).

The influence of vegetation and building morphology on shadow patterns and mean radiant temperatures in urban areas: model development and evaluation

The solar and longwave environmental irradiance geometry (SOLWEIG) model simulates spatial variations of 3-D radiation fluxes and mean radiant temperature (T mrt) as well as shadow patterns in complex urban settings. In this paper, a new vegetation scheme is included in SOLWEIG and evaluated. The new shadow casting algorithm for complex vegetation structures makes it possible to obtain continuous images of shadow patterns and sky view factors taking both buildings and vegetation into account. For the calculation of 3-D radiation fluxes and T mrt, SOLWEIG only requires a limited number of inputs, such as global shortwave radiation, air temperature, relative humidity, geographical information (latitude, longitude and elevation) and urban geometry represented by high-resolution ground and building digital elevation models (DEM). Trees and bushes are represented by separate DEMs. The model is evaluated using 5 days of integral radiation measurements at two sites within a square surrounded by low-rise buildings and vegetation in Göteborg, Sweden (57°N). There is good agreement between modelled and observed values of T mrt, with an overall correspondence of R 2 = 0.91 (p < 0.01, RMSE = 3.1 K). A small overestimation of T mrt is found at locations shadowed by vegetation. Given this good performance a number of suggestions for future development are identified for applications which include for human comfort, building design, planning and evaluation of instrument exposure.

Contribution of land use changes to near-surface air temperatures during recent summer extreme heat events in the Phoenix Metropolitan Area

The impact of 1973–2005 land use–land cover (LULC) changes on near-surface air temperatures during four recent summer extreme heat events (EHEs) are investigated for the arid Phoenix, Arizona, metropolitan area using the Weather Research and Forecasting Model (WRF) in conjunction with the Noah Urban Canopy Model. WRF simulations were carried out for each EHE using LULC for the years 1973, 1985, 1998, and 2005. Comparison of measured near-surface air temperatures and wind speeds for 18 surface stations in the region show a good agreement between observed and simulated data for all simulation periods. The results indicate consistent significant contributions of urban development and accompanying LULC changes to extreme temperatures for the four EHEs. Simulations suggest new urban developments caused an intensification and expansion of the area experiencing extreme temperatures but mainly influenced nighttime temperatures with an increase of up to 10 K. Nighttime temperatures in the existing urban core showed changes of up to 2 K with the ongoing LULC changes. Daytime temperatures were not significantly affected where urban development replaced desert land (increase by 1 K); however, maximum temperatures increased by 2–4 K when irrigated agricultural land was converted to suburban development. According to the model simulations, urban landscaping irrigation contributed to cooling by 0.5–1 K in maximum daytime as well as minimum nighttime 2-m air temperatures in most parts of the urban region. Furthermore, urban development led to a reduction of the already relatively weak nighttime winds and therefore a reduction in advection of cooler air into the city.

Simulated changes in the relationship between tropical ocean temperatures and the western African monsoon during the mid-Holocene

Results from nine coupled ocean-atmosphere simulations have been used to investigate changes in the relationship between the variability of monsoon precipitation over western Africa and tropical sea surface temperatures (SSTs) between the mid-Holocene and the present day. Although the influence of tropical SSTs on the African monsoon is generally overestimated in the control simulations, the models reproduce aspects of the observed modes of variability. Thus, most models reproduce the observed negative correlation between western Sahelian precipitation and SST anomalies in the eastern tropical Pacific, and many of them capture the positive correlation between SST anomalies in the eastern tropical Atlantic and precipitation over the Guinea coastal region. Although the response of individual model to the change in orbital forcing between 6 ka and present differs somewhat, eight of the models show that the strength of the teleconnection between SSTs in the eastern tropical Pacific and Sahelian precipitation is weaker in the mid-Holocene. Some of the models imply that this weakening was associated with a shift towards longer time periods (from 3–5 years in the control simulations toward 4–10 years in the mid-Holocene simulations). The simulated reduction in the teleconnection between eastern tropical Pacific SSTs and Sahelian precipitation appears to be primarily related to a reduction in the atmospheric circulation bridge between the Pacific and West Africa but, depending on the model, other mechanisms such as increased importance of other modes of tropical ocean variability or increased local recycling of monsoonal precipitation can also play a role.

The effect of (mis-specified) GARCH filters on the finite sample distribution of the BDS test

This paper considers the effect of using a GARCH filter on the properties of the BDS test statistic as well as a number of other issues relating to the application of the test. It is found that, for certain values of the user-adjustable parameters, the finite sample distribution of the test is far-removed from asymptotic normality. In particular, when data generated from some completely different model class are filtered through a GARCH model, the frequency of rejection of iid falls, often substantially. The implication of this result is that it might be inappropriate to use non-rejection of iid of the standardised residuals of a GARCH model as evidence that the GARCH model ‘fits’ the data.

Technical note: a simple procedure for removing temporal discontinuities in ERA-Interim upper stratospheric temperatures for use in nudged chemistry-climate model simulations

This note describes a simple procedure for removing unphysical temporal discontinuities in ERA-Interim upper stratospheric global mean temperatures in March 1985 and August 1998 that have arisen due to changes in satellite radiance data used in the assimilation. The derived temperature adjustments (offsets) are suitable for use in stratosphere-resolving chemistry-climate models that are nudged (relaxed) to ERA-Interim winds and temperatures. Simulations using a nudged version of the Canadian Middle Atmosphere Model (CMAM) show that the inclusion of the temperature adjustments produces temperature time series that are devoid of the large jumps in 1985 and 1998. Due to its strong temperature dependence, the simulated upper stratospheric ozone is also shown to vary smoothly in time, unlike in a nudged simulation without the adjustments where abrupt changes in ozone occur at the times of the temperature jumps. While the adjustments to the ERA-Interim temperatures remove significant artefacts in the nudged CMAM simulation, spurious transient effects that arise due to water vapour and persist for about 5 yr after the 1979 switch to ERA-Interim data are identified, underlining the need for caution when analysing trends in runs nudged to reanalyses.