Engineering Papers⌕ Search

Engineering topics

Chen, Xingyuan

Publications and source records attributed to Chen, Xingyuan.

At least 19 records

A Field Guide to Corralling the Chaos: A Conceptual Framework for Using Models to Guide Opportunistic Field Studies of Natural Disturbances

Watersheds regulate biogeochemical processes and provide ecosystem services to human societies, but disturbances can fundamentally alter these processes across space and time. Determining when and where to sample to capture disturbance impacts in watersheds remains a central challenge. Manipulation studies and long-term monitoring are often constrained by scope, and opportunistic studies often lack pre-disturbance data needed to statistically determine disturbance impacts. We identify a persistent knowledge gap: the absence of a clear, transferable framework to guide opportunistic disturbance research where pre-disturbance data collection is not a feasible option. To address this gap, we present a conceptual framework that intentionally integrates modeling and empirical observation in an iterative, stepwise model–experiment workflow. We demonstrate its application through two contrasting case studies: wildfire impacts on headwater streams using a pre-disturbance preparedness approach, and saltwater flooding impacts on coastal forests using an ‘ex-post-facto’ approach. From these applications, we assess strengths, limitations, and the critical role of team science for transferability across disturbance types and study designs. Broadly, this framework offers a scalable path towards more rigorous, timely, and actionable disturbance science that can inform watershed management, hazard risk reduction, and ecosystem resilience.

Coastal Biogeochemistry↗

Water Column Respiration in the Yakima River Basin is Explained by Temperature, Nutrients and Suspended Solids

Understanding aquatic ecosystem metabolism involves the study of two key processes: carbon fixation via primary production and organic C mineralization as total ecosystem respiration (ERtot). In streams and rivers, ERtot includes respiration in the water column (ERwc) and in the sediments (ERsed). While literature surveys suggest that ERsed is often a dominant contributor to ERtot, recent studies indicate that the relative influence of sediment-associated processes versus water column processes can fluctuate along the river continuum. Still, a comprehensive understanding of the factors contributing to these shifts within basins and across stream orders is needed. Here we contribute to this need by measuring ERwc and collecting water samples across 47 sites in the Yakima River basin, Washington, USA. We found that ERwc rates varied throughout the basin during baseflow conditions, ranging from –7.38 to 0.36 g O2 m?3 d?1, and encompassed the range of ERwc literature values. Additionally, by comparing to ERtot estimates for rivers across the contiguous United States, we suggest that the contribution of ERwc rates to reach-scale ERtot rates across the Yakima River was likely highly variable, but we did not test this directly. We observed that temperature, nutrient concentrations (dissolved organic carbon, total dissolved nitrogen), and total suspended solids explained 41% of ERwc variability across the basin. Our findings highlight the potential relevance of water column processes in aquatic ecosystem metabolism, with the Yakima River basin serving as an environmentally diverse river network representative of the larger Columbia River basin that spans much of the Pacific Northwest region of the United States. Our results are generally congruent with previous work, suggesting that the observed variability and suite of associated environmental factors influencing ERwc are potentially transferable across basins.

Laan, Maggi M.↗

Interactive effects of salinity, redox, and colloids on greenhouse gas production and carbon mobility in coastal wetland soils

Coastal wetlands, including freshwater systems near large lakes, rapidly bury carbon, but less is known about how they transport carbon either to marine and lake environments or to the atmosphere as greenhouse gases (GHGs) such as carbon dioxide and methane. This study examines how GHG production and organic matter (OM) mobility in coastal wetland soils vary with the availability of oxygen and other terminal electron acceptors. We also evaluated how OM and redox-sensitive species varied across different size fractions: particulates (0.45–1μm), fine colloids (0.1–0.45μm), and nano particulates plus truly soluble (<0.1μm; NP+S) during 21-day aerobic and anaerobic slurry incubations. Soils were collected from the center of a freshwater coastal wetland (FW-C) in Lake Erie, the upland-wetland edge of the same wetland (FW-E), and the center of a saline coastal wetland (SW-C) in the Pacific Northwest, USA. Anaerobic methane production for FW-E soils were 47 and 27,537 times greater than FW-C and SW-C soils, respectively. High Fe 2+ and dissolved sulfate concentrations in FW-C and SW-C soils suggest that iron and/or sulfate reduction inhibited methanogenesis. Aerobic CO 2 production was highest for both freshwater soils, which had a higher proportion of OM in the NP+S fraction (64±28% and 70±10% for FW-C and FW-E, respectively) and organic C:N ratios reflective of microbial detritus (5.3±5.3 and 5.3±7.0 for FW-E and FW-C, respectively) compared to SW-C, which had a higher fraction of particulate (58±9%) and fine colloidal (19±7%) OM and organic C:N ratios reflective of vegetation detritus (11.4 ± 1.7). The variability in GHG production and shifts in OM size fractionation and composition observed across freshwater and saline soils collected within individual and across different sites reinforce the high spatial variability in the processes controlling OM stability, mobility, and bioavailability in coastal wetland soils.

