Modeling calving process of glacier with dilated polyhedral discrete element method

Mass loss caused by glacier calving is one of the direct contributors to global sea level rise. Reliable calving laws are required for accurate modelling of ice sheet mass balance. Both continuous and discontinuous methods have been used for glacial calving simulations. In this study, the discrete element method(DEM) based on dilated polyhedral elements is introduced to simulate the calving process of a tidewater glacier. Dilated polyhedrons can be obtained from the Minkowski sum of a sphere and a core polyhedron. These elements can be utilized to generate a continuum ice material, where the interaction force between adjacent elements is modeled by constructing bonds at the joints of the common faces. A hybrid fracture model considering fracture energy is introduced. The viscous creep behavior of glaciers on long-term scales is not considered. By applying buoyancy and gravity to the modelled glacier, DEM results show that the calving process is caused by cracks which are initialized at the top of the glacier and spread to the bottom. The results demonstrate the feasibility of using the dilated polyhedral DEM method in glacier simulations, additionally allowing the fragment size of the breaking fragments to be counted. The relationship between crack propagation and internal stress in the glacier is analyzed during calving process. Through the analysis of the Mises stress and the normal stress between the elements, it is found that geometric changes caused by the glacier calving lead to the redistribution of the stress. The tensile stress between the elements is the main influencing factor of glacier ice failure. In addition, the element shape,glacier base friction and buoyancy are studied, the results show that the glacier model based on the dilated polyhedral DEM is sensitive to the above conditions.     

Late Quaternary Alluviation Along Intermittent Streams in Kakadu National Park,Northern Territory

Sedimentation along small, intermittent streams on Kapalga Research Station in Kakadu National Park may have responded to increased base levels following post-glacial flooding of the valleys of the Alligator Rivers. Alternatively, regional climate changes may have controlled sedimentation. Using thermoluminescence dating, we determined that sediments from two streams at Kapalga date from 21.5 ± 4.0 ka. On a third stream, sediments dated from 7.6 ± 1.1 ka, with younger sediments occurring downstream. We interpreted the pre-Holocene dates and the lack of evidence of upstream progradation to indicate that climate variation was more important to sedimentation than base levels. Predicted increases in rainfall variability and in the frequency of high-intensity rainfall under enhanced greenhouse conditions may cause renewed sediment mobilisation. At the outflow of one stream on to the South Alligator flood plain, we found 15 m of sandy alluvia underlying 3-5 m of estuarine muds deposited as a result of sea-level rise. These sandy alluvia dated from about 77 ka at 4 m to more than 300 ka at 19 m depth. These ages are consistent with those recorded on the Magela Creek system, 50 km to the east.     

Late-Variscan antimony mineralisation in the Rheinisches Schiefergebirge, NW Germany: evidence for stibnite precipitation by drastic cooling of high-temperature fluid systems

Antimony-rich vein mineralisation is widespread in the German part of the Variscan orogenic belt. Mineralogical investigation of a representative suite of these deposits, coupled with fluid inclusion characterisation and microthermometry, permits a reconstruction of their genetic evolution. Two structural settings host antimony mineralisation: the cores or flanks of anticlinal zones and major lithological contrasts. Channelled migration of geothermal fluids through permeable rock sequences and later stagnation of fluids in cap-rock situations inside the anticlinal zones led to mineral deposition. The mineralising event is interpreted as relating to input of deep-sourced fluids during late-orogenic exhumation at the transitional stage between collision tectonics and the late-Variscan extensional regime. Fluid inclusion data, chlorite geothermometry and the presence of meneghinite as a characteristic Pb-Sb-sulfosalt mineral in a number of vein systems allows constraints on model P-T conditions at the onset of mineralisation to be made. These are as high as 390 to 440?°C at 0.6–1.0?kbar for the Saarsegen, Apollo and Schöne Freundschaft deposits, with lower temperatures of 320–340?°C being obtained for the Spes deposit. The fluid inclusion data indicate drastic fluid cooling during the mineralising event; minimum temperatures of approximately 150–220?°C are obtained for all deposits at the end of vein quartz formation, which coincided with deposition of stibnite and most of the Pb-Sb sulfosalts. Besides the formation of extensional quartz-stibnite-Pb-sulfosalt veins, the mineralising, low-salinity NaCl-KCl-rich high-temperature tectonic brines have overprinted sulfide assemblages within earlier siderite-(Cu)-Pb-Zn veins. This has led to replacement reaction textures and remobilisation of sulfide components within the vein systems. In contrast with the earlier siderite-(Cu)-Pb-Zn veins, neither the quartz-stibnite-sulfosalt nor the (Cu)-Pb-Sb sulfosalt assemblages were affected by Variscan deformation. Rather, they display characteristic extensional features crosscutting all earlier structures and can thus be assigned to a later phase of mineralisation. Fluid composition characteristics and structural criteria indicate formation in the latest part of the Variscan mineralisation cycle; a post-Variscan genesis being rejected on grounds of conspicuously diverging fluid characteristics. A comparison of antimony deposits in the Rheinisches Schiefergebirge with other late-orogenic deposits elsewhere in the European Variscan belt indicates a significant number of shared features, enabling them to be placed into a common model related to the onset of late-Variscan brittle extensional tectonics.     

