Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Statistical calibration”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Measurement Uncertainty in One-Of-A-Kind Event Data Analysis

A golden standard in science is to repeat an experiment a statistically significant number of times, recording data using the same set of detectors and the same data analysis methodology. In such case experimental error includes both the range of true values generated by repetitions of the experiment, and measurement uncertainty caused by the detector. They are independent. It is a huge and too frequently used simplification, to assume that one can measure multiple repetitions of an identical experiment, resulting in identical true experimental value. Repetitions, as similar is it is experimentally achievable, have unavoidable built-in differences resulting in a range of the true values rather than in a single value. When modern, very sensitive and well calibrated measurement systems are used, this range is not negligible, and sometimes dominates over the measurement uncertainty. Range of true values depends on built-in differences in physics of the experiment. Stochastic physical processes result typically in a broader range of true values than non-stochastic processes do. Measurement uncertainty depends on a measurement method (properties of the detector not of the experiment). Modern measurement methods, including digital ones, frequently make the measurement uncertainty very small. When data from one–of –a kind experiment are analyzed, only the measurement uncertainty is reported. It provides no information about the range of true experimental values, neither about reliability of a reported data point. Reliability of a data point is in general independent from its measurement uncertainty. However, in practice reliable measurement methods frequently have high measurement uncertainty, while low reliability methods are applied to limit measurement uncertainty. Comparison of reliable data with high measurement uncertainty to not so reliable data measured with low uncertainty is discussed – in different scenarios different data analysis methods are applicable. Methods for data analysis from an experiment repeated statistically significant number of times are very well developed. They do not require a detailed expertise in physics of an experiment, nor in the properties of the measurement system used, and meaning of the reported uncertainty is well understood in any scientific community. It all changes when data from one-of-a-kind experiment is analyzed. Analyst’s expertise is required both in the physics of the experiment and in all aspects of the measurement system, all possible malfunctions. Data users must remember that only measurement uncertainty is reported from any one-of-a-kind experiment. Theory with simulations may provide estimation of expected built-in differences in the experiment, and by this of expected range of true values for a given experiment; yet measurement uncertainty can never be used in place of the range of true experimental values.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Measurement of $\nu_\mu$ CC Interactions With Two-Proton Final State in MINERvA

This dissertation presents a measurement of charged–current (CC) muon–neutrino interactions with exactly two protons and no pions in the final state (CC~$2p\,0\pi$), using data collected by the MINERvA detector in the NuMI medium–energy beam at Fermilab. Such two–proton topologies are a sensitive probe of nuclear dynamics in the few–GeV regime, including multi–nucleon correlations (npnh, notably $2p2h$) and intranuclear final–state interactions (FSI) such as pion absorption and nucleon rescattering. A precise experimental characterization of these processes is essential both for neutrino–interaction theory and for reducing systematic uncertainties in oscillation experiments that rely on accurate modeling of neutrino–nucleus interactions. Events are selected by requiring a $\nu_\mu$ CC interaction with a reconstructed $\mu^-$ and two proton tracks originating from a common vertex in MINERvA’s finely segmented scintillator tracker, with no reconstructed mesons. Muon charge and momentum are constrained by matching to the MINOS Near Detector, while proton identification exploits energy–loss profiles and stopping–proton features. Backgrounds from pion–producing channels that enter the signal region through FSI or reconstruction effects are constrained with data–driven sidebands (Michel–electron and isolated–cluster “blob” samples) and tuned via a simultaneous fit across signal and sideband regions. To correct detector resolution and acceptance effects, the analysis employs iterative Bayesian unfolding with extensive validation: statistical pseudo–experiments, and robustness checks against generator systematic “universes” and additional strong shape warps. Single–differential cross sections are reported for three observables tailored to the two–proton final state: the opening–angle cosine $\cos\!\left(\theta_{pp}\right)$, the leading–proton momentum, and the subleading–proton momentum. Systematic uncertainties include contributions from neutrino flux, interaction modeling (e.g., npnh and resonance parameters, pion FSI), and detector response (calibration, reconstruction efficiencies). The resulting distributions provide targeted constraints on the interplay of multi–nucleon dynamics and FSI that shape CC~$2p\,0\pi$ final states on hydrocarbon. Comparisons to modern GENIE–based simulations highlight kinematic regions where model components require refinement. These measurements thus inform generator tuning and improve the reliability of neutrino–energy reconstruction strategies for current and future long–baseline oscillation programs.

