Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Extreme events”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 235 records · Page 13

Life Cycle of Major Sudden Stratospheric Warmings in the Southern Hemisphere from a Multimillennial GCM Simulation

Abstract Sudden stratospheric warmings (SSWs) are rare in the Southern Hemisphere (SH), making it difficult to study possible precursors or subsequent impacts. Using a multimillennial coupled climate model simulation producing 161 SSWs in the SH, we present a detailed study of their life cycle. We show that SH SSWs are predominantly displacement events forced by wave-1 planetary waves, and that a surface signature similar to the negative phase of the Southern Annular Mode (SAM) is detectable up to two months before the onset date, but there is a tendency for a transition from wave 1 before to zonally symmetric anomalies after onset. We identify a strong weakening of the Amundsen Sea low as one of the most prominent precursors, which weakens the climatological wave-2 and wave-3 stationary waves and strengthens wave-1 forcing. Compared to their northern counterparts, SH SSWs generally have a longer time scale, and while there is evidence of pre-onset forcing related to tropical sea surface temperatures, the Indian Ocean dipole is more important than El Niño–Southern Oscillation. Significance Statement Sudden stratospheric warmings (SSWs) are extreme events where the winter polar stratosphere warms within a few days to temperatures usually only experienced in summer. These events are rare in the Southern Hemisphere. Therefore, both the observational record and standard climate model simulations are not enough to understand how SSWs develop, or how they might change surface weather. Here we use very long global climate simulations that produce a large number of SSWs in the Southern Hemisphere to study the development and impact of these events. This includes possible precursors as well as the influence they have on surface weather after they occur.

Meteorology & Atmospheric Sciences↗

Data-driven analysis of relight variability of jet fuels induced by turbulence

For safety purposes, reliable reignition of aircraft engines in the event of flame blow-out is a critical requirement. Typically, an external ignition source in the form of a spark is used to achieve a stable flame in the combustor. However, such forced turbulent ignition may not always successfully relight the combustor, mainly because the state of the combustor cannot be precisely determined. Uncertainty in the turbulent flow inside the combustor, inflow conditions, and spark discharge characteristics can lead to variability in sparking outcomes even for nominally identical operating conditions. Prior studies have shown that of all the uncertain parameters, turbulence is often dominant and can drastically alter ignition behavior. For instance, even when different fuels have similar ignition delay times, their ignition behavior in practical systems can be completely different. In practical operating conditions, it is challenging to understand why ignition fails and how much variation in outcomes can be expected. The focus of this work is to understand relight variability induced by turbulence for two different aircraft fuels, namely Jet-A and a variant named C1. A detailed, previously developed simulation approach is used to generate a large number of successful and failed ignition events. Using this data, the cause of misfire is evaluated based on a discriminant analysis that delineates the difference between turbulent initial conditions that lead to ignition or failure. From the discriminant analysis, a compressed sensing algorithm is then applied to help pinpoint the locations of relevant turbulent features. Findings from the discriminant analysis are confirmed with the time history of near kernel properties. Next, a clustering strategy is used to identify ignition and misfire modes. With this approach, it was determined that the cause of ignition failure is different for the two fuels. While it was found that Jet-A is influenced by fuel entrainment, C1 was found to be more sensitive to small scale turbulence features. Finally, a larger variability is found in the ignition modes of C1, which can be subject to extreme events induced by kernel breakdown.

42 ENGINEERING↗

Stochastic Microgrid Scheduling With Chance‐Constrained Resilience Consideration

