Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “global refinement”

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 721 records · Page 40

Seasonal and static Gravity Field of Mars from MGS, Mars Odyssey and MRO Radio Science

We present a spherical harmonic solution of the static gravity field of Mars to degree and order 120, GMM-3, that has been calculated using the Deep Space Network tracking data of the NASA Mars missions, Mars Global Surveyor (MGS), Mars Odyssey (ODY), and the Mars Reconnaissance Orbiter (MRO). We have also jointly determined spherical harmonic solutions for the static and time-variable gravity field of Mars, and the Mars k(sub 2) Love numbers, exclusive of the gravity contribution of the atmosphere. Consequently, the retrieved time-varying gravity coefficients and the Love number k(sub 2) solely yield seasonal variations in the mass of the polar caps and the solid tides of Mars, respectively. We obtain a Mars Love number k(sub 2) of 0.1697 +/- 0.0027 (3- sigma). The inclusion of MRO tracking data results in improved seasonal gravity field coefficients C(sub 30) and, for the first time, C 50. Refinements of the atmospheric model in our orbit determination program have allowed us to monitor the odd zonal harmonic C(sub 30) for approximately 1.5 solar cycles (16 years). This gravity model shows improved correlations with MOLA topography up to 15% larger at higher harmonics ( l = 60-80) than previous solutions.

Genova, Antonio↗

Streaming Solutions for Fine-Grained Network Traffic Measurements and Analysis

Online network traffic measurements and analysis is critical for detecting and preventing any real-time anomalies in the network. We propose, implement, and evaluate an online, adaptive measurement platform, which utilizes real-time traffic analysis results to refine subsequent traffic measurements. Central to our solution is the concept of Multi-Resolution Tiling (MRT), a heuristic approach that performs sequential analysis of traffic data to zoom into traffic subregions of interest. However, MRT is sensitive to transient traffic spikes. In this paper, we propose three novel traffic streaming algorithms that overcome the limitations of MRT and can cater to varying degrees of computational and storage budgets, detection latency, and accuracy of query response. We evaluate our streaming algorithms on a highly parallel and programmable hardware as well as a traditional software-based platforms. The algorithms demonstrate significant accuracy improvement over MRT in detecting anomalies consisting of synthetic hard-to-track elephant flows and global icebergs. Our proposed algorithms maintain the worst-case complexities of the MRT while incurring only a moderate increase in average resource utilization.

97 MATHEMATICS AND COMPUTING↗

ELM-Wet: Inclusion of a wetland landunit with sub-grid representation of eco-hydrological patches and hydrological forcing in E3SM Land Model (ELM)

Wetlands emit the most biogenic methane (CH4) and present the greatest uncertainty in the global CH4 budget. Modeling these emissions is challenging due to the temporal and spatial variability in wetland structure and CH4 flux rates, along with complex interactions among hydrological, ecological, meteorological, and microbial processes that govern CH4 dynamics. To address these issues, we aim to enhance the accuracy of wetland representation in the U.S. Department of Energy’s Exascale Earth System Model (E3SM) Land Model, ELM. This effort led to the development of ELM-Wet, which incorporates a dedicated wetland landunit with subgrid representation of eco-hydrological patch types. We implement wetland-specific hydrology by imposing site-specific constraints on surface water levels, thereby allowing different patches to sustain varying depths of inundation. Additionally, we refined the calculation of aerenchyma transport diffusivity based on observed conductance across different vegetation types. We validated these enhancements through site-specific simulations of a coastal freshwater wetland, the Salvador WMA Freshwater Marsh (Ameriflux, site ID US-LA2), located in the coast of Louisiana (29.85N,90.29W) at an elevation of 0 m. The site was simulated with ELM-Wet and the default version ELMv1. We use Bayesian Optimization to parameterize CO2 and CH4 fluxes. Eddy covariance observations of CO2 and CH4 fluxes from 2012-2013 were used to train the model and data from 2021 were used for validation. In this repository, we include all the output of all the simulations performed using different versions of ELMv1 and ELM-Wet, all input required to the model, and the Matlab scripts we used to processed the output data. ELM-Wet_flmd.csv includes a detailed description of the datasets files.

