Search NASA⌕ Search

SEARCH · Search NASA

Results for “Tracking error”

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 397 records · Page 22

Automatic signal decoding and sensor stability of a 3-electrode mixed-potential sensor for NO x /NH 3 quantification

Sensors to detect mixtures of NO x /NH 3 are needed to monitor emissions of diesel automobiles where a selective catalytic reduction system uses an NH 3 mediated reaction to reduce NO x . We report on the application of a three electrode La 0.8 Sr 0.2 CrO 3 , Au 0.5 Pd 0.5 , Pt mixed potential sensor using yttria-stabilized-zirconia (YSZ) as a solid electrolyte to NO x /NH 3 sensing. Artificial neural networks were used to automatically decode the concentrations of NO x /NH 3 and errors of less than 15% are achieved. The optimal architecture for ANN decoding and the maximum density of training data points are also determined. The stability of the sensor was monitored by electrochemical impedance spectroscopy. The impedance associated with YSZ oxygen ion conduction and the electrochemical reactions at the three-phase interface are tracked for a period of over 100 days.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Supervisory Control and Data Acquisition for Electrochemical Separation Experimentation

The Python-based program is a laboratory automation tool designed to control and monitor electrochemical systems. The tool was developed for capacitive deionization (CDI) experiments, but it can be used for any system that requires controlled voltage or current segments and multi-parameter monitoring. The program integrates hardware components to run user-defined experimental parameters, providing operational control of a programmable power supply, peristaltic pump, and data acquisition devices. Currently, the program is structured with a workflow that includes an initialization (or pre-run) phase, a main loop, and a post-experiment stabilization (or post-run) phase. The initialization phase prepares and stabilizes the cell, ensuring that the electrodes and solution reach a baseline state before the experiment begins. The main loop consists of multiple voltage segments that repeat, controlling the experiment while recording key parameters such as time, voltage, current, pH, and conductivity. Finally, the post-experiment stabilization phase allows the system to stabilize after the experiment, returning the cell and solution to equilibrium conditions before ending the sequence. The program is designed with four variations, each tailored to different experimental needs. All variations include both the initialization and post-experiment stabilization stages, which run for a set amount of time, voltage, current, and flow rate before and after the main experiment block. The main loop runs for a set number of cycles, as defined by the user input, and each cycle is composed of 2 or 4 segments. The 4 program variations are described as follows: Program 1: The main program includes 2 segments. Each segment is defined to have a set duration, flow rate, voltage, and current. This program measures conductivity, flow rate, voltage, and current. Program 2: The main program expands Program 1 to include 4 segments. Each segment has a specified duration, flow rate, voltage, and current. Like Program 1, it measures conductivity, flow rate, voltage, and current. Program 3: The main program consists of 2 segments, each defined by time, flow rate, voltage, and current. In addition to conductivity, flow rate, voltage, and current, Program 3 collects pH and temperature data through a 4-channel data acquisition device. Program 4: This program independently controls two channels of a multi-channel power supply simultaneously. While conductivity can only be measured for one cell at a time, the dual-channel control makes it possible to operate two cells simultaneously under different voltage/current conditions. The main program includes 2 segments.For each program, all measurements are automatically logged and integrated into a single Excel output file. Data are displayed in numerical format and plotted, both in real time, to track system performance. A key feature of the program is its ability to synchronize all outputs so that every measurement shares a single timestamp, ensuring accurate alignment of voltage, current, pH, conductivity, and pH data.By combining hardware control, real-time monitoring, and unified data collection, this program significantly reduces manual workload and minimizes errors, making it a reliable platform for researchers, engineers, and laboratory technicians conducting CDI experiments, among other electrochemical tests.

