Radiative heating of carbonaceous near-Earth objects as a cause of thermal metamorphism for CK chondrites

Metamorphic CK carbonaceous chondrites display matrix textures that are best explained by a transient thermal event with temperatures in the 550–950 K range and durations in the order of days to years, longer than what is commonly admitted for shock events but shorter than what is required for nuclide decay. We propose that radiative heating of small carbonaceous meteoroids with perihelia close to the Sun could account for the petrological features observed in CK chondrites. Numerical thermal modeling, using favorable known NEOs orbital parameters (perihelion distances between 0.07 and 0.15 AU) and physical properties of CV and CK chondrites (albedo in the range 0.01–0.1, 25% porosity, thermal diffusivity of 0.5–1.5 W m^?1 K^?1), shows that radiative heating can heat carbonaceous meteoroids in the meter size range to core temperatures up to 1050 K, consistent with the metamorphic temperatures estimated for CK chondrites. Sizes of known CV and CK chondrites indicate that all these objects were small meteoroids (radii from a few cm to 2.5 m) prior to their atmospheric entry. Simulations of dynamic orbits for NEO objects suggest that there are numerous such bodies with suitable orbits and properties, even if they are only a small percentage of all NEOs. Radiative heating would be a secondary process (superimposed on parent-body processes) affecting meteoroids formed by the disruption of an initially homogeneous CV3-type parent body. Different petrologic types can be accounted for depending on the sizes and heliocentric distances of the objects in such a swarm.     

MODFLOW as a Configurable Multi-Model Hydrologic Simulator

MODFLOW 6 is the latest in a line of six “core” versions of MODFLOW released by the U.S. Geological Survey. The MODFLOW 6 architecture supports incorporation of additional hydrologic processes, in addition to groundwater flow, and allows interaction between processes. The architecture supports multiple model instances and multiple types of models within a single simulation, a flexible approach to formulating and solving the equations that represent hydrologic processes, and recent advances in interoperability, which allow MODFLOW to be accessed and controlled by external programs. The present version of MODFLOW 6 consolidates popular capabilities available in MODFLOW variants, such as the unstructured grid support in MODFLOW-USG, the Newton-Raphson formulation in MODFLOW-NWT, and the support for partitioned stress boundaries in MODFLOW-CDSS. The flexible multi-model capability allows users to configure MODFLOW 6 simulations to represent the local-grid refinement (LGR) capabilities available in MODFLOW-LGR, the multi-species transport capabilities in MT3DMS, and the coupled variable-density capabilities available in SEAWAT. This paper provides a new, holistic and integrated overview of simulation capabilities made possible by the MODFLOW 6 architecture, and describes how ongoing and future development can take advantage of the program architecture to integrate new capabilities in a way that is minimally invasive and automatically compatible with the existing MODFLOW 6 code.     

Hydraulic-Head Formulation for Density-Dependent Flow and Transport

Density-dependent flow and transport solutions for coastal saltwater intrusion investigations, analyses of fluid injection into deep brines, and studies of convective fingering and instabilities of denser fluids moving through less dense fluids typically formulate the groundwater flow equation in terms of pressure or equivalent freshwater head. A formulation of the flow equation in terms of hydraulic head is presented here as an alternative. The hydraulic-head formulation can facilitate adaptation of existing constant-density groundwater flow codes to include density-driven flow by avoiding the need to convert between freshwater head and hydraulic head within the code and by incorporating density-dependent terms as a compartmentalized “correction” to constant-density calculations already performed by the code. The hydraulic-head formulation also accommodates complexities such as unconfined groundwater flow and Newton-Raphson solution schemes more readily than the freshwater-head formulation. Simulation results are presented for four example problems solved using an implementation of the hydraulic-head formulation in MODFLOW.     

Environmental heterogeneity affects input, storage, and transformation of coarse particulate organic matter in a floodplain mosaic

