Search NASA⌕ Search

SEARCH · Search NASA

Results for “simulation output analysis”

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 1,351 records · Page 75

TARDIS. I. A Constrained Reconstruction Approach to Modeling the z ~ 2.5 Cosmic Web Probed by Lyα Forest Tomography

Recent Lyα forest tomography measurements of the intergalactic medium (IGM) have revealed a wealth of cosmic structures at high redshift (z ~ 2.5). In this work, we present the Tomographic Absorption Reconstruction and Density Inference Scheme (TARDIS), a new chronocosmographic analysis tool for understanding the formation and evolution of these observed structures. Here, we use maximum likelihood techniques with a fast nonlinear gravitational model to reconstruct the initial density field of the observed regions. We find that TARDIS allows accurate reconstruction of smaller-scale structures than standard Wiener-filtering techniques. Applying this technique to mock Lyα forest data sets that simulate ongoing and future surveys such as CLAMATO, Subaru PFS, or the ELTs, we are able to infer the underlying matter density field at observed redshift and classify the cosmic web structures. We find good agreement with the underlying truth in both the characteristic eigenvalues and eigenvectors of the pseudo-deformation tensor, with the eigenvalues inferred from 30 m class telescopes correlated at r = 0.95 relative to the truth. As an output of this method, we are able to further evolve the inferred structures to late time (z = 0) and also track the trajectories of coeval z = 2.5 galaxies to their z = 0 cosmic web environments.

79 ASTRONOMY AND ASTROPHYSICS↗

Aeroelastic analysis for helicopter rotors with blade appended pendulum vibration absorbers. Mathematical derivations and program user's manual

Mathematical development is presented for the expanded capabilities of the United Technologies Research Center (UTRC) G400 Rotor Aeroelastic Analysis. This expanded analysis, G400PA, simulates the dynamics of teetered rotors, blade pendulum vibration absorbers and the higher harmonic excitations resulting from prescribed vibratory hub motions and higher harmonic blade pitch control. Formulations are also presented for calculating the rotor impedance matrix appropriate to these higher harmonic blade excitations. This impedance matrix and the associated vibratory hub loads are intended as the rotor blade characteristics elements for use in the Simplified Coupled Rotor/Fuselage Vibration Analysis (SIMVIB). Sections are included presenting updates to the development of the original G400 theory, and material appropriate to the user of the G400PA computer program. This material includes: (1) a general descriptionof the tructuring of the G400PA FORTRAN coding, (2) a detaild description of the required input data and other useful information for successfully running the program, and (3) a detailed description of the output results.

Bielawa, R. L.↗

Performance analysis of a fault inferring nonlinear detection system algorithm with integrated avionics flight data

This paper presents the performance analysis results of a fault inferring nonlinear detection system (FINDS) using integrated avionics sensor flight data for the NASA ATOPS B-737 aircraft in a Microwave Landing System (MLS) environment. First, an overview of the FINDS algorithm structure is given. Then, aircraft state estimate time histories and statistics for the flight data sensors are discussed. This is followed by an explanation of modifications made to the detection and decision functions in FINDS to improve false alarm and failure detection performance. Next, the failure detection and false alarm performance of the FINDS algorithm are analyzed by injecting bias failures into fourteen sensor outputs over six repetitive runs of the five minutes of flight data. Results indicate that the detection speed, failure level estimation, and false alarm performance show a marked improvement over the previously reported simulation runs. In agreement with earlier results, detection speed is faster for filter measurement sensors such as MLS than for filter input sensors such as flight control accelerometers. Finally, the progress in modifications of the FINDS algorithm design to accommodate flight computer constraints is discussed.

Caglayan, A. K.↗

Analysis and use of VAS satellite data

Four interrelated investigations have examined the analysis and use of VAS satellite data. A case study of VAS-derived mesoscale stability parameters suggested that they would have been a useful supplement to conventional data in the forecasting of thunderstorms on the day of interest. A second investigation examined the roles of first guess and VAS radiometric data in producing sounding retrievals. Broad-scale patterns of the first guess, radiances, and retrievals frequently were similar, whereas small-scale retrieval features, especially in the dew points, were often of uncertain origin. Two research tasks considered 6.7 micron middle tropospheric water vapor imagery. The first utilized radiosonde data to examine causes for two areas of warm brightness temperature. Subsidence associated with a translating jet streak was important. The second task involving water vapor imagery investigated simulated imagery created from LAMPS output and a radiative transfer algorithm. Simulated image patterns were found to compare favorably with those actually observed by VAS. Furthermore, the mass/momentum fields from LAMPS were powerful tools for understanding causes for the image configurations.

