Evaluation of a numerical model of the British–Irish ice sheet using relative sea‐level data: implications for the interpretation of trimline observations

The glacial isostatic adjustment (GIA) of the British Isles is complex due to the interplay between local and non‐local signals. A number of recent studies have modelled the GIA response of the British Isles using relative sea‐level data. This study extends these previous analyses by using output from a numerical glaciological model as input to a GIA model. This is a necessary step towards more realistic GIA models, and although there have been similar studies for the major late Pleistocene ice sheets, this is the first study to do so for the British Isles. We test three reconstructions, classed as ‘minimal’, ‘median’ and ‘maximal’ in terms of their volume at maximum extent, and find it is possible to obtain good data–model fits. The minimal reconstruction is clearly preferred by the sea‐level data. The ice reconstructions tested were not constrained by geomorphological information of past ice extent (lateral and vertical). As a consequence, the reconstructions extend further than much of this information suggests, particularly in terms of ice thickness. It is notable, however, that the reconstructions produce good fits to many sea‐level data from central, mountainous regions (e.g. Scottish highlands), which lends support to the suggestion that trimlines, often used as an constraint on the palaeo ice surface location, are in fact features formed at the transition from warm‐ to cold‐based ice and so mark a minimum constraint on the ice surface altitude. Based on data–model misfits, suggestions are made for improving the ice model reconstructions. However, in many locations, the cause of the misfit could be due to non‐local, predominantly Fennoscandian ice and so interpretation is not straightforward. As a result, we suggest that future analyses of this type consider models and observations for both Fennoscandia and the British Isles. Copyright © 2012 John Wiley & Sons, Ltd.     

Understanding soil‐pipe flow and its role in ephemeral gully erosion

The role of soil pipeflow in ephemeral gully erosion is not well understood. Experiments were conducted on continuous soil pipes to better understand the role of internal erosion of soil pipes and its relation to ephemeral gully development. Soil beds of 140 cm length, 100 cm width and 20 cm depth had a single soil pipe of different initial sizes (2, 4, 6, 8, and 10 mm diameter) extend from a water reservoir to the outlet. Experiments were run on Providence silt loam and Smithdale loam soils under a constant head of 15 cm established for 30 min. Either the tunnel collapsed or the head could not be maintained. Soil pipes that were initially 2 and 4 mm clogged instantaneously at their mouth and did not exhibit flow, whereas, pipes initially ≥ 6 mm enlarged by 268, 397, and 699% on average for the 6, 8, and 10 mm diameters, respectively. Critical shear stress values were found to be essentially zero, and erodibility values gave erosion indexes that were extremely high. The rapid internal erosion resulted in erratic flow and sediment concentrations with periods of no flow as pipes were temporarily clogged followed by surges of high flow and high sediment concentrations. Tensiometers within 6 cm of the soil pipes did not exhibit pressure increases typically associated with pipe clogging. Flow through 10 mm diameter soil pipes exhibited tunnel collapse for both soils tested. Tunnel collapse typically occurred within minutes of flow establishment suggesting that ephemeral gullies could be misinterpreted as being caused by convergent surface flow if observations were made after the runoff event instead of when flow is first established through soil pipes. Copyright © 2011 John Wiley & Sons, Ltd.     

The ultraluminous and hyperluminous infrared galaxies in the Sloan Digital Sky Survey, 2dF Galaxy Redshift Survey and 6dF Galaxy Survey

We present a result of cross-correlating the Infrared Astronomical Satellite Faint Source Catalogue with the spectroscopic catalogues of galaxies in the Fourth Data Release of the Sloan Digital Sky Survey, the Final Data Release of the 2dF Galaxy Redshift Survey (2dFGRS) and the Second Data Release of the 6dF Galaxy Survey. We have identified 324 ultraluminous infrared galaxies (ULIRGs) including 190 newly discovered ULIRGs, and two hyperluminous infrared galaxies. Adding these new ULIRGs, we increase the number of known ULIRGs by about 30 per cent. The reliability of the cross-correlation is estimated using the likelihood ratio method. The incompleteness of our sample introduced by the identification procedure in this study is estimated to be about 5 per cent. Our sample covers the redshift range of   z = 0.037–0.517  with a median redshift of     , which is larger than that     of the sample of previously known ULIRGs.     

Postglacial relative sea‐level observations from Ireland and their role in glacial rebound modelling

The British Isles have been the focus of a number of recent modelling studies owing to the existence of a high‐quality sea‐level dataset for this region and the suitability of these data for constraining shallow earth viscosity structure, local to regional ice sheet histories and the magnitude/timing of global meltwater signals. Until recently, the paucity of both glaciological and relative sea‐level (RSL) data from Ireland has meant that the majority of these glacial isostatic adjustment (GIA) modelling studies of the British Isles region have tended to concentrate on reconstructing ice cover over Britain. However, the recent development of a sea‐level database for Ireland along with emergence of new glaciological data on the spatial extent, thickness and deglacial chronology of the Irish Ice Sheet means it is now possible to revisit this region of the British Isles. Here, we employ these new data to constrain the evolution of the Irish Ice Sheet. We find that in order to reconcile differences between model predictions and RSL evidence, a thick, spatially extensive ice sheet of ～600–700 m over much of north and central Ireland is required at the LGM with very rapid deglaciation after 21 k cal. yr BP. Copyright © 2007 John Wiley & Sons, Ltd.     

