Impacts of saturation-dependent anisotropy on the shrinkage behavior of clay rocks

Geomaterials such as soils and rocks can exhibit inherent anisotropy due to the preferred orientation of mineral grains and/or cracks. They can also be partially saturated with multiple types of fluids occupying the pore space. The anisotropic and unsaturated behaviors of geomaterials can be highly interdependent. Experimental studies have shown that the elastic parameters of rocks evolve with saturation. The effect of saturation has also been shown to differ between directions in transversely isotropic clay rock. This gives rise to saturation-dependent stiffness anisotropy. Similarly, permeability anisotropy can also be saturation-dependent. In this study, constitutive equations accommodating saturation-dependent stiffness and hydraulic anisotropy are presented. A linear function is used to describe the relationship between the elastic parameters and saturation, while the relative permeability–saturation relationship is characterized with a log-linear function. These equations are implemented into a hydromechanical framework to investigate the effects of saturation-dependent properties on the shrinkage behavior of clay rocks. Numerical simulations are presented to demonstrate the role of saturation-dependent stiffness and hydraulic anisotropy in shrinkage behavior. The results highlight that strain anisotropy and time evolution of pore pressures are substantially influenced by saturation-dependent stiffness and hydraulic anisotropy.

     

Evolution of anisotropy with saturation and its implications for the elastoplastic responses of clay rocks

Many clay rocks have distinct bedding planes. Experimental studies have shown that their mechanical properties evolve with the degree of saturation (DOS), often with higher stiffness and strength after drying. For transversely isotropic rocks, the effects of saturation can differ between the bed-normal (BN) and bed-parallel (BP) directions, which gives rise to saturation-dependent stiffness and strength anisotropy. Accurate prediction of the mechanical behavior of clay rocks under partially saturated conditions requires numerical models that can capture the evolving elastic and plastic anisotropy with DOS. In this study, we present an anisotropy framework for coupled solid deformation-fluid flow in unsaturated elastoplastic media. We incorporate saturation-dependent strength anisotropy into an anisotropic modified Cam-Clay (MCC) model and consider the evolving anisotropy in both the elastic and plastic responses. The model was calibrated using experimental data from triaxial tests to demonstrate its capability in capturing strength anisotropy at various levels of saturation. Through numerical simulations, we demonstrate the role of evolving stiffness and strength anisotropy in the mechanical behavior of clay rocks. Plane strain simulations of triaxial compression tests were also conducted to demonstrate the impacts of material anisotropy and DOS on the mechanical and fluid flow responses.     

Multitemporal terrestrial laser scanning point clouds for thaw subsidence observation at Arctic permafrost monitoring sites

This paper investigates different methods for quantifying thaw subsidence using terrestrial laser scanning (TLS) point clouds. Thaw subsidence is a slow (millimetre to centimetre per year) vertical displacement of the ground surface common in ice-rich permafrost-underlain landscapes. It is difficult to quantify thaw subsidence in tundra areas as they often lack stable reference frames. Also, there is no solid ground surface to serve as a basis for elevation measurements, due to a continuous moss–lichen cover. We investigate how an expert-driven method improves the accuracy of benchmark measurements at discrete locations within two sites using multitemporal TLS data of a 1-year period. Our method aggregates multiple experts’ determination of the ground surface in 3D point clouds, collected in a web-based tool. We then compare this to the performance of a fully automated ground surface determination method. Lastly, we quantify ground surface displacement by directly computing multitemporal point cloud distances, thereby extending thaw subsidence observation to an area-based assessment. Using the expert-driven quantification as reference, we validate the other methods, including in-situ benchmark measurements from a conventional field survey. This study demonstrates that quantifying the ground surface using 3D point clouds is more accurate than the field survey method. The expert-driven method achieves an accuracy of 0.1 ± 0.1 cm. Compared to this, in-situ benchmark measurements by single surveyors yield an accuracy of 0.4 ± 1.5 cm. This difference between the two methods is important, considering an observed displacement of 1.4 cm at the sites. Thaw subsidence quantification with the fully automatic benchmark-based method achieves an accuracy of 0.2 ± 0.5 cm and direct point cloud distance computation an accuracy of 0.2 ± 0.9 cm. The range in accuracy is largely influenced by properties of vegetation structure at locations within the sites. The developed methods enable a link of automated quantification and expert judgement for transparent long-term monitoring of permafrost subsidence. © 2020 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd     

