RcppArmadillo : accelerating R with high-performance C++ linear algebra

The R statistical environment and language has demonstrated particular strengths for interactive development of statistical algorithms, as well as data modelling and visualisation. Its current implementation has an interpreter at its core which may result in a performance penalty in comparison to directly executing user algorithms in the native machine code of the host CPU. In contrast, the C++ language has no built-in visualisation capabilities, handling of linear algebra or even basic statistical algorithms; however, user programs are converted to high-performance machine code, ahead of execution. A new method avoids possible speed penalties in R by using the Rcpp extension package in conjunction with the Armadillo C++ matrix library. In addition to the inherent performance advantages of compiled code, Armadillo provides an easy-to-use template-based meta-programming framework, allowing the automatic pooling of several linear algebra operations into one, which in turn can lead to further speedups. With the aid of Rcpp and Armadillo, conversion of linear algebra centered algorithms from R to C++ becomes straightforward. The algorithms retains the overall structure as well as readability, all while maintaining a bidirectional link with the host R environment. Empirical timing comparisons of R and C++ implementations of a Kalman filtering algorithm indicate a speedup of several orders of magnitude.

Adaptive radiotherapy for virally mediated head and neck cancer

Purpose: Virally mediated head and neck cancers (VMHNC) often present with nodal involvement, and are generally considered radioresponsive, resulting in the need for a re-planning CT during radiotherapy (RT) in a subset of patients. We sought to identify a high-risk group based on nodal size to be evaluated in a future prospective adaptive RT trial. Methodology: Between 2005-2010, 121 patients with virally-mediated, node positive nasopharyngeal (EBV positive) or oropharyngeal (HPV positive) cancers, receiving curative intent RT were reviewed. Patients were analysed based on maximum size of the dominant node with a view to grouping them in varying risk categories for the need of re-planning. The frequency and timing of the re-planning scans were also evaluated. Results: Sixteen nasopharyngeal and 105 oropharyngeal tumours were reviewed. Twenty-five (21%) patients underwent a re-planning CT at a median of 22 (range, 0-29) fractions with 1 patient requiring re-planning prior to the commencement of treatment. Based on the analysis, patients were subsequently placed into 3 groups; ≤35mm (Group 1), 36-45mm (Group 2), ≥46mm (Group 3). Re-planning CT’s were performed in Group 1- 8/68 (11.8%), Group 2- 4/28 (14.3%), Group 3- 13/25 (52%). Sample size did not allow statistical analysis to detect a significant difference or exclusion of a lack of difference between the 3 groups. Conclusion: In this series, patients with VMHNC and nodal size > 46mm appear to be a high-risk group for the need of re-planning during a course of definitive radiotherapy. This finding will now be tested in a prospective adaptive RT study.

Adaptive radiotherapy for virally mediated head and neck cancer

Purpose: Virally mediated head and neck cancers (VMHNC) often present with nodal involvement, and are generally considered radioresponsive, resulting in the need for a re-planning CT during radiotherapy (RT) in a subset of patients. We sought to identify a high-risk group based on nodal size to be evaluated in a future prospective adaptive RT trial. Methodology: Between 2005-2010, 121 patients with virally-mediated, node positive nasopharyngeal (EBV positive) or oropharyngeal (HPV positive) cancers, receiving curative intent RT were reviewed. Patients were analysed based on maximum size of the dominant node with a view to grouping them in varying risk categories for the need of re-planning. The frequency and timing of the re-planning scans were also evaluated. Results: Sixteen nasopharyngeal and 105 oropharyngeal tumours were reviewed. Twenty-five (21%) patients underwent a re-planning CT at a median of 22 (range, 0-29) fractions with 1 patient requiring re-planning prior to the commencement of treatment. Based on the analysis, patients were subsequently placed into 3 groups; ≤35mm (Group 1), 36-45mm (Group 2), ≥46mm (Group 3). Re-planning CT’s were performed in Group 1- 8/68 (11.8%), Group 2- 4/28 (14.3%), Group 3- 13/25 (52%). Sample size did not allow statistical analysis to detect a significant difference or exclusion of a lack of difference between the 3 groups. Conclusion: In this series, patients with VMHNC and nodal size > 46mm appear to be a high-risk group for the need of re-planning during a course of definitive radiotherapy. This finding will now be tested in a prospective adaptive RT study.

