Surrogate Modeling of Spatially Heterogeneous Source Terms for Probabilistic Assessment of Repository Performance.
Abstract not provided.
SEARCH · Search NASA
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.
Abstract not provided.
This paper presents the steps followed to GPU-offload parts of the core solver of EFIT-AI, an equilibrium reconstruction code suitable for tokamak experiments and burning plasmas. For this work, we will focus on the fitting procedure that consists of a Grad–Shafranov (GS) equation inverse solver that calculates equilibrium reconstructions on a grid. We will show profiling results of the original code (CPU-baseline), as well as the directives used to GPU-offload the most time-consuming function, initially to compare OpenACC and OpenMP on NVIDIA and AMD GPUs and later on to assess OpenMP performance portability on NVIDIA, AMD and Intel GPUs. We will make a performance comparison for different spatial grid sizes and show the speedup achieved on NVIDIA A100 (Perlmutter-NERSC), AMD MI250X (Frontier-OLCF) and Intel PVC GPUs (Sunspot-ALCF). Finally, we will draw some conclusions and recommendations to achieve high-performance portability for an equilibrium reconstruction code on the new HPC architectures
The efficacy of vegetation dynamics simulations in offline land surface models (LSMs) largely depends on the quality and spatial resolution of meteorological forcing data. In this study, the Princeton Global Meteorological Forcing Data (PMFD) and the high spatial resolution and upscaled China Meteorological Forcing Data (CMFD) were used to drive the Simplified Simple Biosphere model version 4/Top-down Representation of Interactive Foliage and Flora Including Dynamics (SSiB4/TRIFFID) and investigate how meteorological forcing datasets with different spatial resolutions affect simulations over the Tibetan Plateau (TP), a region with complex topography and sparse observations. By comparing the monthly Leaf Area Index (LAI) and Gross Primary Production (GPP) against observations, we found that SSiB4/TRIFFID driven by upscaled CMFD improved the performance in simulating the spatial distributions of LAI and GPP over the TP, reducing RMSEs by 24.3% and 20.5%, respectively. The multi-year averaged GPP decreased from 364.68 gC m –2 yr –1 to 241.21 gC m –2 yr –1 with the percentage bias dropping from 50.2% to –1.7%. When using the high spatial resolution CMFD, the RMSEs of the spatial distributions of LAI and GPP simulations were further reduced by 7.5% and 9.5%, respectively. This study highlights the importance of more realistic and high-resolution forcing data in simulating vegetation growth and carbon exchange between the atmosphere and biosphere over the TP.
Capacity expansion models for the power sector are used to project future decisions over the coming decades by simulating investment and operation decisions for the use of electricity. Due to model performance constraints, these models typically do not explicitly simulate every hour within a year, but instead simulate representative time segments (groups of hours). This paper evaluates different approaches for selecting time segments across three methods: sequential, categorical, and clustering, across a wide range of time-segment quantities, for a total of 204 temporal profiles. To measure the performance of each profile's ability to accurately represent data, the root-mean-square-error of each profile's time segments are compared to the data's original hourly data. The temporal alignment across regions is also measured (i.e., how often windy days align across regions). Different spatial resolutions were applied for a subset of the temporal selection methods to investigate the impact spatial resolution has on performance. This paper provides a framework for measuring the value of different temporal selection methods and of adding more granular data to energy system models. Overall, multi-criteria clustering yields the lowest root-mean-square-error across all datasets evaluated and provides a holistic view of the intertwined relationships between renewable generation and electricity demand.
A computational fluid dynamics study was performed for a cross-flow marine hydrokinetic turbine. The analysis was done in three dimensions and used the unsteady Reynolds-averaged Navier-Stokes solver in the commercial code STAR-CCM+. The base turbine configuration is the RivGen(R) Turbine, designed by the Ocean Renewable Power Company. A convergence and uncertainty analysis was performed for both the spatial and temporal discretization; this was done using the base configuration, which features support struts and helical foils. Both struts and helical blades introduce three-dimensional flow effects, influencing the complex flow phenomenon of dynamic stall. The study compares the relative impact of struts on power performance and blade loading for both helical and straight blades, and found that for this turbine the relative loss in power from struts was lower with helical blades.
Thermal neutron detection in neutron scattering science is a challenging endeavour due to a limited number of technologies which are sensitive to these weakly interacting particles. While many improvements to conventional detectors are underway at various facilities, there is a great opportunity to make a leap in performance by combining the spatial resolution benefits of imaging and the temporal resolution and background discrimination of event-driven radiation detectors. This idea has been realized by using a neutron sensitive scintillator read out by a single-photon sensitive camera based on a Timepix3 ASIC. We demonstrate how such data-driven imaging sensors can enable unprecedented performance in neutron reflectometry using the ASTERIX instrument at the Los Alamos Neutron Scattering Center. Several samples were measured with both the new and a conventional 3 He detector systems. The results from this work demonstrate that these imaging based systems can satisfy performance parameters for the future QIKR reflectometer to be built at the Second Target Station at Oak Ridge National Laboratory. Further improvements to the detector are already underway which will allow streamlined and expedited experiments. We demonstrate at least a two order of magnitude increase in detection rate at an acceptable dead time and introduce a new way of tuning the detector efficiency using light collecting optics to accommodate highly intense direct beams which cannot be measured with any current detectors without severe attenuation. This will allow measurements of complete reflectometry profiles and using a single sample measurement combined with fewer direct beam calibration measurements on QIKR and potentially other reflectometers.
Laser-produced plasma coupled with optical emission spectroscopy (OES) is a promising technique for detecting certain isotopes, with unique capabilities such as standoff and rapid detection and minimal to no sample preparation requirements. The key figure-of-merit for isotopic analysis using optical spectroscopy tools is the linewidth relative to the isotope shift. Although the isotopes of hydrogen (1H, 2H, and 3H) possess large isotopic shifts (1H–2H ≈ 180 pm, 1H–3H ≈ 240 pm), being a light element, the H transitions are susceptible to various broadening mechanisms in the plasma environment. One of the critical parameters that influence the linewidth of a transition in an LPP is the incident laser energy. In the present study, we evaluated the role of laser energy on plume expansion dynamics, deuterium emission intensity, and linewidth in a nanosecond laser-produced Zircaloy-4 plasma. The changes in 2Hα emission intensity and linewidth were investigated for varying laser fluence and time after plasma onset. Spatially resolved and spatially integrated OES were performed and compared to investigate the emission spectral features and linewidth of 2Hα. Monochromatic two-dimensional time-resolved imaging was also performed to understand the morphology of the deuterium and protium emission relative to all species in the plume. Our results showed that 1Hα and 2Hα emissions predominantly occur closer to the target. Measurements of 2Hα linewidth approached similar values at later times of plasma evolution regardless of the laser energy. The linewidths of the 2Hα transition showed insignificant differences between spatially resolved and spatially integrated measurements.
Smart transportation infrastructure relies on networks of heterogeneous sensors - cameras, radars, and lidars - continuously monitoring traffic conditions. However, executing the initial spatial calibration of multiple sensors and the subsequent health monitoring presents significant operational challenges. Environmental factors, mechanical vibrations, and gradual drift cause spatial misalignment, degrading fusion performance and tracking accuracy. Traditional calibration approaches require manual intervention with specialized targets or survey equipment, resulting in service interruptions and high maintenance costs. This work presents an automated framework for initial calibration and continuous health monitoring without human intervention or service disruption. Our approach addresses two critical problems: (1) detecting when sensors become miscalibrated during operation, and (2) automatically re-establishing spatial alignment using only operational traffic data. The health monitoring component analyzes measurement innovations - differences between sensor observations and predicted object states - to detect systematic biases indicative of calibration drift. By computing bias magnitude, directional consistency, and rejection rates, the system identifies miscalibrations as small as 0.5 meters. Unlike traditional methods requiring known calibration targets, our diagnostic operates continuously on live traffic observations, enabling early detection before fusion quality degrades. The automatic recalibration algorithm leverages overlapping sensor fields-of-view and temporal correlation of vehicle observations. Using graph-based optimization, the system automatically discovers which sensor pairs observe common regions, estimates pairwise spatial transformations using RANSAC-based robust estimation, and jointly optimizes all sensor poses through bundle adjustment. The framework handles practical deployment challenges, including different sensor sampling rates (1-10 Hz), varying installation positions, unknown orientations, and limited overlap regions (>10%). When approximate sensor positions are available from installation surveys (+/-1m accuracy), the algorithm additionally estimates sensor orientations, refining both position and rotation to sub-meter and sub-degree accuracy. We validate the framework on multi-hour traffic datasets from six heterogeneous sensors with sampling rates ranging from 1 Hz to 10 Hz. Results demonstrate successful calibration even with sparse overlap (<20%) and automatic detection of miscalibrations exceeding 0.8 meters. This work enables a "deploy-and-forget" sensor infrastructure that maintains calibration autonomously, reducing maintenance costs while improving tracking accuracy. The techniques generalize beyond transportation to any multi-sensor monitoring application requiring robust spatial alignment, including smart cities, industrial monitoring, and surveillance systems.
As an optical processor, a diffractive deep neural network (D2NN) utilizes engineered diffractive surfaces designed through machine learning to perform all-optical information processing, completing its tasks at the speed of light propagation through thin optical layers. With sufficient degrees of freedom, D2NNs can perform arbitrary complex-valued linear transformations using spatially coherent light. Similarly, D2NNs can also perform arbitrary linear intensity transformations with spatially incoherent illumination; however, under spatially incoherent light, these transformations are nonnegative, acting on diffraction-limited optical intensity patterns at the input field of view. Here, we expand the use of spatially incoherent D2NNs to complex-valued information processing for executing arbitrary complex-valued linear transformations using spatially incoherent light. Through simulations, we show that as the number of optimized diffractive features increases beyond a threshold dictated by the multiplication of the input and output space-bandwidth products, a spatially incoherent diffractive visual processor can approximate any complex-valued linear transformation and be used for all-optical image encryption using incoherent illumination. The findings are important for the all-optical processing of information under natural light using various forms of diffractive surface-based optical processors.
The US Department of Energy (DOE) conducts climate simulations on some of the world’s largest supercomputers. These exascale machines use heterogeneous architectures with both CPUs and GPUs, and scientific codes must adapt to make full use of this computing power. Los Alamos National Lab is developing Omega: The Ocean Model for E3SM Global Applications, which is specifically designed for modern exascale computers. It uses external libraries that have been optimized for a variety of architectures to run on different supercomputers. Omega is an unstructured-mesh ocean model based on TRiSK numerical methods. It will be the new ocean component of the DOE’s Energy Exascale Earth System Model (E3SM). The algorithms in Omega follow those of the current ocean component, MPAS-Ocean, but it will be written in C++ rather than Fortran to take advantage of the Kokkos performance portability library. Omega spatial operators are written as Kokkos kernels to run efficiently on both CPUs and GPUs. Work on Omega began in 2023 with a new C++ framework for unstructured mesh partitioning, halo exchanges, parallel IO, and Kokkos interfaces. The current version, Omega-0, is being developed to solve the shallow water equations and at present includes all of the tendency terms but not time stepping. Here we share the results of Omega-0 verification and performance testing. Verification includes unit tests implemented with CTest as well as convergence tests in Polaris, an in-house python package with a large suite of test problems. Performance tests compare simulations conducted on CPUs versus GPUs and across different architectures: tests are run on Frontier, which has AMD “Optimized 3rd Gen EPYC” CPUs and AMD MI250X GPUs, as well as Perlmutter, which is composed of AMD EPYC 7763 CPUs and NVIDIA A100 GPUs.
Systems involving advection-dominated transport through heterogeneous porous and fractured media are ubiquitous in subsurface engineering applications. However, upscaling such systems continues to challenge rigorous modeling efforts, particularly when advection is stronger than diffusion at fine spatial scales (i.e., when the Péclet number is greater than one at length scales that characterize a system’s unit-cells, representative elementary volumes, or averaging regions). Here, in this work, we propose and validate a strategy for extending the Method of Finite Averages (MoFA), a rigorous upscaling methodology for heterogeneous porous media, to upscale transport systems experiencing stronger advection than diffusion at fine scales (i.e., fine-scale Péclet numbers greater than one). We detail the strategy, the physical conditions under which it can be applied while retaining a priori modeling error guarantees, and implement the strategy to obtain a MoFA model for advective-diffusive transport that accommodates advective physics at fine spatial scales. We then perform two numerical experiments considering systems with system-scale Péclet numbers of 300 and 1000 — which correspond to fine-scale Péclet numbers of 30 and 100, respectively — to verify that the error guarantees are met under the strategy. After, we conduct a numerical study to demonstrate the strategy’s advantages over the original MoFA methodology. The results suggest that rigorously-upscaled transport models for heterogeneous porous media experiencing advective physics at finer spatial scales can be derived through MoFA and resolved orders of magnitude faster than their pore-scale counterparts. The results also suggest that the presented strategy is limited to modeling shallow concentration gradients when there are large differences between the time scales related to advection and a system’s temporally-varying boundary conditions. This limitation hinders the strategy’s practicality in modeling more advective systems, and as such, opportunity exists for developing additional strategies that accommodate rapidly-varying boundary conditions — and consequentially, steeper concentration gradients — while modeling advective systems with MoFA.
Two-dimensional (2D) semiconductors with narrow bandgaps are promising candidates for near- and far-infrared (IR) photodetection, particularly in the telecommunication spectral window. However, current low-bandgap IR photodetectors face significant challenges due to their high dark current, increased carrier recombination, and thermally generated noise. Here, in this work, a hybrid phototransistor is demonstrated by integrating direct, contact-free palladium diselenide (PdSe 2 ) as a highly responsive IR detection layer with a non-IR-absorbing molybdenum diselenide (MoSe 2 ) field-effect transistor (FET), using a near-IR source at a wavelength of λ = 1650 nm. Exfoliated PdSe 2 flakes integrated into a back-gated FET architecture exhibit ambipolar transport behavior, with extracted hole and electron mobilities of 24.8 cm 2 V –1 s –1 and 58.4 cm 2 V –1 s –1 , respectively. The devices show a clear photocurrent generation under the illumination of a λ = 1650 nm laser source, achieving a notable responsivity of ∼300 mA W -1 at an applied gate voltage of 15 V, which highlights the suitability of PdSe 2 as a narrow-bandgap material for photodetection. Photoresponsivity saturates and does not have any effect above an applied gate voltage of 15 V. To further tune the photoresponsivity performance continuously with the applied gate voltage, we construct a van der Waals heterostructure phototransistor, where few layers of PdSe 2 are directly transferred onto the 2D channel region of a MoSe 2 FET, while avoiding any contact with the metal electrodes. In this heterostructure, PdSe 2 works as the primary active IR-absorbing layer, while MoSe 2 provides high-performance FET characteristics. This spatial separation of absorption and transport facilitates efficient interlayer charge transfer and charge separation, resulting in high responsivities of up to 972 mA W –1 at near-IR wavelengths and a low power density of 1.5 mW/mm 2 . The responsivity of our photodetector is comparable to that of some state-of-the-art commercially available NIR photodetectors, highlighting the potential of PdSe 2 -based heterostructures as scalable, CMOS-compatible platforms for high-performance near-IR detection.
An accurate characterization of the wind resource available at hub-height is required for an efficient and bankable wind farm project. However, direct measurement of wind speed at the constantly increasing height of the hub of commercial wind turbines is oftentimes challenging and expensive, so that it is common practice to vertically extrapolate the wind resource from lower and more easily accessible levels. Conventional techniques for wind speed vertical extrapolation include the use of a power law and a logarithmic profile. While simple, the limits in accuracy of these methods have been shown in various studies. Recently, machine learning has been proposed as a new method to vertically extrapolate winds. All the published studies on the topic assess the performance of machine learning techniques in vertically extrapolating the wind resource at the same location where the algorithm has been trained. However, in real-world applications, the wind resource is measured at the instrument location, but it then needs to be extrapolated at hub height at the location of the wind turbines within the find farm. To be able to fully recommend the use of machine learning techniques over the simple power law and logarithmic law, the spatial variability of the performance improvements of the machine learning approaches needs to be assessed. Here, we propose a round-robin validation of a machine learning-based method for wind speed extrapolation. We use 20 months of observations at four locations spanning a 100 km wide region at the Southern Great Plains (SGP) atmospheric observatory, in north-central Oklahoma. At each location, we train a random forest to predict 30-min average wind speed at 143 m AGL. We use as input features lidar wind speed at 65 m AGL, time of day, sonic anemometer wind speed at 4 m AGL, turbulent kinetic energy, and Obukhov length. First, we perform a same-site comparison of the performance of the proposed random forest against the conventional techniques for wind speed extrapolation (namely power law and logarithmic profile, with widely accepted stability corrections). We find that the random forest outperforms the power law in vertically extrapolating wind speed in all the considered stability regimes, with a 33% reduction in MAE for stable conditions, and a 31% reduction in unstable conditions. Similar results are found when comparing predictions of extrapolated winds from the logarithmic profile and the random forest with the observed values. Next, we propose a round-robin validation, to use the random forest trained at each site to extrapolate wind speed at the remaining three sites. We find that the performance of the random forest approach degrades when the algorithm is tested at a site different than the training one. However, even under those circumstances, the machine learning-based approach still outperforms the conventional techniques for wind speed extrapolation, with, on average, a reduction in mean absolute error between 15 and 20% over the conventional methods, with the largest benefits obtained under stable conditions.
The nuclear industry is pursuing microreactors that can be factory assembled and deployed to remote regions for reliable power generation. One class of microreactors uses a monolithic metal core block coupled to heat pipes for heat rejection, which results in significant thermal stresses in the monolithic structures. This work describes the initial characterization and test plan for evaluating stainless steel test articles fabricated with embedded sensors for measuring heat pipe performance limits, as well as spatially distributed temperatures and strains during electrically heated thermal testing. The electrically heated testing will be performed in the non-nuclear Microreactor Agile Non-Nuclear Testbed and Single Primary Heat Extraction and Removal Emulator facilities located at Idaho National Laboratory. The goals of these tests are to (1) accurately monitor temperature and strain distributions that result from differential thermal expansion in the test articles and (2) quantify heat rejection limits of heat pipes as a function of operating temperature and working fluid during steady-state and transient operations. More generally, the ability to monitor component and system health during microreactor operation is attractive for providing a high sensor density to inform a limited number of microreactor operators to ultimately reduce operation and maintenance costs and move toward semi-autonomous operation. This report discusses the characterization of embedded thermocouples and fiber optic sensors in relevant test articles, including cylindrical pipes and hexagonal monolithic test articles for heat pipe-based reactors. The sensors were embedded by placing them in machined channels and then building up additional material by using ultrasonic additive manufacturing (UAM). UAM is a solid-state welding process that uses downward pressure and a lateral scrubbing motion to bond thin metal foils to a base material layer by layer. The ultrasonic welding process relies on the plastic deformation of the metal—as opposed to typical melting and solidification—to break oxide scales and bond the metal layers. The characterization of these embedded sensors included evaluating fiber optic signal attenuation, observing residual strain in the fibers, investigating microstructural and mechanical aspects, and demonstrating the sensors under various thermal loads and acoustic vibrations. Post-embedding characterization showed a fine grain structure (<1 μm) near the interfaces of the bonded foils as a result of severe deformation from the welding process. A large increase in hardness was observed at the foil interfaces and the fiber/matrix interface compared with the bulk matrix. Even when compared with the SS304 interfaces, the higher hardness observed around the embedded fiber suggests a higher degree of deformation due to the soft metal coating around the silica fiber core. The distributed fiber-optic temperature sensors and embedded thermocouples reliably measured temperature distributions during steady-state and transient thermal testing. The embedded fiber-optic sensors reliably measured strain during both transient and steady-state testing and properly identified resonant frequencies during acoustic testing.
Given its versatility in drawing power from many sources in the natural world, piezoelectric energy harvesting (PEH) has become increasingly popular. However, its energy harvesting capacities could be enhanced further. Here, a mathematical model that accurately simulates the dynamic behavior and energy harvested can facilitate further improvements in the performance of piezoelectric devices. One of the goals of this study is to create a dependable reduced-order model of a multi-purpose gyroscope. This model will make it possible to compute the harvested voltage and electrical power in a semi-analytical manner. The harvested voltage is often modeled as an average value across the whole electrode surface in piezoelectric devices. We propose a model which provides practical insights toward optimizing the performance of the system by considering a spatially varying electric field across the electrode surface length. Our framework allows investigation of the limits of applicability of the modeling assumptions across a range of load resistances. The differential quadrature method (DQM) provides the basis for the suggested numerical solution. The model is also employed to examine energy harvesting under various resistance loads. The newly developed spatially varying model is evaluated for open- and closed-circuit conditions and is proved to be accurate for various values of load resistance that have not previously been considered. The results show that using a spatially varying model is more versatile when modeling the performance of the piezoelectric multifunctional energy harvester. The performance may be accurately captured by the model for load resistances ranging between 103 Ω and 108 Ω. At optimum load resistance and near 65 KHz, the maximum power output predicted by the spatially varying (SV) model is 1.3 mV, 1.5 mV for the open-circuit (OC) model, and 2.1 mV for the closed circuit (CE) model. At a high-load resistance, the SV and OC models all predict the maximum power output to be 1.9 mV while the CE model predicted the maximum voltage to be 3 mV.
Abstract. Soil moisture (SM) datasets are critical to understanding the global water, energy, and biogeochemical cycles and benefit extensive societal applications. However, individual sources of SM data (e.g., in situ and satellite observations, reanalysis, offline land surface model simulations, Earth system model – ESM – simulations) have source-specific limitations and biases related to the spatiotemporal continuity, resolutions, and modeling and retrieval assumptions. Here, we developed seven global, gap-free, long-term (1970–2016), multilayer (0–10, 10–30, 30–50, and 50–100 cm) SM products at monthly 0.5∘ resolution (available at https://doi.org/10.6084/m9.figshare.13661312.v1; Wang and Mao, 2021) by synthesizing a wide range of SM datasets using three statistical methods (unweighted averaging, optimal linear combination, and emergent constraint). The merged products outperformed their source datasets when evaluated with in situ observations (mean bias from −0.044 to 0.033 m3 m−3, root mean square errors from 0.076 to 0.104 m3 m−3, Pearson correlations from 0.35 to 0.67) and multiple gridded datasets that did not enter merging because of insufficient spatial, temporal, or soil layer coverage. Three of the new SM products, which were produced by applying any of the three merging methods to the source datasets excluding the ESMs, had lower bias and root mean square errors and higher correlations than the ESM-dependent merged products. The ESM-independent products also showed a better ability to capture historical large-scale drought events than the ESM-dependent products. The merged products generally showed reasonable temporal homogeneity and physically plausible global sensitivities to observed meteorological factors, except that the ESM-dependent products underestimated the low-frequency temporal variability in SM and overestimated the high-frequency variability for the 50–100 cm depth. Based on these evaluation results, the three ESM-independent products were finally recommended for future applications because of their better performances than the ESM-dependent ones. Despite uncertainties in the raw SM datasets and fusion methods, these hybrid products create added value over existing SM datasets because of the performance improvement and harmonized spatial, temporal, and vertical coverages, and they provide a new foundation for scientific investigation and resource management.
Offshore aquaculture and marine renewable energy (energy from waves, tides, currents, and ocean gradients) are two developing ocean-based industries. Aquaculture, an industry that has typically relied on diesel for power, is expected to grow globally, presenting an opportunity to reduce greenhouse gas emissions by switching to renewable sources as it expands. As the industry moves further offshore and is situated in more energetic environments, the prospect to co-locate offshore aquaculture with wave energy increases. To improve understanding of this potential, a feasibility assessment was completed to estimate the energy needs and wave resource required to power offshore finfish aquaculture operations. The study found it is possible to power offshore aquaculture operations entirely with wave energy. A spatial analysis was then performed to assess the suitability of co-locating offshore finfish aquaculture and wave energy off California and Hawaii. Suitable locations were identified offshore of O’ahu, Hawaii, and northern California. Southern California was also assessed, using a lower wave resource, based on study areas evaluated by the National Oceanic and Atmospheric Administration to identify Aquaculture Opportunity Areas, and while limited there are suitable locations that may warrant further evaluation. This study presents an analysis into the potential to pair wave energy with offshore aquaculture, and how various factors can help determine suitable areas for co-location. The analyses developed in this study will support future identification of potential sites for development and decision-making to optimize the success of co-locating wave energy resources and offshore finfish aquaculture.
Microbially-explicit models may improve understanding and projections of carbon dynamics in response to future climate change, but their fidelity in simulating global-scale soil heterotrophic respiration (RH), a stringent test for soil biogeochemical models, has never been evaluated. We used statistical global RH products, as well as 7,821 daily site-scale RH measurements, to evaluate the spatio-temporal performance of one first-order decay model (CASA-CNP) and two microbially-explicit biogeochemical models (CORPSE and MIMICS) that were forced by two different climate datasets. CORPSE and MIMICS did not provide any measurable performance improvement; instead, the models were highly sensitive to the meteorological input data used to drive them. Spatial RH variability was generally well simulated except in the northern middle latitudes (~50°N) and arid regions; models captured the seasonal variability of RH well, but showed more divergence in tropic and arctic regions. Our results demonstrate that the next generation of biogeochemical models shows promise, but also needs to be improved for realistic spatio-temporal variability of RH. Finally, we emphasize the importance of net primary production, soil moisture, and soil temperature inputs, and that jointly evaluating soil models for their spatial (global scale) and temporal (site scale) performance provides crucial benchmarks for improving biogeochemical models.