54 ENVIRONMENTAL SCIENCES↗

Quantifying Streambed Grain Size, Uncertainty, and Hydrobiogeochemical Parameters Using Machine Learning Model YOLO

Abstract Streambed grain sizes control river hydro‐biogeochemical (HBGC) processes and functions. However, measuring their quantities, distributions, and uncertainties is challenging due to the diversity and heterogeneity of natural streams. This work presents a photo‐driven, artificial intelligence (AI)‐enabled, and theory‐based workflow for extracting the quantities, distributions, and uncertainties of streambed grain sizes from photos. Specifically, we first trained You Only Look Once, an object detection AI, using 11,977 grain labels from 36 photos collected from nine different stream environments. We demonstrated its accuracy with a coefficient of determination of 0.98, a Nash–Sutcliffe efficiency of 0.98, and a mean absolute relative error of 6.65% in predicting the median grain size of 20 ground‐truth photos representing nine typical stream environments. The AI is then used to extract the grain size distributions and determine their characteristic grain sizes, including the 10th, 50th, 60th, and 84th percentiles, for 1,999 photos taken at 66 sites within a watershed in the Northwest US. The results indicate that the 10th, median, 60th, and 84th percentiles of the grain sizes follow log‐normal distributions, with most likely values of 2.49, 6.62, 7.68, and 10.78 cm, respectively. The average uncertainties associated with these values are 9.70%, 7.33%, 9.27%, and 11.11%, respectively. These data allow for the computation of the quantities, distributions, and uncertainties of streambed HBGC parameters, including Manning's coefficient, Darcy‐Weisbach friction factor, top layer interstitial velocity magnitude, and nitrate uptake velocity. Additionally, major sources of uncertainty in grain sizes and their impact on HBGC parameters are examined.

58 GEOSCIENCES↗

Bridging Hydrological Ensemble Simulation and Learning Using Deep Neural Operators

Ensemble-based simulation and learning (ESnL) has long been used in hydrology for parameter inference, but computational demands of process-based ESnL can be quite high. To address this issue, we propose a deep neural operator learning approach. Neural operators are generic machine learning algorithms that can learn functional mappings between infinite-dimensional spaces, providing a highly flexible tool for scientific machine learning. Our approach is built upon DeepONet, a specific deep neural operator, and is designed to address several common problems in hydrology, namely, model parameter estimation, prediction at ungaged locations, and uncertainty quantification. Here we demonstrate the effectiveness of our DeepONet-based workflow using an existing large model ensemble created for an eastern U.S. watershed that is instrumented with 10 streamflow gages. Results suggest DeepONet achieves high efficiency in learning an ML surrogate model from the model ensemble, with the modified Kling-Gupta Efficiency exceeding 0.9 on holdout test sets. Parameter inference, carried out using the trained DeepONet surrogate model and genetic algorithm, also yields robust results. Additionally, we formulate and train a separate DeepONet model for physics-informed, seq-to-seq streamflow forecasting, which further reduces biases in the pre-trained DeepONet surrogate model. While this study focuses primarily on a single watershed, our approach is general and may be extended to enable learning from model ensembles across multiple basins or models. Thus, this research represents a significant contribution to the application of hybrid machine learning in hydrology.

54 ENVIRONMENTAL SCIENCES↗

Mechanistic definition and prediction of the mass exchange coefficient between rivers and hyporheic zones: The $α$ of two $Ω$s

