Using Genetics to Assess Differentiation Among Raccoons in an Area with Variable Rabies Status in Alabama

The western spread of raccoon rabies in Alabama has been slow and even appears to regress eastward periodically. While the disease has been present in the state for over 30 years, areas in northwest Alabama are devoid of raccoon rabies. This variation resulting in an enzootic area of raccoon rabies primarily in southeastern Alabama may be due to landscape features that hinder the movement of raccoons (i.e., gene flow) among different locations. We used 11 raccoon-specific microsatellite markers to obtain individual genotypes to examine gene flow among areas that were rabies free, enzootic with rabies, or had only sporadic reports of the disease. Samples from 70 individuals were collected from 5 sampling localities in 3 counties. The landscape feature data were collected from geographic information system (GIS) data. We inferred gene flow by estimating FST and by using Bayesian tests to identify genetic clusters. Estimates of pairwise FST indicated genetic differentiation and restricted gene flow between some sites, and an uneven distribution of genetic clusters was observed. Of the landscape features examined (i.e., land cover, elevation, slope, roads, and hydrology), only land cover had an association with genetic differentiation, suggesting this landscape variable may affect gene flow among raccoon populations and thus the spread of raccoon variant of rabies in Alabama.

Studying the influence of a secondary variable in Collocated Cokriging estimates

In this paper the influence of a secondary variable as a function of the correlation with the primary variable for collocated cokriging is examined. For this study five exhaustive data sets were generated in computer, from which samples with 60 and 104 data points were drawn using the stratified random sampling method. These exhaustive data sets were generated departing from a pair of primary and secondary variables showing a good correlation. Then successive sets were generated by adding an amount of white noise in such a way that the correlation gets poorer. Using these samples, it was possible to find out how primary and secondary information is used to estimate an unsampled location according to the correlation level.

The assessment of change in rural landscape and the analysis of the driving forces. Proposal of a spatial model for the rural built environment

The presented study carried out an analysis on rural landscape changes. In particular the study focuses on the understanding of driving forces acting on the rural built environment using a statistical spatial model implemented through GIS techniques. It is well known that the study of landscape changes is essential for a conscious decision making in land planning. From a bibliography review results a general lack of studies dealing with the modeling of rural built environment and hence a theoretical modelling approach for such purpose is needed. The advancement in technology and modernity in building construction and agriculture have gradually changed the rural built environment. In addition, the phenomenon of urbanization of a determined the construction of new volumes that occurred beside abandoned or derelict rural buildings. Consequently there are two types of transformation dynamics affecting mainly the rural built environment that can be observed: the conversion of rural buildings and the increasing of building numbers. It is the specific aim of the presented study to propose a methodology for the development of a spatial model that allows the identification of driving forces that acted on the behaviours of the building allocation. In fact one of the most concerning dynamic nowadays is related to an irrational expansion of buildings sprawl across landscape. The proposed methodology is composed by some conceptual steps that cover different aspects related to the development of a spatial model: the selection of a response variable that better describe the phenomenon under study, the identification of possible driving forces, the sampling methodology concerning the collection of data, the most suitable algorithm to be adopted in relation to statistical theory and method used, the calibration process and evaluation of the model. A different combination of factors in various parts of the territory generated favourable or less favourable conditions for the building allocation and the existence of buildings represents the evidence of such optimum. Conversely the absence of buildings expresses a combination of agents which is not suitable for building allocation. Presence or absence of buildings can be adopted as indicators of such driving conditions, since they represent the expression of the action of driving forces in the land suitability sorting process. The existence of correlation between site selection and hypothetical driving forces, evaluated by means of modeling techniques, provides an evidence of which driving forces are involved in the allocation dynamic and an insight on their level of influence into the process. GIS software by means of spatial analysis tools allows to associate the concept of presence and absence with point futures generating a point process. Presence or absence of buildings at some site locations represent the expression of these driving factors interaction. In case of presences, points represent locations of real existing buildings, conversely absences represent locations were buildings are not existent and so they are generated by a stochastic mechanism. Possible driving forces are selected and the existence of a causal relationship with building allocations is assessed through a spatial model. The adoption of empirical statistical models provides a mechanism for the explanatory variable analysis and for the identification of key driving variables behind the site selection process for new building allocation. The model developed by following the methodology is applied to a case study to test the validity of the methodology. In particular the study area for the testing of the methodology is represented by the New District of Imola characterized by a prevailing agricultural production vocation and were transformation dynamic intensively occurred. The development of the model involved the identification of predictive variables (related to geomorphologic, socio-economic, structural and infrastructural systems of landscape) capable of representing the driving forces responsible for landscape changes.. The calibration of the model is carried out referring to spatial data regarding the periurban and rural area of the study area within the 1975-2005 time period by means of Generalised linear model. The resulting output from the model fit is continuous grid surface where cells assume values ranged from 0 to 1 of probability of building occurrences along the rural and periurban area of the study area. Hence the response variable assesses the changes in the rural built environment occurred in such time interval and is correlated to the selected explanatory variables by means of a generalized linear model using logistic regression. Comparing the probability map obtained from the model to the actual rural building distribution in 2005, the interpretation capability of the model can be evaluated. The proposed model can be also applied to the interpretation of trends which occurred in other study areas, and also referring to different time intervals, depending on the availability of data. The use of suitable data in terms of time, information, and spatial resolution and the costs related to data acquisition, pre-processing, and survey are among the most critical aspects of model implementation. Future in-depth studies can focus on using the proposed model to predict short/medium-range future scenarios for the rural built environment distribution in the study area. In order to predict future scenarios it is necessary to assume that the driving forces do not change and that their levels of influence within the model are not far from those assessed for the time interval used for the calibration.

