A new technique to automatically quantify microstructures of fine grained carbonate mylonites: two-step etching combined with SEM imaging and image analysis

A two-step etching technique for fine-grained calcite mylonites using 0.37% hydrochloric and 0.1% acetic acid produces a topographic relief which reflects the grain boundary geometry. With this technique, calcite grain boundaries become more intensely dissolved than their grain interiors but second phase minerals like dolomite, quartz, feldspars, apatite, hematite and pyrite are not affected by the acid and therefore form topographic peaks. Based on digital backscatter electron images and element distribution maps acquired on a scanning electron microscope, the geometry of calcite and the second phase minerals can be automatically quantified using image analysis software. For research on fine-grained carbonate rocks (e.g. dolomite calcite mixtures), this low-cost approach is an attractive alternative to the generation of manual grain boundary maps based on photographs from ultra-thin sections or orientation contrast images.

Extensional orogenic collapse captured by strike-slip tectonics:

The late Paleozoic collision between Gondwana and Laurussia resulted in the polyphase deformation and magmatism that characterizes the Iberian Massif of the Variscan orogen. In the Central Iberian Zone, initial con- tinental thickening (D1; folding and thrusting) was followed by extensional orogenic collapse (D2) responsible for the exhumation of high-grade rocks coeval to the emplacement of granitoids. This study presents a tectonometamorphic analysis of the Trancoso-Pinhel region (Central Iberian Zone) to ex- plain the processes in place during the transition froman extension-dominated state (D2) to a compression-dom- inated one (D3).Wereveal the existence of low-dipping D2 extensional structures later affected by several pulses of subhorizontal shortening, each of them typified by upright folds and strike-slip shearing (D3, D4 and D5, as identified by superimposition of structures). The D2 Pinhel extensional shear zone separates a low-grade domain from an underlying high-grade domain, and it contributed to the thermal reequilibration of the orogen by facil- itating heat advection from lower parts of the crust, crustal thinning, decompression melting, and magma intru- sion. Progressive lessening of the gravitational disequilibrium carried out by this D2 shear zone led to a switch from subhorizontal extension to compression and the eventual cessation and capture of the Pinhel shear zone by strike-slip tectonics during renewed crustal shortening. High-grade domains of the Pinhel shear zone were folded together with low-grade domains to define the current upright folded structure of the Trancoso-Pinhel re- gion, the D3 Tamames-Marofa-Sátão synform. Newdating of syn-orogenic granitoids (SHRIMP U\\Pb zircon dat- ing) intruding the Pinhel shear zone, together with the already published ages of early extensional fabrics constrain the functioning of this shear zone to ca. 331–311 Ma, with maximum tectonomagmatic activity at ca. 321–317 Ma. The capture and apparent cessation of movement of the Pinhel shear zone occurred at ca. 317– 311 Ma.

Basin analysis and economic geology of the Northern Mount Isa Basin

