Search NASA⌕ Search

SEARCH · Search NASA

Results for “coastal model”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

A Multi-variable, Multi-value Hydrodynamic Modeling Approach to Support Marine Energy and Coastal Resilience Applications

Funded by the Marine Energy Seedling Laboratory Program, the study reported on here aimed to develop a multi-variable, multi-value hydrodynamic modeling approach based on the Finite-Volume Community Ocean Model (FVCOM) to support a wide range of marine renewable energy and coastal resilience applications. We migrated the marine hydrokinetic energy module to simulate tidal energy extraction by tidal turbines under wave-current interactions. We explored the particle tracking module in FVCOM and applied it to Sequim Bay, Washington, to illustrate the surface current field and mimic the pollutant transport near the surface. A tracer-water age module was added to FVCOM to examine water exchange through the entrance. Lastly, a user-specified vertical coordinate option was added to FVCOM to improve simulation of the effect of floating turbines on surface currents and scalar transport. The above-described model development has been tested using a simplified model domain and demonstrated in real-world coastal systems. The outcome of this project will provide a useful modeling code and approach for marine renewable energy and coastal resilience research and applications.

16 TIDAL AND WAVE POWER↗
DOE OSTI2021

Assessment of Model Hub Height Wind Speed Performance Using DOE Lidar Buoy Data

As the offshore wind industry continues to gain momentum both in the U.S. and globally, reanalysis products remain essential tools for the lifecycle of a wind farm. From project siting and production estimates to determining construction and maintenance windows, reanalysis products play an important role. These regional and global long-term models also serve as inputs to higher resolution models. With such wide-ranging impacts on the offshore wind industry, the need for validation of reanalysis products is strong, particularly at typical turbine hub heights. Pacific Northwest National Laboratory (PNNL) operates two AXYS WindSentinelTM buoys for the U.S. Department of Energy (DOE) in order to collect meteorological and oceanographic data, including hub height wind speed, in areas of interest for offshore wind development. This work uses DOE buoy observations off the coasts of New Jersey and Virginia to evaluate the performance of Modern Era Retrospective Analysis for Research and Applications-2 (MERRA-2), the North American Regional Reanalysis (NARR), the European Center for Medium-Range Weather Forecasts (ECMWF) Reanalysis 5th generation (ERA5), and the National Oceanic and Atmospheric Administration (NOAA) assimilation system Rapid Refresh (RAP). Biases and degrees of correspondence are determined for each reanalysis product in order to provide insights on the performance and uncertainty for long-term wind resource characterization. With the baseline metrics determined, this work proceeds with an investigation into the sources of large deviations between the modeled and observed hub height wind speeds during the East Coast lidar buoy deployments. Consistent sources of reanalysis model error include stably stratified flow conditions, high wind shear, flow parallel to the coastline, and tropical and winter storms. For the near shore New Jersey location, offshore winds (i.e. blowing from land to water) tend to be associated with model overestimation of observed hub height wind speed, while winds parallel to the coastline are correlated with model underestimation of wind speed. Similarly, near the Virginia buoy, during southerly winds (parallel to the coastline) models underestimate the observed wind speed due to coastal upwelling. Reanalysis model wind speeds underestimate the hub height winds with increasing stable atmospheric stratification. Reanalysis model biases are seen when the observed wind speeds approach or exceed typical turbine cut-out speeds (around 20 m s-1). Large model errors during winter and tropical storm events, high wind speed events, and ramp events are often the result of wind speed timing and magnitude offsets. Overall, ERA5 provides the most successful representation of observed offshore hub height wind speeds at the U.S. East Coast buoy locations and is therefore best suited as an input boundary condition for model case studies. ERA5 performs well for a variety of atmospheric phenomena, including storms and sea breezes, however RAP, the next most successful reanalysis model, performs best in capturing the frequency of ramp events.

17 WIND ENERGY↗
DOE OSTI2024

