Stream water temperature is considered a “master variable” in environmental processes and human activities. Existing process-based models have difficulties with defining true equation parameters, and sometimes simplifications like assuming constant values influence the accuracy of results. Machine learning models are a highly successful tool for simulating stream temperature, but it is challenging to learn about processes and dynamics from their success. Here we integrate process-based modeling (SNTEMP model) and machine learning by building on a recently developed framework for parameter learning. With this framework, we used a deep neural network to map raw information (like catchment attributes and meteorological forcings) to parameters, and then inspected and fed the results into SNTEMP equations which we implemented in a deep learning platform. We trained the deep neural network across many basins in the conterminous United States in order to maximize the capturing of physical relationships and avoid overfitting. The presented framework has the ability of providing dynamic parameters based on the response of basins to meteorological conditions. The goal of this framework is to minimize the differences between stream temperature observations and SNTEMP outputs in the new platform. Parameter learning allows us to learn model parameters on large scales, providing benefits in efficiency, performance, and generalizability through applying global constraints. This method has also been shown to provide more physically-sensible parameters due to applying a global constraint. This model improves our understanding of how to parameterize the physical processes related to water temperature.

As the sea level rises, it is alarming that the threat from flooding induced by tropical cyclones would cause more severe damages to coastal regions worldwide. In order to address this threat, optimizing coastal protective or mitigation strategies is necessary, given limited resources. The optimization methodology must incorporate feedback from stakeholders for practical use. Multiple interviews were conducted by engineering model developers and social scientists with stakeholders who have first-hand knowledge and varied backgrounds in New York. The protective strategies have been tuned to the critical infrastructure's particular and contextual risks due to flood hazards by engaging and integrating stakeholders' knowledge on the interdependency of the infrastructures and other aspects after the first interview. The second interview was conducted for further model improvement.

Coastal regions are continuously under the threat of flooding induced by tropical cyclones worldwide. These threats continue to increase due to the effects of climate change such as sea-level rise. A number of available protective or mitigation strategies have been examined to address this threat and protect coastlines around the world. However, identifying the most effective strategy given limited resources is a complex question. Optimization methodologies as we have proposed integrate physical analysis and stakeholder feedback to come to a set of best mitigation strategies. This study examines physical and socio-economical aspects of flooding impacts to optimize these strategies. These are then examined including seawalls, elevated promenades, and strategic retreats.

Knowing the heterogeneous crustal structure is essential for understanding the ice dynamics, glacial isostatic adjustment (GIA) and tectonic history in Antarctica. For example, geothermal heat flux (GHF) is a major boundary condition for ice dynamics and the crust thickness and its composition (mafic or felsic) are important factors in GHF. Meanwhile, the GIA signal and its gravity response are essential for detecting mass-balance change and predicting future sea-level change. Errors in the density model used, which may be over 10%, will propagate into the gravity calculations. In this study, we use gravity inversion constrained by seismic depth estimation to recover the heterogeneous crustal structure of Antarctica, and estimate its uncertainties. Specifically, we modify by inversion the density of the uppermost mantle, the crustal density, the Moho depth, and the sedimentary cover thickness with an ensemble model with different density/geometry variation constraints. The output models indicate the most representative model of Antarctic crustal structure within the capacity of the method and current data constraints. Our preliminary results show that crustal density varies between 2.75 to 2.95 g/cm3 while the Moho depth varies between 22 km in Ross Ice Shelf and 54 km in Gamburtsev Subglacial Mountains. Low-density sedimentary basins are modelled at up to 10 km thickness beneath the ice shelf, and 3 km inland of Antarctica. Model also shows mantle density varies from 3.25 to 3.35 g/cm3. These density and thickness variations indicate likely substantial differences in crustal heat production, crustal rheology, and the expected GIA response of Antarctica’s crust.

The UN identified water as a human right in 2010; however, the bodies analyzing and designing water management often do not represent a diverse spectrum of the populations they serve. Villanova Center for Resilient Water Systems (VCRWS) research focuses on supplying resilient and sustainable solutions to global water and stormwater challenges. A diverse team of faculty, staff and students created a group that discusses diversity, equity, and inclusion (DEI) in the context of water resources engineering and associated design and solution considerations. The group has discussed a wide array of topics through the biweekly meetings. Some benefits of this group included the opportunity for personal growth and development, education on incorporating DEI into engineering thought, and re-education on DEI issues. This group served as a mechanism for members to share their grief and concerns over the political, social, and environmental issues that impede positive change. Outputs from this group include (1) competing in the Schiller Challenge through writing a paper on DEI in water engineering as a group, (2) drawing an action plan to further incorporate DEI in the activities of VCRWS through research, symposiums, conferences, public outreach programs, articles and papers, and engineering education, and (3) creating a guideline for creating DEI groups in other departments, research centers and institutions. This presentation will discuss lessons learned from creating and implementing a successful DEI focused working group within the water resource research center and will include some of the future work planned by the VCRWS group, such as incorporating DEI into the Water and Environmental Engineering curriculum through curating audiovisual materials on DEI concepts, and case studies of successes stories of DEI.