The Mount Isa Basin is a new concept used to describe the area of Palaeo- to Mesoproterozoic rocks south of the Murphy Inlier and inappropriately described presently as the Mount Isa Inlier. The new basin concept presented in this thesis allows for the characterisation of basin-wide structural deformation, correlation of mineralisation with particular lithostratigraphic and seismic stratigraphic packages, and the recognition of areas with petroleum exploration potential. The northern depositional margin of the Mount Isa Basin is the metamorphic, intrusive and volcanic complex here referred to as the Murphy Inlier (not the "Murphy Tectonic Ridge"). The eastern, southern and western boundaries of the basin are obscured by younger basins (Carpentaria, Eromanga and Georgina Basins). The Murphy Inlier rocks comprise the seismic basement to the Mount Isa Basin sequence. Evidence for the continuity of the Mount Isa Basin with the McArthur Basin to the northwest and the Willyama Block (Basin) at Broken Hill to the south is presented. These areas combined with several other areas of similar age are believed to have comprised the Carpentarian Superbasin (new term). The application of seismic exploration within Authority to Prospect (ATP) 423P at the northern margin of the basin was critical to the recognition and definition of the Mount Isa Basin. The Mount Isa Basin is structurally analogous to the Palaeozoic Arkoma Basin of Illinois and Arkansas in southern USA but, as with all basins it contains unique characteristics, a function of its individual development history. The Mount Isa Basin evolved in a manner similar to many well described, Phanerozoic plate tectonic driven basins. A full Wilson Cycle is recognised and a plate tectonic model proposed. The northern Mount Isa Basin is defined as the Proterozoic basin area northwest of the Mount Gordon Fault. Deposition in the northern Mount Isa Basin began with a rift sequence of volcaniclastic sediments followed by a passive margin drift phase comprising mostly carbonate rocks. Following the rift and drift phases, major north-south compression produced east-west thrusting in the south of the basin inverting the older sequences. This compression produced an asymmetric epi- or intra-cratonic clastic dominated peripheral foreland basin provenanced in the south and thinning markedly to a stable platform area (the Murphy Inlier) in the north. The fmal major deformation comprised east-west compression producing north-south aligned faults that are particularly prominent at Mount Isa. Potential field studies of the northern Mount Isa Basin, principally using magnetic data (and to a lesser extent gravity data, satellite images and aerial photographs) exhibit remarkable correlation with the reflection seismic data. The potential field data contributed significantly to the unravelling of the northern Mount Isa Basin architecture and deformation. Structurally, the Mount Isa Basin consists of three distinct regions. From the north to the south they are the Bowthorn Block, the Riversleigh Fold Zone and the Cloncurry Orogen (new names). The Bowthom Block, which is located between the Elizabeth Creek Thrust Zone and the Murphy Inlier, consists of an asymmetric wedge of volcanic, carbonate and clastic rocks. It ranges from over 10 000 m stratigraphic thickness in the south to less than 2000 min the north. The Bowthorn Block is relatively undeformed: however, it contains a series of reverse faults trending east-west that are interpreted from seismic data to be down-to-the-north normal faults that have been reactivated as thrusts. The Riversleigh Fold Zone is a folded and faulted region south of the Bowthorn Block, comprising much of the area formerly referred to as the Lawn Hill Platform. The Cloncurry Orogen consists of the area and sequences equivalent to the former Mount Isa Orogen. The name Cloncurry Orogen clearly distinguishes this area from the wider concept of the Mount Isa Basin. The South Nicholson Group and its probable correlatives, the Pilpah Sandstone and Quamby Conglomerate, comprise a later phase of now largely eroded deposits within the Mount Isa Basin. The name South Nicholson Basin is now outmoded as this terminology only applied to the South Nicholson Group unlike the original broader definition in Brown et al. (1968). Cored slimhole stratigraphic and mineral wells drilled by Amoco, Esso, Elf Aquitaine and Carpentaria Exploration prior to 1986, penetrated much of the stratigraphy and intersected both minor oil and gas shows plus excellent potential source rocks. The raw data were reinterpreted and augmented with seismic stratigraphy and source rock data from resampled mineral and petroleum stratigraphic exploration wells for this study. Since 1986, Comalco Aluminium Limited, as operator of a joint venture with Monument Resources Australia Limited and Bridge Oil Limited, recorded approximately 1000 km of reflection seismic data within the basin and drilled one conventional stratigraphic petroleum well, Beamesbrook-1. This work was the first reflection seismic and first conventional petroleum test of the northern Mount Isa Basin. When incorporated into the newly developed foreland basin and maturity models, a grass roots petroleum exploration play was recognised and this led to the present thesis. The Mount Isa Basin was seen to contain excellent source rocks coupled with potential reservoirs and all of the other essential aspects of a conventional petroleum exploration play. This play, although high risk, was commensurate with the enormous and totally untested petroleum potential of the basin. The basin was assessed for hydrocarbons in 1992 with three conventional exploration wells, Desert Creek-1, Argyle Creek-1 and Egilabria-1. These wells also tested and confrrmed the proposed basin model. No commercially viable oil or gas was encountered although evidence of its former existence was found. In addition to the petroleum exploration, indeed as a consequence of it, the association of the extensive base metal and other mineralisation in the Mount Isa Basin with hydrocarbons could not be overlooked. A comprehensive analysis of the available data suggests a link between the migration and possible generation or destruction of hydrocarbons and metal bearing fluids. Consequently, base metal exploration based on hydrocarbon exploration concepts is probably. the most effective technique in such basins. The metal-hydrocarbon-sedimentary basin-plate tectonic association (analogous to Phanerozoic models) is a compelling outcome of this work on the Palaeo- to Mesoproterozoic Mount lsa Basin. Petroleum within the Bowthom Block was apparently destroyed by hot brines that produced many ore deposits elsewhere in the basin.

Structural and textural variations in carbonaceous chondrite matrices

Recent studies of C2 carbonaceous chondrite matrices using high resolu­tion transmission electron microscopy (HRTEM)have shown that structural details of the matrix minerals can be imaged [1-4]. The Murchison and Mighei matrices contain minerals having ordered and disordered mixed-layer structures [1,3,4] in addition to chrysotile- and lizardite-type structures [2].