Traditionally, it is assumed that microgrids transition seamlessly from grid‐connected operation to islanded mode in the event of sudden main grid outages. In reality, the islanding process, especially unintentional islanding, is rarely seamless. Instead, it is subject to voltage and frequency fluctuations caused by the instantaneous disconnection of the point of common coupling (PCC) switch, variations in loads and renewable generation output and even the protection tripping of distributed energy resources (DERs). To mitigate these fluctuations and facilitate a smooth islanding process, we propose a stochastic microgrid scheduling model that incorporates chance‐constrained resilience measures. Specifically, the resilience measure is defined as the probability of successful islanding (PSI), that is, the probability that a microgrid can mitigate the generation‐demand imbalance caused by the disconnection of the PCC switch, variations in load and renewable generation and DER tripping. This measure is modelled using chance constraints. Unlike existing reliability and resilience indices, which typically neglect the possibility of microgrid/DER failure under extreme events and assume their survival while primarily focussing on reducing impact duration or magnitude, the proposed PSI‐based framework explicitly addresses microgrid and DER survival during the islanding transition. The formulated nonlinear chance constraints are approximated using a multiinterval approach and equivalently represented as a mixed‐integer linear programming (MILP) formulation. Case study results validate the proposed method, showing that the PSI estimation error is reduced to less than 8%, compared to approximately 28% with existing methods. Various sensitivity analyses on the DER tripping rate and PSI settings were performed to validate the robustness of the proposed method. In particular, the necessity of accounting for DER tripping in the PSI calculation was demonstrated.

chance constrained optimization↗

Continental United States climate projections based on thermodynamic modification of historical weather

Abstract Regional climate models can be used to examine how past weather events might unfold under different climate conditions by simulating analogue versions of those events with modified thermodynamic conditions (i.e., warming signals). Here, we apply this approach by dynamically downscaling a 40-year sequence of past weather from 1980–2019 driven by atmospheric re-analysis, and then repeating this 40-year sequence a total of 8 times using a range of time-evolving thermodynamic warming signals that follow 4 80-year future warming trajectories from 2020–2099. Warming signals follow two emission scenarios (SSP585 and SSP245) and are derived from two groups of global climate models based on whether they exhibit relatively high or low climate sensitivity. The resulting dataset, which contains 25 hourly and over 200 3-hourly variables at 12 km spatial resolution, can be used to examine a plausible range of future climate conditions in direct reference to previously observed weather and enables a systematic exploration of the ways in which thermodynamic change influences the characteristics of historical extreme events.

54 ENVIRONMENTAL SCIENCES↗

Collaborative Research: Advancing Arctic Climate Projection Capability at Seasonal to Decadal Scales (Final Technical Report)

