Mapping Curvilinear Structures with Local Anisotropy Kriging

Interpolating geo-data with curvilinear structures using geostatistics is often disappointing. Channels, for example, become disconnected sets of lakes when interpolated from point data. In order to improve the interpolation of geological structures (e.g., curvilinear structures), we present a new form of kriging, local anisotropy kriging (LAK). Local anisotropy kriging combines a gradient algorithm from image analysis with kriging in an iterative way. After an initial standard kriging interpolation, the gradient algorithm determines the local anisotropy for each cell in the grid using a search area around the cell. Subsequently, kriging is carried out with the spatially varying anisotropy. The anisotropy calculation and subsequent kriging steps will then succeed until the result is satisfactory in the way of reproducing the curvilinear structures. Depending on the size of the search area more or less detail in the geological structures can be reproduced with LAK. Using test examples we show that LAK interpolates data with curvilinear structures more realistically than standard kriging. In a real world case, using bathymetric data of the Oosterschelde estuary, LAK also proves to be quantitatively superior to standard kriging. Absolute interpolation errors are decreased by 23%. Local anisotropy kriging only uses information from point data, which makes the method very objective, it only presents “what the data can tell.”     

A study of the chemistry of pore fluids and authigenic carbonates in methane seep environments: Kodiak Trench, Hydrate Ridge, Monterey Bay, and Eel River Basin

Analyses of the chemical and isotopic composition of carbonates rocks recovered from methane seepage areas of the Kodiak Trench, Hydrate Ridge, Monterey Bay Clam Flats, and the Eel River Basin, coupled with the studies of the chemistry of the pore fluids, have shown that these carbonates have grown within the sediment column. Geochemical profiles of pore fluids show that, in deep water seeps (Kodiak Trench—4450 m; Monterey Bay—1000 m; Hydrate Ridge—650 m), δ¹³C (DIC) values are low (isotopically light), whereas in the Eel River area ( 350–500 m), δ¹³C (DIC) values are much higher (isotopically heavier). In all cases, the δ¹³C values indicate that processes of methane oxidation, associated with sulfate reduction, are dominant in the shallow sediments. Data on the isotopic composition of authigenic carbonates found at sites in Kodiak Trench, Eel River Basin South, and Eel River Basin North indicate a variable composition and origin in different geochemical environments. Some of the authigenic carbonates from the study sites show a trend in their δ¹³C values similar to those of the pore fluids obtained in their vicinity, suggesting formation at relatively shallow depths, but others indicate formation at greater sediment depths. The latter usually consist of high magnesium calcite or dolomite, which, from their high values of δ¹³C (up to 23‰;) and δ¹⁸O (up to 7.5‰), suggest formation in the deeper horizons of the sediments, in the zone of methanogenesis. These observations are in agreement with observations by other workers at Hydrate Ridge, in Monterey Bay, and in the Eel River Basin.     

The last British Ice Sheet: A review of the evidence utilised in the compilation of the Glacial Map of Britain

