Classification of CM chondrite breccias—Implications for the evaluation of samples from the OSIRIS‐REx and Hayabusa 2 missions

CM chondrites are complex impact (mostly regolith) breccias, in which lithic clasts show various degrees of aqueous alteration. Here, we investigated the degree of alteration of individual clasts within 19 different CM chondrites and CM‐like clasts in three achondrites by chemical analysis of the tochilinite‐cronstedtite‐intergrowths (TCIs; formerly named “poorly characterized phases”). To identify TCIs in various chondritic lithologies, we used backscattered electron (BSE) overview images of polished thin sections, after which appropriate samples underwent electron microprobe measurements. Thus, 75 lithic clasts were classified. In general, the excellent work and specific criteria of Rubin et al. (2007) were used and considered to classify CM breccias in a similar way as ordinary chondrite breccias (e.g., CM2.2‐2.7). In BSE images, TCIs in strongly altered fragments in CM chondrites (CM2.0‐CM2.2) appear dark grayish and show a low contrast to the surrounding material (typically clastic matrix), and can be distinguished from TCIs in moderately (CM2.4‐CM2.6) or less altered fragments (CM2.7‐CM2.9); the latter are bright and have high contrast to the surroundings. We found that an accurate subclassification can be obtained by considering only the “FeO”/SiO₂ ratio of the TCI chemistry. One could also consider the TCIs’ S/SiO₂ ratio and the metal abundance, but these were not used for classification due to several disadvantages. Most of the CM chondrites are finds that have suffered terrestrial weathering in hot and cold deserts. Thus, the observed abundance of metal is susceptible to weathering and may not be a reliable indicator of subtype classification. This study proposes an extended classification scheme based on Rubin’s scale from subtypes CM2.0‐CM2.9 that takes the brecciation into account and includes the minimum to maximum degree of alteration of individual clasts. The range of aqueous alteration in CM chondrites and small spatial scale of mixing of clasts with different alteration histories will be important for interpreting returned samples from the OSIRIS‐REx and Hayabusa 2 missions in the future.     

Combining pixel and object based image analysis of ultra-high resolution multibeam bathymetry and backscatter for habitat mapping in shallow marine waters

Habitat mapping data are increasingly being recognised for their importance in underpinning marine spatial planning. The ability to collect ultra-high resolution (cm) multibeam echosounder (MBES) data in shallow waters has facilitated understanding of the fine-scale distribution of benthic habitats in these areas that are often prone to human disturbance. Developing quantitative and objective approaches to integrate MBES data with ground observations for predictive modelling is essential for ensuring repeatability and providing confidence measures for habitat mapping products. Whilst supervised classification approaches are becoming more common, users are often faced with a decision whether to implement a pixel based (PB) or an object based (OB) image analysis approach, with often limited understanding of the potential influence of that decision on final map products and relative importance of data inputs to patterns observed. In this study, we apply an ensemble learning approach capable of integrating PB and OB Image Analysis from ultra-high resolution MBES bathymetry and backscatter data for mapping benthic habitats in Refuge Cove, a temperate coastal embayment in south-east Australia. We demonstrate the relative importance of PB and OB seafloor derivatives for the five broad benthic habitats that dominate the site. We found that OB and PB approaches performed well with differences in classification accuracy but not discernible statistically. However, a model incorporating elements of both approaches proved to be significantly more accurate than OB or PB methods alone and demonstrate the benefits of using MBES bathymetry and backscatter combined for class discrimination.     

The connection between Labrador Sea buoyancy loss, deep western boundary current strength, and Gulf Stream path in an ocean circulation model