FIAs(Foliation Intersection/Inflection Axes) Preserved in Porphyroblasts,the DNA of Deformation:A Solution to the Puzzle of Deformation and Metamorphism in the Colorado,Rocky Mountains,USA

Abstract: FIAs have been used extensively for more than a decade to unravel deformation and metamorphic puzzles. Orogenic processes developing early during the history or orogenesis challenge scientists because compositional layering in rocks always reactivates where multiple deformations have occurred, leaving little evidence of the history of foliation development preserved in the matrix. The foothills of the Rocky Mountains in Colorado, USA contain a succession of four FIA sets (trends) that would not have been distinguishable if porphyroblasts had not grown during the multiple deformation events that affected these rocks or if they had rotated as these events took place. They reveal that both the partitioning of deformation and the location of isograds changed significantly as the deformation proceeded.     

The Triassic of Timor: Lithostratigraphy, chronostratigraphy and palaeogeography

The palaeontologically rich and lithologically diverse Triassic successions of Timor provide a key stratigraphic and palaeontological link between northwestern Australia and other terranes of former eastern Gondwana (present-day Southeast Asia). Timor is now located in the zone of collision between the northern margin of the Australian continent and island arc terranes bordering the Eurasian plate, with the Triassic successions exposed in a fold-and-thrust belt and an extensive mélange complex. Three formal lithostratigraphic units have been defined previously within the main Triassic succession in Timor (Niof, Aitutu and Babulu formations), with a fourth, the Wai Luli Formation, primarily Jurassic in age but extending down into the Triassic. The Niof Formation (Anisian to Ladinian, possibly also Early Triassic) is a fine-grained deepwater succession, succeeded conformably by the Aitutu and Babulu formations (Ladinian to Norian/Rhaetian), which were deposited contemporaneously, with the Aitutu Formation continuing locally into the Lower Jurassic. The Aitutu Formation consists of deep shelf limestones interbedded with shales and marls, while the Babulu Formation is a deltaic to turbiditic siliciclastic succession. The Late Triassic to Jurassic Wai Luli Formation is characterised by marine shales and marls.Informal stratigraphic units include the Cephalopod Limestone Facies, a Rosso Ammonitico-type deposit, which contains an extremely rich fossil fauna (particularly ammonoids) and ranges through the entire Triassic; and the Fatu Limestone and Pualaca Facies which consists of shallow to marginal marine carbonates (mud mounds, oolitic limestones and reefs) restricted to the Late Triassic. Facies diversity was low during the Early Triassic and Anisian, but became more pronounced from the Ladinian and continuing through the Late Triassic, probably as a consequence of renewed tectonic extension. Triassic extension was not associated with major volcanism, unlike a previous phase of extension in the Early Permian.The Cablac Limestone Formation, originally defined as a Miocene stratigraphic element, is now recognised to be at least partly Late Triassic–Early Jurassic in age, with lithologies comparable to parts of the Fatu Limestone. The stratigraphy of these shallow marine carbonate sequences is clearly in need of rigorous revision, but it is not yet possible to suggest appropriate redefined formations.     

Reassessing the management of groundwater use from sandy aquifers: acidification and base cation depletion exacerbated by drought and groundwater withdrawal on the Gnangara Mound, Western Australia

The combined effects of low rainfall, groundwater withdrawal in excess of 300 GL/year and reduced recharge in areas covered by pine plantations has caused the water table in a sandy unconfined aquifer on the Gnangara Mound in Western Australia to drop by up to 5 m and aquifer storage to decline by about 500 GL over the last 20 years. Groundwater has become acidic in areas of high drawdown, with pH values typically being less than 5.0 at the water table, and elevated concentrations of SO₄ ^2?, Al, Fe, Zn, Cu, Ni and Pb. Trends of increasing acidity and base cation concentrations in deep water supply wells in the Mirrabooka wellfield indicate that about 0.7 keq/ha/year of base cations are being leached from soil within cones of depression of pumping wells. These results indicate that the assessment of the sustainable yields of aquifers under conditions of low rainfall needs to consider geochemical interactions between groundwater, aquifer sediments, soils and vegetation, and not be just based on aquifer hydraulics and water-balance changes.     