Quantifying spatial and temporal dynamics of organic matter (OM) is critical both for understanding ecosystem functioning and for predicting impacts of landscape change. To determine the influence of different habitats and coarse particulate OM (CPOM) types upon floodplain OM dynamics, we quantified aerial input, lateral surface transfer, and surface storage of CPOM over an annual cycle on the near-natural floodplain of the River Tagliamento in NE-Italy. Using these data, we modelled floodplain leaf dynamics, taking account of the spatial distribution and hydrologic connectivity of habitats, and using leaf storage as a response variable. Mean aerial CPOM input to the floodplain was similar from riparian forest and islands, but surface transfer was greater from islands, supporting the suggestion that these habitats act as “islands of fertility” along braided rivers. Leaves were the lateral conveyor of energy to more open parts of the floodplain, whereas CPOM was mainly stored as small wood in vegetated islands and riparian forest. Simulating the loss of habitat diversity (islands, ponds) decreased leaf storage on the whole floodplain, on exposed gravel and in large wood accumulations. In contrast, damming (loss of islands, ponds and floods plus floodplain overgrowth) greatly increased storage on exposed gravel. A random shuffle of habitats led to a storage increase on exposed gravel, while that in large wood accumulations and ponds declined. These results disentangle some of the complexities of CPOM dynamics in floodplain ecosystems, illustrate the value of models in understanding ecosystem functioning at a landscape level, and directly inform river management practice.     

Observations of estuarine circulation and solitary internal waves in a highly energetic tidal channel

Despite vigorous tidal and wind mixing, observations in an estuarine tidal inlet in the Wadden Sea show that during part of the tidal cycle, vertical stratification and internal waves may still develop. Acoustic Doppler current profiler (ADCP) and conductivity, temperature, depth observations, collected over the past 6 years at 13 h anchor stations (ASs), reveal that these occur especially during slack tide, when there is little wind and large freshwater discharge from nearby Lake IJssel. Measurements with a moored ADCP show that in the same tidal phase, strong cross-channel circulation develops, which may suddenly reverse circulation sense due to passing density fronts. In the vertically stratified phase that follows after the front passage, propagating mode-one solitary internal waves are observed. These are resonantly generated during decelerating tidal ebb currents when the (shear) flow passes a transcritical regime (Froude number equal to 1). A combination of photographs (including one from the International Space Station), bathymetric data, and ASs data leads to the discovery of yet another source of internal waves in this area, produced during slackening tide by propagating lee waves that develop over a deep trench. We suggest that both the cross-channel circulation as well as the (solitary) internal waves may locally be of importance for the (re)distribution and transport of sediments and nutrients and may influence tidally averaged transports.     

Submarine hydrothermal activity and gold-rich mineralization at Brothers Volcano, Kermadec Arc, New Zealand

