Search NASA⌕ Search

SEARCH · Search NASA

Results for “SIMULATED ALTITUDE”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Topanga: A kinetic ion plasma code for large-scale ionospheric simulations on magnetohydrodynamic timescales

Topanga is a kinetic ion code developed for simulating large-scale plasma phenomena in the Earth's ionosphere on magnetohydrodynamic timescales. It is a domain-decomposed parallel code that runs on high-performance computing platforms. Features of Topanga include spherical geometry for simplified boundary conditions and computational efficiency; a hybrid plasma model with inertia-less fluid electrons, kinetic ions, and an electric field specified via an Ohm's law; a Maxwell-FDTD (finite difference time domain) plasma model which retains the displacement current in Maxwell's equations and models electron currents in the ionosphere with a tensor conductivity; sponge-layer boundary conditions for absorption of electromagnetic and plasma waves incident on the domain boundaries; and a novel mixed-implicit algorithm for evolving the EM fields inside the Maxwell-FDTD region that is stable over many orders of magnitude in the electron–ion collision frequency. We verify the numerical methods used in Topanga on a pair of test problems. The first test involves modeling a three-dimensional collisionless shock using the hybrid set of equations. The second test involves modeling a spherical TEM mode in vacuum using the Maxwell-FDTD set of equations. Finally, we demonstrate how using the combined set of hybrid and Maxwell-FDTD equations to model the Starfish Prime high-altitude nuclear test recovers a “missing” EM signal on the ground that is not present when using only the hybrid set of equations. The magnitude of this signal in the simulation containing the Maxwell-FDTD region agrees well with the E3a portion of the magnetohydrodynamic electromagnetic pulse from Starfish Prime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Assessing Secondary Ice Production in Continental Clouds Based on AMF Synergistic Remote Sensing Observations

Ice crystals play an important role in radiation and precipitation formation. However, the prediction of ice number concentration has proven to be problematic, because the predicted concentration of ice nucleating particles is often smaller than the observed by several orders of magnitude. This difference between the observed and the predicted number concentration has been explained by so-called secondary ice production (SIP), which includes hypothesized mechanisms like rime splintering, shattering of large frozen drops, and ice-ice collision breakup. To better observe SIP and understand its frequency, production rate, and impact on cloud evolution, we developed a new method for retrieving microphysical properties of concurrent pristine ice and snow aggregates from X-band polarimetric radar observations. Compared to aircraft cloud probe measurements, our retrieved total number concentration and ice water content are overestimated by 98% and 44%, respectively, which outperforms other empirical relationships. The capability of quantifying pristine ice number concentration, which is the most essential observable in SIP studies, opens a new opportunity to build a long-term record of SIP from surface-based radar measurements. We have also conducted simulations using a cloud-resolving model to evaluate the representativeness of the parameterizations of three SIP mechanisms and for studying the impact of SIP on a frontal system. We found a surprisingly good agreement between the observed and retrieved pristine ice number concentrations. The good performance is attributed to the rime splintering process in the temperature zone between–3° and –8°C, and to the ice-ice collision breakup process at higher altitudes. The impact of the shattering of large frozen drops appears minimal. Apart from the enhancement of ice number concentration and increased horizontal coverage of pristine ice, the addition of three SIP processes does not change cloud evolutions dramatically in our simulations.

58 GEOSCIENCES↗
DOE OSTI2021

Upper Troposphere Smoke Injection From Large Areal Fires

Abstract Large areal fires, such as those ignited following a nuclear detonation, can inject smoke into the upper troposphere and lower stratosphere. Detailed fire simulations allow for assessment of how local weather interacts with these fires and affects smoke lofting. In this study, we employ the fire simulation package in the Weather Research and Forecasting (WRF‐Fire) model, Version 4.0.1, to explore how smoke lofting from a fire burning a homogeneous fuel bed changes with varying local winds, relative humidity, and atmospheric boundary‐layer stability for two different‐sized areal fires. The presence of moisture has the greatest influence on the results by raising the altitude of lofting, while faster wind speeds dampen lofting and lower the injection height. Stably stratified conditions inhibit plume propagation compared with neutrally stratified conditions, although the impact of stability is not as strong as that of moisture and winds. These findings highlight the importance of using an appropriate atmospheric profile when simulating large fires, as the local weather can have a meaningful influence on smoke lofting.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