Variable Na(v)1.5 protein expression from the wild-type allele correlates with the penetrance of cardiac conduction disease in the Scn5a(+/-) mouse model

BACKGROUND: Loss-of-function mutations in SCN5A, the gene encoding Na(v)1.5 Na+ channel, are associated with inherited cardiac conduction defects and Brugada syndrome, which both exhibit variable phenotypic penetrance of conduction defects. We investigated the mechanisms of this heterogeneity in a mouse model with heterozygous targeted disruption of Scn5a (Scn5a(+/-) mice) and compared our results to those obtained in patients with loss-of-function mutations in SCN5A. METHODOLOGY/PRINCIPAL FINDINGS: Based on ECG, 10-week-old Scn5a(+/-) mice were divided into 2 subgroups, one displaying severe ventricular conduction defects (QRS interval>18 ms) and one a mild phenotype (QRS< or = 18 ms; QRS in wild-type littermates: 10-18 ms). Phenotypic difference persisted with aging. At 10 weeks, the Na+ channel blocker ajmaline prolonged QRS interval similarly in both groups of Scn5a(+/-) mice. In contrast, in old mice (>53 weeks), ajmaline effect was larger in the severely affected subgroup. These data matched the clinical observations on patients with SCN5A loss-of-function mutations with either severe or mild conduction defects. Ventricular tachycardia developed in 5/10 old severely affected Scn5a(+/-) mice but not in mildly affected ones. Correspondingly, symptomatic SCN5A-mutated Brugada patients had more severe conduction defects than asymptomatic patients. Old severely affected Scn5a(+/-) mice but not mildly affected ones showed extensive cardiac fibrosis. Mildly affected Scn5a(+/-) mice had similar Na(v)1.5 mRNA but higher Na(v)1.5 protein expression, and moderately larger I(Na) current than severely affected Scn5a(+/-) mice. As a consequence, action potential upstroke velocity was more decreased in severely affected Scn5a(+/-) mice than in mildly affected ones. CONCLUSIONS: Scn5a(+/-) mice show similar phenotypic heterogeneity as SCN5A-mutated patients. In Scn5a(+/-) mice, phenotype severity correlates with wild-type Na(v)1.5 protein expression.

Estimation of the postmortem interval by means of ¹H MRS of decomposing brain tissue: influence of ambient temperature

Standard methods for the estimation of the postmortem interval (PMI, time since death), based on the cooling of the corpse, are limited to about 48 h after death. As an alternative, noninvasive postmortem observation of alterations of brain metabolites by means of (1)H MRS has been suggested for an estimation of the PMI at room temperature, so far without including the effect of other ambient temperatures. In order to study the temperature effect, localized (1)H MRS was used to follow brain decomposition in a sheep brain model at four different temperatures between 4 and 26°C with repeated measurements up to 2100 h postmortem. The simultaneous determination of 25 different biochemical compounds at each measurement allowed the time courses of concentration changes to be followed. A sudden and almost simultaneous change of the concentrations of seven compounds was observed after a time span that decreased exponentially from 700 h at 4°C to 30 h at 26°C ambient temperature. As this represents, most probably, the onset of highly variable bacterial decomposition, and thus defines the upper limit for a reliable PMI estimation, data were analyzed only up to this start of bacterial decomposition. As 13 compounds showed unequivocal, reproducible concentration changes during this period while eight showed a linear increase with a slope that was unambiguously related to ambient temperature. Therefore, a single analytical function with PMI and temperature as variables can describe the time courses of metabolite concentrations. Using the inverse of this function, metabolite concentrations determined from a single MR spectrum can be used, together with known ambient temperatures, to calculate the PMI of a corpse. It is concluded that the effect of ambient temperature can be reliably included in the PMI determination by (1)H MRS.