Modeling IR spectral energy distributions: a?pilot study of?starburst parameters and silicate absorption curves for some GOALS galaxies

This paper describes a pilot study into the spectral energy distribution (SED) fitting and the derivation of physical parameters for 19 galaxies observed as part of the Great Observatories All-sky LIRG Survey (GOALS) survey as observed with the Spitzer Space Telescope. For this we have used the pan-spectral fitting tools developed in a series of papers by Dopita and his co-workers. We show that the standard Lee and Draine ??astronomical silicate?? model cannot provide a good fit to the silicate absorption features as observed in the heavily dust-extinguished (A _V??50 mag.) starbursts. We have derived an empirical fit to the ??starburst silicate?? absorption in these objects. This absorption curve is consistent with the silicate grains being systematically larger in starburst environments than in the local Galactic interstellar medium. We demonstrate the sensitivity of the SED fitting to each of the fitted parameters, and derive these parameters for those galaxies which do not have an embedded AGN. This technique is simple and provides reasonably robust and uniform parameters for the starburst, especially as far as the star formation rate, population of old stars, compactness of the starburst region and total foreground extinction are concerned. However, the chemical abundances and the optical extinction cannot be reliably determined by this analysis, and optical SEDs will also be required to provide a complete characterization of the starburst region and of the surrounding galaxy.     

Modeling Ocean Currents Through Complex Random Fields Indexed in Time

Surface ocean currents are often of interest in environmental monitoring. These vectorial data can be reasonably treated as a finite realization of a complex-valued random field, where the decomposition in modulus (current speed) and direction (current direction) of the current field is natural. Moreover, when observations are also available for different time points (other than at several locations), it is useful to evaluate the evolution of their complex correlation over time (rather than in space) and the corresponding modeling which is required for estimation purposes. This paper illustrates a first approach where the temporal profile of surface ocean currents is considered. After introducing the fundamental aspects of the complex formalism of a random field indexed in time, a new class of models suitable for including the temporal component is proposed and applied to describe the time-varying complex covariance function of current data. The analysis concerns ocean current observations, taken hourly on 30 April 2016 through high frequency radar systems at some stations located in the Northeastern Caribbean Sea. The selected complex covariance model indexed in time is used for estimation purposes and its reliability is confirmed by a numerical analysis.

     

Paleoenvironmental dynamics in central-eastern Brazil during the last 23 000 years: tropical peatland record in the Cerrado biome

The Cerrado biome is the second largest in Brazil, but the evolution of the Cerrado during the late Quaternary is not yet fully known. This study identifies paleoenvironmental changes during the last 23 000 years, based on a tropical mountain peatland record, in the Serra do Espinhaço Meridional in central-eastern Brazil. A multi-proxy approach was used that involved palynological analysis, stable isotopes (δ¹³C, δ¹⁵N), geochemistry, radiocarbon dating and multivariate statistics derived from a peatland core from Rio Preto (Minas Gerais state). The study reveals a very humid and cold climate during the late Pleistocene, with an increase in temperature and decrease in humidity at the Pleistocene–Holocene transition. During this period there was strong instability in the landscape (episodes of erosion). At the beginning of the Holocene there was a reduction in humidity with greater landscape stability. The current sub-humid climatic conditions seem to have been established in the mid-/late Holocene, with periods of landscape instability. Our findings agree with other Cerrado records that contradict previously established hypotheses, such as the Amazonian Refuge and the Pleistocene Arc.     

Relationships among water column toxins,cell abundance and chlorophyll concentrations during Karenia brevis blooms

The assumptions that Karenia brevis cell abundance and brevetoxin concentrations are proportional and that cell abundance and chlorophyll are related were tested in a 3-year field study off the west coast of Florida. The relationship between K. brevis cell abundance and brevetoxins (PbTx-2+PbTx-3) in whole water samples was strong (R²=0.92). There was no significant difference between the brevetoxin concentrations in whole water and the >0.7 μm particulate fraction. Only 7% of the total brevetoxin concentration was measured in the <0.7 μm (cell free) filtrate. The relationship of K. brevis cell abundance >5000 cells L⁻¹ with chlorophyll for all cruises and at all depths was robust (R²=0.78). These data substantiate the use of chlorophyll as a proxy for K. brevis cell abundance and K. brevis cell abundance as a proxy for brevetoxins during blooms. The ratios of the brevetoxins, PbTx-2:PbTx-3, was significantly higher in surface water than in bottom water. This information in conjunction with K. brevis growth rates may provide a useful indicator for determining the physiological state of the bloom over time.     

