Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “tidal modeling”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

An improved model of the Earth's gravitational field: GEM-T1

Goddard Earth Model T1 (GEM-T1), which was developed from an analysis of direct satellite tracking observations, is the first in a new series of such models. GEM-T1 is complete to degree and order 36. It was developed using consistent reference parameters and extensive earth and ocean tidal models. It was simultaneously solved for gravitational and tidal terms, earth orientation parameters, and the orbital parameters of 580 individual satellite arcs. The solution used only satellite tracking data acquired on 17 different satellites and is predominantly based upon the precise laser data taken by third generation systems. In all, 800,000 observations were used. A major improvement in field accuracy was obtained. For marine geodetic applications, long wavelength geoidal modeling is twice as good as in earlier satellite-only GEM models. Orbit determination accuracy has also been substantially advanced over a wide range of satellites that have been tested.

Marsh, J. G.↗
NASA NTRS1993

Mars dust-driven tides and their impact on the thermosphere

It has been known since the early Mariner 6, 7, and 9 missions that dust loading of the lower atmosphere and the subsequent aerosol heating during dusty periods impacts the martian middle and upper atmospheres. A quantitative measure of this lower atmosphere forcing was obtained by the Viking 1 and 2 landers, from which observed amplitudes of semidiurnal surface-pressure oscillations were correlated with normal-incidence dust optical depths. It appears that the dominant semidiurnal mode is a good indicator of global dust content or mean dust optical depth, especially during dust storm events. A classical tidal model that reproduces the surface pressure oscillations measured by these Viking landers in 1977 was used to calculate tidal amplitudes and phases up to approximately 43 km. These tidal characteristics were calculated for various dust optical depth conditions ranging from typical dusty periods to global dust storm times. Reasonable extrapolations can be made to higher altitudes if one assumes that the vertically propagating tidal modes continue to grow without dissipation or breaking. It is very likely that gravity waves also play an important role in the structure and dynamics of the middle atmosphere of Mars, since the large topographical relief should produce vigorous gravity wave fluxes. Semidiurnal tidal modes, significantly enhanced by lower atmosphere dust-induced heating, may indeed propagate to the Mars thermosphere (approximately less than 100 km) before breaking and generating turbulence. The preferential enhancement of the semidiurnal tides during dust storm onset is primarily due to the elevation of the tidal heating source in a very dusty atmosphere. The (2,2) semidiurnal tidal tidal mode was shown to have the largest variation with dust optical depth, as measured by Viking lander instruments. Also, the (2,2) mode has the largest vertical wavelength of all the semidiurnal tidal modes, and thus is most likely to penetrate into the thermosphere before breaking and to modify the largely in situ solar-driven behavior otherwise expected. The tides may also be partially responsible for determining the height of the martian homopause (approximately 125 km).

Bougher, Stephen W.↗
NASA NTRS2014

The Interior and Orbital Evolution of Charon as Preserved in Its Geologic Record

Pluto and its largest satellite, Charon, currently orbit in a mutually synchronous state; both bodies continuously show the same face to one another. This orbital configuration is a natural end-state for bodies that have undergone tidal dissipation. In order to achieve this state, both bodies would have experienced tidal heating and stress, with the extent of tidal activity controlled by the orbital evolution of Pluto and Charon and by the interior structure and rheology of each body. As the secondary, Charon would have experienced a larger tidal response than Pluto, which may have manifested as observable tectonism. Unfortunately, there are few constraints on the interiors of Pluto and Charon. In addition, the pathway by which Charon came to occupy its present orbital state is uncertain. If Charon's orbit experienced a high-eccentricity phase, as suggested by some orbital evolution models, tidal effects would have likely been more significant. Therefore, we determine the conditions under which Charon could have experienced tidally-driven geologic activity and the extent to which upcoming New Horizons spacecraft observations could be used to constrain Charon's internal structure and orbital evolution. Using plausible interior structure models that include an ocean layer, we find that tidally-driven tensile fractures would likely have formed on Charon if its eccentricity were on the order of 0.01, especially if Charon were orbiting closer to Pluto than at present. Such fractures could display a variety of azimuths near the equator and near the poles, with the range of azimuths in a given region dependent on longitude; east-west-trending fractures should dominate at mid-latitudes. The fracture patterns we predict indicate that Charon's surface geology could provide constraints on the thickness and viscosity of Charon's ice shell at the time of fracture formation.