Fuelberg, Henry E.↗

Sensitivity of UO 2 fuel performance to microstructural evolutions driven by dilute additives

Use of dilute additives to nuclear fuel is being considered to increase the security of commercial fuel management through traceability of fabricated fuel elements. Taggants, as additives are denoted when included for traceability purposes, may also improve fuel performance, as demonstrated in Cr-containing uranium dioxide as described in the literature, and they may also improve fuel safety. In fact, studies have shown that some additives affect fuel material properties such as grain size and density after sintering. Given the possible range of elements that could be used as additives, the impact of such fuel property variations on the fuel’s thermomechanical behavior becomes relevant. These effects can be evaluated through a sensitivity study of standard fuel models to analyze changes in these properties using a fuel performance code. In this work, the BISON code is being used to investigate these effects through a 2D axisymmetric model of smeared UO 2 fuel pellets and ZIRLO® cladding under realistic pressurized water reactor core irradiation conditions. Here, randomly sampled densities and grain sizes within specified ranges are used as input parameters in the simulations, and several fuel model-related outputs are evaluated. The thermomechanical response of the cladding is also addressed in this study. The simultaneous variation of both input parameters offers a more comprehensive path to identify key sensitivities. Outputs explored include temperature, fission gas release, creep, and radial stress. Results show that although most of these outputs are sensitive to grain size to a certain extent, density mainly affects fuel temperature and elastic strain. Furthermore, sensitivities can vary depending on the radial position within the fuel pellet.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

The impact of low-mode symmetry on inertial fusion energy output in the burning plasma state

Indirect Drive Inertial Confinement Fusion Experiments on the National Ignition Facility (NIF) have achieved a burning plasma state with neutron yields exceeding 170 kJ, roughly 3 times the prior record and a necessary stage for igniting plasmas. The results are achieved despite multiple sources of degradations that lead to high variability in performance. Results shown here, for the first time, include an empirical correction factor for mode-2 asymmetry in the burning plasma regime in addition to previously determined corrections for radiative mix and mode-1. Analysis shows that including these three corrections alone accounts for the measured fusion performance variability in the two highest performing experimental campaigns on the NIF to within error. Here we quantify the performance sensitivity to mode-2 symmetry in the burning plasma regime and apply the results, in the form of an empirical correction to a 1D performance model. Furthermore, we find the sensitivity to mode-2 determined through a series of integrated 2D radiation hydrodynamic simulations to be consistent with the experimentally determined sensitivity only when including alpha-heating.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

2D axial-azimuthal particle-in-cell benchmark for low-temperature partially magnetized plasmas

The increasing need to demonstrate the correctness of computer simulations has highlighted the importance of benchmarks. We define in this paper a representative simulation case to study low-temperature partially-magnetized plasmas. Seven independently developed particle-in-cell codes have simulated this benchmark case, with the same specified conditions. The characteristics of the codes used, such as implementation details or computing times and resources, are given. First, we compare at steady-state the time-averaged axial profiles of three main discharge parameters (axial electric field, ion density and electron temperature). We show that the results obtained exhibit a very good agreement within 5% between all the codes. As E x B discharges are known to cause instabilities propagating in the direction of electron drift, an analysis of these instabilities is then performed and a similar behaviour is retrieved between all the codes. A particular attention has been paid to the numerical convergence by varying the number of macroparticles per cell and we show that the chosen benchmark case displays a good convergence. Detailed outputs are given in the supplementary data, to be used by other similar codes in the perspective of code verification.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Energy-Efficient and Secure S-Box circuit using Symmetric Pass Gate Adiabatic Logic

Differential Power Analysis (DPA) attack is considered to be a main threat while designing cryptographic processors. In cryptographic algorithms like DES and AES, S-Box is used to indeterminate the relationship between the keys and the cipher texts. However, S-box is prone to DPA attack due to its high power consumption. In this paper, we are implementing an energy-efficient 8-bit S-Box circuit using our proposed Symmetric Pass Gate Adiabatic Logic (SPGAL). SPGAL is energy-efficient as compared to the existing DPAresistant adiabatic and non-adiabatic logic families. SPGAL is energy-efficient due to reduction of non-adiabatic loss during the evaluate phase of the outputs. Further, the S-Box circuit implemented using SPGAL is resistant to DPA attacks. The results are verified through SPICE simulations in 180nm technology. SPICE simulations show that the SPGAL based S-Box circuit saves upto 92% and 67% of energy as compared to the conventional CMOS and Secured Quasi-Adiabatic Logic (SQAL) based S-Box circuit. From the simulation results, it is evident that the SPGAL based circuits are energy-efficient as compared to the existing DPAresistant adiabatic and non-adiabatic logic families. In nutshell, SPGAL based gates can be used to build secure hardware for lowpower portable electronic devices and Internet-of-Things (IoT) based electronic devices.