A Vertically Resolved MSE Framework Highlights the Role of the Boundary Layer in Convective Self-Aggregation

Abstract Convective self-aggregation refers to a phenomenon in which random convection can self-organize into large-scale clusters over an ocean surface with uniform temperature in cloud-resolving models. Previous literature studies convective aggregation primarily by analyzing vertically integrated (VI) moist static energy (MSE) variance. That is the global MSE variance, including both the local MSE variance at a given altitude and the covariance of MSE anomalies between different altitudes. Here we present a vertically resolved (VR) MSE framework that focuses on the local MSE variance to study convective self-aggregation. Using a cloud-resolving simulation, we show that the development of self-aggregation is associated with an increase of local MSE variance, and that the diabatic and adiabatic generation of the MSE variance is mainly dominated by the boundary layer (BL; the lowest 2 km). The results agree with recent numerical simulation results and the available potential energy analyses showing that the BL plays a key role in the development of self-aggregation. Additionally, we find that the lower free troposphere (2–4 km) also generates significant MSE variance in the first 15 days. We further present a detailed comparison between the global and local MSE variance frameworks in their mathematical formulation and diagnostic results, highlighting their differences.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Modeling Framework for Bulk Electric Grid Impacts from HEMP E1 and E3 Effects (Tasks 3.1 Final Report)

This report presents a framework to evaluate the impact of a high-altitude electromagnetic pulse (HEMP) event on a bulk electric power grid. This report limits itself to modeling the impact of EMP E1 and E3 components. The co-simulation of E1 and E3 is presented in detail, and the focus of the paper is on the framework rather than actual results. This approach is highly conservative as E1 and E3 are not maximized with the same event characteristics and may only slightly overlap. The actual results shown in this report are based on a synthetic grid with synthetic data and a limited exemplary EMP model. The framework presented can be leveraged and used to analyze the impact of other threat scenarios, both manmade and natural disasters. This report d escribes a Monte-Carlo based methodology to probabilistically quantify the transient response of the power grid to a HEMP event. The approach uses multiple fundamental steps to characterize the system response to HEMP events, focused on the E1 and E3 components of the event. 1) Obtain component failure data related to HEMP events testing of components and creating component failure models. Use the component failure model to create component failure conditional probability density function (PDF) that is a function of the HEMP induced terminal voltage. 2) Model HEMP scenarios and calculate the E1 coupled voltage profiles seen by all system components. Model the same HEMP scenarios and calculate the transformer reactive power consumption profiles due to E3. 3) Sample each component failure PDF to determine which grid components will fail, due to the E1 voltage spike, for each scenario. 4) Perform dynamic simulations that incorporate the predicted component failures from E1 and reactive power consumption at each transformer affected by E3. These simulations allow for secondary transients to affect the relays/protection remaining in service which can lead to cascading outages. 5) Identify the locations and amount of load lost for each scenario through grid dynamic simulation. This can be an indication of the immediate grid impacts from a HEMP event. In addition, perform more detailed analysis to determine critical nodes and system trends. 6) To help realize the longer-term impacts, a security constrained alternating current optimal power flow (ACOPF) is run to maximize critical load served. This report describes a modeling framework to assess the systemic grid impacts due to a HEMP event. This stochastic simulation framework generates a large amount of data for each Monte Carlo replication, including HEMP location and characteristics, relay and component failures, E3 GIC profiles, cascading dynamics including voltage and frequency over time, and final system state. This data can then be analyzed to identify trends, e.g., unique system behavior modes or critical components whose failure is more likely to cause serious systemic effects. The proposed analysis process is demonstrated on a representative system. In order to draw realistic conclusions of the impact of a HEMP event on the grid, a significant amount of work remains with respect to modeling the impact on various grid components.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Presenting a Model to Predict Changing Snow Albedo for Improving Photovoltaic Performance Simulation