Tectonics↗
NASA NTRS1981

Inference of dust opacities for the 1977 Martian great dust storms from Viking Lander 1 pressure data

The tidal heating components for the dusty Martian atmosphere are computed based on dust optical parameters estimated from Viking Lander imaging data, and used to compute the variation of the tidal surface pressure components at the Viking Lander sites as a function of season and the total vertical extinction optical depth of the atmosphere. An atmospheric tidal model is used which is based on the inviscid, hydrostatic primitive equations linearized about a motionless basic state the temperature of which varies only with height, and the profiles of the tidal forcing components are computed using a delta-Eddington approximation to the radiative transfer equations. Comparison of the model results with the observed variations of surface pressure and overhead dust opacity at the Viking Lander 1 site reveal that the dust opacities and optical parameters derived from imaging data are roughly representative of the global dust haze necessary to reproduce the observed surface pressure amplitudes, with the exception of the model-inferred asymmetry parameter, which is smaller during the onset of a great storm. The observed preferential enhancement of the semidiurnal tide with respect to the diurnal tide during dust storm onset is shown to be due primarily to the elevation of the tidal heating source in a very dusty atmosphere.

Zurek, R. W.↗
NASA NTRS2019

Consequences of Tidal Dissipation in a Putative Venusian Ocean

The solar tide in an ancient Venusian ocean is simulated using a dedicated numerical tidal model. Simulations with varying ocean depth and rotational periods ranging from −243 to 64 sidereal Earth days are used to calculate the tidal dissipation rates and associated tidal torque. The results show that the tidal dissipation could have varied by more than 5 orders of magnitude, from 0.001 to 780 GW, depending on rotational period and ocean depth. The associated tidal torque is about 2 orders of magnitude below the present day Venusian atmospheric torque, and could change the Venusian daylength by up to 72 days per million years depending on rotation rate. Consequently, an ocean tide on ancient Venus could have had significant effects on the rotational history of the planet. These calculations have implications for the rotational periods of similarly close-in exoplanetary worlds and the location of the inner edge of the liquid water habitable zone.

planets and satellites↗
NASA NTRS2022

Tidal Insights into Rocky and Icy Bodies: An Introduction and Overview

Solid body tides provide key information on the interior structure, evolution, and origin of the planetary bodies. Our Solar system harbours a very diverse population of planetary bodies, including those composed of rock, ice, gas, or a mixture of all. While a rich arsenal of geophysical methods has been developed over several years to infer knowledge about the interior of the Earth, the inventory of tools to investigate the interiors of other Solar-system bodies remains limited. With seismic data only available for the Earth, the Moon, and Mars, geodetic measurements, including the observation of the tidal response, have become especially valuable and therefore, has played an important role in understanding the interior and history of several Solar system bodies. To use tidal response measurements as a means to obtain constraints on the interior structure of planetary bodies, appropriate understanding of the viscoelastic reaction of the materials from which the planets are formed is needed. Here, we review the fundamental aspects of the tidal modeling and the information on the present-day interior properties and evolution of several planets and moons based on studying their tidal response. We begin with an outline of the theory of viscoelasticity and tidal response. Next, we proceed by discussing the information on the tidal response and the inferred structure of Mercury, Venus, Mars and its moons, the Moon, and the largest satellites of giant planets, obtained from the analysis of the data that has been provided by space missions. We also summarise the upcoming possibilities offered by the currently planned missions.

Tides↗
DOE OSTI2023

Turbulence modeling to aid tidal energy resource characterization in the Western Passage, Maine, USA

Numerical models combined with field measurements are regularly used to characterize tidal energy resources at potential energetic sites. However, most existing works only focus on the tidal hydrodynamic characteristics, and turbulence parameters are often not reported because of the lack of high-quality turbulence measurements and the limitations of numerical models in resolving turbulent eddies. In this study, we used FVCOM - a hydrostatic primitive equation (HPE) model - to characterize the tidal energy resource in the Western Passage, Maine, USA, by taking care of the essential macro-scale turbulence properties. We observed an excellent model performance using the Mellor-Yamada Level 2.5 Turbulence Model; estimating the spatial and vertical distribution of the turbulent kinetic energy and intensity added a new perspective to the site ranking for tidal energy converter (TEC) deployments. In addition, we also examined the role of channel geometry and bathymetry, such as headlands and underwater sills, in enhancing turbulent eddies around potential TEC siting locations. Ultimately, the detailed analysis of the turbulent flow characteristics has changed the site-ranking results and demonstrated that the regional-scale HPE models could be used for the relative understanding of more or less turbulent sites for a refined resource assessment.