Kumar, Dinesh↗

Bounding the Chemistry Conditions in the HM Solvent Extraction Process Criticality Safety Analyses

An approach new to the H-Canyon facility criticality safety analysis has been performed using highly detailed models of the H-Modified (HM) process mixer settlers. The HM process uses tributyl phosphate (TBP) in organic n-paraffin to extract and purify uranium from used nuclear fuel. This occurs in two cycles of counter-current solvent extraction. The aqueous stream is a weak to moderate nitric acid and the feed stock is an aqueous stream. The first Cycle, consisting of two 16-stage banks (A and B-Banks) and one 12-stage bank (C-Bank), is a three step process that separates uranium from the fuel matrix, fission products, and most of the transuranic actinide content. The first Cycle aqueous product is concentrated, adjusted for acidity and fed to the second Cycle. The second Cycle, consisting of one 16-stage bank (D-Bank) and one 12-stage bank (E-Bank), is a two-step process that further purifies the uranium and removes to a great extent remaining actinides and fission products Each stage of the mixer settler banks is modeled as three unique compositions: the mixer as a combined composition of aqueous and organic components and the settler one aqueous and one organic component. This work addresses what is believed to be normal, credible abnormal, and beyond credible abnormal conditions for a particular postulated series of process parameter upsets. The level of detail and the source for the stage-wise compositions are both new and unique to criticality safety analyses of the HM process. Initial work was presented in Ref. 1 and has since been refined. A predictive chemistry model for flowsheet development has been used to solve the solvent extraction equations for the prescribed conditions and the output of that model was then transformed into compositions for each mixture in each stage. The remaining question upon which to achieve consensus is if the approach was appropriately conservative by concluding the chemistry is appropriately bounded. Computational modeling for criticality calculations was performed using KENO-VI in the SCALE 6.1 code package. Computational modeling of the chemistry was performed using the SEPHIS Modification 4 code.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Uncertainty quantification and sensitivity analysis of a nuclear thermal propulsion reactor startup sequence

The research presented in this article describes progress in applying stochastic methods, uncertainty quantification, parametric studies, and variance-based sensitivity analysis (also known as Sobol sensitivity analysis) to a full-core model of a nuclear thermal propulsion (NTP) system simulated via the radiation transport code Griffin to simulate neutronics. Our goal is to develop a reduced-order (surrogate) model that can be rapidly sampled with perturbations to multiple input parameters. In this NTP system, reactivity and power feedback affect the rotation of control drums (CDs), which is itself controlled by a hybrid proportional-integral-derivative (PID) controller actuated by the power demand and reactivity feedback from the numerical model. This model uses reactor kinetic feedback (mean generation time [Λ] and effective delayed neutron fraction [ β eff ] from a transient Griffin simulation executed via Griffin’s improved quasi-static solver to provide the kinetic parameters) as inputs to functions that control the CD rotation angle. By investigating numerous stochastic approaches, we developed a dual-purpose surrogate model of the NTP system, using polynomial regression in the Multiphysics Object-Oriented Simulation Environment (MOOSE) Stochastic Tools Module (STM). The trained model can be rapidly sampled while simultaneously perturbing various input parameters, such as coefficients on the PID control or temperature (directly affecting the neutron cross section). The surrogate model delivers accurate (within 5%) results at speeds orders of magnitude faster (minutes, not days of computational time) than the base model. Once the surrogate model has been trained, distributions of the uncertain parameters can be changed at will to investigate the effects of perturbing multiple inputs as well as the effects of these inputs on the model output. For example, coefficients used in the PID control system may vary due to some type of physical interference, or uncertainty may exist in the temperature of the neutron cross sections in various regions of the reactor. A distribution can be placed on these parameters, and operational boundaries can be determined. The goal of this work is to support development of an advanced control system for operating CDs in a functioning NTP system. This work is a scoping study of the MOOSE STM.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗

Iterative sampling of expensive simulations for faster deep surrogate training