Syrotenko, Vladyslav S. [Tufts U.]↗
DOE OSTI2020

Measurement Error and Resolution in Quantitative Stable Isotope Probing: Implications for Experimental Design

Quantitative stable isotope probing (qSIP) estimates isotope tracer incorporation into DNA of individual microbes and can link microbial biodiversity and biogeochemistry in complex communities. As with any quantitative estimation technique, qSIP involves measurement error, and a fuller understanding of error, precision, and statistical power benefits qSIP experimental design and data interpretation. We used several qSIP data sets—from soil and seawater microbiomes—to evaluate how variance in isotope incorporation estimates depends on organism abundance and resolution of the density fractionation scheme. We assessed statistical power for replicated qSIP studies, plus sensitivity and specificity for unreplicated designs. As a taxon’s abundance increases, the variance of its weighted mean density declines. Nine fractions appear to be a reasonable trade-off between cost and precision for most qSIP applications. Increasing the number of density fractions beyond that reduces variance, although the magnitude of this benefit declines with additional fractions. Our analysis suggests that, if a taxon has an isotope enrichment of 10 atom% excess, there is a 60% chance that this will be detected as significantly different from zero (with alpha 0.1). With five replicates, isotope enrichment of 5 atom% could be detected with power (0.6) and alpha (0.1). Finally, we illustrate the importance of internal standards, which can help to calibrate per sample conversions of %GC to mean weighted density. These results should benefit researchers designing future SIP experiments and provide a useful reference for metagenomic SIP applications where both financial and computational limitations constrain experimental scope.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Statistical Properties of Rainfall in Tasmania (SPORT) Field Campaign Report

The objective of this project is to complement the U.S. Department of Energy Atmospheric Radiation Measurement (ARM) User Facility second ARM Mobile Facility (AMF2) deployment at kennaook-Cape Grim (CAPE-k) with a micro-rain radar (MRR-PRO) mounted on top of one of the ARM containers. This was done to locate the MRR-PRO close to the ARM cloud radars (Ka-band ARM Zenith Radar [KAZR] and Marine W-band ARM Cloud Radar [MWACR]). This will allow for a continuous calibration check of the ARM cloud radars against the MRR-PRO and Parsivel-2 disdrometers. Scientifically, it will add one radar frequency to the analysis of vertical profiles of cloud and precipitation properties (K-band).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Understanding Biases in Sample Preparation Techniques for Coupled Scanning Electron Microscopy and MAMA PuO 2 Morphological Analysis

In this project, the scanning electron microscopy (SEM) sampling method used during the statistical design study (SDS) was investigated to determine if any sampling biases were present in the analyzed data. Using standard particle size distribution powders from the National Institute of Standards and Technology (NIST 1984 standard reference material) with the origin wet dispersion method, it was determined that a bias to smaller particles was present. This was supported by theoretical calculations using Stokes’ law to determine the settling rate of spherical particles of roughly the same size and mass as those found in the SDS. Based on the theoretical calculations, it was determined that the settling rate for each of the 76 powder sets in the SDS could be unique based on specific particle shape and mass distributions, making a universal correction factor/formula not applicable. Therefore, priority shifted to developing an improved wet dispersion method that significantly reduced the particle settling rate for all particle size and shapes. This was achieved by replacing the original solvent (isopropyl alcohol) with a heavy liquid (lithium heteropolytungstates), which dramatically slowed the settling rate and allowed for the capture of a suitable homogeneous aliquot. SEM imaging and Morphological Analysis for Material Attribution (MAMA) software analysis were conducted on the NIST standard, and the SEM/MAMA data were compared to data captured by a dynamic image analysis particle size analyzer. The resulting data confirmed that the new wet dispersion method does indeed deliver an improved representative aliquot to the SEM stub. For instance, in the NIST certificate, the average particle size is ~17.1 µm ± 2.2 µm with a normal distribution. The initial wet dispersion method resulted in a drastically reduced average particle size of 6.1 µm in addition to a non-representative heavy bi-modal distribution whereas the improved LST wet dispersion method resulting in an average particle size that was much closer to the NIST certificate (12.7 µm) with a similar normal distribution. Although the improved method was still short of the NIST certificate average, atomic force microscopy analysis determined that the resulting ~20-25% reduction in size was due to particles sinking into the carbon sticky tape used for SEM imaging. It is believed that that this bias can be calibrated in a much more predicable manner than the original settling rate bias. In addition, the matching normal distribution curves between the NIST certificate and the heavy liquid method indicate a much-improved representative aliquot has been sampled and imaged. A surrogate CeO 2 powder was used to reflect PuO 2 more accurately and to aid in implementing radiological controls and shielding. The resulting data sets from the SEM/MAMA method and the particle size analyzer give almost identical average particle sizes and particle distribution statistics. Future work will re-analyze several select runs from the SDS to determine if morphological signatures can be found with the improved sampling method.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Measurement of CP asymmetries in ${D}_{(s)}^{+}\to \eta {\pi}^{+}$ and ${D}_{(s)}^{+}\to {\eta}^{\prime }{\pi}^{+}$ decays