Solute transport in interconnected rivers and hyporheic zones is typically modeled through dual-domain models where first-order solute mass transfer between the two domains, Ω R and Ω HZ , is represented by a coefficient α. The transient storage model (TSM) is an example of such an approach. In practice, α is determined by fitting the tails of solute tracer breakthrough curves using a TSM. This approach has led to ambiguity regarding α’s physical meaning and transferability. Here, in this work, we investigated the physical basis for α and tested it with virtual experiments through the fully coupled multiphysics model hyporheicFoam for the Ω R – Ω HZ system. hyporheicFoam explicitly simulated coupled flow and solute transport over a kilometer with centimeter-scale resolution. Model results were analyzed to calculate α following its theoretical definition directly. Using the determined α within a TSM enables accurate reproduction of solute transport, underscoring α’s physical relevance and precision.

54 ENVIRONMENTAL SCIENCES↗

River–aquifer interactions enhancing evapotranspiration in a semiarid riparian zone: A modelling study

The hydrologic flows across the river–aquifer interface play an important role in groundwater dynamics and biogeochemical reactions within the subsurface; however, little is known about the effects of river–aquifer interactions on land surface processes. In this study, we developed a fully coupled three-dimensional (3D) land surface and subsurface model at a high resolution (~1 km) that accounts for high-frequency hydrologic exchange flow conditions to investigate how river–aquifer interactions modulate surface water budgets in the Upper Columbia-Priest Rapids watershed, a typical semiarid watershed located in the northwestern United States where river stage fluctuates in response to reservoir releases changing. Our results show that the spatiotemporal dynamics of river–aquifer interactions are highly heterogeneous, driven mainly by river-stage fluctuations. Adding 6.64 × 10 6 m 3 year –1 of water over the watershed from the river to groundwater owing to the lateral flow, river–aquifer interactions led to an increase in soil evaporation and transpiration supplied by higher soil moisture content, particularly in deeper subsurface. In a hypothetic future scenarios where a 5-m rise in river stage was assumed, the hydrologic flow exchange rates were intensified, resulting in higher surface water over the entire watershed. Overall, lateral flow induced by river–aquifer exchanges leads to an increase in evapotranspiration of ~75% in the historical period and of ~83% in the hypothetical future scenario. Finally, our study demonstrates the potential of coupled model as an effective tool for understanding river–aquifer–land surface interactions, and indicates that river–aquifer interactions fundamentally alter the water balance of the riparian zone for the semiarid watershed and will likely become more frequent and intense in the future under the effects of climate change.

54 ENVIRONMENTAL SCIENCES↗

Integrating a Water Tracer Model Into WRF‐Hydro for Characterizing the Effect of Lateral Flow in Hydrologic Simulations

Abstract Most current land models approximate terrestrial hydrological processes as one‐dimensional vertical flow, neglecting lateral water movement from ridges to valleys. Such lateral flow is fundamental at catchment scales and becomes crucial for finer‐scale land models. To test the effect of incorporating lateral flow toward three‐dimensional representations of hydrological processes in the next generation land models, we integrate a water tracer model into the WRF‐Hydro framework to track water movement from precipitation to discharge and evapotranspiration. This hydrologic‐tracer integrated system allows us to identify the key mechanisms by which lateral flow affects the flow paths and transit times in WRF‐Hydro. By comparing modeling experiments with and without lateral routing in two contrasting catchments, we determine the impacts of lateral flow on the transit times of precipitation event‐water. Results show that with limited hydrologic connectivity, lateral flow extends the transit times by reducing (increasing) event‐water drainage loss (accumulation) in ridges (valleys) and allowing reinfiltration of infiltration‐excess flow, which is missing in most land models. On the contrary with high hydrologic connectivity, lateral flow can effectively accelerate the water release to streams and reduce the transit time. However, the transit times are substantially underestimated by the model compared with isotope‐derived estimates, indicating model limitations in representing flow paths and transit times. This study provides some insights on the fundamental differences in terrestrial hydrology simulated by land models with and without lateral flow representation.

54 ENVIRONMENTAL SCIENCES↗

Learning Constitutive Relations From Soil Moisture Data via Physically Constrained Neural Networks