An evaluation of temporally adaptive transformation approaches for solving Richards' equation

Developing robust and efficient numerical solution methods for Richards' equation (RE) continues to be a challenge for certain problems. We consider such a problem here: infiltration into unsaturated porous media initially at static conditions for uniform and non-uniform pore size media. For ponded boundary conditions, a sharp infiltration front results, which propagates through the media. We evaluate the resultant solution method for robustness and efficiency using combinations of variable transformation and adaptive time-stepping methods. Transformation methods introduce a change of variable that results in a smoother solution, which is more amenable to efficient numerical solution. We use adaptive time-stepping methods to adjust the time-step size, and in some cases the order of the solution method, to meet a constraint on nonlinear solution convergence properties or a solution error criterion. Results for three test problems showed that adaptive time-stepping methods provided robust solutions; in most cases transforming the dependent variable led to more efficient solutions than untransformed approaches, especially as the pore-size uniformity increased; and the higher-order adaptive time integration method was robust and the most efficient method evaluated.     

Revised ab initio models for H2-H2 collision-induced absorption at low temperatures

A revised ab initio calculation of the H₂-H₂ collision-induced absorption results in significant differences compared with the work of J. Borysow et al. [Borysow, J., Trafton, L., Frommhold, L., Birnbaum, G., 1985. Astrophys. J. 296, 644-654] for wavenumbers greater than 600 cm⁻¹ and temperatures below 120 K. The revision has significant influence on the spectra of Uranus and Neptune, and essentially removes the need for models with “super-solar” helium abundances or stratospheric hazes to explain the spectrum of Uranus.     

Performance characteristics of a new electronic snow water equivalent sensor in different climates

The US Army ERDC CRREL and the US Department of Agriculture Natural Resources Conservation Service developed a square electronic snow water equivalent (e‐SWE) sensor as an alternative to using fluid‐filled snow pillows to measure SWE. The sensors consist of a centre panel to measure SWE and eight outer panels to buffer edge stress concentrations. Seven 3 m square e‐SWE sensors were installed in five different climate zones. During the 2011–2012 winter, 1.8 and 1.2 m square e‐SWE sensors were installed and operated in Oregon. With the exception of New York State and Newfoundland, the e‐SWE sensors accurately measured SWE, with R² values between the sensor and manual SWE measurements of between 0.86 and 0.98. The e‐SWE sensor at Hogg Pass, Oregon, accurately measured SWE during the past 8 years of operations. In the thin, icy snow of New York during midwinter 2008–2009, the e‐SWE sensors overmeasured SWE because of edge stress concentrations associated with strong icy layers and a shallow snow cover. The New York e‐SWE sensors' measurement accuracy improved in spring 2009 and further improved during the 2011–2012 winter with operating experience. At Santiam Junction, measured SWE from the 1.8 and 1.2 m square e‐SWE sensors agreed well with the snow pillow, 3 m square e‐SWE sensor, and manual SWE measurements until February 2013, when dust and gravel blew onto the testing area resulting in anomalous measurements. © 2014 The Authors. Hydrological Processes published by John Wiley & Sons Ltd.     

Timeline analysis and wavelet multiscale analysis of the AKARI All-Sky Survey at 90 μm

We present a careful analysis of the point-source detection limit of the AKARI All-Sky Survey in the WIDE-S 90-μm band near the North Ecliptic Pole (NEP). Timeline analysis is used to detect IRAS ( Infrared Astronomy Satellite ) sources and then a conversion factor is derived to transform the peak timeline signal to the interpolated 90-μm flux of a source. Combined with a robust noise measurement, the point-source flux detection limit at signal-to-noise ratio  (S/N) > 5  for a single detector row is  1.1 ± 0.1 Jy  which corresponds to a point-source detection limit of the survey of ∼0.4 Jy.
Wavelet transform offers a multiscale representation of the Time Series Data ( tsd ). We calculate the continuous wavelet transform of the tsd and then search for significant wavelet coefficients considered as potential source detections. To discriminate real sources from spurious or moving objects, only sources with confirmation are selected. In our multiscale analysis, IRAS sources selected above 4σ can be identified as the only real sources at the Point Source Scales. We also investigate the correlation between the non- IRAS sources detected in timeline analysis and cirrus emission using wavelet transform and contour plots of wavelet power spectrum. It is shown that the non- IRAS sources are most likely to be caused by excessive noise over a large range of spatial scales rather than real extended structures such as cirrus clouds.