Effects of temperature-dependent material properties and radioactive heat production on simple basin subsidence models

The effects of applying laboratory-derived material parameters to simple thermal basin subsidence models are examined. The temperature-dependent properties of thermal conductivity, specific heat and the coefficient of thermal expansion are considered, as well as conductivity contrasts between the crust and mantle and lithosphere-scale radiogenic heat production. The effects of conductivity and radioactivity are complementary, and can be predicted from the way in which they influence the initial and final steady-state geothermal gradient. A convex-up initial gradient, such as is generated by the conductivity and radioactivity functions and is also assumed to be the case based on independent evidence, will lead to more initial subsidence than in a simple constant-parameter model. The final subsidence will also be greater in the modified model, due to the fact that the stretching event will result in a lithospheric column which is intrinsically cooler than in the pre-deformation case. Because the coefficient of thermal expansion rises with increasing temperature, including its temperature dependence will result in a model with substantially less initial subsidence than one with a constant value. When these parameters are combined into a single model, the initial subsidence is approximately 15–20% less than in the constant-parameter model if radioactivity is not included, and 10% less if it is, while the final subsidence is about 5% greater without radioactivity and 7–9% greater with radioactivity included. Depending on the magnitude of extension, these effects can translate into differences of tens to hundreds of meters when compared to the constant-parameter model.     

Geology of the epithermal Ag–Au Huevos Verdes vein system and San José district, Deseado massif, Patagonia, Argentina

The San José district is located in the northwest part of the Deseado massif and hosts a number of epithermal Ag–Au quartz veins of intermediate sulfidation style, including the Huevos Verdes vein system. Veins are hosted by andesitic rocks of the Bajo Pobre Formation and locally by rhyodacitic pyroclastic rocks of the Chon Aike Formation. New ⁴⁰Ar/³⁹Ar constraints on the age of host rocks and mineralization define Late Jurassic ages of 151.3 ± 0.7 Ma to 144.7 ± 0.1 Ma for volcanic rocks of the Bajo Pobre Formation and of 147.6 ± 1.1 Ma for the Chon Aike Formation. Illite ages of the Huevos Verdes vein system of 140.8 ± 0.2 and 140.5 ± 0.3 Ma are 4 m.y. younger than the volcanic host rock unit. These age dates are among the youngest reported for Jurassic volcanism in the Deseado massif and correlate well with the regional context of magmatic and hydrothermal activity. The Huevos Verdes vein system has a strike length of 2,000 m, with several ore shoots along strike. The vein consists of a pre-ore stage and three main ore stages. Early barren quartz and chalcedony are followed by a mottled quartz stage of coarse saccharoidal quartz with irregular streaks and discontinuous bands of sulfide-rich material. The banded quartz–sulfide stage consists of sulfide-rich bands alternating with bands of quartz and bands of chlorite ± illite. Late-stage sulfide-rich veinlets are associated with kaolinite gangue. Ore minerals are argentite and electrum, together with pyrite, sphalerite, galena, chalcopyrite, minor bornite, covellite, and ruby silver. Wall rock alteration is characterized by narrow (< 3 m) halos of illite and illite/smectite next to veins, grading outward into propylitic alteration. Gangue minerals are dominantly massive quartz intergrown with minor to accessory adularia. Epidote, illite, illite/smectite, and, preferentially at deeper levels, Fe-chlorite gangue indicate near-neutral pH hydrothermal fluids at temperatures of >220°C. Kaolinite occurring with the late sulfide-rich veinlet stage indicates pH < 4 and a temperature of <200°C. The Huevos Verdes system has an overall strike of 325°, dipping on average 65° NE. The orientations of individual ore shoots are controlled by vein strike and intersecting north-northwest-striking faults. We propose a structural model for the time of mineralization of the San José district, consisting of a conjugate shear pair of sinistral north-northwest- and dextral west-northwest-striking faults that correspond to R and R′ in the Riedel shear model and that are related to master faults (M) of north-northeast-strike. Veins of 315° strike can be interpreted as nearly pure extensional fractures (T). Variations in vein strike predict an induced sinistral shear component for strike directions of >315°, whereas strike directions of <315° are predicted with an induced dextral strike–slip movement. The components of the structural model appear to be present on a regional scale and are not restricted to the San José district.     

Determination of Radiogenic 87Sr/86Sr and Stable δ88/86SrSRM987 Isotope Values of Thirteen Mineral,Vegetal and Animal Reference Materials by DS‐TIMS