As photovoltaic (PV) deployment increases worldwide, PV systems are being installed more frequently in locations that experience snow cover. The higher albedo of snow, relative to the ground, increases the performance of PV systems in northern and high-altitude locations by reflecting more light onto the PV modules. Accurate modeling of the snow’s albedo can improve estimates of PV system production. Typical modeling of snow albedo uses a simple two-value model that sets the albedo high when snow is present, and low when snow is not present. However, snow albedo changes over time as snow settles and melts and a binary model does not account for transitional changes, which can be significant. Here, we present and validate a model for estimating snow albedo as it changes over time. The model is simple enough to only require daily snow depth and hourly average temperature data, but can be improved through the addition of site-specific factors, when available. We validate this model to quantify its ability to more accurately predict snow albedo and compare the model’s performance against satellite imagery-based methods for obtaining historical albedo data. In addition, we perform modeling using the System Advisor Model (SAM) to show the impact of changes in albedo on energy modeling for PV systems. Overall, our albedo model has a significantly improved ability to predict the solar insolation on PV modules in real time, especially on bifacial PV modules where reflected irradiance plays a larger role in energy production.

Pike, Christopher (ORCID:0000000155888033)↗
DOE OSTI2020

Neutrinos from the Galactic Center Hosting a Hypernova Remnant

Similar to star-forming galaxies or starburst galaxies, star-forming regions in our Galaxy can host cosmic-ray (CR) accelerators and rich gas as targets of hadronuclear interaction. By our estimations, the IceCube neutrino observatory might detect muon neutrinos from a CR accelerator associated with a molecular cloud complex in our Galaxy. The associated high-energy gamma-ray emission might be observed by the Cherenkov Telescope Array (CTA), High-Altitude Water Cherenkov Gamma-Ray Observatory (HAWC), and Large High Altitude Air Shower Observatory (LHAASO). Furthermore, taking the Galactic Center (GC) region as an example, we assume that a hypernova exploded in the past in the GC. We simulate the acceleration of CRs in the hypernova remnant (HNR) as well as their confinement and escape. The high-energy protons escape from the HNR, diffuse around the GC, interact with molecular clouds, and then produce gamma-rays and neutrinos. In the optimal cases, the GC would be a promising 100 TeV gamma-ray source for LHAASO’s one-month observation. We propose that neutrino-induced searching for starting track-like and high-energy starting events (HESEs) observed by IceCube, from the GC region with a radius of 1°.8, would help us discover the particle accelerator in the GC or constrain our models. Under the constraint from high-energy gamma-ray observations by the H.E.S.S. telescope, we estimate the exposure time needed to make a significant discovery for the optimal cases. The analysis combining observations of IceCube and ANTARES, starting track-like events and HESEs, future observations by neutrino detectors IceCube-Gen2 and KM3net, and gamma-ray telescopes CTA, HAWC, and LHAASO would help to constrain our models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Aerosol-deep convection interaction based on joint cell-thermal tracking in Large Eddy Simulations during the TRACER campaign

In cumulus clouds, aerosol concentrations control cloud droplet concentrations, modifying cloud radiative properties, precipitation processes, and cloud electrification. However, mechanisms of aerosol-deep convection interactions are not well understood due to complex cloud dynamics and microphysics. We investigate the interaction of aerosols with isolated deep convection using Large Eddy Simulations of two cases during the TRacking Aerosol Convection interactions ExpeRiment (TRACER) near Houston, Texas, using a joint cell-thermal tracking algorithm. Cumulus thermals are droplet generators, since supersaturation and droplet nucleation coincide with thermal centers, where the strongest updrafts occur. Primary ice crystal formation does not take place inside thermals, but at layers where previous thermals detrained moisture. As subsequent thermals containing supercooled droplets penetrate these layers, hail and graupel form at or near these thermals. Higher aerosol concentrations result in higher droplet concentrations that suppress drizzle, delay warm rain processes, and transport more moisture aloft. This increases snow and ice amount, as well as graupel and hail, leading to more lightning. Polluted thermals initiate at slightly higher altitudes, and are slightly larger and faster, suggesting a weak invigoration. We also find more thermals per cell, but fewer isolated cells, since convection is more aggregated and intense, especially near the end of the 24 h simulation. Non-linear mesoscale feedback likely triggered by temperature and moisture responses to aerosol-thermal interactions causes the aggregation. Time-lagged aerosol-reinitialization experiments show that the mesoscale response is the predominant forcing for the invigoration. These changes happen within one day, on a smaller scale than previously suggested.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Simulation-aided characterization of a versatile water-based condensation particle counter for atmospheric airborne research