The Regional Arctic System Model (RASM) at process resolving configurations has been used to (i) advance understanding of physical processes and feedbacks involved in Arctic amplification and (ii) understand and potentially reduce uncertainty in prediction of arctic climate change at seasonal to decadal scales. RASM consists the atmosphere (Weather and Research Forecasting model, WRF), ocean (Parallel Ocean Program, POP), sea ice (CICE), land hydrology (Variable Infiltration Capacity model, VIC), river routing scheme (RVIC), marine biogeochemistry components and the coupling framework (CPL7). Its domain is pan-Arctic, with the atmosphere and land components configured on a 50-km or 25-km grid and four configurations of the ocean and sea ice components: 1/12°(~9.3km) or 1/48°(~2.4km) and 45 or 60 vertical layers. These RASM configurations have been motivated by the emerging exascale capability for high performance computing to improve model fidelity. The dynamical downscaling of reanalysis allows comparison of RASM results with observations in place and time to: (i) advance system level understanding of physical processes and coupling involved in an event, (ii) optimize model parameter space, (iii) diagnose and reduce model biases and (iv) produce realistic and consistent across all the components initial conditions for predictions and predictability studies, which are all unique capabilities not available in global Earth System Models (ESMs). An evaluation of RASM 1.0 (Cassano et al. 2017) revealed that it had a cold bias over the oceans and a warm bias over land areas due largely to cloud and radiation biases in the model, with too little cloud cover simulated over land and too much cloud cover simulated over sub-polar oceans. This study has motivated an upgrade to WRF version 3.7.1 in RASM and allowed for the inclusion of the radiative impact of convective clouds. A variety of atmospheric physics parameterizations were evaluated against observations (e.g. data from the Arctic Clouds in Summer Experiment (ACSE); Sedlar et al. 2020) to identify an optimal suite of WRF physics options in RASM. The RASM with the optimized WRF physics were used to study the impact of strong mesoscale winds over the ocean around the southern tip of Greenland (DuVivier and Cassano 2016) and their impact on oceanic convection (DuVivier et al. 2017a). Data from the PolarWinds field campaign were used to evaluate WRF boundary layer physics and resolution impacts on the simulation of a Greenland barrier wind event (DuVivier et al. 2017b). The RVIC streamflow routing model has been implemented in RASM to realistically represent high-resolution streamflow processes (Hamman et al. 2017) and to couple the land buoyancy fluxes to the ocean. The RASM-RVIC high-resolution data set of all coastal freshwater fluxes for the Arctic drainage basin and surrounding areas for 1979-2014 was published as a separate product (https://doi.org/10.5281/zenodo.293037). The fidelity of atmospheric momentum transfer to and the response of polar marine Ekman layer in RASM and Community Earth System Model (CESM) was investigated by Roberts et al. (2015). The increased frequency of oceanic flux exchange in CESM, following the RASM guidance, caused a considerable increase in the median inertial ice speed across the Southern Ocean and parts of the Arctic. A comprehensive evaluation of the RASM1.0 atmosphere-ocean-sea ice-land interface was completed by Brunke et al. (2018). RASM was also demonstrated for its capability to simulate extreme events in agreement with observations in space and time (Lee et al. submitted). In particular, the development of three open water events, known as polynyas, have been simulated north of Greenland in February of 2011, 2017 and 2018, in agreement with satellite observations for the past four decades. The optimized RASM sea ice results have been favorably evaluated against satellite observations and a subset of eleven CMIP6 models (Watts et al. submitted). In a complementary project, Jin et al. (2018) have shown that RASM with higher-resolution and new sea-ice processes contributed to lower model errors in sea-ice conditions, concentrations of nutrients and ice algae, in comparison to results from the coarse-resolution (1°) CESM. In two other complementary studies, RASM results were used (i) to explain areas of concentrated use by bowhead whales, the seasonal progression in the use, and the physical environment within those areas (Citta et al. 2015) and (ii) for a synthesis of fall bowhead whales distribution and migration in the Bering-Chukchi-Beaufort (BCB) Sea to investigate whale movements and feeding to the local ocean hydrography and currents (Citta et al. 2018). However, the multi-decadal output from the CESM Large Ensemble yielded unrealistic forcing. Instead, the shorter NCEP CFSv2 9-month forecasts were successfully tested and afforded an increased ensemble size (~30) to demonstrate gains of dynamical downscaling at sub-seasonal to intra-annual time scales. The improved model physics and coupling among RASM model components have yielded more realistic representation of the sea ice cover and consistent across all model components initial conditions. Consequently, RASM demonstrates significant gains compared to simulation of sea ice in the NCEP reanalysis. In addition, RASM 6-month ensemble forecasts yield very realistic sea ice distribution, which demonstrates both significant gains of dynamical downscaling and the continued impact of the initial conditions on forecasts out to 6 months (https://nps.edu/web/rasm/predictions). A follow up study using RASM for dynamical downscaling of the more realistic CESM initialized Decadal Prediction Large Ensemble output is currently ongoing as part of the DOE RGMA HiLAT-RASM project.

54 ENVIRONMENTAL SCIENCES↗

Elevated Forest Canopy Loss After Wildfires in Moist and Cool Forests in the Pacific Northwest

Wildfire frequency, intensity, and rate of spread are increasing across the Western U.S, resulting in more severe ecosystem impacts. Significant tree mortality can occur years after fire events, but this has received little attention compared to the immediate tree loss during a fire. We overlapped forest cover loss data with burn severity maps in the U.S. Pacific Northwest and quantified the total and delayed forest canopy loss after fires. We found that wildfires resulted in total canopy loss fraction (CLF) of 84%, 53%, and 22% within 3 years in areas burned at high, moderate, and low severity, respectively. The delayed canopy loss accounted for approximately 1/3, 1/2, and 2/3 of the total canopy loss for high, moderate, and low severity burns. Delayed canopy loss was greater in moist and cool areas than in dry and warm areas, likely because tree species in wetter environments were less adapted to survive when fires did occur. Across all forests, delayed CLF doubled as temperature increased from the climatological mean to a hot anomaly and tripled as vapor pressure deficit increased from a wet anomaly to a dry anomaly. Fire impacts on forest ecosystems are likely to intensify under future climate scenarios as wildfires expand into areas that historically experienced infrequent fires. The impacts can also be exacerbated by more frequent compound extreme events, such as droughts, heatwaves, and fires. These findings highlight the urgent need for targeted forest management strategies, particularly in mesic forests, to mitigate future fire impacts.

54 ENVIRONMENTAL SCIENCES↗

Process-oriented evaluation of quasi-stationary Rossby waves and their impact on surface air temperature extremes in dynamical downscaling over North America

Quasi-stationary Rossby waves are a crucial component of the general circulation and play a significant role in regional water and energy cycles, as well as in extreme events. However, process-oriented evaluation for Rossby waves is rarely performed for dynamical downscaling simulations. To close this gap, we evaluate three classes of dynamical downscaling approaches, with a focus on quasi-stationary Rossby waves and their impact on surface air temperature over North America during Northern Hemisphere summer. The three classes of models differ in the way large-scale forcing is provided: a limited-area model (LAM) constrained only by lateral boundary conditions, represented by RegCM4 from the North American branch of the Coordinated Regional Downscaling Experiment (NA-CORDEX), a LAM with spectral nudging to maintain consistency in large-scale dynamics with the forcing data, represented by the Weather Research and Forecasting (WRF) model simulation in NA-CORDEX, and a global variable-resolution model with smoothly varying grid spacings, represented by the Community Atmosphere Model version 5.4, with the Model for Prediction Across Scales (MPAS) as its dynamical core (CAM-MPAS). With no constraints on the atmospheric dynamics, CAM-MPAS exhibits several mean biases in the upper-level circulations over the Pacific Coast region: a weaker subtropical jet, a northward-shifted mid-latitude jet, and an overestimated southerly flow. With the lateral boundary constraint alone, RegCM4 also exhibits weaker jets and overestimated southerly winds off the West Coast. Rossby ray theory reveals that those wind biases direct incoming Rossby waves northward. The erroneously routed Rossby waves distort the relationship between the accumulation of wave activity over the US West Coast and surface temperature anomalies over the Southern Great Plains, which emerges approximately 4 d after the convergence of wave-activity flux in the ERA-Interim reanalysis. Furthermore, the response of heatwaves to the extreme wave activity flux is not reproduced by the two models, a serious drawback as a dynamical downscaling framework is expected to connect large-scale forcing to local-scale phenomena. The WRF model employing spectral nudging is largely free from the aforementioned problems. A pair of sensitivity simulations suggests that spectral nudging is the key to improving the dynamics of quasi-stationary Rossby waves and their impact on surface air temperature. Our results also demonstrate the effectiveness of Rossby wave diagnostics that allow for realistic background flows for assessing the credibility of dynamical downscaling over North America, where incoming Rossby waves propagate through complex circulation patterns before traveling across the continent.

54 ENVIRONMENTAL SCIENCES↗

Can machine learning accelerate process understanding and decision‐relevant predictions of river water quality?

Abstract The global decline of water quality in rivers and streams has resulted in a pressing need to design new watershed management strategies. Water quality can be affected by multiple stressors including population growth, land use change, global warming, and extreme events, with repercussions on human and ecosystem health. A scientific understanding of factors affecting riverine water quality and predictions at local to regional scales, and at sub‐daily to decadal timescales are needed for optimal management of watersheds and river basins. Here, we discuss how machine learning (ML) can enable development of more accurate, computationally tractable, and scalable models for analysis and predictions of river water quality. We review relevant state‐of‐the art applications of ML for water quality models and discuss opportunities to improve the use of ML with emerging computational and mathematical methods for model selection, hyperparameter optimization, incorporating process knowledge into ML models, improving explainablity, uncertainty quantification, and model‐data integration. We then present considerations for using ML to address water quality problems given their scale and complexity, available data and computational resources, and stakeholder needs. When combined with decades of process understanding, interdisciplinary advances in knowledge‐guided ML, information theory, data integration, and analytics can help address fundamental science questions and enable decision‐relevant predictions of riverine water quality.

54 ENVIRONMENTAL SCIENCES↗

Advancing an integrated understanding of land–ocean connections in shaping the marine ecosystems of coastal temperate rainforest ecoregions

Land and ocean ecosystems are strongly connected and mutually interactive. As climate changes and other anthropogenic stressors intensify, the complex pathways that link these systems will strengthen or weaken in ways that are currently beyond reliable prediction. In this review we offer a framework of land–ocean couplings and their role in shaping marine ecosystems in coastal temperate rainforest (CTR) ecoregions, where high freshwater and materials flux result in particularly strong land–ocean connections. Using the largest contiguous expanse of CTR on Earth—the Northeast Pacific CTR (NPCTR)—as a case study, we integrate current understanding of the spatial and temporal scales of interacting processes across the land–ocean continuum, and examine how these processes structure and are defining features of marine ecosystems from nearshore to offshore domains. We look ahead to the potential effects of climate and other anthropogenic changes on the coupled land–ocean meta-ecosystem. Finally, we review key data gaps and provide research recommendations for an integrated, transdisciplinary approach with the intent to guide future evaluations of and management recommendations for ongoing impacts to marine ecosystems of the NPCTR and other CTRs globally. In the light of extreme events including heatwaves, fire, and flooding, which are occurring almost annually, this integrative agenda is not only necessary but urgent.

54 ENVIRONMENTAL SCIENCES↗

Role of Renewables in Grid Resilience

This chapter focuses on revealing the impact of renewable energies, especially distributed energy resources, on modern power system resilience. It is widely recognized that the increasing adoption of renewable energy is bringing more uncertainty and variability, posing new challenges to grid resilience. However, both supply and demand are becoming more active in grid resilience management. The communities themselves would benefit from the utilization of local resources and flexible management to gain resilience during the extreme events.

grid resilience↗

Projected changes in the terrestrial and oceanic regulators of climate variability across sub-Saharan Africa

Future changes in the sign and intensity of ocean–land–atmosphere interactions have been insufficiently studied, despite implications for regional climate change projections, extreme event statistics, and seasonal climate predictability. In response to this deficiency, the present study focuses on projected responses to the enhanced greenhouse effect in: (1) the mean state of the atmosphere and land surface; (2) oceanic and terrestrial drivers of sub-Saharan climate variability; and (3) total seasonal climate predictability of sub-Saharan Africa, a region known for its pronounced land–atmosphere coupling. Analysis focuses on output from 23 Earth System Models in the Coupled Model Intercomparison Project Phase Five for the late twentieth and twenty-first centuries. It is projected that the greatest warming across sub-Saharan Africa will occur over the Sahel, the monsoon season will become more persistent into late summer and autumn, short rains in the Horn of Africa (HOA) will intensify, and leaf area index will increase across the HOA. Stepwise Generalized Equilibrium Feedback Assessment, i.e. a multivariate statistical approach, is applied to the model output over sub-Saharan Africa in order to explore the oceanic and terrestrial drivers of regional climate. The models indicate that the study region’s climate variability is dominated by oceanic drivers, with secondary contributions from soil moisture and very modest impacts from vegetation. Overall, the general model consensus of future projections indicates a concerning diminished seasonal predictability of sub-Saharan African regional climate based on key oceanic and terrestrial predictors and an elevated role of the land surface (associated with soil moisture anomalies) compared to oceanic drivers in regulating regional climate variability.

54 ENVIRONMENTAL SCIENCES↗

Dynamics and drivers of land use and land cover changes in Bangladesh

Abstract Bangladesh has undergone dramatic land use and land cover changes (LULCC) in recent years, but no quantitative analysis of LULCC drivers at the national scale exists so far. Here, we quantified the drivers of major LULCC in combination with biophysical and socioeconomic observations at the sub-district level. We used Landsat satellite data to interpret LULCC from 2000 to 2010 and employed a Global Surface Water Dataset to account for the influences of water seasonality. The results suggest that major LULCC in Bangladesh occur between agricultural land and waterbodies and between forest and shrubland. Exclusion of seasonal waterbodies can improve the accuracy of our LULCC results and driver analysis. Although the gross gain and loss of agricultural land are large on the local scale, the net change (gross gain minus gross loss) at a country scale is almost negligible. Climate dynamics and extreme events and changes in urban and rural households were driving the changes from forest to shrubland in the southeast region. The conversion from agricultural land to standing waterbodies in the southwest region was mainly driven by urban household dynamics, population growth, distance to cities and major roads, and precipitation dynamics. This study, which is the first effort accounting for water seasonality and quantifying biophysical and socioeconomic drivers of LULCC at the national scale, provides a perspective on overall LULCC and underlying drivers over a decadal time scale and national spatial scale and can serve as a scientific basis for developing land policies in Bangladesh.

54 ENVIRONMENTAL SCIENCES↗

A treatment-effect model to quantify human dimensions of disaster impacts: the case of Hurricane Maria in Puerto Rico

Here, we propose a supervised learning approach using publicly available panel data to statistically quantify the specific manifestations of human impacts of an extreme event, such as changes number of suicides, substance abuse, excess mortality, and unemployment. This allows us to conceptually focus our framework on human impacts and how by attributing them to disaster events along widely accepted psychological, economic, and social dimensions. Our modified treatment-effect model allows counterfactual baseline conditions to be posited for each manifestation from which an aggregated quantitative multi-faceted measure of human impacts can be determined. The developed statistical methodology could be beneficial to policymakers who must allocate scarce resources to those communities in greater need. We illustrate the applicability of our approach using annual and monthly panel data from 2012 to 2018 encompassing the 2017 Hurricane Maria event across various municipalities in Puerto Rico. Our statistical modeling methodology stands apart since (i) it explicitly and more realistically captures the effect of different human-oriented manifestations of an actual event and (ii) it is flexible enough to accommodate individual preferences of various stakeholders in how they assign importance to multiple manifestations of human impacts.

54 ENVIRONMENTAL SCIENCES↗

A reforecasting-based dynamic reserve estimation for variable renewable generation and demand uncertainty

The installed capacity of renewables-based energy sources has been increasing in traditional power systems. In order to accommodate the increased variability and uncertainty associated with the deeper penetration of renewable sources like solar and wind, adjusted amounts of dynamic reserve are needed. Although probabilistic dynamic reserve estimation methods have been previously developed, most of them consider the uncertainty to be represented by parametric density functions that tend to perform poorly under extreme events and, moreover, neglect uncertainty introduced by the forecasting model itself. Toward addressing these limitations, this work presents, for the first time, a dynamic reserve estimation method for flexibility that incorporates nonparametric density estimation and a machine learning based reforecasting to provide a day-ahead prediction of the mean and spread of uncertainty around the base forecast. The prediction is, in turn, used to estimate the up and down reserve relative to the base forecast. Here, the present method takes various endogenous and exogenous features, including the calendar variables, as input to estimate the day-ahead reserve. Using a combination of reforecasting and dynamic reserve estimation techniques, the method is shown to adjust better to the dynamic nature of reserve requirements providing only what is needed to accommodate the expected deviations. Considering California Independent System Operator (CAISO) solar, wind and load data over an 18 month period, up to 67% reduction in the amount of reserve capacity needed for a one day reserve and reserve penalty for solar uncertainty is demonstrated. Additionally, the risk of reserve insufficiency in meeting the net demand is reduced by 20% with the proposed method.

14 SOLAR ENERGY↗

Mercury linked to Deccan Traps volcanism, climate change and the end-Cretaceous mass extinction

Mercury (Hg) anomalies linked to Large Igneous Provinces (LIP) volcanism have been identified in sediments across all five major mass extinctions in Earth's history. This study tests whether Hg in marine sediments is a reliable proxy linking Deccan Traps volcanic eruptions to late Maastrichtian global climate warming and the mass extinction at the Cretaceous-Paleogene boundary (KPB). Our primary test site is the Elles section in Tunisia, the auxiliary Global Stratotype Section and Point (GSSP) to El Kef. Elles has the most complete marine sedimentary record and a high average sedimentation rate of ~4.7 cm/ky. Here, we chose the Hor Hahar section in Israel to corroborate the geographic distribution of Hg fallout from Deccan volcanism. Reliability of the Hg proxy over the last 550 ky of the Maastrichtian to early Danian was evaluated based on high-resolution age control (orbital cyclostratigraphy), stable isotope climate record, Hg concentrations, biotic turnover and mass extinction. These results were correlated with the pulsed Deccan eruptive history constrained previously by U-Pb zircon geochronology. Our results support Hg as robust proxy for Deccan volcanism with large Hg spikes marking “extreme event” (EE) pulsed eruptions correlative with climate warming peaks separated by steady, less intense eruptions. Long-term global climate warming began near ~350 ky pre-KPB, reached maximum warming (3–4 °C) between 285 and 200 ky pre-KPB, followed by gradual cooling and rapid temperature drop between 45 and 25 ky pre-KPB. During the last 25 ky before the KPB, multiple Hg EE eruptions correlate with hyperthermal warming that culminated in the rapid mass extinction at Elles during ≤1000 years of the Cretaceous. These latest Cretaceous Hg peaks may correlate with massive, distal, Deccan-sourced lava flows (>1000 km long) that traversed the Indian subcontinent and flowed into the Bay of Bengal, bracketing the mass extinction. These results support Deccan volcanism as a primary driver of the end-Cretaceous mass extinction.

54 ENVIRONMENTAL SCIENCES↗

Guest Editorial: Special issue on flexible and resilient urban energy systems

This guest editorial summarizes the topics and the papers selected for the Special Issue on Flexible and Resilient Urban Energy Systems. After rigorous reviewing process, 21 papers are accepted for publication. These 21 accepted papers cover various aspects of urban energy systems and are distributed as following: situational awareness of urban energy systems (2 papers), quantification metrics of flexibility and resilience of urban energy systems (3 papers), vulnerability modeling of urban energy systems under various extreme events (3 papers), planning of flexible and resilient urban energy systems (4 papers), robust and resilient operation and control of urban energy systems (4 papers), recovery and restoration strategy of urban energy systems (2 papers), and coordination and interoperability of interconnected energy systems (3 papers). The Guest Editorial Board hopes this Special Issue can provide a valuable information for future research and advancements in the field of flexible and resilient urban energy systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Statistical evaluation of microscale stress conditions leading to void nucleation in the weak shock regime

Here, we investigate the heterogeneity of the stress state driven by anisotropic deformation response at the single crystal level through five statistical volume element (SVE) calculations of polycrystalline BCC tantalum. This work focuses on grain boundaries as a prominent material defect type prone to void nucleation based upon experimental observations of predominantly intergranular void nucleation in this material. The SVEs are constructed to be statistically representative of larger volumes of material and are meshed such that mean and standard deviation of grain size and orientation information is reconstructed. The computational meshes feature hexahedral (brick) elements and smooth conformal grain boundaries where significant stress concentration is known to occur, a tail effect of interest in the extreme events process of dynamic ductile damage. An existing micromechanical crystallographic plasticity model shown to capture the single crystal behavior of BCC tantalum well is used to perform the polycrystal calculations. The model includes representation of the non-Schmid effect of non-planar screw dislocation kinetics in tantalum. A three-dimensional stress state time profile predicted by damage modeling of a flyer plate impact experiment is applied as boundary conditions to each SVE. Resulting grain boundary stress state statistics are strongly non-Gaussian. Significant structural evolution is observed within the compressive hold before unloading into tension in the stress profile. Strong angular dependence of grain boundary traction magnitude with shock direction is observed. Non-Schmid effects continue to suggest their influence on propensity of microstructural defect types to nucleate voids. A general void nucleation criterion is proposed using probability theory. The general framework is specified to polycrystalline BCC tantalum in the weak shock regime to include the SVE calculations and literature molecular dynamics calculations of grain boundary void nucleation strength. Probability density functions (PDFs) are used to describe the interaction between the local stress state heterogeneity and the distributed grain boundary void nucleation strength state. A causation entropy maximization procedure removes the requirement for ad hoc selection of a PDF functional form and provides a rigorous procedure for data-based PDF determination. The resulting physically informed PDF describes the spatial appearance frequency of nucleated voids as a function of applied macroscale pressure. Lower length scale physics are thus packaged in a precise and computationally efficient way to provide computational plasticity insight to macroscale dynamic ductile damage models.

36 MATERIALS SCIENCE↗

A hydrogeophysical framework to assess infiltration during a simulated ecosystem-scale flooding experiment

This study presents a framework to quantify changes in soil saturation in response to flooding caused by extreme hydrologic perturbation on coastal ecosystems at the interfaces and transition between terrestrial and aquatic systems. Subsurface heterogeneity limits the use of in situ measurements to quantify subsurface flow during flooding due to the spatial discontinuity in the measured data. While geophysical methods, including time-lapse electrical resistivity imaging (ERI), are increasingly used to monitor soil hydrological processes, their abilities to parameterize flow models have been underutilized. This study combines background ERI, ground penetrating radar (GPR), time-lapse ERI, soil characterization, and a numerical flow model developed using an Advanced Terrestrial Simulator (ATS) code to quantify the infiltration pathway and describe the hydrological dynamics during a simulated flooding experiment. We assessed the use of two conceptual models developed using [1] ERI and GPR data that described the stratigraphic distribution, and time-lapse ERI that mapped permeability contrast, and [2] information from a national soil database for capturing changes in saturation. Combining the ERI and GPR results with soil core data revealed the stratigraphic heterogeneity at the site with a silty clay layer from 1 to 2 m between an overlying loamy topsoil and an underlying saturated silty sand. This silty clay layer could restrict deep infiltration. The time-lapse ERI showed up to a 35% decrease in resistivity, which correlated with soil moisture data (R 2 value > 0.53) and revealed preferential infiltration zones used to inform the flow model. Numerical simulation results from both the geophysics- and soil database-informed models quantified changes in soil saturation with calculated soil moistures that agreed with field data. The geophysics-informed model captured more of the system’s variability, reflective of shallow subsurface heterogeneities. The framework presented will serve as a precursor for a robust ecohydrological model that can describe the impacts of extreme events induced by climate change on coastal ecosystems.

54 ENVIRONMENTAL SCIENCES↗