The sensitivity of the North Atlantic gyre circulation to high latitude buoyancy forcing is explored in a global, non-eddy resolving ocean general circulation model. Increased buoyancy forcing strengthens the deep western boundary current, the northern recirculation gyre, and the North Atlantic Current, which leads to a more realistic Gulf Stream path. High latitude density fluxes and surface water mass transformation are strongly dependent on the choice of sea ice and salinity restoring boundary conditions. Coupling the ocean model to a prognostic sea ice model results in much greater buoyancy loss in the Labrador Sea compared to simulations in which the ocean is forced by prescribed sea ice boundary conditions. A comparison of bulk flux forced hindcast simulations which differ only in their sea ice and salinity restoring forcings reveals the effects of a mixed thermohaline boundary condition transport feedback whereby small, positive temperature and salinity anomalies in subpolar regions are amplified when the gyre spins up as a result of increased buoyancy loss and convection. The primary buoyancy flux effects of the sea ice which cause the simulations to diverge are ice melt, which is less physical in the diagnostic sea ice model, and insulation of the ocean, which is less physical with the prognostic sea ice model. Increased salinity restoring ensures a more realistic net winter buoyancy loss in the Labrador Sea, but it is found that improvements in the Gulf Stream simulation can only be achieved with the excessive buoyancy loss associated with weak salinity restoring.     

Singular value decomposition and cluster analysis as regression diagnostics tools for geodetic applications

It is well known that high-leverage observations significantly affect the estimation of parameters. In geodetic literature, mainly redundancy numbers are used for the detection of single high-leverage observations or of single redundant observations. In this paper a further objective method for the detection of groups of important and less important (and thus redundant) observations is developed. In addition, the parameters which are predominantly affected by these groups of observations are identified. This method thus complements other diagnostics tools, such as, e.g., multiple row diagnostics methods as described in statistical literature (see, e.g., Belsley et al. in Regression diagnostics: identifying influential data and sources of collinearity. Wiley, New York, 1980). The method proposed in this paper is based on geometric aspects of adjustment theory and uses the singular value decomposition of the design matrix of an adjustment problem together with cluster analysis methods for regression diagnostics. It can be applied to any geodetic adjustment problem and can be used for the detection of (groups of) observations that significantly affect the estimated parameters or that are of negligible impact. One of the advantages of the proposed method is the improvement of the reliability of observation plans and thus the reduction of the impact of individual observations (and outliers) on the estimated parameters. This is of particular importance for the very long baseline interferometry technique which serves as an application example of the regression diagnostics tool.     

Characterization of ultrafine aerosol particles in Mt. Etna emissions

Particulate emissions from Mt. Etna in the fine-size range below 100 nm were studied in June and September 1989. The aerosol particles were characterized by size, concentration and photoelectric activity. These quantities are sensitive to the physical and chemical properties of the magma. Concentrations varied from 10⁴ to 10⁷ cm^-3. The size distributions peak below 20 nm (radius) and are very narrow. The particles are generated mainly by nucleation and condensation of magmatic volatiles in a strong temperature gradient. The photoelectric activity of these particles can indicate high magma levels and increased exsolution of volatiles. It is therefore related to the observed activity of the respective crater and may be helpful as a prediction tool when used in conjunction with other volcano-monitoring techniques.     

Correlation between regional and global seismic activity

The regional seismic activity in Sweden is significantly correlated to the global seismic activity. This result suggests that the regional Swedish activity is linked to the global activity by intraplate tectonics, rather than being caused exclusively by regional effects, such as the land uplift after the last glaciation.     

Earthquake data analysis: An example from Sweden

Methods and results of earthquake data analysis are illustrated by the seismicity of Sweden in the period 1951–1976, summarizing a comprehensive research project for the past five years. Starting from the regional structural properties, methods are given for the evaluation of instrumental and macroseismic observations, including error analysis and source identification. The resulting earthquake distribution is investigated in relation to space (various types of seismicity mapping), time (periodicities) and magnitude (frequency-magnitude and energy-magnitude relations). The tectonophysical explanations of the earthquake origin are still subject to considerable uncertainty and can be treated only by way of discussions. The engineering applications, stimulated recently by the nuclear power industry, can be placed on much more reliable ground. They include calculations of seismic risk and fracture risk as well as investigations of certain phenomena of practical consequences, especially rockbursts.This review is primarily addressed to those who have use for seismological data in their work, especially geologists, tectonophysicists and engineers. It is important for them to realize how the seismological data are obtained, their reliability and their range of application. Sweden is here mainly used as an example, but the developed methods have a general applicability to any earthquake region, with proper modifications.In short, our purpose is two-fold: (1) to summarize recent research on the seismicity of Sweden: (2) to introduce all users of seismological data into methods and results.     