Mycobacterium tuberculosis transmission in a country with low tuberculosis incidence: role of immigration and HIV infection

Immigrants from high-burden countries and HIV-coinfected individuals are risk groups for tuberculosis (TB) in countries with low TB incidence. Therefore, we studied their role in transmission of Mycobacterium tuberculosis in Switzerland. We included all TB patients from the Swiss HIV Cohort and a sample of patients from the national TB registry. We identified molecular clusters by spoligotyping and mycobacterial interspersed repetitive-unit-variable-number tandem-repeat (MIRU-VNTR) analysis and used weighted logistic regression adjusted for age and sex to identify risk factors for clustering, taking sampling proportions into account. In total, we analyzed 520 TB cases diagnosed between 2000 and 2008; 401 were foreign born, and 113 were HIV coinfected. The Euro-American M. tuberculosis lineage dominated throughout the study period (378 strains; 72.7%), with no evidence for another lineage, such as the Beijing genotype, emerging. We identified 35 molecular clusters with 90 patients, indicating recent transmission; 31 clusters involved foreign-born patients, and 15 involved HIV-infected patients. Birth origin was not associated with clustering (adjusted odds ratio [aOR], 1.58; 95% confidence interval [CI], 0.73 to 3.43; P = 0.25, comparing Swiss-born with foreign-born patients), but clustering was reduced in HIV-infected patients (aOR, 0.49; 95% CI, 0.26 to 0.93; P = 0.030). Cavitary disease, male sex, and younger age were all associated with molecular clustering. In conclusion, most TB patients in Switzerland were foreign born, but transmission of M. tuberculosis was not more common among immigrants and was reduced in HIV-infected patients followed up in the national HIV cohort study. Continued access to health services and clinical follow-up will be essential to control TB in this population.

Estimation with Interval Censored Data in Longitudinal Studies

In biostatistical applications interest often focuses on the estimation of the distribution of a time-until-event variable T. If one observes whether or not T exceeds an observed monitoring time at a random number of monitoring times, then the data structure is called interval censored data. We extend this data structure by allowing the presence of a possibly time-dependent covariate process that is observed until end of follow up. If one only assumes that the censoring mechanism satisfies coarsening at random, then, by the curve of dimensionality, typically no regular estimators will exist. To fight the curse of dimensionality we follow the approach of Robins and Rotnitzky (1992) by modeling parameters of the censoring mechanism. We model the right-censoring mechanism by modeling the hazard of the follow up time, conditional on T and the covariate process. For the monitoring mechanism we avoid modeling the joint distribution of the monitoring times by only modeling a univariate hazard of the pooled monitoring times, conditional on the follow up time, T, and the covariates process, which can be estimated by treating the pooled sample of monitoring times as i.i.d. In particular, it is assumed that the monitoring times and the right-censoring times only depend on T through the observed covariate process. We introduce inverse probability of censoring weighted (IPCW) estimator of the distribution of T and of smooth functionals thereof which are guaranteed to be consistent and asymptotically normal if we have available correctly specified semiparametric models for the two hazards of the censoring process. Furthermore, given such correctly specified models for these hazards of the censoring process, we propose a one-step estimator which will improve on the IPCW estimator if we correctly specify a lower-dimensional working model for the conditional distribution of T, given the covariate process, that remains consistent and asymptotically normal if this latter working model is misspecified. It is shown that the one-step estimator is efficient if each subject is at most monitored once and the working model contains the truth. In general, it is shown that the one-step estimator optimally uses the surrogate information if the working model contains the truth. It is not optimal in using the interval information provided by the current status indicators at the monitoring times, but simulations in Peterson, van der Laan (1997) show that the efficiency loss is small.

Precise pose recovery of distal locking holes from single calibrated fluoroscopic image via a novel variable decomposition approach

This paper presents a novel variable decomposition approach for pose recovery of the distal locking holes using single calibrated fluoroscopic image. The problem is formulated as a model-based optimal fitting process, where the control variables are decomposed into two sets: (a) the angle between the nail axis and its projection on the imaging plane, and (b) the translation and rotation of the geometrical model of the distal locking hole around the nail axis. By using an iterative algorithm to find the optimal values of the latter set of variables for any given value of the former variable, we reduce the multiple-dimensional model-based optimal fitting problem to a one-dimensional search along a finite interval. We report the results of our in vitro experiments, which demonstrate that the accuracy of our approach is adequate for successful distal locking of intramedullary nails.

