Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “precipitation, modeling, diagnostics”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS2021

Preconditioning of overcast-to-broken cloud transitions by riming in marine cold air outbreaks

Marine cold air outbreaks (CAOs) commonly form overcast cloud decks that transition into broken cloud fields downwind, dramatically altering the local radiation budget. In this study, we investigate the impact of frozen hydrometeors on these transitions. We focus on a CAO case in the NW Atlantic, the location of the multi-year flight campaign ACTIVATE (Aerosol Cloud meTeorology Interactions oVer the western ATlantic Experiment). We use MERRA-2 (Modern-Era Retrospective analysis for Research and Applications, version 2) reanalysis fields to drive large eddy simulations with mixed-phase two-moment microphysics in a Lagrangian framework. We find that transitions are triggered by substantial rain (rainwater paths >25 g m−2), and only simulations that allow for aerosol depletion result in sustained breakups, as observed. Using a range of diagnostic ice nucleating particle concentrations, Ninp, we find that increasing ice progressively accelerates transitions, thus abbreviating the overcast state. Ice particles affect the cloud-topped boundary layer evolution, primarily through riming-related processes prior to substantial rain, leading to (1) a reduction in cloud liquid water, (2) early consumption of cloud condensation nuclei, and (3) early and light precipitation cooling and moistening below cloud. We refer to these three effects collectively as “preconditioning by riming”. Greater boundary layer aerosol concentrations available as cloud condensation nuclei (CCN) delay the onset of substantial rain. However, cloud breakup and low CCN concentration final stages are found to be inevitable in this case, due, primarily, to liquid water path buildup. An ice-modulated cloud transition speed suggests the possibility of a negative cloud–climate feedback. To address prevailing uncertainties in the model representation of mixed-phase processes, the magnitude of ice formation and riming impacts and, thereby, the strength of an associated negative cloud–climate feedback process, requires further observational evaluation by targeting riming hot spots with in situ imaging probes that allow for both the characterization of ice particles and abundance of supercooled droplets.

Marine cold air outbreaks↗
NASA NTRSDate not supplied

Auroral ovals versus flare ribbons: Optical footprints of explosive magnetic field reconfiguration in Earth’s magnetosphere and the Sun’s corona

Despite drastically different plasma conditions in the two systems, eruptive energy release events in solar active regions and in Earth’s magnetosphere have a lot in common. Both solar active region eruptions and the substorm activity in the nightside magnetosphere are known to be associated with an abrupt large-scale reorganization of the supporting magnetic systems enabling the conversion and dissipation of the stored free magnetic energy. The rate of energy accumulation in either system tends to be by orders of magnitude slower than the energy release rate, leading to the formation of an impulsive loading-unloading cycle involving a relatively long pre-reconnection phase (substorm growth phase, pre-flare stressing of an active region magnetic field), a short and explosive reconnection phase (substorm breakup and expansion, eruptive flares and CMEs), and a prolonged post-reconnection phase (substorm recovery phase and the post-flare active region reconfiguration). In Earth’s magnetosphere, this loading-unloading cycle is diagnosed by monitoring characteristic changes in the particle precipitation pattern in the northern and southern auroral ovals marking the boundary between the open and close magnetic flux. The position, shape, structure, and dynamics of the auroral borealis were shown to communicate critically important information about an upcoming substorm, the location and timing of the initial reconnection onset, and the subsequent return of the magnetosphere to a lower-energy state. On the Sun, a similar role is played by the chromospheric flare ribbons which enable a detailed diagnostics of the overlying reconnecting magnetic field. In this talk, we’ll discuss similarities and differences in using auroral and flare ribbon observations for monitoring eruptive energy release in respectively Earth’s magnetosphere and solar corona. Special emphasis will be placed on the quantitative interpretation of the two types of optical signatures for empirical diagnostics of the underlying reconnection activity and for constraining its theoretical models. We’ll also compare forecasting capabilities of the two types of optical signatures in the context of space weather prediction objectives.

