Improved numerical dissipation for time integration algorithms in structural dynamics

A new family of unconditionally stable one-step methods for the direct integration of the equations of structural dynamics is introduced and is shown to possess improved algorithmic damping properties which can be continuously controlled. The new methods are compared with members of the Newmark family, and the Houbolt and Wilson methods.     

Solar Flares and Coronal Mass Ejections: A Statistically Determined Flare Flux – CME Mass Correlation

In an effort to examine the relationship between flare flux and corresponding CME mass, we temporally and spatially correlate all X-ray flares and CMEs in the LASCO and GOES archives from 1996 to 2006. We cross-reference 6733 CMEs having well-measured masses against 12 050 X-ray flares having position information as determined from their optical counterparts. For a given flare, we search in time for CMEs which occur 10 – 80 minutes afterward, and we further require the flare and CME to occur within ± 45° in position angle on the solar disk. There are 826 CME/flare pairs which fit these criteria. Comparing the flare fluxes with CME masses of these paired events, we find CME mass increases with flare flux, following an approximately log-linear, broken relationship: in the limit of lower flare fluxes, log (CME mass)∝0.68×log (flare flux), and in the limit of higher flare fluxes, log (CME mass)∝0.33×log (flare flux). We show that this broken power-law, and in particular the flatter slope at higher flare fluxes, may be due to an observational bias against CMEs associated with the most energetic flares: halo CMEs. Correcting for this bias yields a single power-law relationship of the form log (CME mass)∝0.70×log (flare flux). This function describes the relationship between CME mass and flare flux over at least 3 dex in flare flux, from ≈ 10⁻⁷ – 10⁻⁴ W m⁻².     

Gravitational microlensing: A parallel,large-data implementation

Gravitational lensing allows us to probe the structure of matter on a broad range of astronomical scales, and as light from a distant source traverses an intervening galaxy, compact matter such as planets, stars, and black holes act as individual lenses. The magnification from such microlensing results in rapid brightness fluctuations which reveal not only the properties of the lensing masses, but also the surface brightness distribution in the source. However, while the combination of deflections due to individual stars is linear, the resulting magnifications are highly non-linear, leading to significant computational challenges which currently limit the range of problems which can be tackled. This paper presents a new and novel implementation of a numerical approach to gravitational microlensing, increasing the scale of the problems that can be tackled by more than two orders of magnitude, opening up a new regime of astrophysically interesting problems.     

Super-resolution of THEMIS thermal infrared data: Compositional relationships of surface units below the 100 meter scale on Mars

The Thermal Emission Imaging System (THEMIS) has provided the highest spatial resolution (100 m/pixel) thermal infrared (TIR) data of the surface of Mars to date. These data have enabled the discovery of many small-scale compositional units and helped to better constrain surface processes operating at these scales. However, with higher-resolution visible (VIS) instruments revealing smaller-scale surficial differences, there exists a need to detect and map compositional variability using TIR data at scales below 100 m. Because it is unlikely there will be a higher-resolution TIR instrument sent to Mars in the near future, creative image processing techniques commonly classified under the umbrella of “super-resolution” can be employed to improve or enhance the spatial resolution of the THEMIS TIR data. These approaches typically integrate another higher-resolution dataset and can either be qualitative for visual appeal, quantitative for data accuracy, or some combination of both. The super-resolution approach presented here produces enhanced TIR images that are radiometrically accurate, but also visually appealing. For the technique to be successfully applied, multi- to hyper-spectral data from two different spectral regions are required (e.g., the THEMIS TIR and VIS datasets). The focus here is to introduce this new super-resolution methodology and demonstrate its ability using existing THEMIS IR and VIS data. The quartzofeldspathic deposit in northern Syrtis Major was selected because of the spectral variability detected using the original IR resolution data and to better constrain the relationship between the small-scale surface morphology and areal extent of the deposit as well as its formation process. Despite being associated with the central peaks of two craters, the results here show no positive correlation between the small rocky outcrops and the quartzofeldspathic unit. A gradational contact exists between the unit and basaltic sands within the intercrater eolian material. The super-resolution approach offers an alternative approach to traditional sub-pixel deconvolution identification and provides a higher-resolution IR dataset for thermophysical and spectral analysis on Mars.     