Impacts of Sea‐Level Rise on Coastal Groundwater Table Simulated by an Earth System Model With a Land‐Ocean Coupling Scheme

Abstract Sea‐level rise (SLR) poses a severe threat to the coastal environment through seawater intrusion into freshwater aquifers. The rising groundwater table also exacerbates the risk of pluvial, fluvial, and groundwater flooding in coastal regions. However, current Earth system models (ESMs) commonly ignore the exchanges of water at the land‐ocean interface. To address this gap, we developed a novel land‐ocean hydrologic coupling scheme in a state‐of‐the‐science ESM, the Energy Exascale Earth System Model version 2 (E3SMv2). The new scheme includes the lateral exchange between seawater and groundwater and the vertical infiltration of seawater driven by the SLR‐induced inundation. Simulations were performed with the updated E3SMv2 for the global land‐ocean interface to assess the impacts of SLR on coastal groundwater under a high CO 2 emission scenario. By the middle of this century, seawater infiltration on the inundated areas will be the dominant component in the land‐ocean coupling process, while the lateral subsurface flow exchange will be much smaller. The SLR‐induced seawater infiltration will raise the groundwater levels, enhance evapotranspiration, and increase runoff with distinct spatial patterns globally in the future. Although the coupling process is induced by SLR, we found topography and warming temperature have more control on the coupling impacts, probably due to the relatively modest magnitude of SLR during the selected future period. Overall, our study suggests significant groundwater and seawater exchange at the land‐ocean interface, which needs to be considered in ESMs.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2000

Modeling of Dense Water Production and Salt Transport from Alaskan Coastal Polynyas

The main significance of this paper is that a realistic, three-dimensional, high-resolution primitive equation model has been developed to study the effects of dense water formation in Arctic coastal polynyas. The model includes realistic ambient stratification, realistic bottom topography, and is forced by time-variant surface heat flux, surface salt flux, and time-dependent coastal flow. The salt and heat fluxes, and the surface ice drift, are derived from satellite observations (SSM/I and NSCAT sensors). The model is used to study the stratification, salt transport, and circulation in the vicinity of Barrow Canyon during the 1996/97 winter season. The coastal flow (Alaska coastal current), which is an extension of the Bering Sea throughflow, is formulated in the model using the wind-transport regression. The results show that for the 1996/97 winter the northeastward coastal current exports 13% to 26% of the salt produced by coastal polynyas upstream of Barrow Canyon in 20 to 30 days. The salt export occurs more rapidly during less persistent polynyas. The inclusion of ice-water stress in the model makes the coastal current slightly weaker and much wider due to the combined effects of surface drag and offshore Ekman transport.

Signorini, Sergio R.↗
DOE OSTI2025

Modeling Offshore Wind Farm Performance in Coastal Low-Level Jets Using Coupled Mesoscale-Microscale Large Eddy Simulations

Accurately predicting wind farm reliability under complex offshore atmospheric conditions remains a key challenge, particularly during noncanonical meteorological events such as coastal low-level jets (LLJs). LLJs, characterized by strong nonmonotonic vertical shear and directional veer, depart significantly from the simplified inflow assumptions embedded in conventional design standards, low-fidelity engineering models, and microscale large eddy simulations of the atmospheric boundary layer. In this work, we use the virtual wind farm framework—an exascale, graphics processing unit–accelerated large eddy simulation platform coupled with high-fidelity aeroservoelastic turbine models and advanced mesoscale-microscale coupling via the ExaWind software stack—to investigate turbine responses under realistic LLJ forcing. Simulations are performed over the U.S. North Atlantic offshore domain with the use of meteorological inputs from New York State Energy Research and Development Authority buoy data, focusing on a representative LLJ case impacting the International Energy Agency 15 MW reference turbine. Our results show that LLJs can cause up to 50% power deficits in downstream turbine rows and significantly amplify low-speed shaft and tower loads through nonlinear coupling between complex inflow characteristics and turbine structural dynamics. Two primary mechanisms drive these load amplifications: (1) unique LLJ inflow features—including veer and vertical/lateral shear—and (2) the downstream evolution of the flow under stable thermal stratification, which suppresses turbulence mixing and alters wake recovery. These mechanisms produce streamwise variations in turbine loading not captured by standard hub height–based metrics or existing design load case (DLC) definitions. This study highlights the critical role of rotor-scale flow gradients in driving fatigue and system-level aeroelastic responses, challenging current DLC and control strategies. We advocate the integration of full-flow field, environment-aware wind inputs into load modeling and control algorithms. By leveraging exascale computing to resolve mesoscale-microscale coupling, this work lays the groundwork for next-generation offshore wind turbine design and operation in meteorologically complex marine environments.