Adaptive radiotherapy for virally mediated head and neck cancer

Purpose: Virally mediated head and neck cancers (VMHNC) often present with nodal involvement, and are generally considered radioresponsive, resulting in the need for plan adaptation during radiotherapy in a subset of patients. We sought to identify a high-risk group based on pre-treatment nodal size to be evaluated in a future prospective adaptive radiotherapy trial. Methodology: Between 2005-2010, 121 patients with virally-mediated, node positive nasopharyngeal or oropharyngeal cancers, receiving definitive radiotherapy were reviewed. Patients were analysed based on maximum size of the dominant node at diagnosis with a view to grouping them in varying risk categories for the need of re-planning. The frequency and timing of the re-planning scans were also evaluated. Results: Sixteen nasopharyngeal and 105 oropharyngeal tumours were reviewed. Twenty-five (21%) patients underwent a re-planning CT at a median of 22 (range, 0-29) fractions with 1 patient requiring re-planning prior to the commencement of treatment. Based on the analysis, patients were subsequently placed into 3 groups defined by pre-treatment nodal size; ≤ 35mm (Group 1), 36-45mm (Group 2), ≥ 46mm (Group 3). Applying these groups to the patient cohort, re-planning CT’s were performed in Group 1- 8/68 (11.8%), Group 2- 4/28 (14.3%), Group 3- 13/25 (52%). Conclusion: In this series, patients with VMHNC and nodal size > 46mm appear to be a high-risk group for the need of plan adaptation during a course of definitive radiotherapy. This finding will now be tested in a prospective adaptive radiotherapy study.

Nitrous oxide emissions from irrigated wheat in Australia : impact of irrigation management

Background and Aims: Irrigation management affects soil water dynamics as well as the soil microbial carbon and nitrogen turnover and potentially the biosphere-atmosphere exchange of greenhouse gasses (GHG). We present a study on the effect of three irrigation treatments on the emissions of nitrous oxide (N2O) from irrigated wheat on black vertisols in South-Eastern Queensland, Australia. Methods: Soil N2O fluxes from wheat were monitored over one season with a fully automated system that measured emissions on a sub-daily basis. Measurements were taken from 3 subplots for each treatment within a randomized split-plot design. Results: Highest N2O emissions occurred after rainfall or irrigation and the amount of irrigation water applied was found to influence the magnitude of these “emission pulses”. Daily N2O emissions varied from -0.74 to 20.46 g N2O-N ha-1 day-1 resulting in seasonal losses ranging from 0.43 to 0.75 kg N2O N ha-1 season -1 for the different irrigation treatments. Emission factors (EF = proportion of N fertilizer emitted as N2O) over the wheat cropping season, uncorrected for background emissions, ranged from 0.2 to 0.4% of total N applied for the different treatments. Highest seasonal N2O emissions were observed in the treatment with the highest irrigation intensity; however, the N2O intensity (N2O emission per crop yield) was highest in the treatment with the lowest irrigation intensity. Conclusions: Our data suggest that timing and amount of irrigation can effectively be used to reduce N2O losses from irrigated agricultural systems; however, in order to develop sustainable mitigation strategies the N2O intensity of a cropping system is an important concept that needs to be taken into account.

Characterization of in vitro Chlamydia muridarum persistence and utilization in an in vivo mouse model of Chlamydia vaccine

Problem: Chlamydia trachomatis genital tract infections are easily treated with antibiotics, however the majority of infections are asymptomatic and therefore untreated, highlighting the need for a vaccine. Because most infections are asymptomatic, vaccination could potentially be administered to individuals who may have an acute infection at that time. In such individuals the effect of vaccination on the existing infection is unknown; however one potential outcome could be the development of a persistent infection. In vitro chlamydial persistence has been well characterized in various strains, however there have been no reported studies in C. muridarum. Method of Study: We performed ultrastructural characterization, and transcriptome analysis of selected genes. We then used the transcriptional profiles of the selected genes to examine whether intranasal immunization of mice during an active genital infection would induce persistence in the upper reproductive tract of female mice. Results and Conclusions: We found that persistence developed in the oviducts of mice as a result of immunization. This is a significant finding, not only because it is the first time that C. muridarum persistence has been characterized in vitro, but also due to the fact that there is minimal characterization of in vivo persistence of any chlamydial species. This highlights the importance of the timing of vaccination in individuals.