Modeling and classifying variable width riparian zones utilizing digital elevation models, flood height data, digital soil data and national wetlands inventory : a new approach for riparian zone delineation

Riparian zones are dynamic, transitional ecosystems between aquatic and terrestrial ecosystems with well defined vegetation and soil characteristics. Development of an all-encompassing definition for riparian ecotones, because of their high variability, is challenging. However, there are two primary factors that all riparian ecotones are dependent on: the watercourse and its associated floodplain. Previous approaches to riparian boundary delineation have utilized fixed width buffers, but this methodology has proven to be inadequate as it only takes the watercourse into consideration and ignores critical geomorphology, associated vegetation and soil characteristics. Our approach offers advantages over other previously used methods by utilizing: the geospatial modeling capabilities of ArcMap GIS; a better sampling technique along the water course that can distinguish the 50-year flood plain, which is the optimal hydrologic descriptor of riparian ecotones; the Soil Survey Database (SSURGO) and National Wetland Inventory (NWI) databases to distinguish contiguous areas beyond the 50-year plain; and land use/cover characteristics associated with the delineated riparian zones. The model utilizes spatial data readily available from Federal and State agencies and geospatial clearinghouses. An accuracy assessment was performed to assess the impact of varying the 50-year flood height, changing the DEM spatial resolution (1, 3, 5 and 10m), and positional inaccuracies with the National Hydrography Dataset (NHD) streams layer on the boundary placement of the delineated variable width riparian ecotones area. The result of this study is a robust and automated GIS based model attached to ESRI ArcMap software to delineate and classify variable-width riparian ecotones.

Variable impact on mortality of AIDS-defining events diagnosed during combination antiretroviral therapy: not all AIDS-defining conditions are created equal

BACKGROUND: The extent to which mortality differs following individual acquired immunodeficiency syndrome (AIDS)-defining events (ADEs) has not been assessed among patients initiating combination antiretroviral therapy. METHODS: We analyzed data from 31,620 patients with no prior ADEs who started combination antiretroviral therapy. Cox proportional hazards models were used to estimate mortality hazard ratios for each ADE that occurred in >50 patients, after stratification by cohort and adjustment for sex, HIV transmission group, number of antiretroviral drugs initiated, regimen, age, date of starting combination antiretroviral therapy, and CD4+ cell count and HIV RNA load at initiation of combination antiretroviral therapy. ADEs that occurred in <50 patients were grouped together to form a "rare ADEs" category. RESULTS: During a median follow-up period of 43 months (interquartile range, 19-70 months), 2880 ADEs were diagnosed in 2262 patients; 1146 patients died. The most common ADEs were esophageal candidiasis (in 360 patients), Pneumocystis jiroveci pneumonia (320 patients), and Kaposi sarcoma (308 patients). The greatest mortality hazard ratio was associated with non-Hodgkin's lymphoma (hazard ratio, 17.59; 95% confidence interval, 13.84-22.35) and progressive multifocal leukoencephalopathy (hazard ratio, 10.0; 95% confidence interval, 6.70-14.92). Three groups of ADEs were identified on the basis of the ranked hazard ratios with bootstrapped confidence intervals: severe (non-Hodgkin's lymphoma and progressive multifocal leukoencephalopathy [hazard ratio, 7.26; 95% confidence interval, 5.55-9.48]), moderate (cryptococcosis, cerebral toxoplasmosis, AIDS dementia complex, disseminated Mycobacterium avium complex, and rare ADEs [hazard ratio, 2.35; 95% confidence interval, 1.76-3.13]), and mild (all other ADEs [hazard ratio, 1.47; 95% confidence interval, 1.08-2.00]). CONCLUSIONS: In the combination antiretroviral therapy era, mortality rates subsequent to an ADE depend on the specific diagnosis. The proposed classification of ADEs may be useful in clinical end point trials, prognostic studies, and patient management.

Pharmacogenetics-based population pharmacokinetic analysis of etravirine in HIV-1 infected individuals