Rain gardens are green stormwater infrastructure that are designed to leverage natural processes to mitigate the impacts of urban stormwater through capturing, infiltrating, and filtering run off. Overtime these systems have the potential to buildup fines and nutrients, impacting their sustainable function. A rain garden’s performance depends on its ability to infiltrate runoff which can be reduced by clogging. Another concern is the potential transport of contaminants from rain gardens to groundwater through deep drainage. This study analyses the spatial and temporal distribution of fines and nutrients in three rain gardens through comprehensive field tests, laboratory testing, and computation analysis. Geomorphic studies were performed by integrating the digital elevation models, derived from Lidar surveys, with the FastMech solver within International River Interface Cooperative (iRIC) software, to model shear stress distribution and sediment transport relative to spatial observations of soil texture and nutrient concentrations within the rain garden. The soil properties were also used in creating models of water infiltration and nutrient sorption using Hydrus 1D. Results show that shear stresses in localized sections of each rain garden can be correlated with fines and nutrient distributions, allowing for prioritizing locations for maintenance. To conclude, LiDAR scans, flow and shear stress models, infiltration and nutrient transport models, field and laboratory soil tests can help us understand the surface dynamics and soil attributes, and gradually gain insight into the GSI performance with time.

The MHD with embedded PIC (MHD-EPIC) model makes it feasible to incorporate kinetic physics into a global simulation. Still, this requires a large enough box-shaped PIC domain to accommodate the movement and changes of the magnetic reconnection regions over time. This wastes computational resources on simulating regions with the expensive PIC model where MHD would be sufficient to describe the physics. We have developed a new MHD with Adaptively Embedded PIC (MHD-AEPIC) algorithm that couples the BATS-R-US MHD model with the new FLexible Exascale Kinetic Simulator (FLEKS) PIC code. In the new coupled model the PIC domains can move with the magnetic reconnection regions and adapt to them with an arbitrary shape. In this work, we will first introduce the algorithms for selecting the reconnection regions in the MHD model that need to be resolved with the kinetic PIC model. Then we will compare simulations obtained with MHD-EPIC using fixed PIC regions versus MHD-AEPIC employing adaptive PIC regions to verify that the new model generates reliable results. Finally, we will apply the MHD-AEPIC model to a global magnetic storm simulation and demonstrate the improved efficiency.