Abstract The constitutive relations of the Richardson‐Richards equation encode the macroscopic properties of soil water retention and conductivity. These soil hydraulic functions are commonly represented by models with a handful of parameters. The limited degrees of freedom of such soil hydraulic models constrain our ability to extract soil hydraulic properties from soil moisture data via inverse modeling. We present a new free‐form approach to learning the constitutive relations using physically constrained neural networks. We implemented the inverse modeling framework in a differentiable modeling framework, JAX, to ensure scalability and extensibility. For efficient gradient computations, we implemented implicit differentiation through a nonlinear solver for the Richardson‐Richards equation. We tested the framework against synthetic noisy data and demonstrated its robustness against varying magnitudes of noise and degrees of freedom of the neural networks. We applied the framework to soil moisture data from an upward infiltration experiment and demonstrated that the neural network‐based approach was better fitted to the experimental data than a parametric model and that the framework can learn the constitutive relations.

54 ENVIRONMENTAL SCIENCES↗

Integrated Effects of Site Hydrology and Vegetation on Exchange Fluxes and Nutrient Cycling at a Coastal Terrestrial‐Aquatic Interface

Abstract The complex interactions among soil, vegetation, and site hydrologic conditions driven by precipitation and tidal cycles control the biogeochemical transformations and bi‐directional exchange of carbon and nutrients across the terrestrial–aquatic interfaces (TAIs) in coastal regions. This study uses a highly mechanistic model, Advanced Terrestrial Simulator (ATS)‐PFLOTRAN, to explore how these interactions affect exchanges of materials and carbon and nitrogen cycling. We used a transect in the Chesapeake Bay region that spans zones of open water, coastal wetland, transition, and upland forest. We designed several simulation scenarios to parse the effects of the individual controlling factors and the sensitivity of carbon cycling to reaction rate parameters derived from laboratory experiments. Our simulations reveal an active zone for carbon cycling under the transition zones between the wetland and the upland. Evapotranspiration is found to enhance the exchange fluxes between the surface and subsurface domains, resulting in a higher dissolved oxygen concentration in the TAIs. The transport of organic carbon derived from plant leaves and roots provide an additional source of organic carbon needed for the aerobic respiration and denitrification processes in the TAIs. The variability in reaction rate parameters associated with microbial activities is also found to play a dominant role in controlling the heterogeneity and dynamics of the simulated redox conditions. This modeling‐focused exploratory study enabled us to better understand the complex interactions among soil, water and microbes that govern the hydro‐biogeochemical processes at the TAIs, which is an important step toward representing coastal ecosystems in larger‐scale Earth system models.

54 ENVIRONMENTAL SCIENCES↗

Sinuosity‐Driven Hyporheic Exchange: Hydrodynamics and Biogeochemical Potentials

Abstract Hydrologic exchange processes are critical for ecosystem services along river corridors. Meandering contributes to this exchange by driving channel water, solutes, and energy through the surrounding alluvium, a process called sinuosity‐driven hyporheic exchange. This exchange is embedded within and modulated by the regional groundwater flow (RGF), which compresses the hyporheic zone and potentially diminishes its overall impact. Quantifying the role of sinuosity‐driven hyporheic exchange at the reach‐to‐watershed scale requires a mechanistic understanding of the interplay between drivers (meander planform) and modulators (RGF) and its implications for biogeochemical transformations. Here, we use a 2D, vertically integrated numerical model for flow, transport, and reaction to analyze sinuosity‐driven hyporheic exchange systematically. Using this model, we propose a dimensionless framework to explore the role of meander planform and RGF in hydrodynamics and how they constrain nitrogen cycling. Our results highlight the importance of meander topology for water flow and age. We demonstrate how the meander neck induces a shielding effect that protects the hyporheic zone against RGF, imposing a physical constraint on biogeochemical transformations. Furthermore, we explore the conditions when a meander acts as a net nitrogen source or sink. This transition in the net biogeochemical potential is described by a handful of dimensionless physical and biogeochemical parameters that can be measured or constrained from literature and remote sensing. This work provides a new physically based model that quantifies sinuosity‐driven hyporheic exchange and biogeochemical reactions, a critical step toward their representation in water quality models and the design and assessment of river restoration strategies.

54 ENVIRONMENTAL SCIENCES↗

Data and Scripts Associated with the Manuscript “Water Column Respiration in the Yakima River Basin is Explained by Temperature, Nutrients and Suspended Solids”