Staying cool in a changing climate: Reaching vulnerable populations during heat events

The frequency and intensity of hot weather events are expected to increase globally, threatening human health, especially among the elderly, poor, and chronically ill. Current literature indicates that emergency preparedness plans, heat health warning systems, and related interventions may not be reaching or supporting behavior change among those most vulnerable in heat events. Using a qualitative multiple case study design, we comprehensively examined practices of these populations to stay cool during hot weather (“cooling behaviors”) in four U.S. cities with documented racial/ethnic and socio-economic disparities and diverse heat preparedness strategies: Phoenix, Arizona; Detroit, Michigan; New York City, New York; and Philadelphia, Pennsylvania. Based on semi-structured in-depth interviews we conducted with 173 community members and organizational leaders during 2009–2010, we assessed why vulnerable populations do or do not participate in health-promoting behaviors at home or in their community during heat events, inquiring about perceptions of heat-related threats and vulnerability and the role of social support. While vulnerable populations often recognize heat's potential health threats, many overlook or disassociate from risk factors or rely on experiences living in or visiting warmer climates as a protective factor. Many adopt basic cooling behaviors, but unknowingly harmful behaviors such as improper use of fans and heating and cooling systems are also adopted. Decision-making related to commonly promoted behaviors such as air conditioner use and cooling center attendance is complex, and these resources are often inaccessible financially, physically, or culturally. Interviewees expressed how interpersonal, intergenerational relationships are generally but not always protective, where peer relationships are a valuable mechanism for facilitating cooling behaviors among the elderly during heat events. To prevent disparities in heat morbidity and mortality in an increasingly changing climate, we note the implications of local context, and we broadly inform heat preparedness plans, interventions, and messages by sharing the perspectives and words of community members representing vulnerable populations and leaders who work most closely with them.     

Deep sea diapirs and bottom simulating reflector in Fairway Basin (SW Pacific)

A recent swath-bathymetry and geophysical survey of the R/V L'Atalante in the Fairway Basin between Australia and New Caledonia allowed to confirm the Cretaceous age of the creation of the basin by continental stretching. This first stage of opening of the Fairway Basin is associated with the deposition of a continuous salt/mud layer feeding today numerous diapirs, some of them piercing the 3 to 4 km thick sedimentary cover and reaching the seafloor. In close link with this salt/mud layer a Bottom Simulating Reflector indicator of gas hydrates level occupies a 70000 km² surface at about 500 to 600 m-depth beneath the sea floor. The coexistence of both BSR and diapirs suggests a thermogenic better than biogenic origin for the gas hydrates horizon.     

Mechanisms of ringwoodite formation in shocked meteorites: Evidence from L5 chondrite Dhofar 1970

The formation of the high‐pressure compositional equivalents of olivine and pyroxene has been well‐documented within and surrounding shock‐induced veins in chondritic meteorites, formed by crystallization from a liquid‐ or solid‐state phase transformation. Typically polycrystalline ringwoodite grains have a narrow range of compositions that overlap with those of their olivine precursors, whereas the formation of iron‐enriched ringwoodite has been documented from only a handful of meteorites. Here, we report backscattered electron images, quantitative wavelength‐dispersive spectrometry (WDS) analyses, qualitative WDS elemental X‐ray maps, and micro‐Raman spectra that reveal the presence of Fe‐rich ringwoodite (Fa_44‐63) as fine‐grained (500 nm), polycrystalline rims on olivine (Fa_24‐25) wall rock and as clasts engulfed by shock melt in a previously unstudied L5 chondrite, Dhofar 1970. Crystallization of majorite + magnesiowüstite in the vein interior and metastable mineral assemblages within 35 μm of the vein margin attest to rapid crystallization of a superheated shock melt (>2300 K) from 20─25 GPa to ambient pressure and temperature. The texture and composition of bright polycrystalline ringwoodite rims (Fa_44‐63; MnO 0.01─0.08 wt%) surrounding dark polycrystalline olivine (Fa_8‐14; MnO 0.56─0.65 wt%) implies a solid‐state transformation mechanism in which Fe was preferentially partitioned to ringwoodite. The spatial association between ringwoodite and shock melt suggests that the rapidly fluctuating thermal regimes experienced by chondritic minerals in contact with shock melt are necessary to both drive phase transformation but also to prevent back‐transformation.