16 TIDAL AND WAVE POWER↗
NASA NTRS2019

Consequences of Tidal Dissipation in a Putative Venusian Ocean

The solar tide in an ancient Venusian ocean is simulated using a dedicated numerical tidal model. Simulations with varying ocean depth and rotational periods ranging from minus 243 to 64 sidereal Earth days are used to calculate the tidal dissipation rates and associated tidal torque. The results show that the tidal dissipation could have varied by more than 5 orders of magnitude, from 0.001 to 780 gigawatts (GW), depending on rotational period and ocean depth. The associated tidal torque is about 2 orders of magnitude below the present day Venusian atmospheric torque, and could change the Venusian daylength by up to 72 days per million years depending on rotation rate. Consequently, an ocean tide on ancient Venus could have had significant effects on the rotational history of the planet. These calculations have implications for the rotational periods of similarly close-in exoplanetary worlds and the location of the inner edge of the liquid water habitable zone.

Dynamical evolution and stability↗
NASA NTRS2002

Dynamical Studies of the Middle Atmosphere Using High Resolution Doppler Imager Observations

This report summarizes the activities of NASA grant NAG5-11068, "Dynamicai Studies of the Middle Atmosphere Using High Resolution Doppler Imager Observations." The High Resolution Doppler Imager (HRDI) on the Upper Atmosphere Research Satellite (UARS) has been providing direct measurements of the Earth's horizontal wind field in the stratosphere, mesosphere and lower thermosphere. Mesospheric temperatures, ozone, and O((sup 1) D) densities, and stratospheric aerosol extinctions coefficients, are also retrieved. The goal of HRDI is to measure the vector winds in the stratosphere (10-40 km), mesosphere, and lower thermosphere (approximately 50-120 km) during the day, and the lower thermosphere at night (approximately 95 km) to an accuracy of 5 m/s. The horizontal wind vector is measured by observing the Doppler shift of rotational lines of molecular oxygen along two lines of sight. In addition to winds, temperatures and volume emission rates are determined in the mesosphere and lower thermosphere, from which ozone and O((sup 1) D) concentrations can be derived, and aerosol scattering coefficients are determined in the stratosphere. UARS was launched on September 12, 1991, into a 585-km circular orbit inclined 57 degrees to the equator HRDI was activated September 28, 1991 and following a period of checkout and adjustment of the instrument parameters, scientific observations began November 2, 199 1. HRDI operated nearly continuously from launch until April 1995. At that time the UARS solar array drive failed, forcing the instruments to time-share the available power. From July 1995 to July 1996 HRDI operated approximately 50% of the time. At that point, one of the three spacecraft batteries failed and from then until September 1998 the duty cycle was less than 20% per month, At that time it was determined that HRDI could operate during each daytime pass, which increased the daytime duty cycle to close to l00%, while nighttime operations were limited to about a week per month. In the fall of 1999, the second tape recorder failed requiring a real time contact with a TRDSS satellite to retrieve that data. This resulted in about 60% data collection efficiency. Finally, in the summer of 2000, the second star sensor failed requiring the spacecraft attitude to be controlled by a three axis magnetometer and sun sensor. This resulted in a loss of attitude knowledge but operations continue with the anticipation of correcting the attitude. A new method for determining the tide and mean structure from satellite data in conjunction with a new tidal model has been devised. For brevity, it shall be referred to as the TMAT or Tide-Mean Assimilation Technique. Most previous methods of tidal analysis are based on various ways of slicing the data set.

Skinner, Wilbert↗
NASA NTRS2022

Vertical Coupling by Solar Semidiurnal Tides in the Thermosphere From ICON/MIGHTI Measurements