Searches for CP violation in the decays ${D}_{(s)}^{+}\to \eta {\pi}^{+}$ and ${D}_{(s)}^{+}\to {\eta}^{\prime }{\pi}^{+}$ are performed using pp collision data corresponding to 6 fb –1 of integrated luminosity collected by the LHCb experiment. The calibration channels ${D}_{(s)}^{+}\to \phi {\pi}^{+}$ are used to remove production and detection asymmetries. The resulting CP-violating asymmetries are $\mathcal{A}^{CP} = (D^+ \to ηπ^+) = (0.34 ± 0.66 ± 0.16 ± 0.05)$ %, $\mathcal{A}^{CP} = (D^+_s \to ηπ^+) = (0.32 ± 0.51 ± 0.12)$ %, $\mathcal{A}^{CP} = (D^+ \to η'π^+) = (0.49 ± 0.18 ± 0.06 ± 0.05)$ %, $\mathcal{A}^{CP} = (D^+_s \to η'π^+) = (0.01 ± 0.12 ± 0.08)$ %, where the first uncertainty is statistical, the second is systematic and the third, relevant for the D + channels, is due to the uncertainty on ${\mathcal{A}}^{CP}=\left({D}^{+}\to \phi {\pi}^{+}\right)$. These measurements, currently the most precise for three of the four channels considered, are consistent with the absence of CP violation. A combination of these results with previous LHCb measurements is presented.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Code Release for “Unlocking Extreme Space Weather through Advanced Modeling of Legacy Vela Spacecraft” ER

The software being developed for this project has two mains aims. First, a Bayesian Model Calibration (BMC) procedure is being developed to calibrate a spallation model that simulates protons hitting a spacecraft orbiting earth to real data. The procedure will be developed for general data (there is no data release requested as part of this code release). Second, an inverse physics modeling task is being undertaken to map the number of resulting neutrons observed from this process to the expected number of protons that hit the model. This second task is of statistical interest; to publish on it, the code will need to be open source.

Murph, Alexander↗
DOE OSTI2024

Evaluating the Consistency of Cosmological Distances Using Supernova Siblings in the Near-infrared

Abstract The study of supernova (SN) siblings, supernovae with the same host galaxy, is an important avenue for understanding and measuring the properties of Type Ia SN Ia light curves (LCs). Thus far, sibling analyses have mainly focused on optical LC data. Considering that LCs in the near-infrared (NIR) are expected to be better standard candles than those in the optical, we carry out the first analysis compiling SN siblings with only NIR data. We perform an extensive literature search of all SN siblings and find six sets of siblings with published NIR photometry. We calibrate each set of siblings ensuring they are on homogeneous photometric systems, fit the LCs with the SALT3-NIR and SNooPy models, and find median absolute differences inμvalues between siblings of 0.248 and 0.186 mag, respectively. To evaluate the significance of these differences beyond measurement noise, we run simulations that mimic these LCs and provide an estimate for uncertainty on these median absolute differences of ∼0.052 mag, and we find that, statistically, our analysis rules out the nonexistence of intrinsic scatter in the NIR at the 99% level. When comparing the same sets of SN siblings, we observe a median absolute difference inμvalues between siblings of 0.177 mag when using optical data alone as compared to 0.186 mag when using NIR data alone. It is unclear if these results may be due to limited statistics or poor quality NIR data, all of which will be improved with the Nancy Grace Roman Space Telescope.