17 WIND ENERGY↗
DOE OSTI2025

Disentangling atmospheric, hydrological, and coupling uncertainties in compound flood modeling within a coupled Earth system model

Compound riverine and coastal flooding is usually driven by complex interactions among meteorological, hydrological, and ocean extremes. However, existing efforts to model this phenomenon often do not integrate hydrological processes across atmosphere–land–river–ocean systems, leading to substantial uncertainties that have not been fully examined. To bridge this gap, we leverage the new capabilities of the Energy Exascale Earth System Model (E3SM) that enable a multi-component framework that integrates coastally refined atmospheric, terrestrial, and oceanic components. We evaluate compound uncertainties arising from two-way land–river–ocean coupling in E3SM and track the cascading meteorological and hydrological uncertainties through ensemble simulations over the Delaware River basin and estuary during Hurricane Irene (2011). Our findings highlight the importance of two-way river–ocean coupling to compound flood modeling and demonstrate E3SM's capability in capturing compound flood extent near the coast, with a hit rate over 0.75. Our study shows the growing uncertainties that transition from atmospheric forcings to flood distribution and severity. Furthermore, an analysis based on artificial neural networks is used to assess the roles of hydrological drivers, such as infiltration and soil moisture, in the generation of compound flooding. The response of compound floods to tropical cyclones (TCs) is found to be susceptible to these often overlooked drivers. For instance, the flooded area could increase more than 2-fold (∼2.4) if Hurricane Irene were preceded by an extreme antecedent soil moisture condition (AMC). The results not only support the use of a multi-component framework for interactive flooding processes, but also underscore the necessity of broader definitions of compound flooding that encompass the simultaneous occurrence of intense precipitation, storm surge, and high AMC during TCs.