Wind measurements from the Michelson Interferometer for Global High-resolution Ther-23mospheric Imaging (MIGHTI) instrument on the Ionospheric CONnections (ICON) mis-24sion provide new insights into the semidiurnal tidal spectrum in the thermosphere, cov-25ering latitudes 9◦S-39◦N and altitudes 100-280 km altitude throughout 2020. Latitude26versus day of year (DOY) variability of solar semidiurnal tides SE2, S0, SW1, SW2, SW327and SW4 at 250 km are presented for the first time, and evaluated relative to similar re-28sults at 106 km. Using daytime-only data, height versus latitude and height versus DOY29variability of SE2, S0, SW1. SW3 and SW4 amplitudes and phases are depicted for the30first time, revealing the effects of a dissipative thermosphere on the vertical evolutions31of these tidal structures. SW2 is absent from these depictions due to potential aliasing32by zonal mean winds. The above results are considered in light of the Climatological Tidal33Model of the Thermosphere (CTMT), which is based on fits to tidal winds and temper-34atures from the Thermosphere-Ionosphere-Mesosphere Energetics and Dynamics (TIMED)35mission between 80 and 120 km during 2002-2008, and extrapolated to an altitude of 40036km based on modeled tidal structures propagating in a dissipative thermosphere, but with-37out in-situ sources of excitation due to tide-tide or tide-ion drag nonlinear interactions.38On the basis of comparisons with the CTMT and other characteristics revealed in the39MIGHTI tidal structures, it is concluded that in-situ sources exist for S0, SW1, SW2 and40SW3 in the thermosphere above about 200 km.

Jeffrey M Forbes↗
NASA NTRS1995

A numerical model of nonmigrating diurnal tides between the surface and 65 km

Observations of surface pressure and middle atmosphere temperatures and winds indicate that a substantial nonmigrating component is present in the diurnal tide. The nonmigrating tides, which propagate with a zonal phase speed that is different from the earth's rotation, are attributed to the diurnal heating of geographically fixed sources. In this study we utilize a classical tidal model to examine the propagation characteristics of diurnal tides. The global fields of tropospheric sensible, radiative, and latent heating used to drive the model are supplied from summer and winter diurnal climatologies of the National Center for Atmospheric Research (NCAR) Community Climate Model (CCM2). A novel aspect of this study is the focus on the relative importance of the nonmigrating components. The classical model successfully reproduces many observed features of the low-latitude diurnal surface pressure tides. In the middle atmosphere, the simulated migrating (or sun-synchronous) tide shows qualitative agreement with November-March Limb Infrared Monitor of the Stratosphere (LIMS) observations. Tropospheric solar heating is clearly the dominant driving force for the migrating tide, with secondary contributions from boundary-layer sensible heating and tropospheric latent heat release. The leading modes of the zonal mean tide are also driven chiefly by tropospheric solar heating. The higher-order modes of the zonal mean and eastward propagating tides may be attributed to the joint effects of tropospheric solar heating, sensible heating, and latent heat release. The LIMS and other data reveal features that cannot be explained or examined within the context of the classical model used in the present study. These include upward phase propagation, vertical attenuation, and temporal variations in the migrating diurnal tide.

Lieberman, Ruth S.↗
NASA NTRS1979

Tides and the midnight temperature anomally in the thermosphere

Spectral analysis of the midnight temperature anomaly observed on Atmosphere Explorer-E indicates that the phenomenon can be adequately represented with wave numbers between 1 and 5. Results from a multiconstituent circulation model suggest that this feature and associated variations in composition are the result of a complex interaction between the semidiurnal tide generated in the lower atmosphere and ion-neutral momentum coupling associated with the diurnal variation of the ion density which, in turn, also generates higher order tidal models.

Mayr, H. G.↗
DOE OSTI2025

Marine Toxicity, Biodegradability, and Rolling-Sliding Tribological Performance of Ionic Liquid-Enhanced Environmentally-Acceptable Lubricants for Tidal Turbomachinery

Environmentally acceptable lubricants (EALs) are increasingly being recognized in many fields including waterpower, hydraulics, water transport, agricultural machinery, offshore wind turbines, etc. Specifically, high-performance EALs are demanded for tidal turbomachinery to ensure high efficiency and reliability and avoid the significant risk of direct contamination of the marine ecosystem upon leaks. Here we report a new development of ionic liquid (IL)-enhanced EALs for tidal energy. One short-chain phosphonium phosphate and one short-chain ammonium phosphate ILs were used as the candidate additives and the IL-containing EALs demonstrated significantly improved lubricating performance, much lower toxicity, and increased biodegradability compared with commercial baselines. Specifically, in a rolling-sliding test simulating the operation of a model tidal turbine gearbox bearing, an EAL containing the ILs at a 0.5 wt % concentration demonstrated 40% lower friction, 45% less wear loss, substantially reduced rolling contact fatigue-induced surface damage, and one order of magnitude lower vibration noise compared with a commercial gear oil. In an EPA standard toxicity test, 90 and 70% survival of marine biota was observed when exposed to an EAL containing 5 wt % of the short-chain phosphonium phosphate and ammonium phosphate ILs, respectively, while the selected commercial gear oil and bioderived additive killed all marine biota. In a standard biodegradability test, 2 wt % addition of the phosphonium phosphate IL not only retained the EAL’s ready biodegradability but further boosted the oil decomposition from a range of 60–80% to a higher level of 80–95%. Conversely, adding the commercial bioderived additive downgraded the EAL from readily to inherently biodegradable. Furthermore, this work offers scientific insights for development of ILs as potential EAL additives for marine energy and broader applications.