This paper reviews the evidence presently available (as at December 2003) for the compilation of the Glacial Map of Britain (see [Clark C.D., Evans D.J.A., Khatwa A., Bradwell T., Jordan C.J., Marsh S.H., Mitchell W.A., Bateman, M.D. , 2004. Map and GIS database of glacial landforms and features related to the last British Ice Sheet. Boreas 33, 359–375] and http://www.shef.ac.uk/geography/staff/clark_chris/britice.html) in an effort to stimulate further research on the last British Ice Sheet and promote a reconstruction of ice sheet behaviour based on glacial geology and geomorphology. The wide range of evidence that has been scrutinized for inclusion on the glacial map is assessed with respect to the variability of its quality and quantity and the existing controversies in ice sheet reconstructions. Landforms interpreted as being of unequivocal ice-marginal origin (moraines, ice-contact glacifluvial landforms and lateral meltwater channels) and till sheet margins are used in conjunction with available chronological control to locate former glacier and ice-sheet margins throughout the last glacial cycle. Subglacial landforms (drumlins, flutings and eskers) have been used to demarcate former flow patterns within the ice sheet. The compilation of evidence in a regional map is crucial to any future reconstructions of palaeo-ice sheet dynamics and will provide a clearer understanding of ice sheet configuration, ice divide migration and ice thickness and coverage for the British Ice Sheet as it evolved through the last glacial cycle.     

Bar-top hollows: A new element in the architecture of sandy braided rivers

Discrete hollows in the bar tops of the South Saskatchewan River are described that form a newly-recognized morphological element of sandy braided rivers. These bar-top hollows, which are up to 1.7 m deep and may extend for 10–30 m down and across flow, have a circular to ovoid planform and are shown, through use of ground penetrating radar, to be filled by a series of distinct, often angle-of-repose, foresets. The hollows form by both erosion and bar-top deposition and may be generated by bar-tail accretion, cross-bar channel cutoff and subsequent fill or lateral accretion at the bar-head. Bar-top hollows occur in the upper part of the bar depositional sequence and may thus prove useful indicators for braid bar reconstruction in ancient sediments, and should not be confused with channel scour.     

Fluid inclusion and stable isotope (O, H, C, and S) constraints on the genesis of the Serrinha gold deposit, Gurupi Belt, northern Brazil

The Serrinha gold deposit of the Gurupi Belt, northern Brazil, belongs to the class of orogenic gold deposits. The deposit is hosted in highly strained graphitic schist belonging to a Paleoproterozoic (∼2,160 Ma) metavolcano-sedimentary sequence. The ore-zones are up to 11 m thick, parallel to the regional NW–SE schistosity, and characterized by quartz-carbonate-sulfide veinlets and minor disseminations. Textural and structural data indicate that mineralization was syn- to late-tectonic and postmetamorphic. Fluid inclusion studies identified early CO₂ (CH₄-N₂) and CO₂ (CH₄-N₂)-H₂O-NaCl inclusions that show highly variable phase ratios, CO₂ homogenization, and total homogenization temperatures both to liquid and vapor, interpreted as the product of fluid immiscibility under fluctuating pressure conditions, more or less associated with postentrapment modifications. The ore-bearing fluid typically has 18–33mol% of CO₂, up to 4mol% of N₂, and less than 2mol% of CH₄ and displays moderate to high densities with salinity around 4.5wt% NaCl equiv. Mineralization occurred around 310 to 335°C and 1.3 to 3.0 kbar, based on fluid inclusion homogenization temperatures and oxygen isotope thermometry with estimated oxygen fugacity indicating relatively reduced conditions. Stable isotope data on quartz, carbonate, and fluid inclusions suggest that veins formed from fluids with δ¹⁸O_H2O and δD_H2O (310–335°C) values of +6.2 to +8.4‰ and −19 to −80‰, respectively, which might be metamorphic and/or magmatic and/or mantle-derived. The carbon isotope composition (δ¹³C) varies from −14.2 to −15.7‰ in carbonates; it is −17.6‰ in fluid inclusion CO₂ and −23.6‰ in graphite from the host rock. The δ³⁴S values of pyrite are −2.6 to −7.9‰. The strongly to moderately negative carbon isotope composition of the carbonates and inclusion fluid CO₂ reflects variable contribution of organic carbon to an originally heavier fluid (magmatic, metamorphic, or mantle-derived) at the site of deposition and sulfur isotopes indicate some oxidation of the originally reduced fluid. The deposition of gold is interpreted to have occurred mainly in response to phase separation and fluid-rock interactions such as CO₂ removal and desulfidation reactions that provoked variations in the fluid pH and redox conditions.     

Fractionation of the noble metals by physical processes

During partial melting in the earth’s mantle, the noble metals become fractionated. Os, Ir, Ru, and Rh tend to remain in the mantle residue whereas Pt, Pd, and Re behave mildly incompatible and are sequestered to the silicate melt. There is consensus that sulfide plays a role in the fractionation process; the major noble metal repository in the mantle is sulfide, and most primitive mantle melts are sulfide-saturated when they leave their mantle sources. However, with sulfide–silicate partitioning, the fractionation cannot be modeled properly. All sulfide–silicate partition coefficients are so extremely high that a silicate melt segregating from a mantle source with residual sulfide should be largely platinum-group elements free. We offer a physical alternative to sulfide–silicate chemical partitioning and provide a mechanism of generating a noble metal-rich melt from a sulfide-saturated source: Because sulfide is at least partially molten at asthenospheric temperature, it will behave physically incompatible during melt segregation, and a silicate melt segregating from a mantle residue will entrain molten residual sulfide in suspension and incorporate it in the basaltic pool melt. The noble metal abundances of a basalt then become independent of sulfide–silicate chemical partitioning. They reflect the noble metal abundances in the drained sulfide fraction as well as the total amount of sulfide entrained. Contrary to convention, we suggest that a fertile, sulfide-rich mantle source has more potential to generate a noble metal-enriched basaltic melt than a refractory mantle source depleted by previous partial melting events.     

历史地震损失重建研究

历史地震损失重建是地震损失预测的一种方法,也是用于区域地震风险分析模型验证的一种手段。在介绍了历史地震损失重建一般方法的基础上,对方法中等效全损住宅比例估计方法以及住宅结构地震易损性指数的估计方法进行了重点分析,给出了这两个关键参数的可选择的多种估计方法。     

Pit crater formation on Kilauea volcano, Hawaii.

Most of the known pit craters in Hawaii occur along the East and Southwest Rift Zones of Kilauea volcano. The pit craters typically are either astride a single rift zone fracture or between a pair of rift zone fractures. These fractures are prominent in the pit crater walls. The pit craters are elliptical in plan view, with their major diameters ranging from 8 to 1140 m. They range in depth from 6 m to 186 m. They typically develop with initially steep, locally overhanging walls, but as the walls collapse, the craters fill with talus and become shaped like inverted elliptical cones. None of the craters apparently formed as eruptive vents, although some have been subsequently filled by lava. Devil's Throat is the best-exposed pit crater along the East Rift Zone. It is sited at a `waist' between two east-striking zones of ground cracks; the spacing between the crack zones decreases towards Devil's Throat. East-striking fractures are also prominent in the pit crater walls. Pit craters along the Southwest Rift Zone typically are elongate in plan view along the direction of the rift, have large caves at their bases along the long axes of the craters, and are smaller than those of the East Rift Zone. Some closely spaced pits there have coalesced to form a trough. Based on our observations and mechanical considerations, we infer that pit craters form by stoping over an underlying large-aperture rift zone fracture, and not by piston-like collapse over broad magma bodies or voids. Flow of magma along the underlying fracture may remove stoped blocks and prevent the fracture from being choked with debris. This mechanism is consistent with pit crater location, ground crack patterns, the preferred orientation of fractures in pit crater walls, and pit crater geometry (both in map view and cross-section). The mechanism also fits with observations of stoping into a gaping rift fracture that conducted lava from Kilauea caldera during the 1920s. Additionally, the ratio of pit crater width to depth of 0.5 to 2 is consistent with pit craters forming over a nearly vertical opening mode fracture.