This data package is associated with the publication “Water Column Respiration in the Yakima River Basin is Explained by Temperature, Nutrients and Suspended Solids” published in EGU Biogeochemistry (Laan et al. 2025). In this research, water column respiration (ERwc) data, surface water chemistry data, organic matter (OM) chemistry data, and publicly available geospatial data were used in analysis to evaluate the variability in ERwc at 47 sites across the Yakima River basin in Washington, USA. In addition to this readme, this data package also includes a file-level metadata (FLMD) file that describes each file and a data dictionary (DD) that describes all column/row headers and variable definitions. The data package includes the data inputs, and outputs, and R scripts to reproduce all the analyses performed in the manuscript and create manuscript figures. The data package is comprised of three main folders (Code, Data, and Figures). The Code folder is comprised of four scripts and three analysis-specific subfolders that contain the R scripts to perform the analyses described in the publication and create publication figures. The Data folder is comprised of two “.csv” files and four subfolders that contain data input and output files. The Published_Data folder contains a readme that directs the user to download the appropriate files and add to this folder when using scripts. The Figures folder includes figures from the manuscript in “.pdf” and “.png” formats and a folder with intermediate figure files. This data package is associated with a GitHub repository which can be found at https://github.com/river-corridors-sfa/rcsfa-RC2-SPS-ERwc. We acknowledge the Yakama Nation as owners and caretakers of the lands where we collected some of these data. We thank the Confederated Tribes and Bands of the Yakama Nation Tribal Council and Yakama Nation Fisheries for working with us to facilitate sample collection and optimization of data usage according to their values and worldview.

54 ENVIRONMENTAL SCIENCES↗

Data and scripts associated with the manuscript "Estimating soil moisture and salinity response to simulated coastal flooding using time-lapse electrical resistivity and induced polarization monitoring”

This package contains geophysics datasets generated during the simulated ecosystem-flooding experiment – TEMPEST (Terrestrial Ecosystem Manipulation to Probe the Effects of Storm Treatment). These geophysics data consist of Electrical Resistivity Imaging (ERI) and Induced Polarization (IP) datasets collected to estimate the soil moisture and salinity response to the simulated coastal flooding experiment. In addition to using geophysics datasets to estimate soil moisture and salinity response, petrophysical relationship between the measured resistivity and moisture content and salinity using the multisalinity. This package consists of two folders: “MultiSalinity Experiment” folder and “TEMPEST Experiment and petrophysical conversion” folders. The “MultiSalinity Experiment” folder contains the datasets generated from the multisalinity experiment used to develop the petrophysical reslationship. The “TEMPEST Experiment and petrophysical conversion” folder contains ERI, IP, Ground Penetrating Radar (GPR) and Electromagnetic Induction (EMI) datasets collected during the TEMPEST experiment. This study shows the ability to used geophysical datasets in estimating moisture content and salinity at scale will be useful in calibrating Earth System Models.

54 ENVIRONMENTAL SCIENCES↗

Data and scripts associated with the manuscript "A Hydrogeophysical Framework to Assess Infiltration During a Simulated Ecosystem-Scale Flooding Experiment”

This package contains the data and scripts used in "A hydrogeophysical framework to assess infiltration during a simulated ecosystem-scale flooding experiment" (Adebayo et al., 2023). The Geophysics folder contains two subfolders containing Ground Penetrating Radar (GPR) and Electrical Resistivity Imaging (ERI) datasets acquired between August 23 and August 26, 2021. The GPR subfolder contains raw data (Tempest_Project.GPZ) and processed results using the Sensors and Software EKKO Project software. The ERI subfolder contains two additional subfolders named Static and Timelapse. The Static subfolder contains background resistivity datasets acquired before infiltration, both raw and inverted results using the ResIpy code. The Timelapse subfolder contains repeated resistivity datasets acquired during infiltration on two profiles, along with raw and inverted results using the ResIpy code. The Flow model folder contains ATS model simulation files. The Geophysics-informed subfolder contains simulation datasets for two profiles generated using information from the geophysics data. The National_database-informed subfolder contains simulation datasets for two profiles generated using information from a national soil database. Each subfolder contains initial simulation results, multiple ensembles, and final simulation files. The Sensitivity_analysis folder contains the observation files exported from the flow model, notebooks used in running the code, and sensitivity analysis results showing the sensitivity of different parameters over a 48-hour period. We used mutual information-based global sensitivity analysis for the sensitivity analysis, following Jiang et al., 2022. This study shows a novel approach that integrates geophysical data with hydrological modeling to understand subsurface flow during coastal flooding.