Brothers volcano, of the Kermadec intraoceanic arc, is host to a hydrothermal system unique among seafloor hydrothermal systems known anywhere in the world. It has two distinct vent fields, known as the NW Caldera and Cone sites, whose geology, permeability, vent fluid compositions, mineralogy, and ore-forming conditions are in stark contrast to each other. The NW Caldera site strikes for ??600?m in a SW?CNE direction with chimneys occurring over a ??145-m depth interval, between ??1,690 and 1,545?m. At least 100 dead and active sulfide chimney spires occur in this field and are typically 2?C3?m in height, with some reaching 6?C7?m. Their ages (at time of sampling) fall broadly into three groups: <4, 23, and 35?years old. The chimneys typically occur near the base of individual fault-controlled benches on the caldera wall, striking in lines orthogonal to the slopes. Rarer are massive sulfide crusts 2?C3?m thick. Two main types of chimney predominate: Cu-rich (up to 28.5?wt.% Cu) and, more commonly, Zn-rich (up to 43.8?wt.% Zn). Geochemical results show that Mo, Bi, Co, Se, Sn, and Au (up to 91?ppm) are correlated with the Cu mineralization, whereas Cd, Hg, Sb, Ag, and As are associated with the dominant Zn-rich mineralization. The Cone site comprises the Upper Cone site atop the summit of the recent (main) dacite cone and the Lower Cone site that straddles the summit of an older, smaller, more degraded dacite cone on the NE flank of the main cone. Huge volumes of diffuse venting are seen at the Lower Cone site, in contrast to venting at both the Upper Cone and NW Caldera sites. Individual vents are marked by low-relief (??0.5?m) mounds comprising predominately native sulfur with bacterial mats. Vent fluids of the NW Caldera field are focused, hot (??300°C), acidic (pH????2.8), metal-rich, and gas-poor. Calculated end-member fluids from NW Caldera vents indicate that phase separation has occurred, with Cl values ranging from 93% to 137% of seawater values. By contrast, vent fluids at the Cone site are diffuse, noticeably cooler (??122°C), more acidic (pH?1.9), metal-poor, and gas-rich. Higher-than-seawater values of SO₄ and Mg in the Cone vent fluids show that these ions are being added to the hydrothermal fluid and are not being depleted via normal water/rock interactions. Iron oxide crusts 3?years in age cover the main cone summit and appear to have formed from Fe-rich brines. Evidence for magmatic contributions to the hydrothermal system at Brothers includes: high concentrations of dissolved CO₂ (e.g., 206?mM/kg at the Cone site); high CO₂/³He; negative ??D and ??¹⁸O_H2O for vent fluids; negative ??³⁴S for sulfides (to ?4.6??), sulfur (to ?10.2??), and ??¹⁵N₂ (to ?3.5??); vent fluid pH values to 1.9; and mineral assemblages common to high-sulfidation systems. Changing physicochemical conditions at the Brothers hydrothermal system, and especially the Cone site, occur over periods of months to hundreds of years, as shown by interlayered Cu?+?Au- and Zn-rich zones in chimneys, variable fluid and isotopic compositions, similar shifts in ³He/⁴He values for both Cone and NW Caldera sites, and overprinting of ??magmatic?? mineral assemblages by water/rock-dominated assemblages. Metals, especially Cu and possibly Au, may be entering the hydrothermal system via the dissolution of metal-rich glasses. They are then transported rapidly up into the system via magmatic volatiles utilizing vertical (??2.5?km long), narrow (??300-m diameter) ??pipes,?? consistent with evidence of vent fluids forming at relatively shallow depths. The NW Caldera and Cone sites are considered to represent stages along a continuum between water/rock- and magmatic/hydrothermal-dominated end-members.     

Influence of bottom topography roughness on the jet and inertial recirculation of a mid-latitude gyre

Several numerical experiments are conducted to examine the influence of mesoscale, bottom topography roughness on the inertial circulation of a wind-driven, mid-latitude ocean gyre. The ocean model is based on the quasi-geostrophic formulation, and is eddy-resolving as it features high vertical and horizontal resolutions (six layers and a 10 km grid). An antisymmetrical double-gyre wind stress curl forces the baroclinic modes and generates a strong surface jet. In the case of a flat bottom, inertia and inverse energy cascade force the barotropic mode, and the resulting circulation features strong, barotropic, inertial gyres. The sea-floor roughness inhibits the inertial circulation in the deep layers; the barotropic component of the flow is then forced by eddy-topography interactions, and its energy concentrates at the scales of the topography. As a result, the baroclinicity of the flow is intesified: the barotropic mode is reduced with regard to the baroclinic modes, and the bottom flow (constrained by the mesoscale sea-floor roughness) is decoupled from the surface flow (forced by the gyre-scale wind). Rectified, mesoscale bottom circulation induces an interfacial form stress at the thermocline, which enhances horizontal shear instability and opposes the eastward penetration of the jet. The mean jet is consequently shortened, but the instantaneous jet remains very turbulent, with meanders of large meridional extent. The sea-floor roughness modifies the energy pathways, and the eddies have an even more important role in the establishment of the mean circulation: below the thermocline, rectification processes are dominant, and eddies transfer energy toward permanent mesoscale circulations strongly correlated with topography, whereas above the thermocline mean flow and eddy generation are influenced by the mean bottom circulation through interfacial stress. The topography modifies the vorticity of the barotropic and highest baroclinic modes. Vorticity accumulates at the small topographic scales, and the vorticity content of the highest modes, which is very weak in the flat-bottom case, increases significantly. Few changes occur in surface-intensified modes. In the deep layers of the model, the inverse correlation between relative vorticity and topography at small scales ensures the homogenization of the potential vorticity, which mainly retains the largest scales of the bottom flow and the scale of β.     