additives↗
NASA NTRS2012

Simulation Study of a Follow-on Gravity Mission to GRACE

The gravity recovery and climate experiment (GRACE) has been providing monthly estimates of the Earth's time-variable gravity field since its launch in March 2002. The GRACE gravity estimates are used to study temporal mass variations on global and regional scales, which are largely caused by a redistribution of water mass in the Earth system. The accuracy of the GRACE gravity fields are primarily limited by the satellite-to-satellite range-rate measurement noise, accelerometer errors, attitude errors, orbit errors, and temporal aliasing caused by unmodeled high-frequency variations in the gravity signal. Recent work by Ball Aerospace and Technologies Corp., Boulder, CO has resulted in the successful development of an interferometric laser ranging system to specifically address the limitations of the K-band microwave ranging system that provides the satellite-to-satellite measurements for the GRACE mission. Full numerical simulations are performed for several possible configurations of a GRACE Follow-On (GFO) mission to determine if a future satellite gravity recovery mission equipped with a laser ranging system will provide better estimates of time-variable gravity, thus benefiting many areas of Earth systems research. The laser ranging system improves the range-rate measurement precision to approximately 0.6 nm/s as compared to approx. 0.2 micro-seconds for the GRACE K-band microwave ranging instrument. Four different mission scenarios are simulated to investigate the effect of the better instrument at two different altitudes. The first pair of simulated missions is flown at GRACE altitude (approx. 480 km) assuming on-board accelerometers with the same noise characteristics as those currently used for GRACE. The second pair of missions is flown at an altitude of approx. 250 km which requires a drag-free system to prevent satellite re-entry. In addition to allowing a lower satellite altitude, the drag-free system also reduces the errors associated with the accelerometer. All simulated mission scenarios assume a two satellite co-orbiting pair similar to GRACE in a near-polar, near-circular orbit. A method for local time variable gravity recovery through mass concentration blocks (mascons) is used to form simulated gravity estimates for Greenland and the Amazon region for three GFO configurations and GRACE. Simulation results show that the increased precision of the laser does not improve gravity estimation when flown with on-board accelerometers at the same altitude and spacecraft separation as GRACE, even when time-varying background models are not included. This study also shows that only modest improvement is realized for the best-case scenario (laser, low-altitude, drag-free) as compared to GRACE due to temporal aliasing errors. These errors are caused by high-frequency variations in the hydrology signal and imperfections in the atmospheric, oceanographic, and tidal models which are used to remove unwanted signal. This work concludes that applying the updated technologies alone will not immediately advance the accuracy of the gravity estimates. If the scientific objectives of a GFO mission require more accurate gravity estimates, then future work should focus on improvements in the geophysical models, and ways in which the mission design or data processing could reduce the effects of temporal aliasing.

Loomis, Bryant D.↗
NASA NTRS1995

Atmospheric response in aurora experiment: Observations of E and F region neutral winds in a region of postmidnight diffuse aurora

The goal of the Atmospheric Response in Aurora (ARIA) experiment carried out at Poker Flat, Alaska, on March 3, 1992, was to determine the response of the neutral atmosphere to the long-lived, large-scale forcing that is characteristic of the diffuse aurora in the post midnight sector. A combination of chemical release rocket wind measurements, instrumented rocket composition measurements, and ground-based optical measurements were used to characterize the response of the neutral atmosphere. The rocket measurements were made at the end of a 90-min period of strong Joule heating. We focus on the neutral wind measurements made with the rocket. The forcing was determined by running the assimilated mapping of ionospheric electrodynamics (AMIE) analysis procedure developed at the National Center for Atmospheric Research. The winds expected at the latitude and longitude of the experiment were calculated using the spectral thermospheric general circulation model developed at the Danish Meteorological Institute. Comparisons of the observations and the model suggest that the neutral winds responded strongly in two height ranges. An eastward wind perturbation of approximately 100 m/s developed between 140 and 200 km altitude with a peak near 160 km. A southwestward wind with peak magnitude of approximately 150 m/s developed near 115 km altitude. The large amplitude winds at the lower altitude are particularly surprising. They appear to be associated with the upward propagating semidiurnal tide. However, the amplitude is much larger than predicted by any of the tidal models, and the shear found just below the peak in the winds was nominally unstable with a Richardson number of approximately 0.08.

