Computing Fresnel integrals via modified trapezium rules

In this paper we propose methods for computing Fresnel integrals based on truncated trapezium rule approximations to integrals on the real line, these trapezium rules modified to take into account poles of the integrand near the real axis. Our starting point is a method for computation of the error function of complex argument due to Matta and Reichel (J Math Phys 34:298–307, 1956) and Hunter and Regan (Math Comp 26:539–541, 1972). We construct approximations which we prove are exponentially convergent as a function of N , the number of quadrature points, obtaining explicit error bounds which show that accuracies of 10−15 uniformly on the real line are achieved with N=12 , this confirmed by computations. The approximations we obtain are attractive, additionally, in that they maintain small relative errors for small and large argument, are analytic on the real axis (echoing the analyticity of the Fresnel integrals), and are straightforward to implement.

The ocean's role in polar climate change: asymmetric Arctic and Antarctic responses to greenhouse gas and ozone forcing

In recent decades, the Arctic has been warming and sea ice disappearing. By contrast, the Southern Ocean around Antarctica has been (mainly) cooling and sea-ice extent growing. We argue here that interhemispheric asymmetries in the mean ocean circulation, with sinking in the northern North Atlantic and upwelling around Antarctica, strongly influence the sea-surface temperature (SST) response to anthropogenic greenhouse gas (GHG) forcing, accelerating warming in the Arctic while delaying it in the Antarctic. Furthermore, while the amplitude of GHG forcing has been similar at the poles, significant ozone depletion only occurs over Antarctica. We suggest that the initial response of SST around Antarctica to ozone depletion is one of cooling and only later adds to the GHG-induced warming trend as upwelling of sub-surface warm water associated with stronger surface westerlies impacts surface properties. We organize our discussion around ‘climate response functions’ (CRFs), i.e. the response of the climate to ‘step’ changes in anthropogenic forcing in which GHG and/or ozone-hole forcing is abruptly turned on and the transient response of the climate revealed and studied. Convolutions of known or postulated GHG and ozone-hole forcing functions with their respective CRFs then yield the transient forced SST response (implied by linear response theory), providing a context for discussion of the differing warming/cooling trends in the Arctic and Antarctic. We speculate that the period through which we are now passing may be one in which the delayed warming of SST associated with GHG forcing around Antarctica is largely cancelled by the cooling effects associated with the ozone hole. By mid-century, however, ozone-hole effects may instead be adding to GHG warming around Antarctica but with diminished amplitude as the ozone hole heals. The Arctic, meanwhile, responding to GHG forcing but in a manner amplified by ocean heat transport, may continue to warm at an accelerating rate.

Climate at high-obliquity

The question of climate at high obliquity is raised in the context of both exoplanet studies (e.g. habitability) and paleoclimates studies (evidence for low-latitude glaciation during the Neoproterozoic and the ”Snowball Earth” hypothesis). States of high obliquity, φ, are distinctive in that, for φ ≥54◦, the poles receive more solar radiation in the annual mean than the Equator, opposite to the present day situation. In addition, the seasonal cycle of insolation is extreme, with the poles alternatively “facing” the sun and sheltering in the dark for months. The novelty of our approach is to consider the role of a dynamical ocean in controlling the surface climate at high obliquity, which in turn requires understanding of the surface winds patterns when temperature gradients are reversed. To address these questions, a coupled ocean-atmosphere-sea ice GCM configured on an aquaplanet is employed. Except for the absence of topography and modified obliquity, the set-up is Earth-like. Two large obliquities φ, 54◦ and 90◦, are compared to today’s Earth value, φ=23.5◦. Three key results emerge at high obliquity: 1) despite reversed temper- ature gradients, mid-latitudes surface winds are westerly and trade winds exist at the equator (as for φ=23.5◦) although the westerlies are confined to the summer hemisphere, 2) a habitable planet is possible with mid-latitude temperatures in the range 300-280 K and 3) a stable climate state with an ice cap limited to the equatorial region is unlikely. We clarify the dynamics behind these features (notably by an analysis of the potential vorticity structure and conditions for baroclinic instability of the atmosphere). Interestingly, we find that the absence of a stable partially glaciated state is critically linked to the absence of ocean heat transport during winter, a feature ultimately traced back to the high seasonality of baroclinic instability conditions in the atmosphere.

Characterising the use of urban space: a geochemical case study from Calleva Atrebatum (Silchester, Hampshire, UK) Insula IX during the late first/early second century AD

The geochemical analysis of soil samples from the Roman town of Calleva Atrebatum (Silchester, Hampshire, UK) has been undertaken in order to enhance our understanding of urban occupation during the late first/early second century AD. Samples taken from a variety of occupation deposits within several, contemporary timber buildings, including associated hearths, have been analysed using laboratory-based x-ray fluorescence for a suite of elements (Cu, Zn, Pb, Sr, P and Ca). The patterns of elemental enrichment seen across the site have allowed us to compare and contrast the buildings that were occupied during this time in an attempt to distinguish different uses, such as between domestic and work-space. Two of the buildings stand out as having high concentrations of elements which suggest that they were dirtier work spaces, whilst other buildings appear to be have lower chemical loadings suggesting they were cleaner.