Topographical, structural and geophysical characterization of fracture zones: implications for groundwater flow and vulnerability

The main objective of this study is to evaluate selected geophysical, structural and topographic methods on regional, local, and tunnel and borehole scales, as indicators of the properties of fracture zones or fractures relevant to groundwater flow. Such information serves, for example, groundwater exploration and prediction of the risk of groundwater inflow in underground construction. This study aims to address how the features detected by these methods link to groundwater flow in qualitative and semi-quantitative terms and how well the methods reveal properties of fracturing affecting groundwater flow in the studied sites. The investigated areas are: (1) the Päijänne Tunnel for water-conveyance whose study serves as a verification of structures identified on regional and local scales; (2) the Oitti fuel spill site, to telescope across scales and compare geometries of structural assessment; and (3) Leppävirta, where fracturing and hydrogeological environment have been studied on the scale of a drilled well. The methods applied in this study include: the interpretation of lineaments from topographic data and their comparison with aeromagnetic data; the analysis of geological structures mapped in the Päijänne Tunnel; borehole video surveying; groundwater inflow measurements; groundwater level observations; and information on the tunnel s deterioration as demonstrated by block falls. The study combined geological and geotechnical information on relevant factors governing groundwater inflow into a tunnel and indicators of fracturing, as well as environmental datasets as overlays for spatial analysis using GIS. Geophysical borehole logging and fluid logging were used in Leppävirta to compare the responses of different methods to fracturing and other geological features on the scale of a drilled well. Results from some of the geophysical measurements of boreholes were affected by the large diameter (gamma radiation) or uneven surface (caliper) of these structures. However, different anomalies indicating more fractured upper part of the bedrock traversed by well HN4 in Leppävirta suggest that several methods can be used for detecting fracturing. Fracture trends appear to align similarly on different scales in the zone of the Päijänne Tunnel. For example, similarities of patterns were found between the regional magnetic trends, correlating with orientations of topographic lineaments interpreted as expressions of fracture zones. The same structural orientations as those of the larger structures on local or regional scales were observed in the tunnel, even though a match could not be made in every case. The size and orientation of the observation space (patch of terrain at the surface, tunnel section, or borehole), the characterization method, with its typical sensitivity, and the characteristics of the location, influence the identification of the fracture pattern. Through due consideration of the influence of the sampling geometry and by utilizing complementary fracture characterization methods in tandem, some of the complexities of the relationship between fracturing and groundwater flow can be addressed. The flow connections demonstrated by the response of the groundwater level in monitoring wells to pressure decrease in the tunnel and the transport of MTBE through fractures in bedrock in Oitti, highlight the importance of protecting the tunnel water from a risk of contamination. In general, the largest values of drawdown occurred in monitoring wells closest to the tunnel and/or close to the topographically interpreted fracture zones. It seems that, to some degree, the rate of inflow shows a positive correlation with the level of reinforcement, as both are connected with the fracturing in the bedrock. The following geological features increased the vulnerability of tunnel sections to pollution, especially when several factors affected the same locations: (1) fractured bedrock, particularly with associated groundwater inflow; (2) thin or permeable overburden above fractured rock; (3) a hydraulically conductive layer underneath the surface soil; and (4) a relatively thin bedrock roof above the tunnel. The observed anisotropy of the geological media should ideally be taken into account in the assessment of vulnerability of tunnel sections and eventually for directing protective measures.

Earthquake geology of Myanmar

This thesis describes the active structures of Myanmar and its surrounding regions, and the earthquake geology of the major active structures. Such investigation is needed urgently for this rapidly developing country that has suffered from destructive earthquakes in its long history. To archive a better understanding of the regional active tectonics and the seismic potential in the future, we utilized a global digital elevation model and optical satellite imagery to describe geomorphologic evidence for the principal neotectonic features of the western half of the Southeast Asia mainland. Our investigation shows three distinct active structural systems that accommodate the oblique convergence between the Indian plate and Southeast Asia and the extrusion of Asian territory around the eastern syntaxis of the Himalayan mountain range. Each of these active deformation belts can be further separated into several neotectonic domains, in which structures show distinctive active behaviors from one to another.