Capturing the vertical profiles and horizontal variations of atmospheric aerosols often requires accurate airborne measurements. With the advantage of avoiding health and safety concerns related to the use of butanol or other chemicals, water-based condensation particle counters have emerged to provide measurements under various environments. However, airborne deployments are relatively rare due to the lack of instrument characterization under reduced pressure at flight altitudes. This study investigates the performance of a commercial “versatile” water-based condensation particle counter (vWCPC, model 3789, TSI, Shoreview, MN, USA) under various ambient pressure conditions (500–920 hPa) with a wide range of particle total number concentrations (1500–70 000 cm -3 ). The effect of conditioner temperature on vWCPC 3789 performance at low pressure is examined through numerical simulation and laboratory experiments. We show that the default instrument temperature setting of 30 °C for the conditioner is not suitable for airborne measurement and that the optimal conditioner temperature for low-pressure operation is 27°. Under the optimal conditioner temperature (27°), the 7 nm cut-off size is also maintained. Additionally, we show that insufficient droplet growth becomes more significant under the low-pressure operation. The counting efficiency of the vWCPC 3789 can vary up to 20 % for particles of different chemical compositions (e.g., ammonium sulfate and sucrose particles). However, such variation is independent of pressure.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Impacts of Aerosol and Environmental Conditions on Maritime and Continental Deep Convective Systems Using a Bin Microphysical Model

A series of model simulations were conducted to investigate the effects of cloud condensation nuclei (CCN) loading and convective available potential energy (CAPE) on tropical maritime and midlatitude continental deep convection. Dynamical downscaling from global aerosol reanalysis was used to represent aerosol fields for the two convective regimes. We describe a control run and multiple sensitivity experiments using a limited-area model, employing spectral-bin cloud microphysics. The CCN loading is perturbed between the target maritime and continental conditions, roughly 40–2,000 cm -3 at 850 hPa and 1% supersaturation. Surface precipitation rates monotonically increase with increasing CCN loading for both the maritime and continental situations, while these monotonic increases are disrupted in the simulations with reduced CAPE. The increase in precipitation is in the form of convective precipitation, at the expense of stratiform precipitation. CCN increases promote increases in supercooled cloud water, in agreement with previous modeling studies. However, in the simulations investigated herein, the changes in supercooled water have different impacts on the cloud microphysics in the maritime and continental simulations. Increased supercooled water contents lead to more hail and less graupel in the continental simulation. For the maritime simulation, enhanced supercooled cloud water contents promote an increase in graupel since little or no hail is produced. This distinction is due to the difference in relative magnitudes and peak altitudes of supercooled water and snow amounts, which is further attributable to moisture and dynamical differences in the two cases.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Temperature Trends in the Northwestern Tibetan Plateau Constrained by Ice Core Water Isotopes Over the Past 7,000 Years

Abstract The reasons for the Holocene temperature conundrum, known as the inconsistency between the reconstructed cooling and the inferred warming simulations during the Holocene, remain unclear. Temperature reconstructions from the Tibetan Plateau (TP) provide important insights for understanding the Holocene temperature conundrum due to enhanced sensitivity to climate at high altitudes. Given the significant positive correlation between air temperature and δ 18 O in precipitation over the northern TP, the stable isotopic records of ice cores recovered from this area are widely used for paleotemperature reconstruction. Here we present a new high‐resolution δ 18 O record from the Chongce ice cores to bedrock, dated back to 7 ka BP by the accelerator mass spectrometry (AMS) 14 C dating technique. Our reconstructed temperature record shows a long‐term warming trend until ~2 ka BP, followed by an abrupt change to a relatively cool period until the start of the industrial‐era warming. This record challenges the widely recognized Holocene reconstruction from the neighboring Guliya ice core. It is also different from many previous temperature reconstructions, most of which have summer biases and show a long‐term cooling trend over the past two millennia. In addition, our record shows that temperatures during the recent decades are almost the highest during the past 7 ka BP, highlighting the unusual warming forced by anthropogenic greenhouse gases.