Pulling apart the mid to late Cenozoic magmatic record of the Gulf of California : no room for the Comondú Arc

Understanding the link between tectonic-driven extensional faulting and volcanism is crucial from a hazard perspective in active volcanic environments, while ancient volcanic successions provide records on how volcanic eruption styles, compositions, magnitudes and frequencies can change in response to extension timing, distribution and intensity. Significantly, incorrect tectonic interpretations can be made when the spatial-temporal-compositional trends of, and source contributions to magmatism are not properly considered. This study draws on intimate relationships of volcanism and extension preserved in the Sierra Madre Occidental (SMO) and Gulf of California (GoC) regions of western Mexico. Here, a major Oligocene rhyolitic ignimbrite “flare-up” (>300,000 km3) switched to a dominantly bimodal and mixed effusive-explosive volcanic phase in the Early Miocene (~100,000 km3), associated with distributed extension and opening of numerous grabens. Rhyolitic dome fields were emplaced along graben edges and at intersections of cross-graben and graben-parallel structures during early stages of graben development. Concomitant with this change in rhyolite eruption style was a change in crustal source as revealed by zircon chronochemistry with rapid rates of rhyolite magma generation due to remelting of mid- to upper crustal, highly differentiated igneous rocks emplaced during earlier SMO magmatism. Extension became more focused ~18 Ma resulting in volcanic activity being localised along the site of GoC opening. This localised volcanism (known as the Comondú “arc”) was dominantly effusive and andesite-dacite in composition. This compositional change resulted from increased mixing of basaltic and rhyolitic magmas rather than fluid flux melting of the mantle wedge above the subducting Guadalupe Plate. A poor understanding of space-time relationships of volcanism and extension has thus led to incorrect past tectonic interpretations of Comondú-age volcanism.

Late Oligocene to middle Miocene rifting and syn-extensional magmatism in the southwestern Sierra Madre occidental, Mexico: the beginning of the Gulf of California rift

Although Basin and Range style extension affected several areas of western Mexico since the Late Eocene, extension in the Gulf of California region (the Gulf Extensional Province GEP) is thought to have started as subduction waned and ended at ~14 12.5 Ma. A general consensus also exists in considering the mid Miocene Comondú group as a suprasubduction volcanic arc. Our new integration of the geology of the south east Gulf region, backed by 43 new Ar Ar and U Pb mineral ages and geochemical studies, document a widespread phase of extension in the southern GEP between latest Oligocene and Early Miocene that subsequently focused in the region of the future Gulf in the Middle Miocene. Upper Oligocene to Lower Miocene rocks across the southern Sierra Madre Occidental (SMO)(northern Nayarit and southern Sinaloa) were affected by major ~N S to NNW striking normal faults prior to ~21 Ma. Then, between ~21 and 11 Ma, a system of NNW-SSE high angle extensional faults continued extending the southwestern side of the SMO. Rhyolitic domes, shallow intrusive bodies, and lesser basalts were emplaced along this extensional belt at 20-17 Ma. In northern Sinaloa, large grabens were floored by huge dome complexes at ~21-17 Ma and filled by continental sediments with interlayered basalts dated at 15-14 Ma, a setting and timing very similar to Sonora. Early to Middle Miocene volcanism, including the largely volcaniclastic Comondú strata in Baja California Sur, was thus emplaced in rift basins and was likely associated to decompression melting of upper mantle (inducing crustal partial melting) rather than to fluxing by fluids from the young subducting plate. Along the Nayarit and Sinaloa coast, flatlying basaltic lava flows dated at 11-10 Ma are exposed just above the present sea level. Here, crustal thickness is almost half that in the unextended core of the SMO, implying significant lithosphere stretching before ~11 Ma. Our study shows that rifting began much earlier than Late Miocene and provided a fundamental control on the style and composition of volcanism from at least 30 Ma. We envision a sustained period of lithospheric stretching and magmatism during which the pace and breadth of extension changed at ~20-18 Ma to be narrower and likely more rapid, and again at ~12.5 Ma, when the kinematics of rifting became more oblique.

How do supercontinents break up? — assessing the continental large igneous province (LIP) record of the break-up of Pangea