Worldwide coastal land-margins are prone to many flood hazards such as astronomical tides, tropical cyclones, sea-level rise, and extreme precipitation events. Compound flood events, in which two or more flooding mechanisms occur simultaneously or in close succession (Santiago-Collazo et al., 2019, https://doi.org/10.1016/j.envsoft. 2019.06.002), can exacerbate the inundation impacts due to the highly non-linear interaction of coastal and hydrologic processes. Furthermore, sea-level rise will increase the hazard at low-gradient coastal land-margins when assessing future projections due to its non-linear nuance on the compound flood (Santiago-Collazo et al., 2021, https://doi.org/10.3389/fclim.2021.684035). Therefore, there is an urgent need to develop new technologies capable of comprehensively studying compound flood events and identifying hotspots prone to these inundations. This research aims to develop a technique capable of defining and classifying coastal land- margins based on physically-based criteria due to surface flow hydrodynamics. A one-dimensional (1-D) hydrodynamic model was used to quantify the hydrodynamic response of thousands of different combinations of input parameters (e.g., astronomical tides, storm surge, precipitation, and landscape) that define a coastal land-margin. This 1-D fully-coupled model, based on the shallow water equations, was applied at a national spatial scale, considering several coastal watersheds within the Gulf of Mexico and the US East coast. One of the main goals of this tool is to identify coastal land-margins vulnerable to compound flood hazards over broad spatial scales (e.g., national or global scale). Findings suggest that low-gradient (e.g., slopes less than 0.01 m km-1) coastal land-margins are more susceptible to compound flood impacts than ones with a steeper gradient under most flooding scenarios. Future research will focus on applying this tool on a worldwide basis to test its capabilities at low-resolution, scarce data regions. A worldwide classification of coastal land-margins may help authorities, policy-makers, and professionals converge on better coastal resilience measures, such as comprehensive compound flood analysis to delineate accurate compound flood hazard maps.Full online poster version at agu2021fallmeeting-agu.ipostersessions.com/Default.aspx?s=FA-1F-20-67-21-4E-E7-69-9F-89-1E-33-BB-3D-2D-40

The Interstellar Boundary Explorer (IBEX) mission has shown that variations in the ENA flux from the outer heliosphere are associated with the solar cycle and longer-term variations in the solar wind. In particular, there is a good correlation between the dynamic pressure of the outbound solar wind and variations in the later-observed IBEX ENA flux. The time difference between observations of the outbound solar wind and the heliospheric ENAs with which they correlate ranges from approximately two to six years or more, depending on ENA energy and look direction. This time difference can be used as a means of “sounding” the heliosheath, that is, finding the average distance to the ENA source region in a particular direction. We apply this method to build a three-dimensional map of the heliosphere. We use IBEX ENA data collected over a complete solar cycle, from 2009 through 2019, corrected for survival probability to the inner heliosphere. We divide the data into 56 “macro-pixels” covering the entire sky, and as each point in the sky is sampled once every six months, this gives us a time series of 22 points per macro-pixel on which to time-correlate. Consistent with prior studies and heliospheric models, we find that the shortest distance to the heliopause dHP is slightly south of the nose direction (dHP ~ 110 – 120 au), with a flaring toward the flanks and poles (dHP ~ 160 – 180 au). The heliosphere extends at least ~350 au tailward, which is the distance limit of the technique.

Some of the Earth system data products such as those from NASA airborne and field investigations (a.k.a. campaigns), are highly heterogeneous and cross-disciplinary, making the data extremely challenging to manage. For example, airborne and field campaign measurements tend to be sporadic over a period of time, with large gaps. Data products generated are of various processing levels and utilized for a wide range of inter- and cross-disciplinary research and applications. Data and derived products have been historically stored in a variety of domain-specific standard (and some non-standard) formats and in various locations such as NASA Distributed Active Archive Centers (DAACs), NASA airborne science facilities, field archives, or even individual scientists’ computer hard drives. As a result, airborne and field campaign data products have often been managed and represented differently, making it onerous for data users to find, access, and utilize campaign data. Some difficulties in discovering and accessing the campaign data originate from the incomplete data product and contextual metadata that may contain details relevant to the campaign (e.g. campaign acronym and instrument deployment locations), but tend to lack other significant information needed to understand conditions surrounding the data. Such details can be burdensome to locate after the conclusion of a campaign. Utilizing consistent terminology, essential for improved discovery and reuse, is also challenging due to the variety of involved disciplines. To help address the aforementioned challenges faced by many repositories and data managers handling airborne and field data, this presentation will describe stewardship practices developed by the Airborne Data Management Group (ADMG) within the Interagency Implementation and Advanced Concepts Team (IMPACT) under the NASA’s Earth Science Data systems (ESDS) Program.

Earthquake moment release is localized along a global fault system. This network of branching and anastomosing fractures defines the geometrically complex boundaries of tectonic plates and serves as the locus of contemporary elastic strain energy storage between earthquakes. The slow deformation of the earth’s crust in between earthquakes has been observed geodetically for decades and provides a filtered representation of the underlying earthquake behaviors. Here we describe efforts to model fault system activity at a global scale incorporating both tectonic plate motions and earthquake cycle effects. Interseismic earthquake cycle effects are represented using a first-order quasi-static elastic approximation, and these models yield a unified estimate of slip deficit rates and subduction zone coupling constrained by nominally interseismic geodetic surface velocity estimates. We present key findings from a kinematic global fault system model with 1.6×107 km2 of fault system area including 16 subduction zones and constrained by observations 22,500+ GPS velocities. Further, we describe new approaches to the efficient representation of viscoelastic deformation in large-scale block models and the prospects for high-resolution block scale models that directly image partial fault coupling across the entire global fault system. Because global geodetic observations capture faults behaviors at varying stages throughout the earthquake cycle, consideration of time-dependent deformation including viscous dissipation of coseismically induced stresses is important for accurate imaging of fault coupling. And, because concentrations of fault coupling have been shown to spatially correlate with recent significant earthquakes, being able to estimate partial coupling patterns on a global scale may highlight pending seismicity.

We study how ground frost affects the ambient seismic wavefield recorded by a three-component broadband sensor. By applying machine learning algorithms on continuous seismic data, we can retrieve the seismic signature of the continuous freeze and thaw process at the surface of the ground. The retrieved signature reveals that the presence of ground frost imprints the amplitude of the ambient seismic wavefield, and the energy ratio between horizontal and vertical components (H/V). A regression model can even predict diurnal freeze and thaw patterns based on the seismic data. Thus, we assume that slight changes in the physical properties of the frozen surface, such as the thickness, alter the seismic wavefield. Models of the subsurface with different properties of the ground frost agree with the observations from the field. The penetration depth of the ground frost, the temperature of the frozen ground, and the presence of different modes in the wavefield determine how the seismic wavefield is changing. The findings of this study show the potential of a single seismic station for monitoring frozen bodies near the surface, such as permafrosts.

Recently, the continually increasing availability of seismic data has allowed high-resolution imaging of lithospheric structure beneath the African cratons. In this study, S-wave seismic tomography are combined with high resolution satellite gravity data in an integrated approach to investigate the structure of the cratonic lithosphere of Africa. A new model for the Moho depth and data on the crustal density structure are employed along with global dynamic models to calculate residual topography and mantle gravity residuals. Corrections for thermal effects of an initially juvenile mantle are estimated based on S-wave tomography and mineral physics. Joint inversion of the residuals yields necessary compositional adjustments that allow to recalculate the thermal effects. After several iterations, we obtain a consistent model of upper mantle temperature, thermal and compositional density variations, and Mg# as a measure of depletion, as well as an improved crustal density model. Our results show that thick and cold depleted lithosphere underlies West African, northern to central eastern Congo, and Zimbabwe Cratons. However, for most of these regions, the areal extent of their depleted lithosphere differs from the respective exposed Archean shields. Meanwhile, the lithosphere of Uganda, Tanzania, most of eastern and southern Congo, and the Kaapvaal Craton is thinner, warmer, and shows little or no depletion. Furthermore, the results allow to infer that the lithosphere of the exposed Archean shields of Congo and West African cratons was depleted before the single blocks were merged into their respective cratons.

Mangrove forests with complex root systems contribute to increased coastal protection through drag effects. Previous flume studies proposed a predictive model of drag in Rhizophora mangrove forests based on quadratic drag law. However, its general applicability on mangrove forests in the field has not been tested. To fill this knowledge gap, this study quantified drag in a 17-year-old planted Rhizophora mangrove forest using a comprehensive measurement of hydrodynamics and vegetation morphology. The vegetation projected area density, a, showed an approximate exponential increase towards the bed, mainly due to root branching. This vertical variation led to enhanced vegetation drag per unit water volume relative to velocity with decreasing water depth. Alternatively, the drag per vegetation projected area solely depended on the square of velocity, indicating association with the quadratic drag law. The derived drag coefficient (CD) was 1.0 ± 0.2 for tide-driven currents, consistent with previous flume studies. By using the mean value of derived CD (1.0), it was confirmed that the quadratic drag model expresses well the field-measured drag. We also presented a method for predicting a value for a, another unknown parameter in the drag model, using an empirical Rhizophora root model, and confirmed a successful prediction of a and drag. Therefore, the drag in a Rhizophora mangrove forest can be accurately predicted only using the input parameters of the Rhizophora root model – stem diameter and tree density. This provides insights into effectively implementing the drag model in hydrodynamic models for better representation of mangroves’ coastal protection function.

Floods are the most frequent, costliest natural disasters having devastating consequences on people, infrastructure, and the ecosystem. The accurate and rapid mapping of the flooded areas becomes more crucial when floods strike densely populated cities. During flood events near real-time satellite imagery has been proven to be an efficient management tool for disaster management authorities. However one of the challenges is accurate classification and segmentation of flooded water and permanent water. Binary segmentation using the threshold split-based method is commonly used in this regard, however, the generalization ability of this method is limited due to the effects of backscatter, geographical area, and time of image collection. Recent advancements in deep learning algorithms for image segmentation has demonstrated the excellent potential of Convolutional Neural Networks(CNN) for improving flood detection, although there have been limited studies in this domain due to the lack of large scale labeled flood event dataset. In this project, we present a U-net based deep learning approach by leveraging publicly available Sentinel-1 dataset provided jointly by NASA Interagency Implementation and Advanced Concepts Team and IEEE GRSS Earth Science Informatics Technical Committee. Dataset is composed of 66,810 tiles of 256×256 pixels, distributed respectively across the training, validation and test sets and cover flood events from Nebraska, North Alabama, Bangladesh, Red River North and Florence. Specifically we proposed an Unet architecture based convolutional neural network (CNN) with a backbone of EfficientNetb7, trained against the dataset. We then evaluated the performance of the model with multiple training, testing and validation. Two evaluation methods - Intersection over Union (IOU) and F-Score are adopted to evaluate the model performance. During testing, the model achieved the meanIOU score of 75.06% and F-Score of 74.98%. We hope to further improve the performance of the network by performing hyper-parameter tuning and to develop a model which can be used for near-real-time flood mapping.

Our objective is to test and improve cloud subcolumn generators used for greater realism of scales in the radiation schemes and satellite simulators GCMs. For this purpose, we use as guidance water content fields from active observations by the CloudSat radar (CPR) and the CALIPSO lidar (CALIOP). Cloud products from active sensors while suffering significant sampling and coverage drawbacks have the advantage of resolving both horizontal and vertical variability which is what the generators are designed to produce. Our first order goal is to test the ability of the generators to deliver realistic 2D cloud extinction (cloud optical thickness) fields using, as in GCMs, limited domain-averaged information. Our reference 2D cloud extinction fields fully resolving horizontal (along the track of the satellites) and vertical variability come from combining CloudSat’s 2B-CWC-RVOD (liquid clouds) and CALIPSO-enhanced 2C-ICE (ice clouds) products. The combined fields were improved by introducing a simple scheme to fill liquid cloud extinction values identified as missing by comparing with coincident 2D (phase-specific) cloud masks provided by the CALIPSO-enhanced 2B-CLDCLASS-LIDAR CloudSat product. Our presentation will demonstrate the substantial improvements for low clouds brought by the filling scheme through comparisons with MODIS-Aqua cloud fraction distributions expressed in terms of joint cloud top pressure – cloud optical thickness histograms. Beyond global comparisons, the nature of the improvements become clearer when comparing mean joint histograms segregated by MODIS Cloud Regime (CR): improvement is by design superior for MODIS CRs dominated by low clouds. With the improved 2D extinction fields at hand, we test the skill of two subcolumn generators, one used in the COSP satellite simulator package, and one with more sophisticated cloud overlap implemented in the GEOS global model, to reproduce joint histograms that are statistically similar to the observed counterparts described above (as interpreted by COSP’s MODIS simulator). Our main comparison metrics are the Euclidean distance between observed and generator-produced global or near-global mean joint histograms, and the statistics of Euclidean distances calculated for individual scenes. One full year of data is used to assess whether the more sophisticated cloud generator produces clouds with greater realism in 2D cloud variability.

We propose a new approach to the solution of the wave propagation and full waveform inversions (FWIs) based on a recent advance in deep learning called Physics-Informed Neural Networks (PINNs). In this study, we present an algorithm for PINNs applied to the acoustic wave equation and test the model with both forward wave propagation and FWIs case studies. These synthetic case studies are designed to explore the ability of PINNs to handle varying degrees of structural complexity using both teleseismic plane waves and seismic point sources. PINNs’ meshless formalism allows for a flexible implementation of the wave equation and different types of boundary conditions. For instance, our models demonstrate that PINN automatically satisfies absorbing boundary conditions, a serious computational challenge for common wave propagation solvers. Furthermore, a priori knowledge of the subsurface structure can be seamlessly encoded in PINNs’ formulation. We find that the current state-of-the-art PINNs provide good results for the forward model, even though spectral element or finite difference methods are more efficient and accurate. More importantly, our results demonstrate that PINNs yield excellent results for inversions on all cases considered and with limited computational complexity. Using PINNs as a geophysical inversion solver offers exciting perspectives, not only for the full waveform seismic inversions, but also when dealing with other geophysical datasets (e.g., magnetotellurics, gravity) as well as joint inversions because of its robust framework and simple implementation.

Time series of shipboard observations in the southern Arafura Sea near the Tiwi Islands indicated that the water column dynamics differed between the east and west sides of the islands. On the west side, the water column, characterized by temperature, salinity, and velocity, was barotropic and tidal advection dominated. On the east side, the water column was baroclinic and internal tides were present along with tidal advection. These conditions affected the distribution of the turbidity and fluorescence in the water column. Likewise, the influence of the daily solar radiation cycle reached the bottom on the western side, but was limited to the upper layer above the thermocline on the eastern side. The fluorescence peaks also differed between the east and west sides, with the eastern side dominated by the semidiurnal tides and the western side by the daily solar cycle. Fluorescence integrated over the water column was much higher on the eastern side than the western side. Also on the eastern side, fluorescence was limited to the lower layer, while on the western side, it encompassed the entire water column at times and peaked below the warmer, higher oxygenated water generated by solar radiation and surface mixing. These dynamics have distinct implications for biological productivity and also may affect a proposed tidal power system in the region.