54 ENVIRONMENTAL SCIENCES↗

Earth Global Reference Atmospheric Model 2007 (Earth-GRAM07) Applications for the NASA Constellation Program

Engineering models of the atmosphere are used extensively by the aerospace community for design issues related to vehicle ascent and descent. The Earth Global Reference Atmosphere Model version 2007 (Earth-GRAM07) is the latest in this series and includes a number of new features. Like previous versions, Earth-GRAM07 provides both mean values and perturbations for density, temperature, pressure, and winds, as well as monthly- and geographically-varying trace constituent concentrations. From 0 km to 27 km, thermodynamics and winds are based on the National Oceanic and Atmospheric Administration Global Upper Air Climatic Atlas (GUACA) climatology. For altitudes between 20 km and 120 km, the model uses data from the Middle Atmosphere Program (MAP). Above 120 km, EarthGRAM07 now provides users with a choice of three thermosphere models: the Marshall Engineering Thermosphere (MET-2007) model; the Jacchia-Bowman 2006 thermosphere model (JB2006); and the Naval Research Labs Mass Spectrometer, Incoherent Scatter Radar Extended Model (NRL MSIS E-OO) with the associated Harmonic Wind Model (HWM-93). In place of these datasets, Earth-GRAM07 has the option of using the new 2006 revised Range Reference Atmosphere (RRA) data, the earlier (1983) RRA data, or the user may also provide their own data as an auxiliary profile. Refinements of the perturbation model are also discussed which include wind shears more similar to those observed at the Kennedy Space Center than the previous version Earth-GRAM99.

Leslie, Fred W.↗

Atmospheric electricity/meteorology analysis

This activity focuses on Lightning Imaging Sensor (LIS)/Lightning Mapper Sensor (LMS) algorithm development and applied research. Specifically we are exploring the relationships between (1) global and regional lightning activity and rainfall, and (2) storm electrical development, physics, and the role of the environment. U.S. composite radar-rainfall maps and ground strike lightning maps are used to understand lightning-rainfall relationships at the regional scale. These observations are then compared to SSM/I brightness temperatures to simulate LIS/TRMM multi-sensor algorithm data sets. These data sets are supplied to the WETNET project archive. WSR88-D (NEXRAD) data are also used as it becomes available. The results of this study allow us to examine the information content from lightning imaging sensors in low-earth and geostationary orbits. Analysis of tropical and U.S. data sets continues. A neural network/sensor fusion algorithm is being refined for objectively associating lightning and rainfall with their parent storm systems. Total lightning data from interferometers are being used in conjunction with data from the national lightning network. A 6-year lightning/rainfall climatology has been assembled for LIS sampling studies.

Goodman, Steven J.↗

Ground-based Network and Supersite Measurements for Studying Aerosol Properties and Aerosol-Cloud Interactions