Pang, Hongxi↗
DOE OSTI2023

MODFLOW6 models used to evaluate potential stresses and hydrologic conditions driving water-level fluctuations in well ER-5-3-2, Frenchman Flat, Southern Nevada

The hydrograph for well ER-5-3-2 in Frenchman Flat, southern Nevada, has previously unexplained water-level fluctuations. Four, three-dimensional, groundwater models (MODFLOW 6) were developed to evaluate potential stresses and hydrologic conditions affecting the well ER-5-3-2 hydrograph. Four model scenarios were developed that simulated: (1) wellbore leakage without recharge, (2) wellbore leakage with recharge, (3) shallow (low transmissivity) and deep (high transmissivity) carbonate rocks, and (4) lateral heterogeneity of carbonate rocks. Input and output files for the four model scenarios are in the model and output directories, respectively. Hydraulic conductivity, specific storage, and wellbore-leakage rates (when simulated) were estimated with parameter estimation (PEST) by minimizing a weighted composite, sum-of-squares objective function. The objective function was informed by measurement and Tikhonov regularization observations. Measurement observations included drawdowns from the constant-rate aquifer test and water-level altitudes measured in well ER-5-3-2 from 2001-2021. Tikhonov regularization informed hydraulic conductivity and specific storage parameters that were insensitive to measurement observations, where homogeneity was the preferred relation. Batch files, executables, and MODFLOW 6, PEST, and post-processing utilities are in the ancillary directory. Supplementary data also are included in the ancillary directory, including site information, high-frequency water-level and aquifer-test data, transmissivity estimates, water-chemistry data, and water-temperature analyses. This USGS data release contains data, analyses, and model files for the simulations and analysis results described in U.S. Geological Survey Scientific Investigations Report (https://doi.org/10.3133/sir20225132).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Evaluating Cloud Feedback Components in Observations and Their Representation in Climate Models

This study quantifies the contribution of individual cloud feedbacks to the total short-term cloud feedback in satellite observations over the period 2002–2014 and evaluates how they are represented in climate models. The observed positive total cloud feedback is primarily due to positive high-cloud altitude, extratropical high- and low-cloud optical depth, and land cloud amount feedbacks partially offset by negative tropical marine low-cloud feedback. Seventeen models from the Atmosphere Model Intercomparison Project of the sixth Coupled Model Intercomparison Project are analyzed. The models generally reproduce the observed moderate positive short-term cloud feedback. However, compared to satellite estimates, the models are systematically high-biased in tropical marine low-cloud and land cloud amount feedbacks and systematically low-biased in high-cloud altitude and extratropical high- and low-cloud optical depth feedbacks. Errors in modeled short-term cloud feedback components identified in this analysis highlight the need for improvements in model simulations of the response of high clouds and tropical marine low clouds. Our results suggest that skill in simulating interannual cloud feedback components may not indicate skill in simulating long-term cloud feedback components.

58 GEOSCIENCES↗
DOE OSTI2023

Assessing the vulnerability of solar inverters to EMPs: Port testing, PCI modeling, and protection strategies

Renewable energy sources are becoming an ever-larger contributor to the power grid. These renewable energy sources depend upon the power electronic devices, specifically inverters, being essential for connecting Photovoltaic (PV) generation to the grid. However, the Electromagnetic Pulses (EMPs) caused by the high-altitude nuclear explosions can generate fast broad-band pulses with nanosecond rise time, potentially causing damage or destruction to electronic components. To assess the vulnerability of PV inverters to high-altitude EMPs, the port testing and Pulsed Current Injection (PCI) modeling schemes are proposed based on the port impedance analysis. Wide-band frequency measurements are achieved by fusing impedance results from three vector network analyzers. Then, a PCI model is used to simulate the induced response to EMP, with two typical immunity levels of EC5 and EC8 tested. Here, the experiment successfully excites the induced voltage and current under EMP, where the voltage and current can reach 1500V/40A and 8000V/150Aunder EC5 andEC8, respectively. The port vulnerability analysis results demonstrate that only some ports can survive under EC5. To defend against the impact of EMP, three protection strategies are discussed.

42 ENGINEERING↗
DOE OSTI2020

Direct and semi-direct radiative forcing of biomass-burning aerosols over the southeast Atlantic (SEA) and its sensitivity to absorbing properties: a regional climate modeling study

Simulations are performed for the period 2000-2015 by two different regional climate models, ALADIN and RegCM, to quantify the direct and semi-direct radiative effects of biomass-burning aerosols (BBAs) in the southeast Atlantic (SEA) region. Different simulations have been performed using strongly absorbing BBAs in accordance with recent in situ observations over the SEA. For the July-August-September (JAS) season, the single scattering albedo (SSA) and total aerosol optical depth (AOD) simulated by the ALADIN and RegCM models are consistent with the MACv2 climatology and MERRA-2 and CAMS-RA reanalyses near the biomass-burning emission sources. However, the above-cloud AOD is slightly underestimated compared to satellite (MODIS and POLDER) data during the transport over the SEA. The direct radiative effect exerted at the continental and oceanic surfaces by BBAs is significant in both models and the radiative effects at the top of the atmosphere indicate a remarkable regional contrast over SEA (in all-sky conditions), with a cooling (warming) north (south) of 10°S, which is in agreement with the recent MACv2 climatology. In addition, the two models indicate that BBAs are responsible for an important shortwave radiative heating of ~0.5-1K per day over SEA during JAS with maxima between 2 and 4 km a.m.s.l. (above mean sea level). At these altitudes, BBAs increase air temperature by ~0.2-0.5 K, with the highest values being co-located with low stratocumulus clouds. Vertical changes in air temperature limit the subsidence of air mass over SEA, creating a cyclonic anomaly. The opposite effect is simulated over the continent due to the increase in lower troposphere stability. The BBA semi-direct effect on the lower troposphere circulation is found to be consistent between the two models. Changes in the cloud fraction are moderate in response to the presence of smoke, and the models differ over the Gulf of Guinea. Finally, the results indicate an important sensitivity of the direct and semi-direct effects to the absorbing properties of BBAs. Over the stratocumulus (Sc) region, DRE varies from +0.94 W m -2 (scattering BBAs) to +3.93 W m -2 (most absorbing BBAs).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Methods to Evaluate Subcolumn Profiles Based on Two-Point Diagnostics

In atmospheric models, stochastic generation of subgrid-scale profiles or “subcolumns” has been used for a variety of purposes. Such subcolumns can be generated from subgrid probability density functions (PDFs) at different vertical levels, when such PDFs are available. To do so, the generator needs to decide how strongly points should be correlated in the vertical, that is, how much the values should be overlapped. This is sometimes called “PDF overlap.” To assess vertical correlation in a simplified, observable setting, here the vertical correlation of vertical velocity in subcloud layers is examined. Doppler lidar is used to evaluate the vertical profiles of vertical velocity produced by a large-eddy simulation (LES) model and the Subgrid Importance Latin Hypercube Sampler (SILHS) subcolumn generator. In order to diagnose unrealistic features in subcolumn profiles, various statistical diagnostics are examined here, including the bivariate PDF of vertical velocity at two separated points (i.e., altitudes), the two-point velocity correlation, the integral correlation length, the PDF of two-point velocity differences, and the skewness and kurtosis of two-point velocity differences. The profiles produced by LES match lidar well, except that they are too smooth at small scales. The profiles produced by SILHS exhibit sharp jumps from updraft to downdraft that are not observed in the lidar data. To reduce the generation of these unrealistically sharp jumps, the SILHS sampling method is revised. The diagnostics confirm that the revised sampling method reduces the overprediction of sharp jumps.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Effects of spatial resolution on WRF v3.8.1 simulated meteorology over the central Himalaya

The sensitive ecosystem of the central Himalayan (CH) region, which is experiencing enhanced stress from anthropogenic forcing, requires adequate atmospheric observations and an improved representation of the Himalaya in the models. However, the accuracy of atmospheric models remains limited in this region due to highly complex mountainous topography. This article delineates the effects of spatial resolution on the modeled meteorology and dynamics over the CH by utilizing the Weather Research and Forecasting (WRF) model extensively evaluated against the Ganges Valley Aerosol Experiment (GVAX) observations during the summer monsoon. The WRF simulation is performed over a domain (d01) encompassing northern India at 15 km x 15 km resolution and two nests (d02 at 5 km x 5 km and d03 at 1 km x 1 km) centered over the CH, with boundary conditions from the respective parent domains. WRF simulations reveal higher variability in meteorology, e.g., relative humidity (RH = 70.3%-96.1%) and wind speed (WS = 1.1-4.2m s(-1)), compared to the ERA-Interim reanalysis (RH D 80.0%-85.0%, WS D 1.2-2.3m s(-1)) over northern India owing to the higher resolution. WRF-simulated temporal evolution of meteorological variables is found to agree with balloon-borne measurements, with stronger correlations aloft (r = 0.44-0.92) than those in the lower troposphere (r = 0.18-0.48). The model overestimates temperature (warm bias by 2.8 degrees C) and underestimates RH (dry bias by 6.4 %) at the surface in d01. Model results show a significant improvement in d03 (P = 827.6 hPa, T = 19.8 degrees C, RH = 92.3 %), closer to the GVAX observations (P = 801.4 hPa, T = 19.5 degrees C, RH = 94.7 %). Interpolating the output from the coarser domains (d01, d02) to the altitude of the station reduces the biases in pressure and temperature; however, it suppresses the diurnal variations, highlighting the importance of well-resolved terrain (d03). Temporal variations in near-surface P, T, and RH are also reproduced by WRF in d03 to an extent (r>0:5). A sensitivity simulation incorporating the feedback from the nested domain demonstrates the improvement in simulated P, T, and RH over the CH. Our study shows that the WRF model setup at finer spatial resolution can significantly reduce the biases in simulated meteorology, and such an improved representation of the CH can be adopted through domain feedback into regional-scale simulations. Interestingly, WRF simulates a dominant easterly wind component at 1 km x 1 km resolution (d03), which is missing in the coarse simulations; however, the frequency of southeasterlies remains underestimated. The model simulation implementing a highresolution (3 s) topography input (SRTM) improved the prediction of wind directions; nevertheless, further improvements are required to better reproduce the observed local-scale dynamics over the CH.

Singh, Jaydeep↗
DOE OSTI2020

A New Method for Ice-Ice Aggregation in the Adaptive Habit Model

In this paper, a novel methodology for modeling ice-ice aggregation is presented. This methodology combines a modified hydrodynamic collection algorithm with bulk aggregate characteristic information from an offline simulator that collects ice particles, namely, the Ice Particle and Aggregate Simulator, and has been implemented into the Adaptive Habit Microphysics scheme in the Weather Research and Forecasting Model. Aggregates, or snow, are formed via collection of cloud ice particles, where initial ice characteristics and the resulting geometry determine aggregate characteristics. Upon implementation, idealized squall-line simulations are performed to examine the new methodology in comparison with commonly used bulk microphysics schemes. It is found that the adaptive habit aggregation parameterization develops snow and reduces ice mass and number concentrations compared to other schemes. The development of aggregates through the new methodology cascades into other interesting effects, including enhancements in ice and snow growth, as well as homogeneous freezing. Further microphysical analyses reveal varying sensitivities, where snow processes are most sensitive to the new parameterization, followed by ice, then cloud, rain, and graupel processes. Further, the new scheme results in enhancements in surface precipitation due to the persistence of snow at lower altitudes. This persistence is a result of shape-dependent melting and sublimation, increasing the residence time. Moreover, these low-level enhancements are reflected in increases in radar reflectivity at the surface and its spatial distribution. Finally, the ability to predict snow shape and density allows for the simulation of polarimetric radar quantities, resulting in signature enhancements compared to schemes that do not consider spatial and temporal variations in snow shape and density.

54 ENVIRONMENTAL SCIENCES↗