Astronomy & Astrophysics↗
DOE OSTI2019

Experimental Evaluation of the Fast-Neutron Passive Collar (FNPC) Performance and Comparison with Scintillation Based Collars

This report describes a new method for the verification of fresh low enrichment uranium (LEU) fuel assemblies that can reduce International Atomic Energy Agency (IAEA) inspection time in the field and be completely independent of operator declared burnable poison content. The passive measurement makes unattended mode operation practical, and thus provides the potential for reducing inspector time in the field for future applications. The new simplified passive neutron measurement uses the 238 U spontaneous fission in LEU fresh fuel rods to self-interrogate the 235 U mass in the neighboring rods. The measured response increases linearly with the increase in the linear density (LD) of 235 U. The relatively high efficiency of the passive neutron detector (24%) provides good statistical precision (~ 2%) for the singles, doubles, and triples counting rates in the thermal-neutron mode. The three observables makes possible multiplicity analysis so that both the 235 U LD and the effective burnable poison content can be determined independent of the operator declarations. This closes the potential diversion path of miss-declaration the burnable poison content. This report presents the detector design, the as-built system, the detector parameter measurements, the calibration, and the assay uncertainty estimates. The calibration measurements were performed using the Los Alamos Laboratory mockup PWR fuel assembly. The measurements were made in both the fastneutron mode (with Cd/Gd metal liners) and the thermal-neutron mode (no Cd/Gd liners).

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

On the uncertainty of long-period return values of extreme daily precipitation

Methods for calculating return values of extreme precipitation and their uncertainty are compared using daily precipitation rates over the Western U.S. and Southwestern Canada from a large ensemble of climate model simulations. The roles of return-value estimation procedures and sample size in uncertainty are evaluated for various return periods. We compare two different generalized extreme value (GEV) parameter estimation techniques, namely L-moments and maximum likelihood (MLE), as well as empirical techniques. Even for very large datasets, confidence intervals calculated using GEV techniques are narrower than those calculated using empirical methods. Furthermore, the more efficient L-moments parameter estimation techniques result in narrower confidence intervals than MLE parameter estimation techniques at small sample sizes, but similar best estimates. It should be noted that we do not claim that either parameter fitting technique is better calibrated than the other to estimate long period return values. While a non-stationary MLE methodology is readily available to estimate GEV parameters, it is not for the L-moments method. Comparison of uncertainty quantification methods are found to yield significantly different estimates for small sample sizes but converge to similar results as sample size increases. Finally, practical recommendations about the length and size of climate model ensemble simulations and the choice of statistical methods to robustly estimate long period return values of extreme daily precipitation statistics and quantify their uncertainty.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Characterization of UV optical components for photon detector calibration in liquid argon TPCs

Large liquid argon time projection chambers (LArTPCs)require stable and well-characterized delivery of ultraviolet (UV)light for in situ calibration of photosensors at cryogenictemperatures. This article reports bench-top and cryogenicmeasurements of the optical components used in a UV lightcalibration system, including multi-mode fused-silica fibers,SMA-to-SMA connectors, optical fiber feedthroughs, andlight-diffuser assemblies. Light loss in several fiber types and SMAconnectors was measured across wavelengths from275–970 mm. In addition,light-loss measurements of the tested fibers after severalliquid-nitrogen thermal cycles showed no statistically significantdegradation relative to baseline measurements, and high-rate pulsedexposure (30–90 million pulses from a275 nm LED) likewise showed nomeasurable aging in jacketed fibers. A compact, palm-sized,3D-printed PEEK diffuser housing with stacked UV-grade fused-silicadiffusers yields Lambertian emission and the most uniform angulardistribution. Optical components exhibiting improved UV transmissionwere deployed successfully in multiple DUNE small- and large-scaleprototypes, demonstrating reliable operation of UV light calibrationsystem. These findings inform component selection and calibrationprocedures for achieving reliable, uniform UV light delivery inlarge-scale cryogenic detectors such as DUNE.