From radiometric principles, it is expected that the retrieved properties of extensive aerosols and clouds from reflected/emitted measurements by satellite (and/or aircraft) should be consistent with those retrieved from transmitted/emitted radiance observed at the surface. Although space-borne remote sensing observations contain large spatial domain, they are often plagued by contamination of surface signatures. Thus, ground-based in-situ and remote-sensing measurements, where signals come directly from atmospheric constituents, the sun, and the Earth-atmosphere interactions, provide additional information content for comparisons that confirm quantitatively the usefulness of the integrated surface, aircraft, and satellite datasets. The development and deployment of AERONET (AErosol RObotic NETwork) sunphotometer network and SMART-COMMIT (Surface-sensing Measurements for Atmospheric Radiative Transfer - Chemical, Optical & Microphysical Measurements of In-situ Troposphere) mobile supersite are aimed for the optimal utilization of collocated ground-based observations as constraints to yield higher fidelity satellite retrievals and to determine any sampling bias due to target conditions. To characterize the regional natural and anthropogenic aerosols, AERONET is an internationally federated network of unique sunphotometry that contains more than 250 permanent sites worldwide. Since 1993, there are more than 480 million aerosol optical depth observations and about 15 sites have continuous records longer than 10 years for annual/seasonal trend analyses. To quantify the energetics of the surface-atmosphere system and the atmospheric processes, SMART-COMMIT instrument into three categories: flux radiometer, radiance sensor and in-situ probe. Through participation in many satellite remote-sensing/retrieval and validation projects over eight years, SMART-COMMIT have gradually refine( and been proven vital for field deployment. In this paper, we will demonstrate the capability of AERONET SMART-COMMIT in current Asian Monsoon Year-2008 campaigns that are designed and being executed to study the compelling variability in temporal scale of both anthropogenic and natural aerosols (e.g., airborne dust, smoke, mega-city pollutant). Feedback mechanisms between aerosol radiative effects and monsoon dynamics have been recently proposed, however there is a lack of consensus on whether aerosol forcing would be more likely to enhance or reduce the strength of the monsoon circulation. We envision robust approaches which well-collocated ground-based measurements and space-borne observations will greatly advance our understanding of absorbing aerosols (e.g., "Global Dimming" vs. "Elevated Heat-Pump" effects) on aerosol cloud water cycle interactions.

Tsay, Si-Chee↗

Improving the Use of Surface-Sensitive Radiances in the GMAO Hybrid-4DEnVar System

The planetary boundary layer (PBL) was designated as an incubation-class targeted observable (TO) in the 2018 Decadal Survey. As no single instrument source will provide enough information to constrain the global PBL, the assimilation of a wide range of observations in data assimilation (DA) systems will play a critical role, and GMAO has put in efforts to enhance surface-sensitive radiance assimilation. Although a vast number of microwave radiance observations are used in the GEOS global DA system, very few surface-sensitive radiances are currently used over land due to large uncertainty of land surface emissivity in the CRTM as well as cloud detection issue. Dynamically varying emissivity is retrieved from observations of window channels in the GEOS system for non-scattering FOVs and applied to sounding channel assimilation. Moreover, the original radiance bias correction is found not to work well over land, and shows drifting bias when the original emissivity sensitivity bias predictor with a dynamically varying emissivity retrieval is used. Hence, the radiance bias correction has been modified, and the quality control procedure has been adapted accordingly. Cycled experiment results show improvement in the temperature forecast at lower model levels; further refinement of this system continues.

Yanqiu Zhu↗

Absorption of wireless radiation in the child versus adult brain and eye from cell phone conversation or virtual reality

Highlights: • More cell phone radiation is absorbed by children's inner brain tissues than adults’. • Children's radiofrequency radiation exposure should be reduced. • Further research to evaluate the risks to the eye from use of VR is urgently needed. • It is biologically relevant and feasible to reduce the standards’ averaging volume. • Current methods to determine wireless device compliance should be revised. Children's brains are more susceptible to hazardous exposures, and are thought to absorb higher doses of radiation from cell phones in some regions of the brain. Globally the numbers and applications of wireless devices are increasing rapidly, but since 1997 safety testing has relied on a large, homogenous, adult male head phantom to simulate exposures; the “Standard Anthropomorphic Mannequin” (SAM) is used to estimate only whether tissue temperature will be increased by more than 1 Celsius degree in the periphery. The present work employs anatomically based modeling currently used to set standards for surgical and medical devices, that incorporates heterogeneous characteristics of age and anatomy. Modeling of a cell phone held to the ear, or of virtual reality devices in front of the eyes, reveals that young eyes and brains absorb substantially higher local radiation doses than adults’. Age-specific simulations indicate the need to apply refined methods for regulatory compliance testing; and for public education regarding manufacturers' advice to keep phones off the body, and prudent use to limit exposures, particularly to protect the young.

54 ENVIRONMENTAL SCIENCES↗

Artificial Intelligence-Enhanced CMIP6 Climate Projections Across the Conterminous United States

This dataset comprises high-resolution climate projections at 1/24 degree grid (~4km) over the conterminous United States (CONUS) based on ten Global Climate Models (GCMs) that are part of the Coupled Models Intercomparison Project phase 6 (CMIP6). The CMIP6 GCMs are downscaled using two artificial intelligence (AI) techniques, primarily based on the computer vision approach called super-resolution. We train two separate networks: super-resolution convolutional neural network (SRCNN) and super-resolution generative adversarial network (SRGAN). The networks are trained using Daymet observations, originally available at a 1 km resolution. For training purposes, the Daymet data is interpolated to 1/24 degree (~4km), 0.25 degree and 1 degree, which serve as high, intermediate and low-resolution inputs respectively. For each of the SRCNN and SRGAN network, we use a two-step resolution enhancement, the first step generates 4x refinement from 1 degree to 0.25 degree and the second step generates 6x refinement from 0.25 degree to 1/24 degree (~4km). We downscale daily scale precipitation, maximum temperature and minimum temperature for the six CMIP6 GCMs for 1980 to 2019 in the historical period and 2020 to 2059 in the near-term future under the shared socioeconomic pathway 585 and 245 (SSP585 and SSP245) emission scenarios. We also perform double bias-correction with Daymet observations using a quantile mapping approach, first for GCMs prior to making predictions at 1 degree grid and second after making final predictions at ~4km.

13 HYDRO ENERGY↗

Aerosol Complexity and Implications for Predictability and Short-Term Forecasting

There are clear NWP and climate impacts from including aerosol radiative and cloud interactions. Changes in dynamics and cloud fields affect aerosol lifecycle, plume height, long-range transport, overall forcing of the climate system, etc. Inclusion of aerosols in NWP systems has benefit to surface field biases (e.g., T2m, U10m). Including aerosol affects has impact on analysis increments and can have statistically significant impacts on, e.g., tropical cyclogenesis. Above points are made especially with respect to aerosol radiative interactions, but aerosol-cloud interaction is a bigger signal on the global system. Many of these impacts are realized even in models with relatively simple (bulk) aerosol schemes (approx.10 -20 tracers). Simple schemes though imply simple representation of aerosol absorption and importantly for aerosol-cloud interaction particle-size distribution. Even so, more complex schemes exhibit a lot of diversity between different models, with issues such as size selection both for emitted particles and for modes. Prospects for complex sectional schemes to tune modal (and even bulk) schemes toward better selection of size representation. I think this is a ripe topic for more research -Systematic documentation of benefits of no vs. climatological vs. interactive (direct and then direct+indirect) aerosols. Document aerosol impact on analysis increments, inclusion in NWP data assimilation operator -Further refinement of baseline assumptions in model design (e.g., absorption, particle size distribution). Did not get into model resolution and interplay of other physical processes with aerosols (e.g., moist physics, obviously important), chemistry

Predictability↗

Radar Techniques for the Measurement of Mercury's Obliquity and Librations

Evidence for a global magnetic field at Mercury has raised perplexing questions regarding the existence and nature of the planet's core. The problems related to Mercury's core are of great interest because of their implications on how planets evolve thermally and on how they generate magnetic fields. Peale showed that the measurement of four quantities could place constraints on the size and state of Mercury's core. The required parameters in Peale's experiment are the C(sub 20) and C(sub 22) coefficients in the spherical harmonic expansion of the gravity field, the planet's obliquity theta, and the amplitude phi of the forced libration in longitude. While the gravitational harmonics are best determined from an orbiting spacecraft, Earth-based radar observations may provide refined estimates of the obliquity and amplitude of the libration. We describe several techniques and attempt to quantify their potential for the proposed measurements. Additional information is contained in the original extended abstract.

Margot, J. L.↗

The formation of massive primordial stars in the presence of moderate UV backgrounds

Radiative feedback produced by stellar populations played a vital role in early structure formation. In particular, photons below the Lyman limit can escape the star-forming regions and produce a background ultraviolet (UV) flux, which consequently may influence the pristine halos far away from the radiation sources. These photons can quench the formation of molecular hydrogen by photodetachment of H{sup –}. In this study, we explore the impact of such UV radiation on fragmentation in massive primordial halos of a few times 10{sup 7} M {sub ☉}. To accomplish this goal, we perform high resolution cosmological simulations for two distinct halos and vary the strength of the impinging background UV field in units of J {sub 21} assuming a blackbody radiation spectrum with a characteristic temperature of T {sub rad} = 10{sup 4} K. We further make use of sink particles to follow the evolution for 10,000 yr after reaching the maximum refinement level. No vigorous fragmentation is observed in UV-illuminated halos while the accretion rate changes according to the thermal properties. Our findings show that a few 10{sup 2}-10{sup 4} solar mass protostars are formed when halos are irradiated by J {sub 21} = 10-500 at z > 10 and suggest a strong relation between the strength of the UV flux and mass of a protostar. This mode of star formation is quite different from minihalos, as higher accretion rates of about 0.01-0.1 M {sub ☉} yr{sup –1} are observed by the end of our simulations. The resulting massive stars are potential cradles for the formation of intermediate-mass black holes at earlier cosmic times and contribute to the formation of a global X-ray background.

79 ASTRONOMY AND ASTROPHYSICS↗

Hydrogen and Ammonia Blending With Natural Gas: Ignition Delay Times and Chemical Kinetic Model Validation at Gas Turbine Relevant Conditions

Abstract Ignition delay times from undiluted mixtures of natural gas (NG)/H2/Air and NG/NH3/Air were measured using a high-pressure shock tube at the University of Central Florida. The combustion temperatures were experimentally tested between 1000–1500 K near a constant pressure of 25 bar. As mentioned, mixtures were kept undiluted to replicate the same chemistry pathways seen in gas turbine combustion chambers. Recorded combustion pressures exceeded 200 bar due to the large energy release, hence why these were performed at the high-pressure shock tube facility. The data is compared to the predictions of the NUIGMech 1.1 mechanism for chemical kinetic model validation and refinement. An exceptional agreement was shown for stoichiometric conditions in all cases but strayed at lean and rich equivalence ratios, especially in the lower temperature regime of H2 addition and all temperature ranges of the baseline NG mixture. Hydrogen addition also decreased ignition delay times by nearly 90%, while NH3 fuel addition made no noticeable difference in ignition time. NG/NH3 exhibited similar chemistry to pure NG under the same conditions, which is shown in a sensitivity analysis. The reaction CH3 + O2 = CH3O + O is identified and suggested as a possible modification target to improve model performance. Increasing the robustness of chemical kinetic models via experimental validation will directly aid in designing next-generation combustion chambers for use in gas turbines, which in turn will greatly lower global emissions and reduce greenhouse effects.

Pierro, Michael↗

On principles, methods and recent advances in studies towards a GPS-based control system for geodesy and geodynamics

Although Very Long Baseline Interferometry (VLBI) and Satellite Laser Ranging (SLR) are becoming increasingly important tools for geodynamic studies, their future role may well be fulfilled by using alternative techniques such as those utilizing the signals from the Global Positioning System (GPS). GPS, without the full implementation of the system, already offers a favorable combination of cost and accuracy and has consistently demonstrated the capability to provide high precision densification control in the regional and local areas of the VLBI and SLR networks. This report reviews VLBI and SLR vis-a-vis GPS and outlines the capabilities and limitations of each technique and how their complementary application can be of benefit to geodetic and geodynamic operations. It demonstrates, albeit with a limited data set, that dual-frequency GPS observations and interferometric type analysis techniques make possible the modelling of the GPS orbits for several days with an accuracy of a few meters. The use of VLBI or SLR sites as fiducial stations together with refinements in the orbit determination procedures can greatly reduce the systematic errors in the GPS satellite orbits used to compute the positions of non-fiducial locations. In general, repeatability and comparison with VLBI of the GPS determined locations are of the order of between 2 parts in 10 to the 7th power and 5 parts in 10 to the 8th power for baseline lengths less than 2000 km. This report is mainly a synthesis of problems, assumptions, methods and recent advances in the studies towards the establishment of a GPS-based system for geodesy and geodynamics and is one phase in the continuing effort for the development of such a system. To some, including the author, it seems reasonable to expect within the next few years that more evidence will show GPS to be as a powerful and reliable a tool as mobile VLBI and SLR are today, but largely more economical.

Delikaraoglou, Demitris↗

An Improved Model of the Crustal Structure of Mars

The first reliable model of the structure of the crust and upper mantle of Mars from remote observations was produced using data from the Mars Orbiter Laser Altimeter (MOLA) and the Radio Science investigation of Mars Global Surveyor (MGS). That model assumed a uniform crustal density and solved for the global variations in crustal thickness using a gravity field derived from preliminary MGS tracking. In that study, spherical harmonic potential coefficients were derived to degree and order 80, but crustal structure was interpreted cautiously to degree 60, or 360 km wavelength, owing to the presence of noise. Tracking normal equations have since been generated to degree 75, to degree 80 (supplemented by altimetric crossovers), and recently to degree 90, using new constants for the orientation of the spin pole and the rotation rate of Mars provided by the IAU2000 rotation model. Gravity models now incorporate tracking data coverage from the Primary and Extended MGS missions and the early phases of the Mars Odyssey mission. In the present study we exploit these advances in gravity modeling to present a refined crustal inversion, which we also interpret in the context of Mars' internal structure and thermal evolution.

Zuber, M. T.↗

How Certain are We of the Uncertainties in Recent Ozone Profile Trend Assessments of Merged Limbo Ccultation Records? Challenges and Possible Ways Forward

Most recent assessments of long-term changes in the vertical distribution of ozone (by e.g. WMO and SI2N) rely on data sets that integrate observations by multiple instruments. Several merged satellite ozone profile records have been developed over the past few years; each considers a particular set of instruments and adopts a particular merging strategy. Their intercomparison by Tummon et al. revealed that the current merging schemes are not sufficiently refined to correct for all major differences between the limb/occultation records. This shortcoming introduces uncertainties that need to be known to obtain a sound interpretation of the different satellite-based trend studies. In practice however, producing realistic uncertainty estimates is an intricate task which depends on a sufficiently detailed understanding of the characteristics of each contributing data record and on the subsequent interplay and propagation of these through the merging scheme. Our presentation discusses these challenges in the context of limb/occultation ozone profile records, but they are equally relevant for other instruments and atmospheric measurements. We start by showing how the NDACC and GAW-affiliated ground-based networks of ozonesonde and lidar instruments allowed us to characterize fourteen limb/occultation ozone profile records, together providing a global view over the last three decades. Our prime focus will be on techniques to estimate long-term drift since our results suggest this is the main driver of the major trend differences between the merged data sets. The single-instrument drift estimates are then used for a tentative estimate of the systematic uncertainty in the profile trends from merged data records. We conclude by reflecting on possible further steps needed to improve the merging algorithms and to obtain a better characterization of the uncertainties involved.

Merging Algorithms↗

In-Flight Characterization of the Lunar Orbiter Laser Altimeter Instrument Pointing and Far-Field Pattern

The Lunar Orbiter Laser Altimeter (LOLA) aboard the Lunar Reconnaissance Orbiter (LRO) has collected nearly seven billion measurements of surface height on the Moon with an absolute accuracy of ∼1 m and a precision of ∼10 cm. Converting time-of-flight laser altimeter measurements to topographic elevations requires accurate knowledge of the laser pointing with respect to the spacecraft body-fixed coordinate system. To that end, we have utilized altimetric crossovers from LOLA, as well as bidirectional observations of the LOLA laser and receiver boresight via an Earth-based laser tracking ground station. Based on a sample of ∼780,000 globally distributed crossovers from the circular-orbit phase of LRO's mission (∼27 months), we derive corrections to the LOLA laser boresight. These corrections improve the cross-track and along-track agreement of the crossovers by 24% and 33%, respectively, yielding RMS residuals of ∼10 m. Since early in the LRO mission, the bidirectional laser tracking experiments have confirmed a pointing anomaly when the LOLA instrument is facing toward deep space or the night side of the Moon and have allowed the reconstruction of the laser far-field pattern and receiver telescope pointing. By conducting such experiments shortly after launch and nearly eight years later, we have directly measured changes in the laser characteristics and obtained critical data to understand the laser behavior and refine the instrument pointing model. The methods and results presented here are also relevant to the design, fabrication, and operation of future planetary laser altimeters and their long-term behavior in the space environment.

Barker, M. K.↗

Metrics for Intercomparison of Remapping Algorithms (MIRA) protocol applied to Earth system models

Abstract. Strongly coupled nonlinear phenomena such as those described by Earth system models (ESMs) are composed of multiple component models with independent mesh topologies and scalable numerical solvers. A common operation in ESMs is to remap or interpolate component solution fields defined on their computational mesh to another mesh with a different combinatorial structure and decomposition, e.g., from the atmosphere to the ocean, during the temporal integration of the coupled system. Several remapping schemes are currently in use or available for ESMs. However, a unified approach to compare the properties of these different schemes has not been attempted previously. We present a rigorous methodology for the evaluation and intercomparison of remapping methods through an independently implemented suite of metrics that measure the ability of a method to adhere to constraints such as grid independence, monotonicity, global conservation, and local extrema or feature preservation. A comprehensive set of numerical evaluations is conducted based on a progression of scalar fields from idealized and smooth to more general climate data with strong discontinuities and strict bounds. We examine four remapping algorithms with distinct design approaches, namely ESMF Regrid (Hill et al., 2004), TempestRemap (Ullrich and Taylor, 2015), generalized moving least squares (GMLS) (Trask and Kuberry, 2020) with post-processing filters, and WLS-ENOR (Li et al., 2020). By repeated iterative application of the high-order remapping methods to the test fields, we verify the accuracy of each scheme in terms of their observed convergence order for smooth data and determine the bounded error propagation using challenging, realistic field data on both uniform and regionally refined mesh cases. In addition to retaining high-order accuracy under idealized conditions, the methods also demonstrate robust remapping performance when dealing with non-smooth data. There is a failure to maintain monotonicity in the traditional L2-minimization approaches used in ESMF and TempestRemap, in contrast to stable recovery through nonlinear filters used in both meshless GMLS and hybrid mesh-based WLS-ENOR schemes. Local feature preservation analysis indicates that high-order methods perform better than low-order dissipative schemes for all test cases. The behavior of these remappers remains consistent when applied on regionally refined meshes, indicating mesh-invariant implementations. The MIRA intercomparison protocol proposed in this paper and the detailed comparison of the four algorithms demonstrate that the new schemes, namely GMLS and WLS-ENOR, are competitive compared to standard conservative minimization methods requiring computation of mesh intersections. The work presented in this paper provides a foundation that can be extended to include complex field definitions, realistic mesh topologies, and spectral element discretizations, thereby allowing for a more complete analysis of production-ready remapping packages.

58 GEOSCIENCES↗