LIP emplacement is linked to the timing and evolution of supercontinental break-up. LIP-related break-up produces volcanic rifted margins, new and large (up to 108 km2) ocean basins, and new, smaller continents that undergo dispersal and potentially reassembly (e.g., India). However, not all continental LIPs lead to continental rupture. We analysed the <330 Ma continental LIP record(following final assembly of Pangea) to find relationships between LIP event attributes (e.g., igneous volume, extent, distance from pre-existing continental margin) and ocean basin attributes (e.g., length of new ocean basin/rifted margin) and how these varied during the progressive break up of Pangea. No correlation exists between LIP magnitude and size of the subsequent ocean basin or rifted margin. Our review suggests a three-phased break-up history of Pangea: 1) “Preconditioning” phase (∼330–200 Ma): LIP events (n=7) occurred largely around the supercontinental margin clustering today in Asia, with a low (<20%) rifting success rate. The Panjal Traps at ∼280 Ma may represent the first continental rupturing event of Pangea, resulting in continental ribboning along the Tethyan margin; 2) “Main Break-up” phase (∼200–100 Ma): numerous large LIP events(n=10) in the supercontinent interior, resulting in highly successful fragmentation (90%) and large, new ocean basins(e.g., Central/South Atlantic, Indian, >3000 km long); 3) “Waning” phase (∼100–0 Ma): Declining LIP magnitudes (n=6), greater proximity to continental margins (e.g., Madagascar, North Atlantic, Afro-Arabia, Sierra Madre) producing smaller ocean basins (<2600 km long). How Pangea broke up may thus have implications for earlier supercontinent reconstructions and LIP record.

Consequences of extensional tectonics on volcanic eruption style, compositions and source regions : new insights from the southern Sierra Madre Occidental and Gulf of California regions, western Mexico

Understanding the link between tectonic-driven extensional faulting and volcanism is crucial from a hazard perspective in active volcanic environments, while ancient volcanic successions provide records on how volcanic eruption styles, compositions, magnitudes and frequencies can change in response to extension timing, distribution and intensity. This study draws on intimate relationships of volcanism and extension preserved in the Sierra Madre Occidental (SMO) and Gulf of California (GoC) regions of western Mexico. Here, a major Oligocene rhyolitic ignimbrite “flare-up” (>300,000 km3) switched to a dominantly bimodal and mixed effusive-explosive volcanic phase in the Early Miocene (~100,000 km3), associated with distributed extension and opening of numerous grabens. Rhyolitic dome fields were emplaced along graben edges and at intersections of cross-graben and graben-parallel structures during early stages of graben development. Concomitant with this change in rhyolite eruption style was a change in crustal source as revealed by zircon chronochemistry with rapid rates of rhyolite magma generation due to remelting of mid- to upper crustal, highly differentiated igneous rocks emplaced during earlier SMO magmatism. Extension became more focused ~18 Ma resulting in volcanic activity being localised along the site of GoC opening. This localised volcanism (known as the Comondú “arc”) was dominantly effusive and andesite-dacite in composition. This compositional change resulted from increased mixing of basaltic and rhyolitic magmas rather than fluid flux melting of the mantle wedge above the subducting Guadalupe Plate. A poor understanding of space-time relationships of volcanism and extension has thus led to incorrect past tectonic interpretations of Comondú-age volcanism.

QGECE research on delineation of Australian geothermal resources

The geology/reservoir program of the Queensland Geothermal Energy Centre of Excellence (QGECE) has the mission to improve the existing knowledge and develop new innovative scientific approaches for the identification of geothermal resources in Australia, with a particular focus on Queensland. Specifically, the QGECE geology/reservoir program is currently (1) producing a comprehensive geochemical dataset for high heat producing rocks, (2) conducting detailed mineralogical and geochronological studies of granites and hydrothermal alteration minerals, and ; (3) investigating the Cooper Basin representing a superb natural laboratory for understanding of radiogenic heat enrichment process and possible involvement of mantle heat flow. Seven research projects have been established, which are being conducted largely as PhD studies. In the preliminary studies, high quality and valuable results were obtained to address the research topics of understanding the causes and timing of heat producing element enrichment.

Inter-cycle variability of in-cylinder pressure parameters in an ethanol fumigated common rail diesel engine