Valentino, Lauren [Argonne National Laboratory (AN↗

Technical Note: Intrafractional changes in time lag relationship between anterior–posterior external and superior–inferior internal motion signals in abdominal tumor sites

Purpose: To investigate constancy, within a treatment session, of the time lag relationship between implanted markers in abdominal tumors and an external motion surrogate. Methods: Six gastroesophageal junction and three pancreatic cancer patients (IRB-approved protocol) received two cone-beam CTs (CBCT), one before and one after treatment. Time between scans was less than 30 min. Each patient had at least one implanted fiducial marker near the tumor. In all scans, abdominal displacement (Varian RPM) was recorded as the external motion signal. Purpose-built software tracked fiducials, representing internal signal, in CBCT projection images. Time lag between superior–inferior (SI) internal and anterior–posterior external signals was found by maximizing the correlation coefficient in each breathing cycle and averaging over all cycles. Time-lag-induced discrepancy between internal SI position and that predicted from the external signal (external prediction error) was also calculated. Results: Mean ± standard deviation time lag, over all scans and patients, was 0.10 ± 0.07 s (range 0.01–0.36 s). External signal lagged the internal in 17/18 scans. Change in time lag between pre- and post-treatment CBCT was 0.06 ± 0.07 s (range 0.01–0.22 s), corresponding to 3.1% ± 3.7% (range 0.6%–10.8%) of gate width (range 1.6–3.1 s). In only one patient, change in time lag exceeded 10% of the gate width. External prediction error over all scans of all patients varied from 0.1 ± 0.1 to 1.6 ± 0.4 mm. Conclusions: Time lag between internal motion along SI and external signals is small compared to the treatment gate width of abdominal patients examined in this study. Change in time lag within a treatment session, inferred from pre- to post-treatment measurements is also small, suggesting that a single measurement of time lag at the session start is adequate. These findings require confirmation in a larger number of patients.

60 APPLIED LIFE SCIENCES↗

Applying Particle Swarm Optimization and Extended Kalman Filtering to Model Kaplan Generation Dynamics for Hydropower Systems

Variable renewable generation is increasing the need for hydropower plants to provide fast and flexible grid support, which places new demands on plant-level dynamic models used for monitoring, control, and operational decision-making. This need is especially important for hydroelectric systems, where turbine and generator dynamics are strongly coupled, nonlinear, and time-varying, making accurate real-time representation difficult. To address this problem, this paper develops a digital twin (DT) framework for a synchronous generator–Kaplan turbine system using an explicit separation of slow turbine dynamics and fast generator dynamics. The turbine subsystem is represented by a six-coefficient model, whose parameters are identified offline using particle swarm optimization, while the generator subsystem is updated online through an extended Kalman filter for real-time state and parameter estimation. These models are integrated within a closed-loop simulation that includes a proportional–integral–derivative–double-derivative governor and excitation system, allowing the DT to track plant behavior under realistic operating conditions. Unlike prior studies that treat turbine and generator modeling separately or rely mainly on simulated inputs, the proposed framework is validated using real operational data from a hydropower plant. Results show that the DT reproduces terminal voltage, active power, and reactive power with a normalized root mean square error of approximately 5%. This hybrid offline–online formulation constitutes the main contribution of the work, providing an adaptive and practically deployable DT for hydropower systems with direct relevance to control improvement, performance monitoring, and grid-support applications under high renewable penetration.

13 HYDRO ENERGY↗

MO-FG-CAMPUS-JeP3-04: Feasibility Study of Real-Time Ultrasound Monitoring for Abdominal Stereotactic Body Radiation Therapy

Purpose: Ultrasound is ideal for real-time monitoring in radiotherapy with high soft tissue contrast, non-ionization, portability, and cost effectiveness. Few studies investigated clinical application of real-time ultrasound monitoring for abdominal stereotactic body radiation therapy (SBRT). This study aims to demonstrate the feasibility of real-time monitoring of 3D target motion using 4D ultrasound. Methods: An ultrasound probe holding system was designed to allow clinician to freely move and lock ultrasound probe. For phantom study, an abdominal ultrasound phantom was secured on a 2D programmable respiratory motion stage. One side of the stage was elevated than another side to generate 3D motion. The motion stage made periodic breath-hold movement. Phantom movement tracked by infrared camera was considered as ground truth. For volunteer study three healthy subjects underwent the same setup for abdominal SBRT with active breath control (ABC). 4D ultrasound B-mode images were acquired for both phantom and volunteers for real-time monitoring. 10 breath-hold cycles were monitored for each experiment. For phantom, the target motion tracked by ultrasound was compared with motion tracked by infrared camera. For healthy volunteers, the reproducibility of ABC breath-hold was evaluated. Results: Volunteer study showed the ultrasound system fitted well to the clinical SBRT setup. The reproducibility for 10 breath-holds is less than 2 mm in three directions for all three volunteers. For phantom study the motion between inspiration and expiration captured by camera (ground truth) is 2.35±0.02 mm, 1.28±0.04 mm, 8.85±0.03 mm in LR, AP, SI directly, respectively. The motion monitored by ultrasound is 2.21±0.07 mm, 1.32±0.12mm, 9.10±0.08mm, respectively. The motion monitoring error in any direction is less than 0.5 mm. Conclusion: The volunteer study proved the clinical feasibility of real-time ultrasound monitoring for abdominal SBRT. The phantom and volunteer ABC studies demonstrated sub-millimeter accuracy of 3D motion movement monitoring.

60 APPLIED LIFE SCIENCES↗

Multivariate analysis for the estimation of target localization errors in fiducial marker-based radiotherapy

Purpose: To assess the target localization error (TLE) in terms of the distance between the target and the localization point estimated from the surrogates (|TMD|), the average of respiratory motion for the surrogates and the target (|aRM|), and the number of fiducial markers used for estimating the target (n). Methods: This study enrolled 17 lung cancer patients who subsequently underwent four fractions of real-time tumor tracking irradiation. Four or five fiducial markers were implanted around the lung tumor. The three-dimensional (3D) distance between the tumor and markers was at maximum 58.7 mm. One of the markers was used as the target (P{sub t}), and those markers with a 3D |TMD{sub n}| ≤ 58.7 mm at end-exhalation were then selected. The estimated target position (P{sub e}) was calculated from a localization point consisting of one to three markers except P{sub t}. Respiratory motion for P{sub t} and P{sub e} was defined as the root mean square of each displacement, and |aRM| was calculated from the mean value. TLE was defined as the root mean square of each difference between P{sub t} and P{sub e} during the monitoring of each fraction. These procedures were performed repeatedly using the remaining markers. To provide the best guidance on the answer with n and |TMD|, fiducial markers with a 3D |aRM ≥ 10 mm were selected. Finally, a total of 205, 282, and 76 TLEs that fulfilled the 3D |TMD| and 3D |aRM| criteria were obtained for n = 1, 2, and 3, respectively. Multiple regression analysis (MRA) was used to evaluate TLE as a function of |TMD| and |aRM| in each n. Results: |TMD| for n = 1 was larger than that for n = 3. Moreover, |aRM| was almost constant for all n, indicating a similar scale for the marker’s motion near the lung tumor. MRA showed that |aRM| in the left–right direction was the major cause of TLE; however, the contribution made little difference to the 3D TLE because of the small amount of motion in the left–right direction. The TLE calculated from the MRA ({sup MRA}TLE) increased as |TMD| and |aRM| increased and adversely decreased with each increment of n. The median 3D {sup MRA}TLE was 2.0 mm (range, 0.6–4.3 mm) for n = 1, 1.8 mm (range, 0.4–4.0 mm) for n = 2, and 1.6 mm (range, 0.3–3.7 mm) for n = 3. Although statistical significance between n = 1 and n = 3 was observed in all directions, the absolute average difference and the standard deviation of the {sup MRA}TLE between n = 1 and n = 3 were 0.5 and 0.2 mm, respectively. Conclusions: A large |TMD| and |aRM| increased the differences in TLE between each n; however, the difference in 3D {sup MRA}TLEs was, at most, 0.6 mm. Thus, the authors conclude that it is acceptable to continue fiducial marker-based radiotherapy as long as |TMD| is maintained at ≤58.7 mm for a 3D |aRM| ≥ 10 mm.

60 APPLIED LIFE SCIENCES↗

Ptychographic X-ray speckle tracking with multi-layer Laue lens systems

The ever-increasing brightness of synchrotron radiation sources demands improved X-ray optics to utilize their capability for imaging and probing biological cells, nano-devices and functional matter on the nanometre scale with chemical sensitivity. Hard X-rays are ideal for high-resolution imaging and spectroscopic applications owing to their short wavelength, high penetrating power and chemical sensitivity. The penetrating power that makes X-rays useful for imaging also makes focusing them technologically challenging. Recent developments in layer deposition techniques have enabled the fabrication of a series of highly focusing X-ray lenses, known as wedged multi-layer Laue lenses. Improvements to the lens design and fabrication technique demand an accurate, robust, in situ and at-wavelength characterization method. To this end, a modified form of the speckle tracking wavefront metrology method has been developed. The ptychographic X-ray speckle tracking method is capable of operating with highly divergent wavefields. A useful by-product of this method is that it also provides high-resolution and aberration-free projection images of extended specimens. Three separate experiments using this method are reported, where the ray path angles have been resolved to within 4 nrad with an imaging resolution of 45 nm (full period). This method does not require a high degree of coherence, making it suitable for laboratory-based X-ray sources. Likewise, it is robust to errors in the registered sample positions, making it suitable for X-ray free-electron laser facilities, where beam-pointing fluctuations can be problematic for wavefront metrology.

36 MATERIALS SCIENCE↗

A mesoscopic link-transmission-model able to track individual vehicles

Macroscopic traffic flow is a common choice for large-scale traffic simulations. These models do not provide individual-specific metrics as outputs. However, this treatment is necessary in agent-based-models, as in, for example, assigning routes based on personal characteristics. Here, in this paper, we propose an extension of the link-transmission-model, an efficient and yet accurate discretization of the Lighthill-Whitham-Richards (LWR) model, which allow vehicles to be tracked individually while keeping the main features of the underlying model. The extension comprises modifying the link and node models to ensure that the flow between links is always at discrete levels. Therefore, every unit of flow is associated with one individual vehicle moving from its current to its next link. An upper bound of the discretization error is provided. We show that the proposed model resembles its continuous counterpart on lane drop, merge, and diverge cases. In addition, we apply the model into three different networks to validate its applicability in large networks. Finally, we also confirm the parameter transferability between continuous and discrete models and that both can well reproduce field data.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Topological structure of complex predictions

Abstract Current complex prediction models are the result of fitting deep neural networks, graph convolutional networks or transducers to a set of training data. A key challenge with these models is that they are highly parameterized, which makes describing and interpreting the prediction strategies difficult. We use topological data analysis to transform these complex prediction models into a simplified topological view of the prediction landscape. The result is a map of the predictions that enables inspection of the model results with more specificity than dimensionality-reduction methods such as tSNE and UMAP. The methods scale up to large datasets across different domains. We present a case study of a transformer-based model previously designed to predict expression levels of a piece of DNA in thousands of genomic tracks. When the model is used to study mutations in the BRCA1 gene, our topological analysis shows that it is sensitive to the location of a mutation and the exon structure of BRCA1 in ways that cannot be found with tools based on dimensionality reduction. Moreover, the topological framework offers multiple ways to inspect results, including an error estimate that is more accurate than model uncertainty. Further studies show how these ideas produce useful results in graph-based learning and image classification.

Computer Science↗

Shaded fraction and backtracking in single-axis trackers on rolling terrain

A generalized closed-form equation for the shaded collector fraction in solar arrays on rolling or undulating terrain is provided for single-axis tracking and fixed-tilt systems. The equation accounts for different rotation angles between the shaded and shading trackers, cross-axis slope between the two trackers, and offset between the collector plane and axis of rotation. The validity of the equation is demonstrated through comparison with numerical ray-tracing simulations and remaining minor sources of error are quantified. Additionally, a simple procedure to determine backtracking rotations for each row in an array installed on the rolling terrain (varying in the direction perpendicular to the tracker axes) is provided. The backtracking equation accounts for a desired shaded fraction (including complete shade avoidance) as well as an axis-collector offset. Test cases are provided to facilitate implementation of these equations.

14 SOLAR ENERGY↗

Computational fluid dynamics analysis of char conversion in Sandia’s pressurized entrained flow reactor

Design and analysis of practical reactors utilizing solid feedstocks rely on reaction rate parameters that are typically generated in lab-scale reactors. Evaluation of the reaction rate information often relies on assumptions of uniform temperature, velocity, and species distributions in the reactor, in lieu of detailed measurements that provide local information. This assumption might be a source of substantial error, since reactor designs can impose significant inhomogeneities, leading to data misinterpretation. Spatially resolved reactor simulations help understand the key processes within the reactor and support the identification of severe variations of temperature, velocity, and species distributions. In this work, Sandia’s pressurized entrained flow reactor is modeled to identify inhomogeneities in the reaction zone. Tracer particles are tracked through the reactor to estimate the residence times and burnout ratio of introduced coal char particles in gasifying environments. The results reveal a complex mixing environment for the cool gas and particles entering the reactor along the centerline and the main high-speed hot gas reactor flow. Furthermore, the computational fluid dynamics (CFD) results show that flow asymmetries are introduced through the use of a horizontal gas pre-heating section that connects to the vertical reactor tube. Computed particle temperatures and residence times in the reactor differ substantially from the idealized plug flow conditions typically evoked in interpreting experimental measurements. Furthermore, experimental measurements and CFD analysis of heat flow through porous refractory insulation suggest that for the investigated conditions (1350 °C, <20 atm), the thermal conductivity of the insulation does not increase substantially with increasing pressure.

47 OTHER INSTRUMENTATION↗

Online monitoring of lanthanide species with combined spectroscopy in flowing aqueous aerosol systems

Combined spectroscopic analysis through absorption and laser-induced breakdown spectroscopy (LIBS) was used to monitor Nd and Pr concentrations in a flowing aqueous system. A unique, online sampling approach was employed that allowed for an approximately closed analysis loop of liquids from a reservoir. Absorption spectroscopy was performed with an optical flow cell, and LIBS was performed on an aerosol stream. Multivariate calibrations based on combined absorption and LIBS signals were built for Nd and Pr and then used to monitor concentrations in mixed solutions in a series of spiking tests. In these tests, the concentrations of Nd and Pr in solution were intermittently changed while spectroscopic signals were monitored in real-time. The combined spectroscopic signals and multivariate models were successful in monitoring changing concentrations of lanthanide species with high accuracy and minimal latency. Root-mean squared error of predictions were 0.015 mol L −1 and 0.019 mol L −1 for Nd and Pr respectively, and these lanthanides were able to be monitored to an accuracy of < 0.2 wt%. This work demonstrates both the capabilities of a sampling system for online analysis of liquids and the capabilities of multimodal spectroscopic characterization for real-time, continuous tracking of species in potentially hazardous liquid systems.

absorbance spectroscopy↗

TH-AB-201-07: Filmless Treatment Localization QA for the CyberKnife System

Purpose: Accuray recommends daily evaluation of the treatment localization and delivery systems (TLS/TDS) of the CyberKnife. The vendor-provided solution is a Winston-Lutz-type test that evaluates film shadows from an orthogonal beam pair (known as AQA). Since film-based techniques are inherently inefficient and potentially inconsistent and uncertain, this study explores a method which provides a comparable test with greater efficiency, consistency, and certainty. This test uses the QAStereoChecker (QASC, Standard Imaging, Inc., Middleton, WI), a high-resolution flat-panel detector with coupled fiducial markers for automated alignment. Fiducial tracking is used to achieve high translational and rotational position accuracy. Methods: A plan is generated delivering five circular beams, with varying orientation and angular incidence. Several numeric quantities are calculated for each beam: eccentricity, centroid location, area, major-axis length, minor-axis length, and orientation angle. Baseline values were acquired and repeatability of baselines analyzed. Next, errors were induced in the path calibration of the CK, and the test repeated. A correlative study was performed between the induced errors and quantities measured using the QASC. Based on vendor recommendations, this test should be able to detect a TLS/TDS offset of 0.5mm. Results: Centroid shifts correlated well with induced plane-perpendicular offsets (p < 0.01). Induced vertical shifts correlated best with the absolute average deviation of eccentricities (p < 0.05). The values of these metrics which correlated with the threshold of 0.5mm induced deviation were used as individual pass/fail criteria. These were then used to evaluate induced offsets which shifted the CK in all axes (a clinically-realistic offset), with a total offset of 0.5mm. This test provided high and specificity and sensitivity. Conclusion: From setup to analysis, this filmless TLS/TDS test requires 4 minutes, as opposed to 15–20 minutes for film-based methods. The techniques introduced can potentially isolate errors in individual joints of the CK robot. Spectrum Medical Physics, LLC of Greenville, SC has a consulting contract with Standard Imaging of Middleton, WI.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

In Situ Height and Width Estimation of Sorghum Plants from 2.5d Infrared Images

Plant phenotyping, or the measurement of plant traits such as stem width and plant height, is a critical step in the development and evaluation of higher yield biofuel crops. Phenotyping allows biologists to quantitatively estimate the biomass of plant varieties and therefore their potential for biofuel production. Manual phenotyping is costly, time-consuming, and errorprone, requiring a person to walk through the fields measuring individual plants with a tape measure and notebook. In this work we describe an alternative system consisting of an autonomous robot equipped with two infrared cameras that travels through fields, collecting 2.5D image data of sorghum plants. We develop novel image processing based algorithms to estimate plant height and stem width from the image data. Our proposed method has the advantage of working in situ using images of plants from only one side. This allows phenotypic data to be collected nondestructively throughout the growing cycle, providing biologists with valuable information on crop growth patterns. Our approach first estimates plant heights and stem widths from individual frames. It then uses tracking algorithms to refine these estimates across frames and avoid double counting the same plant in multiple frames. The result is a histogram of stem widths and plant heights for each plot of a particular genetically engineered sorghum variety. In-field testing and comparison with human collected ground truth data demonstrates that our system achieves 13% average absolute error for stem width estimation and 15% average absolute error for plant height estimation.

09 BIOMASS FUELS↗

Computing material volume fractions on a superimposed mesh as applied to Monte Carlo particle transport simulations

Here, we present a newly implemented ray tracing algorithm in OpenMC for efficiently computing material volume fractions on superimposed meshes in complex geometries. By firing rays along each coordinate direction through the geometry, the approach accumulates track-length data in each mesh element, thereby determining the fractional composition of each material. Scaling studies on three different models—a random tetrahedra configuration, the Frascati Neutron Generator ITER dose rate benchmark, and a stellarator design—show excellent parallel performance, with nearly linear speedup on modern multi-threaded and distributed-memory systems. An analysis of the residual error relative to high-resolution reference solutions demonstrated that under optimal conditions it decreases as 1/R, where R is the number of rays fired, making it straightforward to achieve user-prescribed accuracy. This new functionality enables practical, mesh-based approaches for detailed nuclear analyses in production Monte Carlo workflows without resorting to expensive, fully conformal or unstructured meshing.

Monte Carlo↗

Autoencoder Neural Network for chemically reacting systems

Incorporating detailed chemical kinetic models is critical for accurate simulations of reacting flows. However, detailed models involve a large number of thermochemical (TC) state variables. Solving the governing equations to evolve these TC variables becomes impractical for real-world applications. In this work, we propose an autoencoder (AE) neural network (NN)-based reduced model to accelerate such simulations. The AE NN is first trained to find a low-dimensional latent representation of the TC states. Then, the evolving state of a chemical system can be tracked by solving the equations of the latent variables instead of the original TC equations. We demonstrate the reduced model in a syngas CO/H 2 combustion system, using training data collected from canonical perfectly stirred reactors (PSRs). It is found that the AE model can reduce the dimension of the combustion system from 12 to 2 while maintaining low reconstruction error and excellent elemental mass conservation for the test dataset. In the a posteriori test, the combustion states obtained from solving the two latent equations are compared to those from solving the 12 equations of the full model. The AE reduced method is found to be able to capture the diverse combustion states on the top two branches of the S-curve well including the extinction turning point, but with higher prediction errors for states near the ignition turning point.

97 MATHEMATICS AND COMPUTING↗

Numerical Modelling of the Optical Stochastic Cooling Experiment at IOTA

A proof-of-principle optical-stochastic cooling (OSC) experiment is currently in its commissioning phase at the Fermilab’s IOTA ring. In support of this experiment, we recently implemented an OSC element in the ELEGANT tracking program. The model, based on a semi-analytic description of OSC [*], supports the simulation of a large number of macroparticles (10⁴-10⁶) over many turns (10⁶). This paper showcases the simulation capabilities to investigate the beam dynamics in the presence of cooling (or self-interacting radiation field in general) and quantify the impact of various sources of error (e.g. transverse and phase jitter), guide data analysis.

43 PARTICLE ACCELERATORS↗

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↗