V. M. Uritskiy↗
DOE OSTI2022

Sensitivity of Cloud-Radiative Effects to Cloud Fraction Parametrizations in Tropical, Mid-Latitude and Arctic Kilometre-Scale Simulations

The Regional Atmosphere (RA) configuration of the Met Office Unified Model currently requires different cloud fraction parametrizations (CFP) for tropical and midlatitude simulations. To explore the scope for unification of these two RA configurations, this paper presents a detailed evaluation of simulations over tropical, midlatitude and arctic domains, with two different diagnostic CFPs, a prognostic CFP, and no CFP at all. Furthermore, a novel, hybrid approach was used that treats liquid cloud diagnostically and ice cloud prognostically. Using observations from three U.S. Department of Energy Atmospheric Radiation Measurement supersites, it is shown that none of these CFPs stands out as superior over all domains. Over the frequently overcast Arctic, the all-or-nothing approach best captures the cloud radiative properties. Conversely, CFPs are of benefit in regions with frequent partial cloudiness, such as the midlatitudes and the tropics. However, their improved cloud radiative properties often hide an error compensation. All models underestimate overcast, low-base cloud with small water paths in convective environments. In addition, mid-latitude overcast, low-base, optically thick clouds in the morning, possibly associated with overnight convection, are frequently too broken. Diagnostic schemes compensate for these errors by producing spurious, scattered afternoon cloud, which could be due to a correct cloud response to too eager convective initiation. Winter clouds over the midlatitudes are improved when liquid cloud is represented diagnostically with a bimodal saturation-departure PDF, without error compensation. While it is difficult to unify the RA across the globe around a single CFP scheme, the newly proposed Hybrid scheme performs reasonably well for cloud cover across all regions. It also exhibits SW biases that are smaller than most other configurations and is less affected by excessive liquid water paths and compensating errors than fully diagnostic schemes. Surface precipitation is fairly insensitive to the CFP in the simulations shown here.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1986

Results of various studies made with the NCAR Thermospheric General Circulation Model (TGCM) (invited review)

The NCAR thermospheric general circulation model (TGCM) has been used for a variety of thermospheric dynamic studies. It has also been used to compare model predictions with measurements made from various ground-based Fabry-Perot interferometer stations, incoherent scatter radar stations and the Dynamics Explorer satellites. The various input and output features of the model are described. These include the specification of solar EUV fluxes, and descriptions of empirical models to specify auroral particle precipitation, ion drag, and magnetospheric convection. Results are presented for solstice conditions giving the model perturbation temperature and circulation response to solar heating forcing alone and also with the inclusion of magnetospheric convections for two different dawn-dusk potential drops, 20 and 60 kV respectively. Results at two constant pressure levels Z =+1 at 300 km and Z= -4 at 120 km are presented for both the winter and summer polar cap regions. The circulation over the Northern Hemisphere polar cap in both the upper and lower thermosphere are presented along with a figure showing that the circulation is mainly a non-divergent irrotational flow responding to ion drag. The results of a study made on the Southern Hemisphere polar cap during October 1981 where Dynamics Explorer satellite measurements of winds, temperature and composition are compared to TGCM predictions are also presented. A diagnostic package that has been developed to analyze the balance of forces operating in the TGCM is presented next illustrating that in the F-region ion drag and pressure provide the main force balance and in the E-region ion drag, pressure and the coriolis forces provide the main balance. The TGCM prediction for the June 10, 1983 total solar eclipse are next presented showing a thermospheric disturbance following the path of totality. Finally, results are presented giving the global circulation, temperature and composition structure of the thermosphere for solar minimum conditions at equinox with 60 kV magnetospheric convection forcing at high latitudes.

Roble, R. G.↗
NASA NTRS2006

The Strings of Eta Carina: The HST/STIS Spectra and [Ca II]