Feng, Dongyu [Pacific Northwest National Laborator↗
NASA NTRS2018

Enhancing New York City's Resilience to Sea Level Rise and Increased Coastal Flooding

The large population and major economic assets along New York City's extensive waterfront face exposure to sea level rise (SLR) and coastal flooding. The New York City Panel on Climate Change (NPCC), an advisory group of academic and private-sector experts, provides the city with scientific guidance on future climate risks. Here we highlight current NPCC research on sea level rise, coastal flooding, with some of the City's current and planned responses.Twentieth century SLR in New York City (2.8 cm/decade) exceeded the global average (1.2-1.9 cm/decade), underscoring a greater regional risk. In 2015, the NPCC projected a 2080s SLR of 46-99 cm relative to 2000-2004 (25th -75th percentile) at the Battery, with high-end SLR estimates (90th percentile) of 1.9 m by 2100. Growing evidence of potential instability of the West Antarctic Ice Sheet (WAIS) suggests the possibility of even higher future sea levels. We therefore present a new low probability, yet high impact SLR scenario for New York City, which incorporates these ice sheet instabilities. The new SLR projections will be combined with coastal flood return period curves for the 100-year storm flood levels. Related ongoing research examines changes in storm characteristics, surge--sea level rise interactions, and mapping of neighborhoods at risk.Guided by NPCC findings, New York City undertakes programs to strengthen coastal defenses, tailored to specific neighborhood needs. NPCC in collaboration with the city, continues to investigate vulnerability to extreme climate events, such as inland floods and coastal storms, and to develop stronger links with community-based stakeholder groups. New York City's plans to enhance coastal urban resiliency stand as a model for other urban coastal centers as they prepare for climate change.

Gornitz, Vivien↗
NASA NTRS1996

Demonstration of a Real Time Capability to Produce Tidal Heights and Currents for Naval Operational Use: A Cast Study for the West Coast of Africa (Liberia)

This report documents an existing capability to produce operationally relevant products on sea level and currents from a tides/storm surge model for any coastal region around the world within 48 hours from the time of the request. The model is ready for transition to the Naval Oceanographic Office (NAVOCEANO) for potential contingency use anywhere around the world. A recent application to naval operations offshore Liberia illustrates this. Mississippi State University, in collaboration with the University of Colorado and NAVOCEANO, successfully deployed the Colorado University Rapidly Relocatable Nestable Tides and Storm Surge (CURReNTSS) model that predicts sea surface height, tidal currents and storm surge, and provided operational products on tidal sea level and currents in the littoral region off south-western coast of Africa. This report summarizes the results of this collaborative effort in an actual contingency use of the relocatable model, summarizes the lessons learned, and provides recommendations for further evaluation and transition of this modeling capability to operational use.

CASE STUDIES↗
NASA NTRS1992

Numerical simulation of topographically forced mesoscale variability in a well-mixed marine layer

A simple nonlinear numerical model of a well-mixed marine layer is employed to examine topographically forced mesoscale variability off coastal California. This model is utilized to simulate a persistent wind maximum observed near Point Conception during northwesterly winds. The model solves for the two horizontal components of the marine-layer height and the marine-layer wind.

Eddington, Lee W.↗
NASA NTRS1980

Planetary boundary-layer wind model evaluation at a mid-Atlantic coastal site

Detailed measurements of the mean flow and turbulence were made with the use of a micrometeorological facility consisting of an instrumented 76-m tall tower located within a 100-m distance from the Atlantic Ocean at Wallops Island, Virginia. Under moderately strong wind conditions, the popular neutral boundary layer flow model fails to provide an adequate description of the actual flow. In addition to detailed flow information for all wind directions, averages of the important flow parameters used for design such as vertical distribution of mean velocity, turbulence intensities and turbulence integral scales were presented for wind direction sectors with near uniform upstream terrain. Power spectra of the three velocity components for the prevailing northwesterly and southerly winds are discussed.

Tieleman, H. W.↗
NASA NTRSDate not supplied

Exploring Flooded Fraction Prediction through Machine Learning Models Focusing on Medical Infrastructure in the Southeast U.S. Coastal Areas

Rising sea levels due to climate change increasingly threaten medical infrastructure through flooding. This study develops machine learning models to predict flood exposure for 11,508 medical facilities in the southeastern coastal regions of the United States by integrating datasets including meteorological, hydrological, topographic, and geological data, the Natural Risk Index, and historical flood records from NASA, HIFLD, and FEMA. Six regression models, namely Linear Regression, Support Vector Regression, Random Forest, k-Nearest Neighbors, XGBoost, and Artificial Neural Networks, are trained using 16 explanatory variables identified through literature review and correlation analysis. Data preprocessing employs the SMOGN for class imbalance and Winsorization for outliers. Model performance is evaluated using MAE, MSE, and RMSE, with Random Forest and XGBoost models achieving the highest performance (MSE of 2.58e-5 and 3.69e-5, respectively). This multifactorial approach allows the models to capture complex flood-influencing relationships, enhancing adaptability and performance across geographic regions. Future work focuses on expanding across the U.S. and developing a near real-time flood monitoring system.

Jihoon Chung↗
NASA NTRS2009

Combined MODIS/AMSR-E SST Composites for Regional Weather Applications

Recent applications of a high resolution MODlS composite SST product have clearly shown the importance of developing high-resolution SST data sets for coastal applications and modeling. In general, coupling between the oceans and atmospheres has been closely linked to SST gradients and fronts, indicating a need for high resolution SSTs, specifically in the areas of large gradients associated with coastal regions. Thus an accurate determination of SST gradients has become critical for determining the appropriate air-sea coupling and the influence on ocean modeling. Recent research is focused on improving the accuracy and spatial coverage of the current operational MODIS SST composite product provided by the Short-term Prediction Research and Transition (SPORT) project and distributed to the community. GHRSST-PP MODlS data and microwave AMSR-E data are being combined to produce composite data sets for both the West Coast and East Coast of the United States, including the Gulf of Mexico. The use of 1 km MODIS data has explicit advantages over other SST products including its global coverage and high resolution. The AMSR-E data will reduce the latency of the composites. A strategy for utilizing the error characteristics contained in the GHRSST data has been developed. This strategy will include using the error characteristics directly to calculate weights in the SST composites, uncertainty maps based on the composite biases and RMS errors, and latency products calculated in the compositing process. Recent accomplishments include the development of an enhanced compositing approach based on the error-weighted combination of recent clear MODIS SST values, where the error contributions come from measurement error, potential cloud contamination, and data latency sources. Future plans call for the inclusion of AMSR-E SST values with appropriate weights based upon measurement accuracy, MODIS-AMSR-E SST bias, and latency.

Jedlovec, Gary↗
DOE OSTI2023

Multi-decade high-resolution regional hindcasts for wave energy resource characterization in U.S. coastal waters

Long-term, high-resolution, regional wave hindcast datasets were generated using unstructured-grid Simulating WAves Nearshore (SWAN) models for the U.S. coastal waters to support nearshore wave energy development in the U.S. including those bordering U.S. territorial islands. The model domains resolved the entire U.S. exclusive economic zones, with a spatial resolution of approximately 200 m nearshore. The regional SWAN models were driven by the global WAVEWATCH III® model outputs and run for a 42-year period from 1979 to 2020. Extensive model validations were performed using buoy observations and altimeter data. Regional resource characterization was performed based on hindcast data points at 2 km from shore and along the 100 m isobath. Aggregations of wave resource parameters were produced, and spatial and seasonal variations were analyzed for all the regions. Wave resource metrics recommended by international standards, including a 3-hour time series of six resource parameters, hourly frequency- and directionally resolved wave spectra at selected “virtual buoy” locations, and average-annual values of omni-directional wave power, significant wave height, and energy period are publicly disseminated through an Amazon Web Service and a Marine Energy Atlas web application tool to facilitate wave energy research and a wide range of coastal ocean applications.

16 TIDAL AND WAVE POWER↗
DOE OSTI2025

Performance of wind assessment datasets in United States coastal areas

The atmospheric dynamics that occur near the intersection of land and water offer exciting and challenging opportunities for wind energy deployment in coastal locations. New models and tools are continually being developed in support of wind resource assessment, and three recent products are explored in this work for their performance in representing characteristics of the wind resource at coastal locations: the Global Wind Atlas 3 (GWA3), the 2023 National Offshore Wind dataset (NOW-23), and the wind climate simulations that are a component of the Wind Integration National Dataset (WIND) Toolkit Long-Term Ensemble Dataset (WTK-LED Climate). These relatively new products are freely available and user-friendly so that anyone – from a utility-scale developer to a resident or business owner – can evaluate the potential for wind energy generation at their location of interest. The validations in this work provide guidance on the accuracy of wind resource assessments for coastal customers interested in installing small or midsize wind turbines (≤ 1 MW in capacity) to support energy needs at the residential, business, or community scale, such as the island and remotely located participants of the U.S. Department of Energy's Energy Transitions Initiative Partnership Project. At 23 coastal locations across the United States, dataset performance varies according to different evaluation metrics. All three recent datasets tend to overestimate the observed coastal wind resource. GWA3 produces the smallest annual average wind speed relative errors, whereas WTK-LED Climate is in best agreement in terms of representing diurnal wind speed cycles. NOW-23 is the highest performing of the datasets for representing seasonal and interannual trends in the coastal wind resource. While GWA3 and WTK-LED Climate are relatively insensitive to the dataset output heights selected for wind resource assessment at small and midsize wind turbine hub heights (20–60 m), significant variation in the NOW-23 representation of wind shear across the wind profile in the lowest 100 m of the atmosphere leads to notable differences in wind speed estimates according to the dataset output heights selected for evaluation. GWA3 exhibits challenges in the representation of observed wind speed diurnal cycles at small and midsize turbine hub heights, likely due to the dataset's consistent treatment of hourly wind speed trends regardless of altitude.

17 WIND ENERGY↗
DOE OSTI2026

Improving coastal water level estimation by merging nadir-only satellite altimetry data into a hydrodynamic model

Providing robust real time flood warnings is of paramount importance to coastal communities. Although state-of-the-art hydrodynamic models are capable of robustly predicting Coastal Water Levels (CWL), unresolved drivers affecting level fluctuations are often not represented by the model governing equations. This work evaluates a novel method to improve the performance of the ADvanced CIRCulation (ADCIRC) hydrodynamic model by assimilating observations from four nadir-only satellite altimetry missions against a set of National Oceanic and Atmospheric Administration (NOAA) gauge stations located across the entire U.S. East Coast. Two different types of simulations were performed – Open Loop (OL) and Data Assimilation (DA). Five different simulations were performed where four different satellite altimetry observations were assimilated individually and combined with two different scenarios – with and without considering the data quality flags. Results indicate that, despite their limited spatial coverage, merging nadir-only observations into ADCIRC from the newly launched Surface Water and Ocean Topography (SWOT)’s nadir altimeter can improve the model performance at 76% of the gauge locations, whereas Sentinel-6 improves 73% of the total locations, Jason-3 74%, and SARAL 21%. Furthermore, combining observations from SWOT-nadir, Jason-3, and Sentinel-6 can improve the ADCIRC performance at more than 80% of the gauge locations for 107-day simulation. Nadir-only satellite altimetry observations can be useful for improving the model performance even if flagged as “poor quality” near the coast. When the flagged data are disregarded, SWOT can improve ADCIRC at 78%, Sentinel-6 at 73%, Jason-3 at 53%, and SARAL at 21% of the gauge locations. The ability to improve the model simulations largely depends on the availability of a satellite overpass nearby. Therefore, model performance can be further enhanced if satellite observations are available during a storm surge event, stressing the importance of frequent satellite overpasses.

Aafnan Bhuiyan, Soelem↗
DOE OSTI2026

Compound flood simulations indicate rising public exposure to sewage-contaminated waters in Waikīkī, Hawai‘i

Waikīkī, Hawai‘i serves as a valuable case study for flood vulnerability due to its low elevation, present-day exposure to multiple flood sources, and role as the state’s economic engine. This study applies the high-resolution, open-source WRF-Hydro-CUFA (Coastal Urban Flood Application) model to simulate compound flood interactions among pluvial, fluvial, coastal, storm drain-driven, and subsurface processes, and to evaluate how flood characteristics in Waikīkī evolve under projected sea-level rise (SLR). Simulated water levels in an estuarine drainage canal, serving as the primary outlet for gravity-flow drainage from Waikīkī, show strong agreement with observations (r ≥ 0.935, R 2 ≥ 0.874), accurately reproducing flooding observed during recent 5-year and 50-year storm events. Simulations incorporating SLR demonstrate a shift in dominant flood drivers as sea levels rise. In addition, results identify critical thresholds at which elevated estuarine tailwater conditions cause widespread drainage backflow, eventually occurring in the absence of rainfall. These findings suggest that tidally driven flooding will progressively dominate over rainfall-driven events, increasing public exposure to highly contaminated canal and groundwater sources. The results underscore the urgent need to integrate water quality issues and drainage failure into coastal flood management strategies for low-lying urban environments, particularly those reliant on polluted estuarine waterways for stormwater management.

Compound flooding↗
NASA NTRS2021

Hawai‘i Island Disasters: Using NASA Earth Observations to Assess Coastal Flood Risk with Measures of Land Cover Change, Flood Extent, and Vulnerability for Adaptation and Mitigation Planning on Hawai‘i Island

As the County of Hawai‘i faces an increased risk of extreme flooding events, sea-level rise, and other hazards associated with climate change, the need for building geospatial capacity to make better-informed decisions is critical. The County of Hawai‘i and Arizona State University partnered with NASA DEVELOP to complete a macro-scale risk analysis for the island of Hawai‘i analyzing flooding, land cover, vulnerability, and exposure factors using Earth observations and socio-economic data. The team assessed the variation in urban coastal vulnerability around the entire island of Hawai‘i, using satellite imagery of coastal land cover typology from satellite products such as Landsat 8 Operational Land Imager (OLI), Landsat 7 Enhanced Thematic Mapper Plus (ETM+), and Sentinel-1 Synthetic Aperture Radar (SAR). The team made a sharable geodatabase containing datasets modeling vulnerability to coastal flooding as well as the Hawai‘i Flood Risk Toolbox (HiFloRT) which contains multiple tools for the County to map land cover, extreme rainfall and flood extent across the Island. The end products will allow the County of Hawai‘i to establish a protocol and standard framework for the utilization of Earth observations in future planning.

Garren Kalter↗
DOE OSTI2022

Improving tide-estuary representation in MPAS-Ocean (Final Report)

The overarching goal of this research is to produce a skill-assessed tide-estuary representation, including BGC (biogeochemistry) components of the ecosystem, into the MPAS-O global climate modeling system. Specifically, in the original proposal we proposed one-way nesting, into global MPAS-O, of a high-resolution model for TES (tidal and estuary systems) including the coupling of the BGC components. We have completed the proposed tasks during the past 4 years. Specifically, after extensive consultation with the MPAS-O group during the first year, we have decided to go beyond the original proposal and to two-way nest SCHISM into MPAS-O. As a result, SCHISM-CoSiNE (with CoSiNE being the BGC model) is now fully integrated into MPAS-O as a subroutine (more precisely, an analysis member) and can be invoked as MPAS-O is running with two-way exchange of information between the TES and global models (with proper spatial interpolation done on the fly). As originally proposed, we have also developed the new finite-volume (FV) model MPAS-OI (aka ELCIRC-sub) for TES. We have conducted several tests using the coupled MPAS-SCHISM- CoSiNE model and also MPAS-OI for Columbia River, San Francisco Bay and Chesapeake Bay. Four journal publications have been produced so far as the result of the project and 1 is under review (cf. Publication section). With these achievements, we will be ready to conduct fully coupled simulations of TES and global ocean anywhere, as soon as some missing capabilities in MPAS-O (e.g. tides, improved ability to handle the boundary conditions, and streamlining of the global BGC model) are worked out. The significance of the outcome of this project cannot be overstated. It marks the first time, as far as we know, that a global climate model is tightly coupled to a bona fide nearshore and terrestrial model (SCHISM system); in fact SCHISM is the first bona fide compound flooding model that seamlessly covers the coastal ocean, estuaries, rivers and watersheds (Ye et al. 2020; Zhang et al. 2020). The coupled model will reduce the errors and uncertainties of the MPAS-O global climate model system. It also complements and extends on-going E3SM project, as the latter seeks in the long-term to include components required to simulate impacts of sea-level change and storm surge on coastal regions, including wave models and focusing resolution in coastal and storm-track regions, as an important component of the Cryosphere-ocean system. Note that a spectral wave model (the unstructured-grid version of Wave Watch III, aka Wind Wave Model) is already incorporated inside the SCHISM system and is therefore directly accessible in the coupled system.

58 GEOSCIENCES↗