Climate control on the evolution of Late Pleistocene alluvial‐fan and aeolian sand‐sheet systems in NW Germany

The Late Pleistocene was characterized by rapid climate oscillations with alternation of warm and cold periods that lasted up to several thousand years. Although much work has been carried out on the palaeoclimate reconstruction, a direct correlation of ice‐core, marine and terrestrial records is still difficult. Here we present new data from late Middle Pleniglacial to Lateglacial alluvial‐fan and aeolian sand‐sheet deposits in northwestern Germany. Records of Late Pleniglacial alluvial fans in central Europe are very rare, and OSL dating is used to determine the timing of fan aggradation. In contrast to fluvial systems that commonly show a delay between climate change and incision/aggradation, the small alluvial‐fan systems of the Senne area responded rapidly to climatic changes and therefore act as important terrestrial climate archives for this time span. The onset of alluvial‐fan deposition correlates with the climate change from warm to cold at the end of MIS 3 (29.3±3.2 ka). Strong fan progradation started at 24.4±2.8 ka and may be related to a period of higher humidity. The vertical stacking pattern of sedimentary facies and channel styles indicate a subsequrent overall decrease in water and sediment supply, with less sustained discharges and more sporadic runoffs from the catchment area, corresponding to an increasing aridity in central Europe during the Late Pleniglacial. Major phases of channel incision and fan aggradation may have been controlled by millennial‐scale Dansgaard–Oeschger cycles. The incision of channel systems is attributed to unstable climate phases at cold–warm (dry–wet) or warm–cold (wet–dry) transitions. The alluvial‐fan deposits are bounded by an erosion surface and are overlain by aeolian sand‐sheets that were periodically affected by flash‐floods. This unconformity might be correlated with the Beuningen Gravel Bed, which is an important marker horizon in deposits of the Late Pleniglacial resulting from deflation under polar desert conditions. The deposition of aeolian sand‐sheet systems (19.6±2.1 to 13.1±1.5 ka) indicates a rapid increase in aridity at the end of the Late Pleniglacial. Intercalated flash‐floods deposits and palaeosols (Finow type) point to temporarily wet conditions during the Lateglacial. The formation of an ephemeral channel network probably marks the warm‐cold transition from the Allerød to the Younger Dryas.     

Emplacement of a Debris Avalanche during the 1883 eruption of Krakatau (Sunda Straits,Indonesia)

Conclusion The data collected during the Mentawai cruise help to clarify understanding of the 1883 eruption of Krakatau. We have previously discussed the weaknesses of the interpretation of Williams (1941) and others (Self and Rampino 1981) and emphasized that only a Mount St. Helens-type collapse during the course of the eruption could account for all the characteristics of the eruption and of the related deposits.The discovery on land of deposits attributable to a debris-avalanche, in the stratigraphic position where they were expected, is a strong argument for the validity of our scenario.Marine surveys confirm that the sea bottom around Krakatau is covered by a thick ignimbritic deposit. But the presence of this deposit does not invalidate the presence of a debris-avalanche deposit under the ignimbrites. The hummocky morphology favours this hypothesis.Flank-failure of volcanoes is generally considered as a very efficient mechanism for triggering tsunamis (Kienle et al. 1987; Siebert et al. 1987). However, the majority of the volcanoes where flank-failure has been described are tall and bulky and the collapse of a broad edifice like Krakatau may be surprising. However the geological evidence shows that such a mechanism can act at various scales; for example the flank collapse of Mayu Yama volcano (height 700 m, volume 0,3 km³), a parasitic cone of Unzen volcano (Japan), triggered a debris-avalanche into the sea that was 1 km long, with a characteristic hummocky surface; the resulting tsunami killed 9528 people (Katayama 1974). In the same way, a partial collapse of Iliwerung volcano, Indonesia (50 × 10⁶ m³) in July 1979, triggered a tsunami which killed several hundred people (McClelland et al. 1989). At Krakatau, the main summit was 822 m asl; the collapse took place along the edge of the prehistoric caldera and this structural unconformity probably facilitated the triggering of the process.