Long linear, filamentary ejecta, are found to move at very high velocity external to the Homunculus, the circumstellar hourglass-shaped ejecta surrounding Eta Carinae. The origin of the strings is a puzzle. As an example, the Weigelt Blobs have N at 10X solar and C, O at 0.01X solar abundance, along with He/H significantly enhanced. This abundance pattern is evidence for extreme CNO-processing. Similarly, the Strontium Filament has Ti/Ni at 100X solar, presumably due to the lack of oxygen to form Ti-oxide precipitates onto dust grains. We have obtained 2-D spectra with the HST/STIS of the Strontium Filament and a portion of a string. These deep spectral exposures, at moderate dispersion, span much of the near red spectral region from 5000 to 9000A. We have identified twelve emission lines in these spectra with proper velocities and spatial structure of this string and obtained line ratios for [Ca II] (7293/7325A) and [Fe Ill (7157/8619A) which are useful for determining physical conditions in this nebulosity. In an attempt to use the [Ca II] ratio to determine the physical parameters, and ultimately the abundances in the strings, we have constructed a statistical equilibrium model for Ca II , including radiative and collisional rates. These results incorporate our newly calculated atomic data for levels n = 3,4,5 and 6 configurations of Ca II. The aim is to compute the [Ca II] line ratios and use them as a diagnostic of the physical parameters. Using the [Fe II] ratio we find that for Te=10,000 K, the electron density is Ne approx.10(exp 6)/cu cm. We plan to use the [Ca II] ratio to confirm this result. Then, we will extend the use of this multilevel model Ca II atom to study the physical conditions of the Strontium filament where eight lines of Ca II, both allowed and forbidden, had been identified. With the physical conditions determined, we will be able to derive reliable estimates for the gas phase abundances in the strings.

Melendez, M. B.↗
NASA NTRS2007

Evidence for Interannual to Decadal Variations in Hadley and Walker Circulations and Links to Water and Energy Fluxes

Mass and energy transports associated with the Hadley and Walker circulations are important components of the earth s climate system and are strongly linked to hydrologic processes. Interannual to decadal variation in these flows likely signify a combination of natural climate noise as well as a response to anthropgenic forcing. There remains considerable uncertainty in quantifying variations in these flows. Evidence in the surface pressure record supports a weakening of the Walker circulation over the Pacific in recent decades. Conversely the NCEP / NCAR and ERA 40 reanalyses indicate that the Hadley circulation has increased in strength over the last two decades, though these analyses depict significantly different mass circulation changes. Interestingly, the NCEP - II / DOE reanalysis contains essentially no Hadley circulation changes. Most climate model integrations anticipate a weakening of both tropical circulations associated with stronger static stability. Clearly there is much uncertainty not only with the mass transports, but also how they are linked to water and energy balance of the planet through variations in turbulent heat and radiative fluxes and horizontal exports / imports of energy. Here we examine heat and water budget variations from a number of reanalysis products and focus on the linear and nonlinear response of ENSO warm and cold events as opportunities to study budget variations over the past 15-20 years. Our analysis addresses such questions as To what extent do Hadley and Walker Cell variations compensate each other on mass and energy transport? Do static stability adjustments appear to constrain fractional precipitation response vs. fractional water vapor response? We appeal to constraints offered by GPCP precipitation, SSWI ocean evaporation estimates, and ISCCP-FD radiative fluxes, and other satellite data sets to interpret and confirm reanalysis-based diagnostics. Using our findings we also attempt to place in context the recent findings that tropical ocean evaporation increased by order 5% or more during the 1990s, reconciling this with GPCP precipitation variations.

Robertson, Franklin↗
DOE OSTI2023

On the Role of Macrophysics and Microphysics in Km-Scale Simulations of Mixed-Phase Clouds During Cold Air Outbreaks

Regional atmospheric models struggle to maintain supercooled liquid in mixed-phase clouds during polar cold-air outbreaks (CAOs). Previous studies focused on the parameterization of aerosol, microphysics and turbulence to understand the origin of this widespread model bias. This study investigates the role of sub-grid cloud fraction parameterizations (CFP) in the simulation of mixed-phase cloud. Km-scale simulations are performed for a large number of CAO cases over Norway, for which continuous ground observations were collected at one site over 6 months. We use a novel analysis that attributes the cloud-radiative errors to deficiencies in specific cloud regimes. We show that the CFP matters for cloud-radiative effects in CAOs, but that it is probably not the primary cause of the lack of liquid water in simulated mixed-phase clouds. Of all the CFP sensitivities explored in this study, the prognostic representation of both liquid and ice shows most promise in increasing the liquid water path. A newly proposed hybrid CFP with prognostic frozen and diagnostic liquid cloud fraction reproduces some of the benefits of the prognostic scheme at reduced cost and complexity. All simulations appear to produce too large precipitation particles. As a result, reducing the snow deposition rate increases the liquid water path and improves the precipitation size distributions, while simulations are less sensitive to reduced riming rates.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Weathering Incongruence in Mountainous Mediterranean Climates Recorded by Stream Lithium Isotope Ratios

Lithium isotope ratios (δ 7 Li) of rivers are increasingly serving as a diagnostic of the balance between chemical and physical weathering contributions to overall landscape denudation rates. Here, we show that intermediate weathering intensities and highly enriched stream δ 7 Li values typically associated with lowland floodplains can also describe small upland watersheds subject to cool, wet climates. This behavior is revealed by stream δ 7 Li between +22.4 and +23.5‰ within a Critical Zone observatory located in the Cévennes region of southern France, where dilute stream solute concentrations and significant atmospheric deposition otherwise mask evidence of incongruence. The water-rock reaction pathways underlying this behavior are quantified through a multicomponent, isotope-enabled reactive transport model. Using geochemical characterization of soil profiles, bedrock, and long-term stream samples as constraints, we evolve the simulation from an initially unweathered granite to a steady state weathering profile which reflects the balance between (a) fluid infiltration and drainage and (b) bedrock uplift and soil erosion. Enriched stream δ 7 Li occurs because Li is strongly incorporated into actively precipitating secondary clay phases beyond what prior laboratory experiments have suggested. Chemical weathering incongruence is maintained despite relatively slow reaction rates and moderate clay accumulation due to a combination of two factors. First, reactive primary mineral phases persist across the weathering profile and effectively “shield” the secondary clays from resolubilization due to their greater solubility. Second, the clays accumulating in the near-surface profile are relatively mature weathering byproducts. These factors promote characteristically low total dissolved solute export from the catchment despite significant input of exogenous dust.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2002

Better Weather Prediction and Climate Diagnostics Using Rainfall Measurements from Space

Progress in understanding of the role of water in global weather and climate is currently limited by our knowledge of the spatial and temporal variability of primary hydrological fields such as precipitation and evaporation. The Tropical Rainfall Measuring Mission (TRMM) has recently demonstrated that use of microwave-based rainfall observations from space in data assimilation can provide better climate data sets and improve short-range weather forecasting. At NASA, we have been exploring non-traditional approaches to assimilating TRMM Microwave Imager (TMI) and Special Sensor Microwavehager (SSM/I) surface rain rate and latent heating profile information in global systems. In this talk we show that assimilating microwave rain rates using a continuous variational assimilation scheme based on moisture tendency corrections improves quantitative precipitation estimates (QPE) and related clouds, radiation energy fluxes, and large-scale circulations in the Goddard Earth Observing System (GEOS) reanalyses. Short-range forecasts initialized with these improved analyses also yield better QPE scores and storm track predictions for Hurricanes Bonnie and Floyd. We present a status report on current efforts to assimilate convective and stratiform latent heating profile information within the general variational framework of model parameter estimation to seek further improvements. Within the next 5 years, there will be a gradual increase in microwave rain products available from operational and research satellites, culminating to a target constellation of 9 satellites to provide global rain measurements every 3 hours with the proposed Global Precipitation Measurement (GPM) mission in 2007/2008. Based on what has been learned from TRMM, there is a high degree of confidence that these observations can play a'major role in improving weather forecasts and producing better global datasets for understanding the Earth's water and energy cycle. The key to success is to adopt an integrated approach to retrieval, validation, modeling, and data assimilation in a coordinated end-to-end observation-application program.

Hou, Arthur↗
NASA NTRS2002

A Variable-Resolution Stretched-Grid General Circulation Model and Data Assimilation System with Multiple Areas of Interest: Studying the Anomalous Regional Climate Events of 1998

The new stretched-grid design with multiple (four) areas of interest, one at each global quadrant, is implemented into both a stretched-grid GCM (general circulation model) and a stretched-grid data assimilation system (DAS). The four areas of interest include: the U.S./Northern Mexico, the El Nino area/Central South America, India/China, and the Eastern Indian Ocean/Australia. Both the stretched-grid GCM and DAS annual (November 1997 through December 1998) integrations are performed with 50 km regional resolution. The efficient regional down-scaling to mesoscales is obtained for each of the four areas of interest while the consistent interactions between regional and global scales and the high quality of global circulation, are preserved. This is the advantage of the stretched-grid approach. The global variable resolution DAS incorporating the stretched-grid GCM has been developed and tested as an efficient tool for producing regional analyses and diagnostics with enhanced mesoscale resolution. The anomalous regional climate events of 1998 that occurred over the U.S., Mexico, South America, China, India, African Sahel, and Australia are investigated in both simulation and data assimilation modes. Tree assimilated products are also used, along with gauge precipitation data, for validating the simulation results. The obtained results show that the stretched-grid GCM and DAS are capable of producing realistic high quality simulated and assimilated products at mesoscale resolution for regional climate studies and applications.

Fox-Rabinovitz, Michael S.↗
DOE OSTI2023

Isotopic Constraints on the Nature of Primary Precipitates in Archean–Early Paleoproterozoic Iron Formations from Determinations of the Iron Phonon Density of States of Greenalite and 2L- and 6L-Ferrihydrite

Iron formations (IFs) are chemical sedimentary rocks that were widely deposited before the Great Oxidation Event (GOE) around 2.4- 2.2 Ga. It is generally thought that IFs precipitated as hydrated Fe 3+ oxides (HFOs) such as ferrihydrite following surface oxidation of Fe 2+ -rich, anoxic deep waters. This model often implicates biological oxidation and underpins reconstructions of marine nutrient concentrations. However, nanoscale petrography indicates that an Fe 2+ silicate, greenalite, is a common primary mineral in well-preserved IFs, motivating an alternative depositional model of anoxic ferrous silicate precipitation. It is unclear, however, if Fe 2+ -rich silicates can produce the Fe isotopic variations in IFs that are well explained by Fe 2+ oxidation. To address this question, we constrain the equilibrium Fe isotopic ( 56 Fe/ 54 Fe) fractionation of greenalite and ferrihydrite by determining the iron phonon densities of states for those minerals. Here, we use ab initio density functional theory (DFT + U) calculations and nuclear resonant inelastic X-ray scattering spectroscopy to show that ferrous greenalite should be isotopically lighter than ferrihydrite by similar to 1-1.2 parts per thousand at equilibrium, and fractionation should scale linearly with increasing Fe 3+ content in greenalite. By anchoring ferrihydrite-greenalite mineral pair fractionations to published experimental Fe isotopic fractionations between HFOs and aqueous Fe 2+ , we show that ferrous greenalite may produce all but the heaviest pre-GOE Fe isotopic compositions and mixed valence greenalites can produce the entire record. Our results suggest that heavy Fe isotope enrichments alone are not diagnostic of primary IF mineralogies, and ferrihydrite and partially oxidized or even purely ferrous greenalite are all viable primary IF mineralogies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1992

Diagnosing coupled jet-streak circulations for a northern plains snow band from the operational nested-grid model

On 17 March 1989, moderate to heavy snow developed in a 100- to 200-km-wide band extending from South Dakota to northern Michigan. The 4- to 8-inch snowfall within this band was not associated with major cyclogenesis, and developed 500 to 600 km north of a stationary surface front. A diagnostic analysis based on an application of the General Meteorological Package (GEMPAK 5.0) to a numerical simulation from the operational nested-grid model (NGM) is utilized to show that the development of this snow band is related to the interaction of two upper-tropospheric jet streaks and their associated transverse circulation patterns. The eastward propagation of a jet streak from the West Coast toward the middle United States and to the south of a slower-moving jet along the U.S.-Canadian border led to a merger of the ascent maxima associated with the direct and indirect circulations of the northern and southern jets, respectively. The snow band developed as the ascending branches of the jet-streak circulation patterns merged, with the eastward propagation of the heaviest snow linked to the motion of the coupled circulation pattern. The study also demonstrates the usefulness of the operational NGM for providing the higher-resolution datasets required to relate the evolution of jet-streak circulation patterns to the development of mesoscale precipitation bands.

Hakim, Gregory J.↗
DOE OSTI2024

Estimating Uncertainty in Simulated ENSO Statistics

Abstract Large ensembles of model simulations are frequently used to reduce the impact of internal variability when evaluating climate models and assessing climate change induced trends. However, the optimal number of ensemble members required to distinguish model biases and climate change signals from internal variability varies across models and metrics. Here we analyze the mean, variance and skewness of precipitation and sea surface temperature in the eastern equatorial Pacific region often used to describe the El Niño–Southern Oscillation (ENSO), obtained from large ensembles of Coupled model intercomparison project phase 6 climate simulations. Leveraging established statistical theory, we develop and assess equations to estimate, a priori, the ensemble size or simulation length required to limit sampling‐based uncertainties in ENSO statistics to within a desired tolerance. Our results confirm that the uncertainty of these statistics decreases with the square root of the time series length and/or ensemble size. Moreover, we demonstrate that uncertainties of these statistics are generally comparable when computed using either pre‐industrial control or historical runs. This suggests that pre‐industrial runs can sometimes be used to estimate the expected uncertainty of statistics computed from an existing historical member or ensemble, and the number of simulation years (run duration and/or ensemble size) required to adequately characterize the statistic. This advance allows us to use existing simulations (e.g., control runs that are performed during model development) to design ensembles that can sufficiently limit diagnostic uncertainties arising from simulated internal variability. These results may well be applicable to variables and regions beyond ENSO.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2010

The Impact of the Evolving Satellite Data Record on Reanalysis Water and Energy Fluxes During the Past 30 Years

Retrospective analyses (reanalyses) use a fixed assimilation model to take diverse observations and synthesize consistent, time-dependent fields of state variables and fluxes (e.g. temperature, moisture, momentum, turbulent and radiative fluxes). Because they offer data sets of these quantities at regular space / time intervals, atmospheric reanalyses have become a mainstay of the climate community for diagnostic purposes and for driving offline ocean and land models. Of course, one weakness of these data sets is the susceptibility of the flux products to uncertainties because of shortcomings in parameterized model physics. Another issue, perhaps less appreciated, is the fact that the discreet changes in the evolving observational system, particularly from satellite sensors, may also introduce artifacts in the time series of quantities. In this paper we examine the ability of the NASA MERRA (Modern Era Retrospective Analysis for Research and Applications) and other recent reanalyses to determine variability in the climate system over the satellite record (~ the last 30 years). In particular we highlight the effect on reanalyses of discontinuities at the junctures of the onset of passive microwave imaging (Special Sensor Microwave Imager) in late 1987 as well as improved sounding and imaging with the Advanced Microwave Sounding Unit, AMSU-A, in 1998. We examine these data sets from two perspectives. The first is the ability to capture modes of variability that have coherent spatial structure (e.g. ENSO events and near-decadal coupling to SST changes) and how these modes are contained within trends in near global averages of key quantities. Secondly, we consider diagnostics that measure the consistency in energetic scaling in the hydrologic cycle, particularly the fractional changes in column-integrated water vapor versus precipitation as they are coupled to radiative flux constraints. These results will be discussed in the context of implications for science objectives and priorities of the NASA Energy and Water Cycle Study, NEWS.

Robertson, Franklin R.↗
NASA NTRS2022

Lightning Over Central Canada: Skill Assessment for Various Land-Atmosphere Model Configurations and Lightning Indices Over a Boreal Study Area

Current lightning predictions are uncertain because they rely on empirical diagnostic relationships and often use coarse-scale climate scenario simulations in which deep convection is parameterized. Previous studies demonstrated that simulations with convection permitting resolutions improve lightning predictions compared to coarser-grid simulations using convection parameterizations for different geographical locations but not over the boreal zone. In this study, lightning simulations with the NASA Unified-Weather Research and Forecasting (NU-WRF) model are evaluated over a domain including the Great Slave Lake in Canada, for six lightning seasons. The simulations are performed at convection parameterized (9 km) and convection-permitting (3 km) resolution using the Goddard 4ICE and the Thompson microphysics schemes. Four lightning indices are evaluated against observations from the Canadian Lightning Detection Network, in terms of spatiotemporal frequency distribution, spatial pattern, daily climatology, and an event-based overall skill assessment. The Thompson scheme is, regardless of the spatial resolution, superior to the Goddard 4ICE scheme in predicting daily climatology but worse in predicting the spatial patterns of lightning occurrence. Results indicate that lightning estimation benefits from modeling at convection-permitting resolution, in particular for the ice based lightning indices. In contrast, the product of convective available potential energy and precipitation rate proved to be the most robust index that was largely invariant to varying spatial resolution. Finally, this study reveals issues of the models to reproduce the observed spatial pattern of lightning well, which might be related to an insufficient representation of land surface heterogeneity, including peatlands, in the study area.

Jonas Mortelmans↗
DOE OSTI2025

Yield From Iowa's First Commercial Miscanthus Fields: Implications of Spatial Variability for Productivity and Sustainability Beyond Research Plots

The cultivation of sterile giant miscanthus (Miscanthus × giganteus, M × g) for bioenergy and bioproducts has expanded into grain-cropped land in the United States (US) as local markets developed for this high-yielding perennial grass (10–30 Mg DM ha −1 ). However, the magnitude of spatial and temporal variability in yield within US Corn Belt fields, along with impacts on economic return and sustainable land management, is poorly understood. This study established a diagnostic model relating remote sensing-derived vegetation indices to ground truth data from 105 hand-harvested stem biomass samples, which were strategically selected to represent the full range of vegetation index observations. The high-resolution satellite-sensed vegetation indices captured > 90% of the yield variation measured within fields. This model was then used to predict yield variability and assess economic performance across four of the first commercial M × g fields in the Corn Belt state of Iowa, US. Significant spatial variability in biomass dry matter (DM) yields (9.3–18.1 Mg DM ha −1 ) and net profits ($\$$83 to $\$$1211.5 ha −1 ) was observed. All fields were profitable in all site-years. When low profit occurred, it was explained by limited management experience of the crop in Iowa. The breakeven yield at a selling price of $\$$130 Mg −1 varied from 9.0–12.1 Mg ha −1 at 15% moisture content (7.6–10.3 Mg DM ha −1 ). Breakeven prices ranged from $\$$73 to $\$$122.4 Mg −1 , matching ranges used in the Department of Energy Billion Ton Report (US Department of Energy, 2023). Notably, M × g yield and profits were commensurate with grain crops particularly with favorable precipitation. This study provides insight on the M × g management “learning curve”, performance on marginal land and in drought conditions, and demonstrates that addressing yield gaps, reducing costs, and implementing precision agriculture strategies can enhance profitability. These findings emphasize the value of remote sensing technologies in guiding sustainable and competitive commercial-scale M × g production.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Formation and fate of freshwater on an ice floe in the Central Arctic

The melt of snow and sea ice during the Arctic summer is a significant source of relatively fresh meltwater. The fate of this freshwater, whether in surface melt ponds or thin layers underneath the ice and in leads, impacts atmosphere–ice–ocean interactions and their subsequent coupled evolution. Here, we combine analyses of datasets from the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition (June–July 2020) for a process study on the formation and fate of sea ice freshwater on ice floes in the Central Arctic. Our freshwater budget analyses suggest that a relatively high fraction (58 %) is derived from surface melt. Additionally, the contribution from stored precipitation (snowmelt) outweighs by 5 times the input from in situ summer precipitation (rain). The magnitude and rate of local meltwater production are remarkably similar to those observed on the prior Surface Heat Budget of the Arctic Ocean (SHEBA) campaign, where the cumulative summer freshwater production totaled around 1 m during both. A relatively small fraction (10 %) of freshwater from melt remains in ponds, which is higher on more deformed second-year ice (SYI) compared to first-year ice (FYI) later in the summer. Most meltwater drains laterally and vertically, with vertical drainage enabling storage of freshwater internally in the ice by freshening brine channels. In the upper ocean, freshwater can accumulate in transient meltwater layers on the order of 0.1 to 1 m thick in leads and under the ice. The presence of such layers substantially impacts the coupled system by reducing bottom melt and allowing false bottom growth; reducing heat, nutrient, and gas exchange; and influencing ecosystem productivity. Regardless, the majority fraction of freshwater from melt is inferred to be ultimately incorporated into the upper ocean (75 %) or stored internally in the ice (14 %). Terms such as the annual sea ice freshwater production and meltwater storage in ponds could be used in future work as diagnostics for global climate and process models. For example, the range of values from the CESM2 climate model roughly encapsulate the observed total freshwater production, while storage in melt ponds is underestimated by about 50 %, suggesting pond drainage terms as a key process for investigation.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2010

The Effect of Satellite Observing System Changes on MERRA Water and Energy Fluxes

Because reanalysis data sets offer state variables and fluxes at regular space / time intervals, atmospheric reanalyses have become a mainstay of the climate community for diagnostic purposes and for driving offline ocean and land models. Although one weakness of these data sets is the susceptibility of the flux products to uncertainties because of shortcomings in parameterized model physics, another issue, perhaps less appreciated, is the fact that continual but discreet changes in the evolving observational system, particularly from satellite sensors, may also introduce artifacts in the time series of quantities. In this paper we examine the ability of the NASA MERRA (Modern Era Retrospective Analysis for Research and Applications) and other recent reanalyses to determine variability in the climate system over the satellite record (approximately the last 30 years). In particular we highlight the effect on the reanalysis of discontinuities at the junctures of the onset of passive microwave imaging (Special Sensor Microwave Imager) in late 1987 as well as improved sounding and imaging with the Advanced Microwave Sounding Unit, AMSU-A, in 1998. We first examine MERRA fluxes from the perspective of how physical modes of variability (e.g. ENSO events, Pacific Decadal Variability) are contaminated by artificial step-like trends induced by the onset of new moisture data these two satellite observing systems. Secondly, we show how Redundancy Analysis, a statistical regression methodology, is effective in relating these artifact signals in the moisture and temperature analysis increments to their presence in the physical flux terms (e.g. precipitation, radiation). This procedure is shown to be effective greatly reducing the artificial trends in the flux quantities.

Robertson, Franklin R.↗