The effects of ethanol fumigation on the inter-cycle variability of key in-cylinder pressure parameters in a modern common rail diesel engine have been investigated. Specifically, maximum rate of pressure rise, peak pressure, peak pressure timing and ignition delay were investigated. A new methodology for investigating the start of combustion was also proposed and demonstrated—which is particularly useful with noisy in-cylinder pressure data as it can have a significant effect on the calculation of an accurate net rate of heat release indicator diagram. Inter-cycle variability has been traditionally investigated using the coefficient of variation. However, deeper insight into engine operation is given by presenting the results as kernel density estimates; hence, allowing investigation of otherwise unnoticed phenomena, including: multi-modal and skewed behaviour. This study has found that operation of a common rail diesel engine with high ethanol substitutions (>20% at full load, >30% at three quarter load) results in a significant reduction in ignition delay. Further, this study also concluded that if the engine is operated with absolute air to fuel ratios (mole basis) less than 80, the inter-cycle variability is substantially increased compared to normal operation.

Algorithm for queue estimation with loop detector of time occupancy in off-ramps on signalized motorways

Long traffic queues on off-ramps significantly compromise the safety and throughput of motorways. Obtaining accurate queue information is crucial for countermeasure strategies. However, it is challenging to estimate traffic queues with locally installed inductive loop detectors. This paper deals with the problem of queue estimation with the interpretation of queuing dynamics and the corresponding time-occupancy distribution over motorway off-ramps. A novel algorithm for real-time queue estimation with two detectors is presented and discussed. Results derived from microscopic traffic simulation validated the effectiveness of the algorithm and revealed some of its useful features: (a) long and intermediate traffic queues could be accurately measured, (b) relatively simple detector input (i.e., time occupancy) was required, and (c) the estimation philosophy was independent with signal timing changes and provided the potential to cooperate with advanced strategies for signal control. Some issues concerning field implementation are also discussed.

Optimized strategy for integrated traffic and transit signal control

This paper proposes a unique and innovative approach to integrate transit signal priority control into a traffic adaptive signal control strategy. The proposed strategy was named OSTRAC (Optimized Strategy for integrated TRAffic and TRAnsit signal Control). The cornerstones of OSTRAC include an online microscopic traffic f low prediction model and a Genetic Algorithm (GA) based traffic signal timing module. A sensitivity analysis was conducted to determine the critical GA parameters. The developed traffic f low model demonstrated reliable prediction results through a test. OSTRAC was evaluated by comparing its performance to three other signal control strategies. The evaluation results revealed that OSTRAC efficiently and effectively reduced delay time of general traffic and also transit vehicles.

Morphometry and abnormalities of the feet of Kaimanawa feral horses in New Zealand

Objective: The present study investigated the foot health of the Kaimanawa feral horse population and tested the hypotheses that horses would have a large range of foot morphology and that the incidence of foot abnormality would be significantly high. Procedures: Abnormality was defined as a variation from what the two veterinarian assessors considered as optimal morphology and which was considered to impact negatively on the structure and/or function of the foot. Fifteen morphometric variables were measured on four calibrated photographic views of all four feet of 20 adult Kaimanawa feral horses. Four morphometric variables were measured from the lateromedial radiographs of the left forefoot of each horse. In addition, the study identified the incidence of gross abnormality observed on the photographs and radiographs of all 80 feet. Results: There was a large variation between horses in the morphometric dimensions, indicating an inconsistent foot type. Mean hoof variables were outside the normal range recommended by veterinarians and hoof care providers; 35% of all feet had a long toe conformation and 15% had a mediolateral imbalance. Abnormalities included lateral (85% of horses) and dorsal (90% of horses) wall flares, presence of laminar rings (80% of horses) and bull-nose tip of the distal phalanx (75% of horses). Both hypotheses were therefore accepted. Conclusions: The Kaimanawa feral horse population demonstrated a broad range of foot abnormalities and we propose that one reason for the questionable foot health and conformation is lack of abrasive wearing by the environment. In comparison with other feral horse populations in Australia and America there may be less pressure on the natural selection of the foot of the Kaimanawa horses by the forgiving environment of the Kaimanawa Ranges. Contrary to popular belief, the feral horse foot type should not be used as an ideal model for the domestic horse foot.