Deep neural network (DNN) surrogates of expensive physics simulations are enabling a rapid change in the way that common experimental design and analysis tasks are approached. Surrogate models allow simulations to be performed in parallel and separately from downstream tasks, thereby enabling analyses that would be impossible with the simulation in-the-loop; surrogates based on DNNs can effectively emulate diverse non-scalar data of the types collected in fusion and laboratory-astrophysics experiments. The challenge is in training the surrogate model, for which large ensembles of physics simulations must be run, preferably without wasting computational effort on uninteresting simulations. Here, in this paper, we present an iterative sampling scheme that can preferentially propose simulations in interesting regions of parameter space without neglecting unexplored regions, allowing high-quality and wide-ranging surrogate models to be trained using 2–3 times fewer simulations compare to space-filling designs. Our approach uses an explicit importance function defined on the simulation output space, balanced against a measure of simulation density which serves as a proxy for surrogate accuracy. It is easy to implement and can be tuned to find interesting simulations early in the study, allowing surrogates to be trained quickly and refined as new simulations become available; this represents an important step towards the routine generation of deep surrogate models quickly enough to be truly relevant to experimental work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Modular HF Isolated MV String Inverters Enable a New Paradigm for Large PV Farms

The “Modular HF Isolated MV String Inverters Enable a New Paradigm for Large PV Farms” project focuses on exploring alternative power converter and system-level plant configurations to achieve the lowest cost and highest energy output for a given solar plus storage (e.g., PV + battery) plant, including the use of medium voltage (MV) collection while taking into account detailed models of all elements. To realize this objective, four approaches were utilized (i) employ a novel Medium Voltage String Inverter (MVSI) topology (soft switching solid state transformer – S4T) to convert 1000 Vdc to 4.16 kVac; (ii) plant collection using standard, low-cost overhead MV distribution network; (iii) enable energy storage integration without additional converter cost to achieve dispatchability of the PV resource; and (iv) provide advanced functionality (autonomous operation, track ISO signals for dynamic balancing and ancillary services, and PV farm operation as a virtual grid resource). Subsequently and in alignment with the previously mentioned approaches, the project was structured in five efforts (i) S4T MVSI simulation and design; (ii) system analysis and storage optimization; (iii) financial analysis; (iv) power converter prototype build and test; and (v) regulatory and commercial impact study. The outcomes provided by each effort can be summarized as follows (i) Project explored the use of MV AC distribution architecture for hybrid PV+storage utility-scale PV farms; (ii) Detailed loss and LCOE analysis for AC and DC side BESS architecture, including multiple converter topologies, as well as for proposed MVSI/MDCT systems; (iii) MVSI was built and holds promise but needs lower-cost high-voltage Si-C devices, which does not seem possible in the near term; (iv) MDCT provides a simpler modular building block – validated through HIL and farm level modeling, simulation and experimental validation; (v) 300 kVA MDCT prototype built and tested, technology is being commercialized; and (vi) Regulatory model of utility building PV plants, where PV panels are treated as DC generation (IPP), seems viable and can allow improved grid integration.

14 SOLAR ENERGY↗

DEEPER: An Intergrated Platform for Deeper Roots

Crops with deeper roots would have multiple benefits, including better drought tolerance, reduced requirement for nitrogen fertilizer, and better sequestration of atmospheric CO 2 . DEEPER is an integrated platform of phenomic, genomic, and in silico technologies to generate maize lines with deeper roots. DEEPER is: LEADER (Leaf Elemental Accumulation from Deep Roots) is a breakthrough technology to nondestructively measure rooting depth by using the plant itself as a sensor. LEADER uses handheld X-ray Fluorescence spectrometry to quantify foliar accumulation of elements that are differentially distributed in the soil profile. LEADER is nondestructive and is orders of magnitude cheaper, faster, and more precise than any competing assay of rooting depth in the field. LEADER is able to distinguish deep-rooted from shallow-rooted maize lines in the field without the need for costly and noisy soil coring. RootRobot/DIRT3D, to automatically phenotype root architecture in any field, combining RootRobot, a mechatronics platform to excavate, clean, section, and image mature root crowns, with DIRT3D, software to quantify architectural traits in 3D. Anatomics, a high-throughput platform to phenotype root anatomy, combining LAT 2.0, a technology for 3D imaging of root anatomy and composition, with RootScan3D, software to automatically extract 3D anatomical and cell wall composition metrics from LAT 2.0 output. Using this platform we discovered two novel root traits, parenchyma cell wall thickness and multiseriate cortical sclerenchyma, that improve rooting depth and drought tolerance in maize and wheat. OpenSimRoot/Deep, software to simulate root interaction with hard subsoils. Using this platform we discovered novel concepts regarding how to increase crop rooting depth by modulating how individual root axes respond to hard soil. DeepGenes, a toolkit of genes, parent lines, and genomic selection strategies to enable breeding hybrids with deeper roots. We discovered 3 novel root genes that increase rooting depth in maize and wheat. DEEPER discovered novel root phenotypes for deeper rooting, and delivered validated ideotypes for deeper-rooted maize; novel technologies to rapidly assess root depth, root architecture and anatomy in field-grown plants; novel software tools for root modeling and 3D image analysis of root architecture and anatomy; and validated genes and genomic selection models to deploy traits for deeper rooting in maize breeding. Each DEEPER technology is transformative in its own right, and exceeds existing technologies. They are mutually synergistic, deployable for field-grown plants, and are ready for application. The phenotyping and modeling technologies are readily applicable to many crops, and genetic leads in maize may have utility in other grasses. Taken as a whole they represent a transformative platform to develop deeper-rooted crops, with greater drought tolerance, reduced fertilizer requirement, and greater carbon sequestration.