The effect of fluid composition on the behavior of well cemented, quartz-rich sandstone during faulting

Previous studies [O'Kane, A., Onasch, C.M., Farver, J., 2007. The role of fluids in low-temperature, fault-related deformation of quartz arenite. Journal of Structural Geology 29, 819--836; Cook, J., Dunne, W.M., Onasch, C.M., 2006. Development of a dilatant damage zone along a thrust relay in a low-porosity quartz arenite, Journal of Structural Geology 28, 776–792] found that quartz arenite within two fault zones in the Appalachian foreland thrust belt displays very different structural styles and histories despite deforming at similar pressures and temperatures during the Alleghanian orogeny. A comparison of the grain-scale deformation and fluid histories using transmitted and cathodoluminescence microscopy and fluid inclusion microthermometry, shows that fluid composition was a controlling factor for causing these differences. The Cove fault zone deformed in the presence of aqueous fluids, first a CaCl₂ brine and then an iron-rich fluid. The precipitation of quartz cement from the brine kept pace with brittle deformation, maintaining overall rock cohesion in the fault zone. The introduction of iron-rich fluids caused a switch from precipitation to dissolution of quartz, along with precipitation of goethite. In a damage zone along a backthrust in the Cave Mountain anticline, early deformation occurred in the presence of an aqueous fluid from which quartz was precipitated. The latest deformation, however, occurred in the presence of a methane-rich fluid, which inhibited the precipitation of quartz cement producing porous breccias and open fractures despite deformation at 6 km depth. Fluid composition not only affected cementation in the fault zones, but also the selection of grain-scale deformation mechanisms. Therefore, it is a controlling factor in determining the behavior and strength of these fault zones.     

The thermoinsulation effect of snow cover within a climate model

We use a state of the art climate model (CAM3–CLM3) to investigate the sensitivity of surface climate and land surface processes to treatments of snow thermal conductivity. In the first set of experiments, the thermal conductivity of snow at each grid cell is set to that of the underlying soil (SC-SOIL), effectively eliminating any insulation effect. This scenario is compared against a control run (CTRL), where snow thermal conductivity is determined as a prognostic function of snow density. In the second set of experiments, high (SC-HI) and low (SC-LO) thermal conductivity values for snow are prescribed, based on upper and lower observed limits. These two scenarios are used to envelop model sensitivity to the range of realistic observed thermal conductivities. In both sets of experiments, the high conductivity/low insulation cases show increased heat exchange, with anomalous heat fluxes from the soil to the atmosphere during the winter and from the atmosphere to the soil during the summer. The increase in surface heat exchange leads to soil cooling of up to 20 K in the winter, anomalies that persist (though damped) into the summer season. The heat exchange also drives an asymmetric seasonal response in near-surface air temperatures, with boreal winter anomalies of +6 K and boreal summer anomalies of −2 K. On an annual basis there is a net loss of heat from the soil and increases in ground ice, leading to reductions in infiltration, evapotranspiration, and photosynthesis. Our results show land surface processes and the surface climate within CAM3–CLM3 are sensitive to the treatment of snow thermal conductivity.     

Textural variation in the pyrite-rich ore deposits of the Røros district, Trondheim Region, Norway: implications for pyrite deformation mechanisms

The Røros district is a pyrite-rich polymetallic sulfide orefield in the southeastern part of the Trondheim region, Central Norwegian Caledonides. All of the ore deposits at Røros are hosted within a Cambrian to Silurian succession that was deformed and metamorphosed at lower greenschist to lower amphibolite facies conditions during the Caledonian orogeny. Samples from five individual deposits across the orefield have been analyzed using a combination of reflected light petrographic observation, orientation contrast imaging, and electron backscatter diffraction. Results indicate that, whereas samples from each ore deposit have a variety of different textures, all of them preserve plastic deformation in pyrite grains that occurred at peak metamorphic conditions characterized by the development of internal lattice misorientation within pyrite grains and low-angle (～2°) dislocation walls. These observations indicate that the principal deformation mechanisms at peak metamorphic conditions were dislocation glide and creep. The preservation of brittle fracturing represents later overprinting events.