In order to better understand the behaviors of active structures, we focused on the active characteristics of the right-lateral Sagaing fault and the oblique subducting northern Sunda megathrust in the second part of this thesis. The detailed geomorphic investigations along these two major plate-interface faults revealed the recent slip behavior of these structures, and plausible recurrence intervals of major seismic events. We also documented the ground deformation of the 2011 Tarlay earthquake in remote eastern Myanmar from remote sensing datasets and post-earthquake field investigations. The field observation and the remote sensing measurements of surface ruptures of the Tarlay earthquake are the first study of this kind in the Myanmar region.

Geology of the Pacoima Area, Los Angeles County, California

The Pacoima area is located on an isolated hill in the northeast section of the San Fernando, the northeast portion of the Pacoima Quadrangle, Los Angeles County, California. Within it are exposed more than 2300 feet of Tertiary rocks, which comprise three units of Middle Miocene (?) age, and approximately 950 feet of Jurassic (?) granite basement. The formations are characterized by their mode of occurrence, marine and terrestial origin, diverse lithology, and structural features.

The basement complex is composed of intrusive granite, small masses of granodiorite and a granodiorite gneiss with the development of schistosity in sections. During the long period of erosion of the metamorphics, the granitic rocks were exposed and may have provided clastic constituents for the overlying formations.

As a result of rapid sedimentation in a transitional environment, the Middle Miocene Twin Peaks formation was laid down unconformably on the granite. This formation is essentially a large thinning bed of gray to buff pebble and cobble conglomerate grading to coarse yellow sandstone. The contact of conglomerate and granite is characterized by its faulted and depositional nature.

Beds of extrusive andesite, basalt porphyry, compact vesicular amygdaloidal basalts, andesite breccia, interbedded feldspathic sands and clays of terrestial origin, and mudflow breccia comprise the Pacoima formation which overlies the Twin Peaks formation unconformably. A transgressing shallow sea accompanied settling of the region and initiated deposition of fine clastic sediments.

The marine Topanga (?) formation is composed of brown to gray coarse sandstone grading into interbedded buff sandstones and gray shales. Intrusions of rhyolitedacite and ash beds mark continued but sporatic volcanism during this period.

The area mapped represents an arch in the Tertiary sediments. Forces that produced the uplift of the granite structural high created stresses that were relieved by jointing and faulting. Vertical and horizontal movement along these faults has displaced beds, offset contacts and complicated their structure. Uplift and erosion have exposed the present sequence of beds which dip gently to the northeast. The isolated hill is believed to be in an early stage of maturity.

Geology of the Pacoima Hills

The Pacoima Hills lie between Foothill Boulevard and the San Fernando Road, three miles southeast of San Fernando, California. In this area are exposed Jurassic(?) granodiorite intruded in older gneiss, and a mid Miocene Topango (?) sedimentary section lying in both fault and sedimentary contact with the intrusive complex. Two distinct lava flows and a small laccolith of andesite occur within the Topango (?) formation. The principal structural feature is an anticline plunging steeply northward. An upward acting force is postulated to have produced this anticline; upon cessation of the force, normal faulting occurred with consequent down-dropping of north-south blocks.

Structural analysis of Turtle mountain (Alberta) using digital elevation model: Toward a progressive failure

In 1903, the eastern slope of Turtle Mountain (Alberta) was affected by a 30 M m3-rockslide named Frank Slide that resulted in more than 70 casualties. Assuming that the main discontinuity sets, including bedding, control part of the slope morphology, the structural features of Turtle Mountain were investigated using a digital elevation model (DEM). Using new landscape analysis techniques, we have identified three main joint and fault sets. These results are in agreement with those sets identified through field observations. Landscape analysis techniques, using a DEM, confirm and refine the most recent geology model of the Frank Slide. The rockslide was initiated along bedding and a fault at the base of the slope and propagated up slope by a regressive process following a surface composed of pre-existing discontinuities. The DEM analysis also permits the identification of important geological structures along the 1903 slide scar. Based on the so called Sloping Local Base Level (SLBL) an estimation was made of the present unstable volumes in the main scar delimited by the cracks, and around the south area of the scar (South Peak). The SLBL is a method permitting a geometric interpretation of the failure surface based on a DEM. Finally we propose a failure mechanism permitting the progressive failure of the rock mass that considers gentle dipping wedges (30°). The prisms or wedges defined by two discontinuity sets permit the creation of a failure surface by progressive failure. Such structures are more commonly observed in recent rockslides. This method is efficient and is recommended as a preliminary analysis prior to field investigation.

Structural studies in the Southern Province, South of Sudbury, Ontario