Larsen, M. F.↗
NASA NTRS1989

On the diurnal variation of noctilucent clouds

The relationship between nuctilucent clouds (NLC), which are observed from the ground usually within 2 or 3 hrs of local midnight, and the polar mesospheric clouds (PMC), which are observed by satellites in full daylight, as well as the reason for the differences in their optical properties and their observed heights are investigated. Based on a suggestion that these differences can be attributed to a diurnal variation in the properties of a single type of cloud, two published models of the diurnal and semidiurnal variations of temperature, vertical wind speed, and eddy diffusion coefficient at high latitude to simulate the evolution of ice clouds over a 24-hr period. The results show that the minimum in temperature at about 2000 hours LT causes a sharp maximum in scattered brightness to occur about 1 hour before local midnight, with up to a factor of 7 variation in cloud brightness between noon and midnight. It is noted, however, that considerable uncertainties exist in these tidal models.

Jensen, Eric↗
NASA NTRS1993

A Simulation Atlas of Tidal Features in Galaxies

Detailed simulations of tidally induced structure in disk galaxies have either concentrated on specific systems or consisted of a few encounters with relatively small numbers of particles and no self-gravity. Observers need a 'dictionary' of simulations that covers many encounter parameters with fine morphological resolution and includes effects of self-gravitation. Observers can then search the dictionary for the parameters that best match a particular observed morphology. Alternatively, the dictionary can be used with observational samples for statistical studies of system parameters. To fill this need, we present a survey of model tidal encounters using a self-gravitating, 180,000 particle, two-component ('stars' and 'gas') disk. A wide variety of fascinating morphologies results. There are 86 different encounters that vary orbit tilt, perigalacticon distance, galaxy to companion mass ratio, and the amount of halo dark matter relative to the disk. For morphological comparisons, over 1700 images of the entire survey are available in video form. While there is a rich variety of tidal structure covering much of this parameter space, some general patterns may be remarked. There is a strong orbital inclination dependence of the symmetry of tidal patterns, most symmetric for planar orbits and nearly one-sided for polar encounters. Retrograde encounters produce only broad fanlike global patterns, but rich small-scale internal structure. In both kinds of encounter, our numerical resolution allows us to track internal spiral structure driven by the outer material arms, especially in the lighter halo simulations. We note also that polar encounters generate series of expanding, essentially non-rotating loops resembling shell structures in some respects.

Howard, Sethanne↗
DOE OSTI2023

LANL Institutional Computing Close-out Report for Project t22_ocean_time_step

This year, my team made a very productive use of LANL IC time. Four papers were published that used IC resources, and two more are under review. These publications fall into three categories: 1.) Improving tide modeling in the global ocean by adding self attraction and loading (Barton et. al. 2022), ice shelf cavities (Pal et al. 2023), and local time stepping (Lilly et al. 2023). 2.) A new sea ice numerical formulation (Capodaglio 2023). 3.) Performance comparisons, methods, and test cases for ocean model development, This includes a verification suite for ocean models (Bishnu, submitted) and a comparison between Julia and Fortran (Strauss, submitted). A new time-stepping method was introduced in Calandrini et al. (2022). These publications, and the use of LANL IC resources, go hand-in-hand with our mentoring efforts to train young scientists. Two of the lead authors are graduate students who are conducting their PhD research: Kristin Barton at the University of Michigan and Jeremy Lilly at Oregon State. These are both their first publications, and they both have DOE funding and DOE mentors. In addition, Bishnu is a post-doctoral researcher at LANL; Strauss conducted his research as a senior in high school; and Capodaglio, Calandrini and Pal are all early-career scientists who were converted to staff in 2021 or 2022. Here we highlight two publications on improvements in tidal modeling: Barton et al. 2022 and Lilly et al. 2023.

58 GEOSCIENCES↗