Regional and species specific sexual reproductive patterns of three zooxanthellate scleractinian corals across the Eastern Tropical Pacific

Sexual reproduction of zooxanthellate scleractinian corals in the Eastern Tropical Pacific (ETP) is influenced by the interactive effects of regional and local oceanographic conditions, as well as a variety of other environmental factors. Differences in spatial and temporal gamete development and reproductive patterns of three widespread scleractinian corals of this region—Pocillopora damicornis (branching colony morphology, characterized as hermaphrodite broadcaster), Pavona gigantea (massive colonies, characterized as gonochoric broadcaster/sequential co‐sexually hermaphroditic) and Porites panamensis (encrusting colonies, characterized as gonochoric brooder)—were evaluated at local and regional scales across the ETP. This area extended from the Gulf of California (24°N) to the southern coast of Ecuador (–1°S), including offshore islands, using existing data pooled from prior studies. Predictive models were employed based on environmental variables: sea surface temperature, daylight hours, diffuse attenuation co‐efficient at 490 nm and photosynthetically active radiation. Datasets were extracted from satellite images ( https://oceancolor.gsfc.nasa.gov ) and analysed using WAM_STATIST ver. 6.33® software to obtain monthly average data from each site. The spatial (region, sub‐region and site) and seasonal (wet, dry) variation in reproductive activity (% colonies with gametes) differed among the three species; significant interactions were season × sub‐region for P. damicornis, season × site (sub‐region) for P. gigantea, and season × site for P. panamensis, for which sub‐region was not considered as a factor. The predictive models also suggest that gamete production/maturation of P. damicornis and P. gigantea is influenced by local differences in sea surface temperature and daylight hours. Porites panamensis was not correlated with any environmental variable examined. Variations in local and regional reproductive developmental patterns are likely an acclimatization response by each species imposed by the timing and duration of local stressor events. This analysis has provided insights into the diverse local and regional physical drivers that affect species responses and acclimatization in sexual reproduction across the ETP.     

Using Resin‐Based 3D Printing to Build Geometrically Accurate Proxies of Porous Sedimentary Rocks

Three‐dimensional (3D) printing is capable of transforming intricate digital models into tangible objects, allowing geoscientists to replicate the geometry of 3D pore networks of sedimentary rocks. We provide a refined method for building scalable pore‐network models (“proxies”) using stereolithography 3D printing that can be used in repeated flow experiments (e.g., core flooding, permeametry, porosimetry). Typically, this workflow involves two steps, model design and 3D printing. In this study, we explore how the addition of post‐processing and validation can reduce uncertainty in the 3D‐printed proxy accuracy (difference of proxy geometry from the digital model). Post‐processing is a multi‐step cleaning of porous proxies involving pressurized ethanol flushing and oven drying. Proxies are validated by: (1) helium porosimetry and (2) digital measurements of porosity from thin‐section images of 3D‐printed proxies. 3D printer resolution was determined by measuring the smallest open channel in 3D‐printed “gap test” wafers. This resolution (400 µm) was insufficient to build porosity of Fontainebleau sandstone (～13%) from computed tomography data at the sample's natural scale, so proxies were printed at 15‐, 23‐, and 30‐fold magnifications to validate the workflow. Helium porosities of the 3D‐printed proxies differed from digital calculations by up to 7% points. Results improved after pressurized flushing with ethanol (e.g., porosity difference reduced to ～1% point), though uncertainties remain regarding the nature of sub‐micron “artifact” pores imparted by the 3D printing process. This study shows the benefits of including post‐processing and validation in any workflow to produce porous rock proxies.     

Understanding fracture distribution and its relation to knickpoint evolution in the Rio Icacos watershed (Luquillo Critical Zone Observatory,Puerto Rico) using landscape-scale hydrogeophysics

The Rio Icacos watershed in the Luquillo Mountains (Puerto Rico) is unique due to its extremely rapid weathering rates. The watershed is incised into a quartz diorite that has developed a large knickzone defining the river profile. Regolith thickness within the watershed generally decreases from 20 to 30 m at the ridges to several meters in the quartz diorite-dominated valley to tens of centimeters near the major river knickpoint, as determined from previous studies. Above the knickzone, we observe spheroidal corestones, but below this weathering is much less apparent. Measured erosion rates from previous studies are also high in the knickzone compared with upper elevations within the river profile. A suite of near-surface geophysical methods (i.e. ground penetrating radar and terrain conductivity) capable of fast data acquisition in rugged landscapes, was deployed at kilometer scales to characterize critical zone structure. Concentrations of chaotic ground penetrating radar (GPR) reflections and diffraction hyperbolas with low electrical conductivity were observed in vertical zones that outcrop at the land surface as areas of intense fracturing and spheroidally weathered corestones. The width of these fractured and weathered zones showed an increase with proximity to the knickpoint, and was attributed to dilation of these sub-vertical fractures near the knickpoint, as postulated theoretically by a stress model calculated for the topographic variability across the knickzone in the Rio Icacos, and that shows a release of compressive stress near the knickpoint. We hypothesize that erosion rates increase in the knickzone because of this inferred dilation of fractures. Specifically, opened fractures could enhance access of water and in turn promote spalling, erosion, and spheroidal weathering. This study shows that ground-based hydrogeophysical methods used at the landscape-scale (traditionally applied at smaller scales) can be used to explore critical zone architecture at the scales needed to explain the extreme variability in erosion rates across river profiles. © 2018 John Wiley & Sons, Ltd.     

Postglacial to Holocene landscape evolution and process rates in steep alpine catchments

Climate change and high magnitude mass wasting events pose adverse societal effects and hazards, especially in alpine regions. Quantification of such geomorphic processes and their rates is therefore critical but is often hampered by the lack of appropriate techniques and the various spatiotemporal scales involved in these studies. Here we exploit both in situ cosmogenic beryllium-10 (¹⁰Be) and carbon-14 (¹⁴C) nuclide concentrations for deducing exposure ages and tracing of sediment through small alpine debris flow catchments in central Switzerland. The sediment cascade and modern processes we track from the source areas, through debris flow torrents to their final export out into sink regions with cosmogenic nuclides over an unprecedented five-year time series with seasonal resolution. Data from a seismic survey and a 90 m core revealed a glacially overdeepened basin, filled with glacial and paraglacial sediments. Surface exposure dating of fan boulders and radiocarbon ages constrain the valley fill from the last deglaciation until the Holocene and show that most of the fan existed in early Holocene times already. Current fan processes are controlled by episodic debris flow activity, snow (firn) and rock avalanches. Field investigations, digital elevation models (DEMs) of difference and geomorphic analysis agree with sediment fingerprinting with cosmogenic nuclides, highlighting that the bulk of material exported today at the outlet of the subcatchments derives from the lower fans. Cosmogenic nuclide concentrations steadily decrease from headwater sources to distal fan channels due to the incorporation of material with lower nuclide concentrations. Further downstream the admixture of sediment from catchments with less frequent debris flow activity can dilute the cosmogenic nuclide signals from debris flow dominated catchments but may also reach thresholds where buffering is limited. Consequently, careful assessment of boundary conditions and driving forces is required when apparent denudation rates derived from cosmogenic nuclide analysis are upscaled to larger regions. © 2018 John Wiley & Sons, Ltd.     

Multi‐method (14C, 36Cl, 234U/230Th) age bracketing of the Tschirgant rock avalanche (Eastern Alps): implications for absolute dating of catastrophic mass‐wasting

Correct and precise age determination of prehistorical catastrophic rock‐slope failures prerequisites any hypotheses relating this type of mass wasting to past climatic regimes or palaeo‐seismic records. Despite good exposure, easy accessibility and a long tradition of absolute dating, the age of the 230 million m³ carbonate‐lithic Tschirgant rock avalanche event of the Eastern Alps (Austria) still is relatively poorly constrained. We herein review the age of mass‐wasting based on a total of 17 absolute ages produced with three different methods (¹⁴C, ³⁶Cl, ²³⁴U/²³⁰Th). Chlorine‐36 (³⁶Cl) cosmogenic surface exposure dating of five boulders of the rock avalanche deposit indicates a mean event age of 3.06 ± 0.62 ka. Uranium‐234/thorium‐230 (²³⁴U/²³⁰Th) dating of soda‐straw stalactites formed in microcaves beneath boulders indicate mean precipitation ages of three individual soda straws at 3.20 ± 0.26 ka, 3.04 ± 0.10 ka and 2.81 ± 0.15 ka; notwithstanding potential internal errors, these ages provide an ‘older‐than’ (ante quam) proxy for mass‐wasting. Based on radiocarbon ages (nine sites) only, it was previously suggested that the present rock avalanche deposit represents two successive failures (3.75 ± 0.19 ka bp , 3.15 ± 0.19 ka bp ). There is, however, no evidence for two events neither in surface outcrops nor in LiDAR derived imagery and drill logs. The temporal distribution of all absolute ages (¹⁴C, ³⁶Cl, ²³⁴U/²³⁰Th) also does not necessarily indicate two successive events but suggest that a single catastrophic mass‐wasting took place between 3.4 and 2.4 ka bp . Taking into account the maximum age boundary given by reinterpreted radiocarbon datings and the minimum U/Th‐ages of calcite precipitations within the rock avalanche deposits, a most probable event age of 3.01 ± 0.10 ka bp can be proposed. Our results underscore the difficulty to accurately date catastrophic rock slope failures, but also the potential to increase the accuracy of age determination by combining methods. Copyright © 2016 John Wiley & Sons, Ltd.     

Quantifying glacial erosion on a limestone bed and the relevance for landscape development in the Alps

Glacial erosion is the basic process that has shaped the landscapes of the Alps. Despite intense research over centuries, and the use of various techniques, determination of glacial erosion rates remains challenging. This is not only because the location where the process occurs is almost inaccessible, but also because it is dependent on many different factors, including ice thickness and velocity, glacier thermal regime and lithology. Reported glacial erosion rates range over several orders of magnitude (0.01 to >10 mm a⁻¹). Most studies focus on crystalline bedrock, whereas few researchers have investigated glacial erosion on limestone. Here we analyse glacially polished bedrock surfaces at the recently deglaciated forefield of the Tsanfleuron glacier, Swiss Alps. The nearly horizontally bedded limestone hosts a well-developed karst system. Meltwater from the glacier drains into the subsurface within a few metres of the ice margin. By combining geomorphological mapping, measurement of cosmogenic ³⁶Cl concentrations of glacially eroded bedrock surfaces and a numerical model (MECED), we quantify at each sample location the amount of rock removed during glacier occupation. The glacial erosion rates calculated from these values range from 0 to 0.08 mm a⁻¹. These are orders of magnitude lower than values measured at comparable sites on crystalline bedrock. The high ³⁶Cl concentrations we measured show that the Tsanfleuron glacier was unable to effectively erode the gently dipping, strongly karstified limestone. We suggest that this effect may play a key role in formation and preservation over many glacial cycles of high-elevation, low-relief limestone plateaus in the Alps. © 2020 John Wiley & Sons Ltd.     

Laboratory and theoretical evaluation of impact of packing density,particle shape,and uniformity coefficient on erodibility of coarse-grained soil particles

Sedimentation – including erosion, transport, and deposition of coarse-grained particles – is a primary and growing environmental, engineering, and agricultural issue around the world. Soil erosion occurs when the hydrodynamic force induced by flowing water exceeds the geotechnical resistance of soils, as measured by critical shear stress for initiation of soil-particle motion. Even though various quantitative methods have been suggested with respect to different types of soil, the most widely accepted formula to estimate critical shear stress for coarse-grained soil is a direct function of the median grain size of the soil particles; however, the erosion resistance of soils also varies with other geotechnical properties, such as packing density, particle shape, and uniformity coefficient. Thus, in this study, a combined rolling–lift model for particle detachment was derived based on theoretical analysis. A series of experimental flume tests were conducted with specimens prepared with standard soil types, as well as laboratory-prepared mixtures of coarse-grained soil to validate the theoretical model and determine the effect of other geotechnical properties on the erosion characteristics of coarse grains, coupled with the effect of median particle size. The results indicated that the median grain size is the primary variable determining the resistance of coarse grains, but the critical shear stress also varies with the packing density of the soil matrix. In addition, angular particles show more erosion resistance than rounded particles, and the erosion potential of a soil decreased when the grain is well graded (higher value of uniformity coefficient). Additionally, regression analysis was performed to quantify the effect of each parameter on the critical shear stress of coarse grains. © 2020 John Wiley & Sons, Ltd.     

Bottom water circulation in Cascadia Basin

A combination of beta spiral and minimum length inverse methods, along with a compilation of historical and recent high-resolution CTD data, are used to produce a quantitative estimate of the subthermocline circulation in Cascadia Basin. Flow in the North Pacific Deep Water, from 900-1900 m, is characterized by a basin-scale anticyclonic gyre. Below 2000 m, two water masses are present within the basin interior, distinguished by different potential temperature-salinity lines. These water masses, referred to as Cascadia Basin Bottom Water (CBBW) and Cascadia Basin Deep Water (CBDW), are separated by a transition zone at about 2400 m depth. Below the depth where it freely communicates with the broader North Pacific, Cascadia Basin is renewed by northward flow through deep gaps in the Blanco Fracture Zone that feeds the lower limb of a vertical circulation cell within the CBBW. Lower CBBW gradually warms and returns to the south at lighter density. Isopycnal layer renewal times, based on combined lateral and diapycnal advective fluxes, increase upwards from the bottom. The densest layer, existing in the southeast quadrant of the basin below 2850 m, has an advective flushing time of 0.6 years. The total volume flushing time for the entire CBBW is 2.4 years, corresponding to an average water parcel residence time of 4.7 years. Geothermal heating at the Cascadia Basin seafloor produces a characteristic bottom-intensified temperature anomaly and plays an important role in the conversion of cold bottom water to lighter density within the CBBW. Although covering only about 0.05% of the global seafloor, the combined effects of bottom heat flux and diapycnal mixing within Cascadia Basin provide about 2-3% of the total required global input to the upward branch of the global thermohaline circulation.