Rocks correlated with the Hough Lake and Quirke Lake Groups of the Huronian Supergroup form part of a northeasterly trending corridor that separates 1750 Ma granitic intrusive rocks of the Chief Lake batholith from the 1850 Ma mafic intrusive rocks of the Sudbury Igneous Complex. This corridor is dissected by two major structural features; the Murray Fault Zone (MFZ) and the Long Lake Fault (LLF). Detailed structural mapping and microstructural analysis indicates that the LLF, which has juxtaposed Huronian rocks of different deformation style and metamorphism grade, was a more significant plane of dislocation than the MFZ. The sense of displacement along the LLF is high angle reverse in which rocks to the southeast have been raised relative to those in the northwest. South of the LLF Huronian rocks underwent ductile defonnation at amphibolite facies conditions. The strain was constrictional, defined by a triaxial strain ellipsoid in which X > Y > z. Calculations of a regional k value were approximately 1.3. Penetrative ductile defonnation resulted in the development of a preferred crystallographic orientation in quartz as well as the elongation of quartz grains to fonn a regional southeast-northwest trending, subvertical lineation. Similar lithologies north of the LLF underwent dominantly brittle deformation under greenschist facies conditions. Deformation north of the LLF is characterized by the thrusting of structural blocks to form angular discordances in bedding orientation which were previously interpreted as folds. Ductile deformation occurred between 1750 and 1238 Ma and is correlated with a regional period of south over north reverse faulting that effected much of the southern Sudbury region. Post dating the reverse faulting event was a period of sedimentation as a conglomerate unit was deposited on vertically bedded Huronian rocks. Rocks in the study area were intruded by both mafic and felsic dykes. The 1238 Ma mafic dykes appear to have been offset during a period of dextral strike slip displacement along the major fault'). Indirect evidence indicates that this event occurred after the thrusting at 950 to 1100 Ma associated with the Grenvillian Orogeny.

Geology of the Rankin Inlet area, Northwest Territories /

The Rankin Inlet area, on the west shore of Hudson Bay in the Northwest Territories, is in the Churchill Structural Province. Metamorphosed volcanic and sedimentary rocks, previously mapped as Archean and part of the Kaminak Group, underlie most of the area. The Rankin Inlet Group consists of greywacke, with minor conglomeratic greywacke, quartzite and dolomite, overlain by massive and pillowed basaltic flows. Gabbro sills intrude the sediments near the base of the volcanic sequence and three serpentinite sills outcrop at the base of the volcanic sequence. The sediments are in fault-contact with quartz monzonite to the south and were intruded by granitic rocks to the northwest. Two periods of folding were defined by the mapping. The first generation folds are recumbent isoclinal folds, with northwest-trending and northeast-dipping axial planes, formed through gravitational sliding. The second generation folds are symmetrically disposed about the axis of the granitic intrusion and have east-southeast trending and nearly vertical axial planes. Whole-rock analysis of 64 rock samples indicates that metasomatic alteration accompanied the intrusion of both the granitic rocks and the serpentinite. The volcanic rocks, gabbro and serpentinite were derived from a magma of oceanic tholeiitic affinities. The stratigraphic sequence and chemistry of the volcanic rocks of the Rankin Inlet Group indicate that this assemblage is correlative with the Hurwitz Group rather than the Kaminak Group and is therefore Aphebian in age.

Geology of the Gander Group, Gander River Ultramafic Belt and Davidsville Group in the Jonathan's Pond - Weir's Pond Area, northeast Newfoundland /