Objectives: Etravirine (ETV) is metabolized by cytochrome P450 (CYP) 3A, 2C9, and 2C19. Metabolites are glucuronidated by uridine diphosphate glucuronosyltransferases (UGT). To identify the potential impact of genetic and non-genetic factors involved in ETV metabolism, we carried out a two-step pharmacogenetics-based population pharmacokinetic study in HIV-1 infected individuals. Materials and methods: The study population included 144 individuals contributing 289 ETV plasma concentrations and four individuals contributing 23 ETV plasma concentrations collected in a rich sampling design. Genetic variants [n=125 single-nucleotide polymorphisms (SNPs)] in 34 genes with a predicted role in ETV metabolism were selected. A first step population pharmacokinetic model included non-genetic and known genetic factors (seven SNPs in CYP2C, one SNP in CYP3A5) as covariates. Post-hoc individual ETV clearance (CL) was used in a second (discovery) step, in which the effect of the remaining 98 SNPs in CYP3A, P450 cytochrome oxidoreductase (POR), nuclear receptor genes, and UGTs was investigated. Results: A one-compartment model with zero-order absorption best characterized ETV pharmacokinetics. The average ETV CL was 41 (l/h) (CV 51.1%), the volume of distribution was 1325 l, and the mean absorption time was 1.2 h. The administration of darunavir/ritonavir or tenofovir was the only non-genetic covariate influencing ETV CL significantly, resulting in a 40% [95% confidence interval (CI): 13–69%] and a 42% (95% CI: 17–68%) increase in ETV CL, respectively. Carriers of rs4244285 (CYP2C19*2) had 23% (8–38%) lower ETV CL. Co-administered antiretroviral agents and genetic factors explained 16% of the variance in ETV concentrations. None of the SNPs in the discovery step influenced ETV CL. Conclusion: ETV concentrations are highly variable, and co-administered antiretroviral agents and genetic factors explained only a modest part of the interindividual variability in ETV elimination. Opposing effects of interacting drugs effectively abrogate genetic influences on ETV CL, and vice-versa.

Investigation of Quaternary diatoms in sediment core CRP-1 from the Ross Sea, Antarctica

In the first season of drilling, the Cape Roberts Project (CRP) recovered one drillcore (CRP-l) from Roberts Ridge in western McMurdo Sound, Ross Sea, Antarctica Diatom biostratigraphy places the upper six lithostratigraphic units (Units 1.1, 2.1, 2.2, 2.3, 3.1, and 4.1) of CRP-l (0.0 to 43.15 mbsf) within the Quaternary. Both non-marine and marine Quaternary diatoms occur in variable abundance in the Quaternary interval of CRP- 1 Biostratigraphic data resolve two Quaternary time slices or events within CRP-1. Marine diatom assemblages in Units 4.1 and 3.1 represent sedimentation within the diatom Actinocyclus ingens Zone (1.35 to 0.66 Ma). Further refinement of the age of Unit 3.l places deposition in the interval 1.15 to 0.75 Ma based on the common occurrence of Thalassiosira elliptipora and correlation to the Southern Ocean acme of this taxon The absence of ActiActinocyclus ingens and the presence ot Thalassiosira antarctica in Unit 2.2 require a younger zonal assignment for this interval, within the diatom Thalassiosira lentiginosa Zone (0.66 to 0.0 Ma). A new diatom species. Rouxia leventerae, is described from marine assemblages of Units 2.2, 2.3, 3.1, and 4.l. Lithostratigraphic Unit 3.1 (33.82 to 31.89 mbsf) is a bryozoan-dominated skeletal-carbonate facies. Low abundance of Fragilariopsis curta and Fragilariopsis cylindrus within this unit combined with the relatively high abundance of species associated with open water indicates deposition in waters that remained ice free for much or all of the year Diatom assemblages suggest carbonate deposition in Unit 3.1 is linked to a significant early Pleistocene event in McMurdo Sound, when elevated surface-water temperatures inhibited the formation of sea ice.

Stable isotope record and lithic fragments distribution in sediments of MIS9-14 of the North Atlantic

Stable isotope and ice-rafted debris records from three core sites in the mid-latitude North Atlantic (IODP Site U1313, MD01-2446, MD03-2699) are combined with records of ODP Sites 1056/1058 and 980 to reconstruct hydrographic conditions during the middle Pleistocene spanning Marine Isotope Stages (MIS) 9-14 (300-540 ka). Core MD03-2699 is the first high-resolution mid-Brunhes record from the North Atlantic's eastern boundary upwelling system covering the complete MIS 11c interval and MIS 13. The array of sites reflect western and eastern basin boundary current as well as north to south transect sampling of subpolar and transitional water masses and allow the reconstruction of transport pathways in the upper limb of the North Atlantic's circulation. Hydrographic conditions in the surface and deep ocean during peak interglacial MIS 9 and 11 were similar among all the sites with relative stable conditions and confirm prolonged warmth during MIS 11c also for the mid-latitudes. Sea surface temperature (SST) reconstructions further reveal that in the mid-latitude North Atlantic MIS 11c is associated with two plateaus, the younger one of which is slightly warmer. Enhanced subsurface northward heat transport in the eastern boundary current system, especially during early MIS 11c, is denoted by the presence of tropical planktic foraminifer species and raises the question how strongly it impacted the Portuguese upwelling system. Deep water ventilation at the onset of MIS 11c significantly preceded surface water ventilation. Although MIS 13 was generally colder and more variable than the younger interglacials the surface water circulation scheme was the same. The greatest differences between the sites existed during the glacial inceptions and glacials. Then a north - south trending hydrographic front separated the nearshore and offshore waters off Portugal. While offshore waters originated from the North Atlantic Current as indicated by the similarities between the records of IODP Site U1313, ODP Site 980 and MD01-2446, nearshore waters as recorded in core MD03-2699 derived from the Azores Current and thus the subtropical gyre. Except for MIS 12, Azores Current influence seems to be related to eastern boundary system dynamics and not to changes in the Atlantic overturning circulation.