This work presents new ⁸⁷Sr/⁸⁶Sr and δ^88/86Sr_SRM987 isotopic values of thirteen mineral, vegetal and animal reference materials. Except for UB‐N, all our results are consistent with previously published data. Our results highlight intermediate precisions among the best presently published and a non‐significant systematic shift with the calculated δ^88/86Sr_SRM987 mean values for the three most analysed reference materials in the literature (i.e., IAPSO, BCR‐2 and JCp‐1). By comparison with the literature and between two distinct digestions, a significant bias of δ^88/86Sr_SRM987 values was highlighted for two reference materials (UB‐N and GS‐N). It has also been shown that digestion protocols (nitric and multi‐acid) have a moderate impact on the δ^88/86Sr_SRM987 isotopic values for the Jls‐1 reference materials suggesting that a nitric acid digestion of carbonate can be used without significant bias from partial digestion of non‐carbonate impurities. Different δ^88/86Sr_SRM987 values were measured after two independent Sr/matrix separations, according to the same protocol, for a fat‐rich organic reference material (BCR‐380R) and have been related to a potential post‐digestion heterogeneity. Finally, the δ^88/86Sr_SRM987 value differences measured between animal‐vegetal and between coral‐seawater reference materials agree with the previously published results, highlighting an Sr isotopic fractionation along the trophic chain and during carbonate precipitation.     

Emplacement age and thermal footprint of the diamondiferous Ellendale E9 lamproite pipe, Western Australia

The diamondiferous Ellendale 9 (E9) pipe is a funnel-shaped maar-diatreme volcano consisting of inward-dipping tuff sequences intruded by lamproite plugs and dykes. The host rocks for the E9 pipe are Permian sandstones. The multiple lithological contacts exposed within the mined maar volcano provide a natural laboratory in which to study the effect of volcanic processes on U–Th–Pb–He systematics. Zircon from the regional sandstone and E9 lamproite display a bimodal distribution of ages on (U–Th)/He–U/Pb plots. The zircon U/Pb ages for the E9 pipe (n?=?52) range from 440 to 2,725 Ma, while the cluster of (U–Th)/He ages for the lamproite dyke zircon indicate that dyke emplacement occurred at 20.6?±?2.8 Ma, concordant with a maximum emplacement age of about ≤22 Ma from phlogopite ⁴⁰Ar/³⁹Ar. These ages indicate a xenocrystic origin for the zircon entrained in the E9 dyke. The U/Pb ages of detrital zircon from the regional sandstone host (373–3,248 Ma; n?=?41) are indistinguishable from those of the lamproite zircon xenocrysts, whereas the detrital zircon in the host sandstone yield (U–Th)/He ages from 260 to 1,500 Ma. A thermochronology traverse across the E9 lamproite dyke reveals that the zircon (U–Th)/He ages in the host sandstone have not been significantly thermally reset during dyke emplacement, even at the contact. The capability of the zircon (U–Th)/He method to distinguish deep, mantle source lithologies from upper crustal source lithologies could be used in geochemical exploration for diamonds. Pre-screening of detrital samples using etching and helium assay methods will improve the efficiency and decrease the cost of greenfields exploration.     

Transition from arid to hyper-arid environment in the southern Levant deserts as recorded by early Pleistocene cummulic Aridisols

The time at which deserts established their current arid or hyper-arid conditions remains a fundamental question regarding the history of Earth. Cosmogenic isotope exposure ages of desert pavement and welded, calcic–gypsic–salic Reg soils that developed on relatively flat alluvial surfaces ～2 Ma ago in the Negev Desert indicate long geomorphic stability under extremely dry conditions. Over a short interval during their initial stage of development between 1–2 Ma, these cumulative soils are characterized by calcic soils reaching maximum stage III of carbonate morphology. This interval is the only period when calcic soil horizons formed on stable abandoned alluvial surfaces in the southern Negev Desert. Since ～1 Ma pedogenesis changed toward more arid soil environment and the formation of gypsic–salic soil horizons that were later followed by dust accumulation. The dichotomy of only moderately-developed calcic soil (stages II–III) during a relatively long time interval (10⁵–10⁶ years) indicates an arid environment that does not support continuous development but only occasional calcic soil formation. The very low δ¹⁸O and relatively high δ¹³C values of these early pedogenic carbonates support soil formation under arid climatic conditions. Such an environment was probably characterized by rare and relatively longer duration rainstorms which occasionally allowed deeper infiltration of rainwater and longer retention of soil moisture. This, in turn enabled the growth of sparse vegetation that enhanced deposition of pedogenic carbonate. At ～1 Ma these rare events of slightly wetter conditions ceased and less atmospheric moisture reached the southern Negev Desert leading to deposition of soluble salts and dust deposited in the soils. The combination of long-term hyperaridity, scarcity of vegetation and lack of bioturbation, salts cementation, dust accumulation and tight desert pavement cover, has protected the surfaces from erosion forming one of the most remarkably stable landscapes on Earth, a landscape that essentially has not eroded, but accumulated salt and dust for more than 10⁶ yr.     

Stochastic wobble of accretion discs and jets from turbulent rocket torques