59 BASIC BIOLOGICAL SCIENCES↗

Structure Computation of Quiet Spike[Trademark] Flight-Test Data During Envelope Expansion

System identification or mathematical modeling is used in the aerospace community for development of simulation models for robust control law design. These models are often described as linear time-invariant processes. Nevertheless, it is well known that the underlying process is often nonlinear. The reason for using a linear approach has been due to the lack of a proper set of tools for the identification of nonlinear systems. Over the past several decades, the controls and biomedical communities have made great advances in developing tools for the identification of nonlinear systems. These approaches are robust and readily applicable to aerospace systems. In this paper, we show the application of one such nonlinear system identification technique, structure detection, for the analysis of F-15B Quiet Spike(TradeMark) aeroservoelastic flight-test data. Structure detection is concerned with the selection of a subset of candidate terms that best describe the observed output. This is a necessary procedure to compute an efficient system description that may afford greater insight into the functionality of the system or a simpler controller design. Structure computation as a tool for black-box modeling may be of critical importance for the development of robust parsimonious models for the flight-test community. Moreover, this approach may lead to efficient strategies for rapid envelope expansion, which may save significant development time and costs. The objectives of this study are to demonstrate via analysis of F-15B Quiet Spike aeroservoelastic flight-test data for several flight conditions that 1) linear models are inefficient for modeling aeroservoelastic data, 2) nonlinear identification provides a parsimonious model description while providing a high percent fit for cross-validated data, and 3) the model structure and parameters vary as the flight condition is altered.

Kukreja, Sunil L.↗

Analysis of the Planetary Boundary Layer Height During DISCOVER-AQ Baltimore–Washington, D.C., with Lidar and High-Resolution WRF Modeling

The daytime planetary boundary layer (PBL) was examined for the Deriving Information on Surface Conditions from Column and Vertically Resolved Observations Relevant to Air Quality (DISCOVER-AQ) Baltimore (Maryland)–Washington, D.C., campaign of July 2011 using PBL height (PBLH) retrievals from aerosol backscatter measurements from ground-based micropulse lidar (MPL), the NASA Langley Research Center airborne High Spectral Resolution Lidar-1 (HSRL-1), and the Cloud–Aerosol Lidar with Orthogonal Polarization (CALIOP) on the Cloud–Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) satellite. High-resolution Weather Research and Forecasting (WRF) Model simulations with horizontal grid spacing of 1 km and different combinations of PBL schemes, urban parameterization, and sea surface temperature inputs were evaluated against PBLHs derived from lidars, ozonesondes, and radiosondes. MPL and WRF PBLHs depicted a growing PBL in the morning that reached a peak height by midafternoon. WRF PBLHs calculated from gridded output profiles generally showed more rapid growth and higher peak heights than did the MPLs, and all WRF–lidar differences were dependent on model configuration, PBLH calculation method, and synoptic conditions. At inland locations, WRF simulated an earlier descent of the PBL top in the afternoon relative to the MPL retrievals and radiosonde PBLHs. At Edgewood, Maryland, the influence of the Chesapeake Bay breeze on the PBLH was captured by both the ozonesonde and WRF data but generally not by the MPL PBLH retrievals because of generally weaker gradients in the aerosol backscatter profile and limited normalized relative backscatter data near the top height of the marine layer.