Preliminary results and documentation of sediment cores CRP2 and CRP-2A from the Ross Sea off Cape Roberts, Antarctica

The site for CRP-2, 14 km east of Cape Roberts (77.006°S; 163.719°E), was selected to overlap the early Miocene strata cored in nearby CRP-1, and to sample deeper into the east-dipping strata near the western margin ofe he Victoria Land Basin to investigate Palaeogene climatic and tectonic history. CRP-2 was cored from 5 to 57 mbsf (metres below the sea floor) (core recovery 91 %), with a deviation resulting in CRP-2A being cored at the same site. CRP-2A reached down to 624mbsf (recovery 95%), and to strata with an age of c. 33-35 Ma. Drilling took place from 16 October to 25 November 1998, on 2.0-2.2 m of sea ice and through 178 m of water. Core fractures and other physical properties, such as sonic velocity, density and magnetic susceptibility, were measured throughout the core. Down-hole logs for these and other properties were run from 63 to 167 mbsf and subsequently from 200 to 623 mbsf, although density and velocity data could be obtained only to 440 mbsf because of hole collapse. Sonic velocity averages c. 2.0 km S-1 for the upper part of the hole, but there is an sharp increase to c. 3.0 km s-1 and also a slight angular unconformity, at 306 mbsf, corresponding most likely to the early/late Oligocene boundary (c. 28-30 Ma). Velocity then increases irregularly to around 3.6 km s-1 at the bottom of the hole, which is estimated to lie 120 m above the V4/V5 boundary. The higher velocities below 306 mbsf probably reflect more extensive carbonate and common pyrite cementation, in patches, nodules, bedding-parallel masses and as vein infills. Dip of the strata also increases down-hole from 3° in the upper 300 in to over 10° at the bottom. Temperature gradient is 21° k-1. Over 2 000 fractures were logged through the hole. Borehole televiewer imagery was obtained for the interval from 200 to 440 mbsf to orient the fractures for stress field analysis. Lithostratigraphical descriptions on a scale of 1:20 are presented for the full length of the core, along with core box images, as a 200 page supplement to this issue. The hole initially passed through a layer of muddy gravel to 5.5 mbsf (Lithological Sub-Unit or LSU 1.1), and then into a Quaternary diatom-bearing clast-rich diamicton to 21 mbsf (LSU 2. l), with an interval of alternating compact diamicton and loose sand, and containing a rich Pliocene foraminiferal fauna, to 27 mbsf (LSU 2.2). The unit beneath this (LSU 3.1) has similar physical properties (sonic velocity, porosity, magnetic susceptibility) and includes diamictites of similar character to those of LSU 2.1 and 2.2, but an early Miocene (c. 19 Ma) diatom assemblage at 28 mbsf (top of LSU 3.1) shows that this sub-unit is part of the older section. The strata beneath 27 mbsf, primary target for the project, extend from early Miocene to perhaps latest Eocene age, and are largely cyclic glacimarine nearshore to offshore sediments. They are described as 41 lithological sub-units and interpreted in terms of 12 recurrent lithofacies. These are 1) mudstone, 2) inter-stratified mudstone and sandstone, 3) muddy very fine to coarse sandstone, 4) well-sorted stratified fine sandstone, 5) moderately to well-sorted, medium-grained sandstone, 6) stratified diamictite, 7) massive diamictite, 8) rhythmically inter-stratified sandstone and mudstone, 9) clast-supported conglomerate, 10) matrix-supported conglomerate, 11) mudstone breccia and 12) volcaniclastic sediment. Sequence stratigraphical analysis has identified 22 unconformity-bounded depositional sequences in pre- Pliocene strata. They typically comprise a four-part architecture involving, in ascending order, 1) a sharp-based coarse-grained unit (Facies 6,7,9 or 10), 2) a fining-upward succession of sandstones (Facies 3 and 4), 3) a mudstone interval (Facies l), in some cases coarsening upward to muddy sandstones (Facies 3), and 4) a sharp-based sandstone dominated succession (mainly Facies 4). The cyclicity recorded by the strata is interpreted in terms of a glacier ice margin retreating and advancing from land to the west, and of rises and falls in sea level. Analysis of sequence periodicity awaits afirmer chronology. However, apreliminary spectral analysis of magnetic susceptibility for a deepwater mudstone within one of the sequences (from 339 to 347 mbsf) reveals ratios between hierarchical levels that are similar to those of the three Milankovitch orbital forcing periodicities. The strata contain a wide range of fossils, the most abundant being marine diatoms. These commonly form up to 5% of the sediment, though in places the core is barren (notably between 300 and 412 mbsf). Fifty samples out of 250 reviewed were studied in detail. The assemblages define ten biostratigraphical zones, some of them based on local or as yet undescribed forms. The assemblages are neritic, and largely planktonic, suggesting that the sea floor was mostly below the photic zone throughout deposition of the corcd sequence. Calcareous nannofossils, representing incursions of ocean surface waters, are much less common (72 out of 183 samples examined) and restricted to mudstone intervals a few tens of metres thick, but are important for dating. Foraminifera are also sparse (73 out of 135 samples) and represented only by calcareous benthic species. Changing assemblages indicate a shift from inshore environments in the early Oligocenc to outer shelf in the late Oligocenc, returning to inshore in the early Miocene. Marine palynomorplis yielded large numbers of well-preserved forms from most of the 116 samples examined. The new in situ assemblagc found last year in CRP-1 is extended down into the late Oligocene and a further new assemblage is found in the early Oligoccnc. Many taxa are new, and cannot us yet contribute to an improved understanding of chronology or ecology. Marine invertebrate macrofossils, mostly molluscs and serpulid tubes, are scattered throughout the core. Preservation is good in mudstones but poor in other lithologies. Climate on land is reflected in the content of terrestrial palynomorphs, which are extremely scarce down to c. 300 mbsf. Some forms are reworked, and others represent a low growing sparse tundra with at least one species of Nothofagus. Beneath this level, a significantly greater diversity and abundance suggests a milder climate and a low diversity woody vegetation in the early Oligocene, but still far short of the richness found in known Eocene strata of the region. Sedimentary facies in the oldest strata also suggest a milder climate in the oldest strata cored, with indications of substantial glacial melt-water discharges, but are typical of a coldcr climate in late Oligocene and early Miocene times. Clast analyses from diamictites reveal weak to random fabrics, suggesting either lack of ice-contact deposition or post-depositional modification, but periods when ice grounded at the drill site are inferred from thin zones of in-situ brecciated rock and soft-sediment folding. These are more common above c. 300 mbsf, perhaps reflecting more extensive glacial advances during deposition of those strata. Erosion of the adjacent Transantarctic Mountains through Jurassic basalt and dolerite-intruded Beacon strata into basement rocks beneath is recorded by petrographical studies of clast and sand grain assemblages. Core below 310 mbsf contains a dominance of fine-grained Jurassic dolerite and basalt fragments along with Beacon-derived coal debris and rounded quartz grains, whereas the strata above this level have a much higher proportion of basement derived granitoids, implying that the large areas of the adjacent mountains had been eroded to basement by the end of the early Oligocene. There is little indication of rift-related volcanism below 310 mbsf. Above this, however, basaltic and trachytic tephras are common, especially from 280 to 200 mbsf, from 150 to 46 mbsf, and in Pliocene LSU 2.2 from 21 to 27 mbsf. The largest volcanic eruptions generated layers of coarse (up to 1 cm) trachytic pumice lapilli between 97 and 114 mbsf. The thickest of these (1.2 m at 112 mbsf) may have produced an eruptive column extending tens of km into the stratosphere. A source within a few tens of km of the drill site is considered most likely. Present age estimates for the pre-Pliocene sequence are based mainly on biostratigraphy (using mainly marine diatoms and to a lesser extent calcareous nannofossils), with the age of the tephra from 112 to 114 mbsf (21.44k0.05 Ma from 84 crystals by Ar-Ar) as a key reference point. Although there are varied and well-preserved microfossil assemblages through most of the sequence (notably of diatoms and marine palynomorphs), they comprise largely taxa either known only locally or as yet undescribed. In addition, sequence stratigraphical analysis and features in the core itself indicate numerous disconformities. The present estimate from diatom assemblages is that the interval from 27 to 130 mbsf is early Miocene in age (c. 19 to 23.5 Ma), consistent with the Ar-Ar age from 112 to 114 mbsf. Diatom assemblages also indicate that the late Oligocene epoch extends from c. 130 to 307 mbsf, which is supported by late Oligocene nannofossils from 130 to 185 mbsf. Strata from 307 to 412 mbsf have no age-diagnostic assemblages, but below this early Oligocene diatoms and nannofossils have been recovered. A nannoflora at the bottom of the hole is consistent with an earliest Oligocene or latest Eocene age. Magnetostratigraphical studies based on about 1000 samples, 700 of which have so far undergone demagnetisation treatment, have provided a polarity stratigraphy of 12 pre-Pliocene magnetozones. Samples above 270 mbsf are of consistently high quality. Below this, magnetic behaviour is more variable. A preliminary age-depth plot using the Magnetic Polarity Time Scale (MPTS) and constrained by biostratigraphical data suggests that episodes of relatively rapid sedimentation took place at CRP-2 during Oligocene times (c. 100 m/My), but that more than half of the record was lost in a few major and many minor disconformities. Age estimates from Sr isotopes in shell debris and further tephra dating are expected to lead to a better comparison with the MPTS. CRP-2/2A has recorded a history of subsidence of the Victoria Land Basin margin that is similar to that found in CIROS-170 km to the south, reflecting stability in both basin and the adjacent mountains in late Cenozoic times, but with slow net accumulation in the middle Cenozoic. The climatic indicators from both drill holes show a similar correspondence, indicating polar conditions for the Quaternary but with sub-polar conditions in the early Miocene-late Oligocene and indications of warmer conditions still in the early Oligocene. Correlation between the CRP-2A core and seismic records shows that seismic units V3 and V4, both widespread in the Victoria Land Basin, represent a period of fluctuating ice margins and glacimarine sedimentation. The next drill hole, CRP-3, is expected to core deep into V5 and extend this record of climate and tectonics still further back in time.