Hydrofracturing in response to the development of an overpressurised subglacial meltwater system during drumlin formation: an example from Anglesey,NW Wales

This paper presents the results of a detailed study of a complex hydrofracture system and host diamictons exposed within a longitudinal section through an elongate drumlin located to the west of Cemlyn Bay, Anglesey, NW Wales. This complex, laterally extensive sand, silt and clay filled hydrofracture system was active over a prolonged period and is thought to have developed beneath the Late Devensian (Weichselian) Irish Sea Ice Stream as it overrode this part of NW Anglesey. The sediment-fill to the hydrofracture system is deformed with kinematic indicators (folds, thrusts, augen) recording a SW-directed sense of shear, consistent with the regional ice flow direction across this part of the island. The lack of any geomorphological evidence for active retreat of the Irish Sea ice across Anglesey has led to the conclusion that hydrofracturing at the Cemlyn Bay site occurred within the bed of the Irish Sea Ice Stream whilst this relatively faster flowing corridor of ice was actively overriding the island. Shear imposed by the overriding ice led to the development of a subglacial shear zone which facilitated the propagation of the hydrofracture system with the laterally extensive feeder sills occurring parallel to Y-type Riedel shears. Although a subglacial setting beneath the active Irish Sea Ice Stream can be argued for the Cemlyn Bay hydrofracture system, its relationship to the formation of the ‘host’ drumlin remains uncertain. However, evidence presented here suggests that hydrofracturing may have occurred during the later stages or post landform development in response to the migration of overpressurised meltwater within the bed of the Irish Sea ice; possibly accompanying the local thinning and shutdown of the Irish Sea Ice Stream on Anglesey.     

Whipple telescope observations of LS I +61 303: 2004–2006

In this paper we present the results of the past two years observations on the galactic microquasar LS I +61 303 with the Whipple 10 m gamma-ray telescope. The recent MAGIC detection of the source between 200 GeV and 4 TeV suggests that the source is periodic with very high energy (VHE) gamma-ray emission linked to its orbital cycle. The entire 50-hour data set obtained with Whipple from 2004 to 2006 was analyzed with no reliable detection resulting. The upper limits obtained in the 2005–2006 season covered several of the same epochs as the MAGIC Telescope detections, albeit with lower sensitivity. Upper limits are placed on emission during the orbital phases of 0→0.1 and 0.8→1, phases which are not included in the MAGIC data set.      

An investigation of gravitational lens determinations of H0 in quintessence cosmologies

There is growing evidence that the majority of the energy density of the Universe is not baryonic or dark matter, but rather it resides in an exotic component with negative pressure. The nature of this 'quintessence' influences our view of the Universe, modifying angular diameter and luminosity distances. Here, we examine the influence of a quintessence component upon gravitational lens time-delays. As well as a static quintessence component, an evolving equation of state is also considered. It is found that the equation of state of the quintessence component and its evolution influence the value of the Hubble constant derived from gravitational lenses. However, the differences between evolving and non-evolving cosmologies are relatively small. We undertake a suite of Monte Carlo simulations to examine the potential constraints that can be placed on the universal equation of state from the monitoring of gravitational lens systems, and demonstrate that at least an order of magnitude more lenses than currently known will have to be discovered and analysed to accurately probe any quintessence component.     

The magnitude distribution, perihelion distribution and flux of long-period comets

The mass distribution and perihelion distribution of long-period comets are re-assessed. The mass distribution index is found to be 1.598±0.016 , indicating that the distribution is somewhat steeper than was obtained by previous analyses of an amalgam of all the available historical data. The number of long-period comets that have orbital perihelion distances, q , that fall in a specific q to q +d q range is found to be independent of q . It is also noted that the flux of long-period comets to the inner Solar system has remained constant throughout recorded history.
The number of long-period comets, , per 1-au interval of perihelion distance, per year, brighter than H , entering the inner Solar system is found to be given by log₁₀ =−2.607+0.359 H . It is therefore estimated that, for example, about 0.5, 30 and 2000 long-period comets with absolute magnitudes brighter than 0, 5 and 10 respectively pass the Sun on orbits with perihelion distances less than 2.0 au, every century.