Behera, B. [South Dakota Sch. Mines Tech.] (ORCID:↗
DOE OSTI2022

An Algorithm to Calibrate and Correct the Response to Unpolarized Radiation of the X-Ray Polarimeter Onboard IXPE

The Gas Pixel Detector (GPD) is an X-ray polarimeter to fly onboard IXPE and other missions. To correctly measure the source polarization, the response of IXPE’s GPDs to unpolarized radiation has to be calibrated and corrected. In this paper, we describe the way such response is measured with laboratory sources and the algorithm to apply such correction to the observations of celestial sources. The latter allows to correct the response to polarization of single photons, therefore allowing great flexibility in all the subsequent analysis. Our correction approach is tested against both monochromatic and nonmonochromatic laboratory sources and with simulations, finding that it correctly retrieves the polarization up to the statistical limits of the planned IXPE observations.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗
DOE OSTI2022

Sample-Efficient Adaptive Calibration of Quantum Networks Using Bayesian Optimization

All physical systems employed for quantum information tasks must act as unbiased carriers of encoded quantum states. Ensuring such indistinguishability of information carriers is a major challenge in many quantum information applications, including advanced quantum communication protocols. For photons, the workhorses of quantum communication networks, it is difficult to obtain and maintain their indistinguishability because of environment-induced transformations and loss imparted by communication channels, especially in noisy scenarios. Conventional strategies to mitigate these transformations often require hardware or software overhead that is restrictive (e.g., adding noise), infeasible (e.g., on a satellite), or time-consuming for deployed networks. In this work we propose and develop resource-efficient Bayesian optimization techniques to rapidly and adaptively calibrate the indistinguishability of individual photons for quantum networks using only information derived from their measurement. To experimentally validate our approach, we demonstrate the optimization of Hong-Ou-Mandel interference between two photons-a central task in quantum networking- finding rapid, efficient, and reliable convergence towards maximal photon indistinguishability in the presence of high loss and shot noise. We expect our resource-optimized and experimentally friendly methodology will allow fast and reliable calibration of indistinguishable quanta, a necessary task in distributed quantum computing, communications, and sensing, as well as for fundamental investigations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Partial Least Squares, Experimental Design, and Near-Infrared Spectrophotometry for the Remote Quantification of Nitric Acid Concentration and Temperature

Near-infrared spectrophotometry and partial least squares regression (PLSR) were evaluated to create a pleasantly simple yet effective approach for measuring HNO3 concentration with varying temperature levels. A training set, which covered HNO3 concentrations (0.1–8 M) and temperature (10–40 °C), was selected using a D-optimal design to minimize the number of samples required in the calibration set for PLSR analysis. The top D-optimal-selected PLSR models had root mean squared error of prediction values of 1.4% for HNO 3 and 4.0% for temperature. The PLSR models built from spectra collected on static samples were validated against flow tests including HNO 3 concentration and temperature gradients to test abnormal conditions (e.g., bubbles) and the model performance between sample points in the factor space. Based on cross-validation and prediction modeling statistics, the designed near-infrared absorption approach can provide remote, quantitative analysis of HNO 3 concentration and temperature for production-oriented applications in facilities where laser safety challenges would inhibit the implementation of other optical techniques (e.g., Raman spectroscopy) and in which space, time, and/or resources are constrained. The experimental design approach effectively minimized the number of samples in the training set and maintained or improved PLSR model performance, which makes the described chemometric approach more amenable to nuclear field applications.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2024

Deep learning-based predictive models for laser direct drive at the Omega Laser Facility

The rich and complex physics of inertial confinement fusion provides a unique and challenging space for high-fidelity first-principles modeling. Consequently, simulation codes that are used to design experiments are computationally expensive and lack the predictive capability required for extensive parameter exploration in search of a high-performing design for laser direct drive. In this article, we present two deep-learning-based predictive models intended to address these difficulties. The first model (TL DNN) acts as a fast emulator of simulations as well as experiments at the Omega Laser Facility. This model is trained on a simulation database and subsequently calibrated on experimental data using transfer learning. To facilitate the development of this model, an autoencoder is developed to reduce the dimensionality of the input space by compressing the laser pulse input. The model predicts key experimental scalar observables of Omega experiments with high accuracy and minimal computational cost. This deep neural net enables rapid exploration of a high-dimensional input parameter space for an optimal implosion design. The second model (DNN SM+) aims to extend the statistical modeling work of Lees et al. [Phys. Rev. Lett. 127, 105001 (2021)], by increasing the complexity of the model space and allowing for coupling between degradation terms. Since the model capacity of DNN SM+ is higher than the model of Lees et al., DNN SM+ can potentially provide an improvement in predictive capability, and we use this model to provide insight into complicated degradation dependencies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Adaptive Client Selection in Federated Learning: A Network Anomaly Detection Use Case

Federated Learning (FL) has become a ubiquitous approach for training machine learning models on decentralized data, addressing the myriad privacy concerns inherent in traditional centralized methods. However, the efficiency of FL depends on effective client selection and robust privacy preservation mechanisms. Inadequate client selection may lead to suboptimal model performance, while insufficient privacy measures risk exposing sensitive data. This paper proposes a client selection framework for FL that integrates differential privacy and fault tolerance. Our adaptive approach dynamically adjusts the number of selected clients based on model performance and system constraints, ensuring privacy through calibrated noise addition. We evaluate our method on a network anomaly detection use case using the UNSW-NB15 and ROAD datasets. Results show up to a 7% increase in accuracy and a 25% reduction in training time compared to FedL2P. Moreover, we highlight the trade-offs between privacy budgets and model performance, with higher privacy budgets reducing noise and improving accuracy. Our fault tolerance mechanism, while causing a slight performance drop, enhances robustness to client failures. Statistical validation using Mann-Whitney U tests confirms the significance of these improvements (p < 0.05).

Marfo, William [University of Texas at El Paso,Dep↗
DOE OSTI2025

Spatially explicit terrestrial carbon densities for calibrating the carbon cycle in human-Earth system Models

Soil and vegetation carbon stocks play a critical role in human-Earth system models. These stocks (denominated as densities in MgC/ha) affect variables such as land use change emissions and also influence land use change pathways under climate forcing scenarios where terrestrial carbon is assigned a carbon price. Here we present reharmonized soil and vegetation carbon densities both at the 5-arcmin resolution grid cell level and also aggregated to 235 water sheds for 4 land use types (Cropland, Grazed land, Urban land and unmanaged vegetation) and 15 unmanaged land cover types. Moreover, we use the distribution of carbon within and across pixels to define statistical "states" of carbon, once again differentiated by land type. These statistical states are used to define a range of possible carbon values that can be used for defining initial conditions of soil and vegetation carbon in human-Earth system models. We implement these data in a state-of-the-art multi sector dynamics model, namely the Global Change Analysis Model (GCAM), and show that these new data improve several land use responses, especially when terrestrial carbon is assigned a carbon price.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Coupling localized Noah-MP-Crop model with the WRF model improved dynamic crop growth simulation across Northeast China

Croplands play a critical role in regulating the energy and moisture exchanges between the land surface and atmosphere. However, the interactions between cropland and climate are usually poorly represented due to a lack of detailed representation in crop types and field management. Here, we coupled the Noah-MP-Crop model with the state-of-the-art Weather Research and Forecasting (WRF) model to explore and evaluate the crop growth dynamics in response to climate variations across Northeast China. The default parameters of the crop model were not exactly suitable for the agricultural ecosystems in Northeast China. The detailed cropland distribution, and crop phenology parameters including growing degree days (GDD) and planting (harvesting) date were first created using multi-source remote sensing products and reanalysis data, and was then successfully used to simulate the growth and yield for corn and soybean and associated energy exchanges. We also optimized and calibrated other crop parameters using the time-series of the Moderate Resolution Imaging Spectroradiometer (MODIS) land surface products. The modified crop model substantially improved the simulation of crop growth, plant physiology, and biomass accumulation for both corn and soybean. Coupling the localized dynamic crop model into the WRF led to considerable decreases in the simulated mean-absolute-errors (MAEs) and biases of the leaf area index, evapotranspiration, and gross primary production compared with the MODIS observed values. Compared with the statistical yield from each province, the modified crop model underestimated the corn yield from 11.1% to 48.6%, whereas overestimated the soybean yield from 16.5% to 162.6%.

54 ENVIRONMENTAL SCIENCES↗