The study area is situated in NE Newfoundland between Gander Lake and the north coast and on the boundary between the Gander and Botwood tectonostratigraphic zones (Williams et al., 1974). The area is underlain by three NE trending units; the Gander Group, the Gander River Ultramafic Belt (the GRUB) and the Davidsville Group. The easternmost Gander Group consists of a thick, psammitic unit composed predominantly of psammitic schist and a thinner, mixed unit of semipelitic and pelitic schist with minor psammite. The mixed unit may stratigraphically overlie the psammitic unit or be a lateral facies equivalent of the latter. No fossils have been recovered from the Gander Group. The GRUB is a terrain of mafic and ultramafic plutonic rocks with minor pillow lava and plagiogranite. It is interpreted to be a dismembered ophiolite in thrust contact with the Gander Group. The westernmost Davidsville Group consists of a basal conglomerate, believed deposited unconformably upon the GRUB from which it was derived, and an upper unit of greywacke and slate, mostly of turbidite origin, with minor limestone and calcareous sandstone. The limestone, which lies near the base of the unit, contains Upper Llanvirn to Lower Llandeilo fossils. The Gander and Davidsville Groups display distinctly different sedimentological , structural and metamorphic histories. The Gander Group consists of quartz-rich, relatively mature sediment. It has suffered three pre-Llanvirn deformations, of which the main deformation, Dp produced a major, NE-N-facing recumbent anticline in the southern part of the study area. Middle greenschist conditions existed from D^ to D- with growth of metamorphic minerals during each dynamic and static phase. In contrast, the mineralogically immature Davidsville Group sediment contains abundant mafic and ultramafic detritus which is absent from the Gander Group. The Davidsville Group displays the effects of a single penetrative deformation with localized D_ and D_ features, all of which can be shown to postdate D_ in the Gander Group. Rotation of the flat Gander S- into a subvertical orientation near the contact with the GRUB and the Davidsville Group is believed to be a Davidsville D^ feature. Regional metamorphism in the Davidsville Group is lower greenschist with a single growth phase, MS . These sedimentological, structural and metamorphic differences between the Gander and Davidsville Groups persist even where the GRUB is absent and the two units are in contact, indicating that the tectonic histories of the Gander and Davidsville Groups are distinctly different. Structural features in the GRUB, locally the result of multiple deformations, may be the result of Gander and/or Davidsville deformations. Metamorphism is in the greenschist facies. Geochemical analyses of the pillow lava suggest that these rocks were formed in a back-arc basin. Mafic intrusives in the Gander Group appear to be the result of magraatism separate from that producing the pillow lava. The Gander Group is interpreted to be a continental rise prism deposited on the eastern margin of the Late Precambrian-Lower Paleozoic lapetus Ocean. The GRUB, oceanic crust possibly formed in a marginal basin to the west, is believed to have been thrust eastward over the Gander Group, deforming the latter, during the pre-Llanvirnian, possibly Precambrian, Ganderian Orogeny. The Middle Ordovician and younger Davidsville Group was derived from, and deposited unconformably on, this deformed terrain. Deformation of the Davidsville Group occurred during the Middle Devonian Acadian Orogeny.

Structural analysis of interior layered deposits in north and south Coprates Chasma, Mars

Interior layered deposits within an embayment in the northern as well as near the southern wall of Coprates Chasma in the Valles Marineris, Mars are studied using HRSC, CTX, HiRISE and CRISM data. In the northern embayment, layered deposits outcrop in three separate locations (a western deposit, a central deposit and an eastern deposit). The central layered deposit in the north has a stratigraphic thickness of 2 km. The western layered deposit abuts against the chasma wall appearing to have a relatively un-eroded depositional surface. The eastern deposit is near a landslide scar which appears to have exposed basement layering showing downward displacement. This northern embayment is suggested to have been an ancestral basin. The triangular edged deposit near the southern wall of Coprates Chasma has an elongated mound protruding from the central edge and is suggested to be the outer limits of a fault block which is back rotated 6° south. The rotation may be the result of the Valles Marineris opening.

Modelling structural zeros in compositional data

This analysis was stimulated by the real data analysis problem of household expenditure data. The full dataset contains expenditure data for a sample of 1224 households. The expenditure is broken down at 2 hierarchical levels: 9 major levels (e.g. housing, food, utilities etc.) and 92 minor levels. There are also 5 factors and 5 covariates at the household level. Not surprisingly, there are a small number of zeros at the major level, but many zeros at the minor level. The question is how best to model the zeros. Clearly, models that try to add a small amount to the zero terms are not appropriate in general as at least some of the zeros are clearly structural, e.g. alcohol/tobacco for households that are teetotal. The key question then is how to build suitable conditional models. For example, is the sub-composition of spending excluding alcohol/tobacco similar for teetotal and non-teetotal households? In other words, we are looking for sub-compositional independence. Also, what determines whether a household is teetotal? Can we assume that it is independent of the composition? In general, whether teetotal will clearly depend on the household level variables, so we need to be able to model this dependence. The other tricky question is that with zeros on more than one component, we need to be able to model dependence and independence of zeros on the different components. Lastly, while some zeros are structural, others may not be, for example, for expenditure on durables, it may be chance as to whether a particular household spends money on durables within the sample period. This would clearly be distinguishable if we had longitudinal data, but may still be distinguishable by looking at the distribution, on the assumption that random zeros will usually be for situations where any non-zero expenditure is not small. While this analysis is based on around economic data, the ideas carry over to many other situations, including geological data, where minerals may be missing for structural reasons (similar to alcohol), or missing because they occur only in random regions which may be missed in a sample (similar to the durables)