54 ENVIRONMENTAL SCIENCES↗

Data associated with the manuscript “Radiative impact of record-breaking wildfires from integrated ground-based data” collected in Richland, Washington in September 2020

This data package is associated with the publication “Radiative impact of record-breaking wildfires from integrated ground-based data” submitted to Nature Scientific Reports (Kassianov et al., 2024). Data from ground-based measurements of shortwave and spectrally resolved irradiance and aerosol optical depth (AOD) in the visible and near-infrared spectral ranges were assessed to quantify the radiative impact of the September 2020 wildfires that occurred in the Western United States. Data were collected in September 2020 by several ground-based instruments at the Atmospheric Measurements Laboratory (AML) located in Richland, Washington (46.3451, -119.2792). These data include (1) Aerosol Optical Depth (AOD); (2) spectrally resolved and shortwave (SW) irradiances; (3) backscatter profiles; (4) total sky images; and (5) near-surface ambient air temperatures.The data package consists of five sub-directories: (1) “AML_Ceilometer_”; (2)” AML_CSPHOT_”; (3) “AML_MFRSR_irradiances_”; (4) “AML_SW_irradiances_and_Temp_”; (5) “AML_TSI_images_”; and 6 files stored at the directory level, including the readme, file-level metadata file, and data dictionary. The file-level metadata file (the file ending in “_flmd.csv”) lists all files contained in this data package and descriptions for each. The data dictionary (the file ending in “_dd.csv”) describes each tabular column header’s unit, definition, and structure. Below are descriptions of each sub-directory:“AML_Ceilometer_” includes ceilometer data collected at the AML. These files contain the corresponding narratives of data. Details related to the ceilometer data can be found in Morris (2016). “AML_CSPHOT_” includes ascii files with high-temporal resolution (about 10-15 min) AML CSPHOT data and their daily-averaged counterparts. These two files contain the corresponding narratives of data. Details related to the CSPHOT data can be found in Gregory (2011). “AML_MFRSR_irradiances_” includes ascii files with the AML MFRSR-measured diffuse, normal, and total spectrally resolved irradiance. Details related to the MFRSR data can be found in Hodges and Michalsky (2016) and Koontz et al. (2013). “AML_SW_irradiances_+_Temp_” includes near-surface ambient air temperature and SW irradiances, namely direct normal, diffuse hemispherical, and total hemispheric (global), measured at the AML. These files also incorporate the corresponding narratives of data. Details related to the SW irradiances can be found in Andreas et al. (2018). “AML_TSI_images_” includes Total Sky Images (TSIs) collected at the AML. Details related to the TSI data can be found in Morris (2005).

54 ENVIRONMENTAL SCIENCES↗

Data for: A Synoptic System for Capturing Ecosystem Control Points Across Terrestrial-Aquatic Interfaces

The investigation of how climate change and water level fluctuations impact variable and interconnected ecosystems, like the interfaces between terrestrial and aquatic environments, requires the collection and integration of many data types. We describe an integrative and autonomous environmental monitoring approach that uses environmental sensors and data loggers to monitor surface water, groundwater, soil, and vegetation changes and generate essential data for predictive models. We established the network at seven sites along the Chesapeake Bay and Lake Erie coastlines, including a large-scale flood manipulation experiment, collectively generating over three million observations per month. Such sensor networks hold great promise for tracking and comprehending environmental changes where land and water intersect. The sensor system and overall approach to sensor management that we have designed is intended to be widely accessible for research teams spanning in size from an individual investigator to large multi-institution projects. This dataset shows example data generated by the sensor network described above. Data output for data loggers connected to groundwater water quality sondes measuring dissolved oxygen, pH, oxidative redox potential (ORP), groundwater elevation, groundwater salinity and temperature, and groundwater elevation; replicate soil moisture and conductivity probes installed at 10 and 30 cm below the ground surface; rainfall and solar radiation; and mean sap flow from 8 replicate probes. Sensors are installed at upland forest (UP), wetlands (W), and transitional locations between these ecosystems (TR). Example datasets are from Goodwin Islands in the Chesapeake Bay and Portage River along the Lake Erie coastline.

54 ENVIRONMENTAL SCIENCES↗