Biostratigraphy of sediment core CRP-3 from the Ross Sea, Antarctica

Early Oligocene siliceous microfossils were recovered in the upper c. 193 m of the CRP-3 drillcore. Although abundance and preservation are highly variable through this section, approximately 130 siliceous microfossil taxa were identified, including diatoms, silicoflagellates, ebridians, chrysophycean cysts, and endoskeletal dinoflagellates. Well-preserved and abundant assemblages characterize samples in the upper c. 70 m and indicate deposition in a coastal setting with water depths between 50 and 200 m. Abundance fluctuations over narrow intervals in the upper c. 70 mbsf are interpreted to reflect environmental changes that were either conducive or deleterious to growth and preservation of siliceous microfossils. Only poorly-preserved (dissolved, replaced, and/or fragmented) siliceous microfossils are present from c. 70 to 193 mbsf. Diatom biostratigraphy indicates that the CRP-3 section down to c. 193 mbsf is early Oligocene in age. The lack of significant changes in composition of the siliceous microfossil assemblage suggests that no major hiatuses are present in this interval. The first occurrence (FO) of Cavitatus jouseanus at 48.44 mbsf marks the base of the Cavitatus jouseanus Zone. This datum is inferred to be near the base of Subchron C12n at c. 30.9 Ma. The FO of Rhizosolenia antarctica at 68.60 mbsf marks the base of the Rhizosolenia antarctica Zone. The FO of this taxon is correlated in deep-sea sections to Chron C13 (33.1 to 33.6 Ma). However, the lower range of R. antarctica is interpreted as incomplete in the CRP-3 drillcore, as it is truncated at an underlying interval of poor preservation: therefore, an age of c. 33.1 to 30.9 Ma is inferred for interval between c. 70 and 50 mbsf. The absence of Hemiaulus caracteristicus from diatom-bearing interval of CRP-3 further indicates an age younger than c. 33 Ma (Subchron C13n) for strata above c. 193 mbsf. Siliceous microfossil assemblages in CRP-3 are significantly different from the late Eocene assemblages reported CIROS-1 drillcore. The absence of H. caracteristicus, Stephanopyxis splendidus, and Pterotheca danica, and the ebridians Ebriopsis crenulata, Parebriopsis fallax, and Pseudoammodochium dictyoides in CRP-3 indicates that the upper 200 m of the CRP-3 drillcore is equivalent to part of the stratigraphic interval missing within the unconformity at c. 366 mbsf in CIROS-1.