Chironomid-based paleosalinity records in southern British Columbia, Canada: a comparison of transfer functions

Chironomid remains from Big Lake, British Columbia were analysed and paleosalinities were estimated using a pre-existing transfer function and several developed using new regression methods. A two component partial-least-squares model (PLS-2) had the highest coefficient of determination (R² _(Jackknifed) = 0.75) and lowest root-mean-squared error-of-prediction (RMSEP). As compared to the pre-existing model, it was also less sensitive to the influence of rare taxa. Nevertheless, the marginally larger R² _(Jackknifed) and lower RMSEP do not clearly identify a single best model. The models were applied to Big, Mahoney and Kilpoola lakes, revealing the sensitivity of paleosalinity inferences to model selection. A synopsis of chironomid-based paleosalinities in British Columbia and their correspondence with other paleoclimatic data are presented and discussed.     

Reconstructing spatio-temporal patterns of debris-flow activity using dendrogeomorphological methods

Debris flows are a major threat in many parts of the Alps, where they repeatedly cause severe damage to infrastructure and transportation corridors or even loss of life. Nonetheless, the spatial behavior of past debris-flow activity and the analysis of areas affected during particular events have been widely neglected in reconstructions so far. It was therefore the purpose of this study to reconstruct spatio-temporal patterns of past debris flows on a forested cone in the Swiss Alps (Bruchji torrent, Blatten, Valais). The analysis of past events was based on a detailed geomorphic map (1:1000) of all forms related to debris flows as well as on tree-ring series from 401 heavily affected trees (Larix decidua Mill. and Picea abies (L.) Karst.) growing in or next to deposits. The samples were analyzed and growth disturbances related to debris-flow activity assessed, such as tangential rows of traumatic resin ducts, the onset of reaction wood or abrupt growth suppression or release.In total, 960 growth disturbances were identified in the samples, belonging to 40 different event years between A.D. 1867 and 2005. In addition, the coupling of tree-ring data with the geomorphic map allowed reconstruction of eleven formerly active channels and spatial representation of individual events. Based on our results we believe that before 1935, debris flows preferentially used those channels located in the western part of the cone, whereas the eastern part of the cone remained widely unaffected. The spatial representation of the 40 events also allowed identification of five different spatial patterns for debris flows at the study site.     

Dual models of pore spaces

We present a model for pore spaces that consists of two parts related by duality: (1) a decomposition of an open polyhedral pore space into open contractible pore bodies separated by relatively open interfaces and (2) a pore network that is homotopy equivalent to the pore space. The dual model is unique and free of parameters, but it relies on regularity conditions for the pore space. We show how to approximate any pore space by the interior of a polyhedral complex such that the regularity conditions are fulfilled. Thus, we are able to calculate the dual model from synthetic porous media and images of real porous media. The pore bodies are unions of relatively open Delaunay cells with respect to the corners of the pore boundary, and the pore network consists of certain at most two-dimensional (2D) Voronoi cells with respect to the corners of the pore boundary. The pore network describes the neighborhood relations between the pore bodies. In particular, any relatively open 2D Delaunay face f separating two pore bodies has a unique (relatively open) dual network edge. In our model, f is a pore throat only if it is hit by its dual network edge. Thus, as opposed to widespread intuition, any pore throat is convex, and adjacent pore bodies are not necessarily separated by pore throats. Due to the duality between the pore network and the decomposition of the pore space into pore bodies it is straightforward to store the geometrical properties of the pore bodies [pore throats] as attributes of the dual network vertices [edges]. Such an attributed network is used to perform 2D drainage simulations. The results agree very well with those from a pore-morphology based modeling approach performed directly on the digital image of a porous medium. Contractibility of the pore bodies and homotopy equivalence of the pore space and the pore network is proven using discrete Morse theory and the nerve theorem from combinatorial topology.