Models of accretion discs and their associated outflows often incorporate assumptions of axisymmetry and symmetry across the disc plane. However, for turbulent discs these symmetries only apply to averaged quantities and do not apply locally. The local asymmetries can induce local imbalances in outflow power across the disc mid-plane, which can in turn induce local tilting torques. Here we calculate the effect of the resulting stochastic torques on disc annuli that are a consequence of standard mean field accretion disc models. The torques induce a random walk of the vector perpendicular to the plane of each averaged annulus. This random walk is characterized by a radially dependent diffusion coefficient which we calculate for small angle tilt. We use the coefficient to calculate a radially dependent time-scale for annular tilt and associated jet wobble. The wobble time depends on the square of the wander angle so the age of a given system determines the maximum wobble angle. We apply this to examples of blazars, young stellar objects and binary engines of pre-planetary nebulae and microquasars. It is noteworthy that for an averaging time   t _w∼ 3 d  , we estimate a wobble angle for jets in SS 433 of  θ∼ 0.8°  , not inconsistent with observational data. In general the non-periodic nature of the stochastic wobble could distinguish it from faster periodic jet precession.     

Why chromatic imaging matters

During the last two decades, the first generation of beam combiners at the Very Large Telescope Interferometer has proved the importance of optical interferometry for high-angular resolution astrophysical studies in the near- and mid-infrared. With the advent of 4-beam combiners at the VLTI, the u ? v coverage per pointing increases significantly, providing an opportunity to use reconstructed images as powerful scientific tools. Therefore, interferometric imaging is already a key feature of the new generation of VLTI instruments, as well as for other interferometric facilities like CHARA and JWST. It is thus imperative to account for the current image reconstruction capabilities and their expected evolutions in the coming years. Here, we present a general overview of the current situation of optical interferometric image reconstruction with a focus on new wavelength-dependent information, highlighting its main advantages and limitations. As an Appendix we include several cookbooks describing the usage and installation of several state-of-the art image reconstruction packages. To illustrate the current capabilities of the software available to the community, we recovered chromatic images, from simulated MATISSE data, using the MCMC software SQUEEZE. With these images, we aim at showing the importance of selecting good regularization functions and their impact on the reconstruction.     

Numerical modeling of nonlinear water waves over heterogeneous porous beds

The transformation of nonlinear water waves over porous beds is studied by applying a numerical model based on Chen's [2006. Fully nonlinear Boussinesq-type equations for waves and currents over porous beds. Journal of Engineering Mechanics, 132:2, 220–230] Boussinesq-type equations for highly nonlinear waves on permeable beds. The numerical model uses a high-order time-marching solution and fourth-order finite-difference schemes for discretization of first-order spatial derivatives to obtain a computational accuracy consistent with the model equations. By forcing the wave celerity and spatial porous-damping rate of the linearized model to match the exact linear theory for horizontal porous bed over a prescribed range of relative depths, the values of the model parameters are optimally determined. Numerical simulations of the damped wave propagation over finite-thickness porous layer demonstrate the accuracy of both the numerical model and governing equations, which have been shown by prior theoretical analyses to be accurate for both nominal and thick porous layers. These simulations also elucidate on the significance of the higher-order porous-damping terms and the influence of the hydraulic parameters. Application of the model to the simulation of the wave field around a laboratory-scale submerged porous mound provides a measure of its capability, as well as useful insight into the scaling of the porous-resistance coefficients. For application to heterogeneous porous beds, the assumption of weak spatial variation of the porous resistance is examined using truncated forms of the governing equations. The results indicate that the complete set of Boussinesq-type equations is applicable to porous beds of nonhomogeneous makeup.     

A Leaky-Conduit Model of Transient Flow in Karstic Aquifers

Karst Flow Model (KFM) simulates transient flow in an unconfined karstic aquifer having a well-developed conduit system. KFM treats the springshed as a two-dimensional porous matrix containing a triangulated irregular network of leaky conduits. The number and location of conduits can be specified arbitrarily, perhaps using field information as a guide, or generated automatically. Conduit networks can be tree-like or braided. Rainwater that has infiltrated down from the surface leaks into the conduits from the adjacent porous matrix at a rate dictated by Darcy’s law, then flows turbulently to the spring via the conduits. KFM is calibrated using the known steady state; geometry and recharge determine the steady fluxes in the conduits, and the head distribution determines conduit gradients and sizes. Spring flow can vary with time due to spatially and temporally variable recharge and due to prescribed variations in the elevation of the spring. KFM is illustrated by four examples run on a test aquifer consisting of 27 nodes, 42 elements, and 26 conduits. Three examples (drought, uniform rainstorm, storm-water input to one element) are simulations, while the fourth uses data from a spring-basin flooding event. The qualitative fit between the predicted and observed spring discharge in the fourth example provides support of the hypothesis that the dynamic behavior of a karst conduit system is an emergent property of a self-organized system, largely independent of the locations and properties of individual conduits.