Hegarty, Jennifer D.↗

A Generalized Transformer-Based Pulse Detection Algorithm

Pulse-like signals are ubiquitous in the field of single molecule analysis, e.g., electrical or optical pulses caused by analyte translocations in nanopores. The primary challenge in processing pulse-like signals is to capture the pulses in noisy backgrounds, but current methods are subjectively based on a user-defined threshold for pulse recognition. Here, we propose a generalized machine-learning based method, named pulse detection transformer (PETR), for pulse detection. PETR determines the start and end time points of individual pulses, thereby singling out pulse segments in a time-sequential trace. It is objective without needing to specify any threshold. It provides a generalized interface for downstream algorithms for specific application scenarios. PETR is validated using both simulated and experimental nanopore translocation data. It returns a competitive performance in detecting pulses through assessing them with several standard metrics. Finally, the generalization nature of the PETR output is demonstrated using two representative algorithms for feature extraction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Global warming intensity of biofuel derived from switchgrass grown on marginal land in Michigan

Energy crops for biofuel production, especially switchgrass (Panicum virgatum), are of interest from a climate change perspective. Here, we use outputs from a crop growth model and life cycle assessment (LCA) to examine the global warming intensity (GWI; g CO 2 MJ –1 ) and greenhouse gas (GHG) mitigation potential (Mg CO 2 year –1 ) of biofuel systems based on a spatially explicit analysis of switchgrass grown on marginal land (abandoned former cropland) in Michigan, USA. We find that marginal lands in Michigan can annually produce over 0.57 hm 3 of liquid biofuel derived from nitrogen-fertilized switchgrass, mitigating 1.2–1.5 Tg of CO 2 year –1 . About 96% of these biofuels can meet the Renewable Fuel Standard (60% reduction in lifecycle GHG emissions compared with conventional gasoline; GWI ≤37.2 g CO 2 MJ –1 ). Furthermore, 73%–75% of these biofuels are carbon-negative (GWI less than zero) due to enhanced soil organic carbon (SOC) sequestration. However, simulations indicate that SOC levels would fail to increase and even decrease on the 11% of lands where SOC stocks >>200 Mg C ha –1 , leading to carbon intensities greater than gasoline. Results highlight the strong climate mitigation potential of switchgrass grown on marginal lands as well as the needs to avoid carbon rich soils such as histosols and wetlands and to ensure that productivity will be sufficient to provide net mitigation.

09 BIOMASS FUELS↗

Influence of Lake Ice Biases in Reanalysis Data on Downscaled Climate Simulations over the Great Lakes Region

This data package contains observation-based and model-simulated datasets (all provided in NetCDF format) for evaluating how wintertime lake-ice representation affects regional weather and climate over the Laurentian Great Lakes (freshwater lake ecosystem) during the high–ice-cover winter of 2009. The observational component includes: (1) Stage IV gridded precipitation at 4 km, hourly resolution for January–February 2009 over the Great Lakes region (radar–gauge multisensor precipitation analyses); (2) Great Lakes Surface Environmental Analysis (GLSEA) satellite-derived lake-ice coverage at 1.3 km, daily resolution for the 2009 winter months, providing ice coverage over Lakes Superior, Michigan, Huron, Erie, and Ontario; and (3) in situ measurements at the Standard Rock site on Lake Superior from the Great Lakes Evaporation Network (GLEN) at hourly resolution, including near-surface atmospheric variables and sensible and latent heat fluxes (air–lake exchange) at a fixed point location. The modeling component provides corresponding fields from two simulations, both archived at 4 km, hourly resolution: a standalone Weather Research Forecasting model (WRF) run driven by the European Centre for Medium-Range Weather Forecasts (ECMWF) Reanalysis 5 (ERA5), and a two-way coupled model using WRF and the Finite Volume Community Ocean Model (WRF-FVCOM, a 3-D hydrodynamic lake model). These outputs include variables relevant to air–lake interaction and lake-effect processes (e.g., near-surface temperature, humidity, wind, precipitation, and surface turbulent fluxes), enabling direct comparison with the observational datasets. Users can analyze and visualize these NetCDF files with common tools such as Python (e.g., xarray, netCDF4, numpy, pandas), NCO/CDO, Panoply, or ncview; NetCDF variables can also be converted to other formats (e.g., CSV, GeoTIFF) using these utilities.

EARTH SCIENCE > ATMOSPHERE > ATMOSPHERIC TEMPERATU↗