Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “modeling framework”

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 163 records · Page 9

Model validation - A connection between robust control and identification

The gap between the models used in control synthesis and those obtained from identification experiments is considered by investigating the connection between uncertain models and data. The model validation problem addressed is: given experimental data and a model with both additive noise and norm-bounded perturbations, is it possible that the model could produce the observed input-output data? This problem is studied for the standard H-infinity/mu framework models. A necessary condition for such a model to describe an experimental datum is obtained. For a large class of models in the robust control framework, this condition is computable as the solution of a quadratic optimization problem.

Smith, Roy S.↗

Applications and Improvement of a Coupled, Global and Cloud-Resolving Modeling System

Recently Grabowski (2001) and Khairoutdinov and Randall (2001) have proposed the use of 2D CFWs as a "super parameterization" [or multi-scale modeling framework (MMF)] to represent cloud processes within atmospheric general circulation models (GCMs). In the MMF, a fine-resolution 2D CRM takes the place of the single-column parameterization used in conventional GCMs. A prototype Goddard MMF based on the 2D Goddard Cumulus Ensemble (GCE) model and the Goddard finite volume general circulation model (fvGCM) is now being developed. The prototype includes the fvGCM run at 2.50 x 20 horizontal resolution with 32 vertical layers from the surface to 1 mb and the 2D (x-z) GCE using 64 horizontal and 32 vertical grid points with 4 km horizontal resolution and a cyclic lateral boundary. The time step for the 2D GCE would be 15 seconds, and the fvGCM-GCE coupling frequency would be 30 minutes (i.e. the fvGCM physical time step). We have successfully developed an fvGCM-GCE coupler for this prototype. Because the vertical coordinate of the fvGCM (a terrain-following floating Lagrangian coordinate) is different from that of the GCE (a z coordinate), vertical interpolations between the two coordinates are needed in the coupler. In interpolating fields from the GCE to fvGCM, we use an existing fvGCM finite- volume piecewise parabolic mapping (PPM) algorithm, which conserves the mass, momentum, and total energy. A new finite-volume PPM algorithm, which conserves the mass, momentum and moist static energy in the z coordinate, is being developed for interpolating fields from the fvGCM to the GCE. In the meeting, we will discuss the major differences between the two MMFs (i.e., the CSU MMF and the Goddard MMF). We will also present performance and critical issues related to the MMFs. In addition, we will present multi-dimensional cloud datasets (i.e., a cloud data library) generated by the Goddard MMF that will be provided to the global modeling community to help improve the representation and performance of moist processes in climate models and to improve our understanding of cloud processes globally (the software tools needed to produce cloud statistics and to identify various types of clouds and cloud systems from both high-resolution satellite and model data will be also presented).

Tao, W.-K.↗

The Impact of the Assimilation of AIRS Radiance Measurements on Short-term Weather Forecasts

Advanced spaceborne instruments have the ability to improve the horizontal and vertical characterization of temperature and water vapor in the atmosphere through the explicit use of hyperspectral thermal infrared radiance measurements. The incorporation of these measurements into a data assimilation system provides a means to continuously characterize a three-dimensional, instantaneous atmospheric state necessary for the time integration of numerical weather forecasts. Measurements from the National Aeronautics and Space Administration (NASA) Atmospheric Infrared Sounder (AIRS) are incorporated into the gridpoint statistical interpolation (GSI) three-dimensional variational (3D-Var) assimilation system to provide improved initial conditions for use in a mesoscale modeling framework mimicking that of the operational North American Mesoscale (NAM) model. The methodologies for the incorporation of the measurements into the system are presented. Though the measurements have been shown to have a positive impact in global modeling systems, the measurements are further constrained in this system as the model top is physically lower than the global systems and there is no ozone characterization in the background state. For a study period, the measurements are shown to have positive impact on both the analysis state as well as subsequently spawned short-term (0-48 hr) forecasts, particularly in forecasted geopotential height and precipitation fields. At 48 hr, height anomaly correlations showed an improvement in forecast skill of 2.3 hours relative to a system without the AIRS measurements. Similarly, the equitable threat and bias scores of precipitation forecasts of 25 mm (6 hr)-1 were shown to be improved by 8% and 7%, respectively.

McCarty, Will↗

The Global Methane Budget 2000–2017

Understanding and quantifying the global methane (CH4) budget is important for assessing realistic pathways to mitigate climate change. Atmospheric emissions and concentrations of CH4 continue to increase, making CH4 the second most important human-influenced greenhouse gas in terms of climate forcing, after carbon dioxide (CO2). The relative importance of CH4 compared to CO2 depends on its shorter atmospheric lifetime, stronger warming potential, and variations in atmospheric growth rate over the past decade, the causes of which are still debated. Two major challenges in reducing uncertainties in the atmospheric growth rate arise from the variety of geographically overlapping CH4 sources and from the destruction of CH4 by short-lived hydroxyl radicals (OH). To address these challenges, we have established a consortium of multidisciplinary scientists under the umbrella of the Global Carbon Project to synthesize and stimulate new research aimed at improving and regularly updating the global methane budget. Following Saunois et al. (2016), we present here the second version of the living review paper dedicated to the decadal methane budget, integrating results of top-down studies (atmospheric observations within an atmospheric inverse-modelling framework) and bottom-up estimates (including process-based models for estimating land surface emissions and atmospheric chemistry, inventories of anthropogenic emissions, and data-driven extrapolations). For the 2008–2017 decade, global methane emissions are estimated by atmospheric inversions (a top-down approach) to be 576 Tg CH4/yr (range 550–594, corresponding to the minimum and maximum estimates of the model ensemble). Of this total, 359 Tg CH4/yr or ∼ 60 % is attributed to anthropogenic sources, that is emissions caused by direct human activity (i.e. anthropogenic emissions; range 336–376 Tg CH4/yr or 50 %–65 %). The mean annual total emission for the new decade (2008–2017) is 29 Tg CH4/yr larger than our estimate for the previous decade (2000–2009), and 24 Tg CH4/yr larger than the one reported in the previous budget for 2003–2012 (Saunois et al., 2016). Since 2012, global CH4 emissions have been tracking the warmest scenarios assessed by the Intergovernmental Panel on Climate Change. Bottom-up methods suggest almost 30 % larger global emissions (737 Tg CH4/yr, range 594–881) than top-down inversion methods. Indeed, bottom-up estimates for natural sources such as natural wetlands, other inland water systems, and geological sources are higher than top-down estimates. The atmospheric constraints on the top-down budget suggest that at least some of these bottom-up emissions are overestimated. The latitudinal distribution of atmospheric observation-based emissions indicates a predominance of tropical emissions (∼ 65 % of the global budget, < 30° N) compared to mid-latitudes (∼ 30 %, 30–60° N) and high northern latitudes (∼ 4 %, 60–90° N). The most important source of uncertainty in the methane budget is attributable to natural emissions, especially those from wetlands and other inland waters. Some of our global source estimates are smaller than those in previously published budgets (Saunois et al., 2016; Kirschke et al., 2013). In particular wetland emissions are about 35 Tg CH4/yr lower due to improved partition wetlands and other inland waters. Emissions from geological sources and wild animals are also found to be smaller by 7 Tg CH4/yr by 8 Tg CH4/yr, respectively. However, the overall discrepancy between bottom-up and top-down estimates has been reduced by only 5 % compared to Saunois et al. (2016), due to a higher estimate of emissions from inland waters, highlighting the need for more detailed research on emissions factors. Priorities for improving the methane budget include (i) a global, high-resolution map of water-saturated soils and inundated areas emitting methane based on a robust classification of different types of emitting habitats; (ii) further development of process-based models for inland-water emissions; (iii) intensification of methane observations at local scales (e.g., FLUXNET-CH4 measurements) and urban-scale monitoring to constrain bottom-up land surface models, and at regional scales (surface networks and satellites) to constrain atmospheric inversions; (iv) improvements of transport models and the representation of photochemical sinks in top-down inversions; and (v) development of a 3D variational inversion system using isotopic and/or co-emitted species such as ethane to improve source partitioning.

methane budget↗

Implementation of Charged Particle Behavior in Discrete Element Method (DEM) Simulations

Lunar landers will agitate the surface of the Moon with an exhaust plume during descent which will, due to the particulate nature of the lunar regolith, loosen and eject grains from the surface. This ejection is not only coupled with the charged plume gas, but also results in significant particle-particle interactions. Settling of these grains after plume effects have subsided takes much longer than expected in a ballistic trajectory. The prevailing hypothesis attributes this behavior to the accumulated charge on the particles. We are thus developing a discrete element method (DEM) approach to explore these charged particle interactions on the lunar surface. The Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) Improved for General Granular and Granular Heat Transfer Simulations (LIGGGHTS) software package provides a DEM modeling framework for granular interactions. It includes many complexities such as non-spherical particle shapes, cohesion and frictional forces, and heat transfer, but has no provision for inter-particle electrostatic forces and charge transfer that are important to examine in the lunar environment. In this work, a standard Coulomb potential and a Yukawa potential are integrated into the LIGGGHTS framework to provide a basis for particle-particle electrostatic interactions, as well as a gravitational potential to enable inter-grain gravitational interactions. A preliminary approach to charge transfer between grains incorporating properties such as work function and electrical conductivity to the library of available material characteristics will be presented. Several scenarios have been simulated that include charged particle interactions within a diffuse granular gas, settling of charged grains into a regolith bed, sliding of granular material along an incline, and vibration of settled grains to produce a behavior similar to fluidization. There are numerous challenges to incorporate realistic interactions between complex lunar particles. Currently, grains are modeled to behave as if the entirety of the charge acts at the center of mass, such as conductors with spherical symmetry and insulators with homogeneously distributed charge. We are developing improvements that will include the use of non-spherical particle geometries, as well as reasonable approximations of insulating/dielectric materials that have non-uniform charge distributions. The cases simulated thus far will be examined in a relevant environment within a vacuum chamber to validate the simulations. These simulations will be bounded by experiments utilizing high-speed camera observations of the motion for validation. The grains in the experiment will exchange charge during their motion and this can be quantified by collection within a charge measurement device such as a Faraday cup. Such a device may be modeled within the software by defining an integration region and computing the contained charge as a function of simulation time, allowing for side-by-side comparison of simulated and measured bulk charging results. Any differences will be reconciled by updating the mathematical mechanisms described within the simulation suite. Successfully combining results from experiments within a relevant environment into the LIGGGHTS framework will improve modeling of the charged grain dynamics experienced on the Moon to provide insights into dust behavior for future lunar exploration missions.

Electrostatics↗

Computational Aspects of Data Assimilation and the ESMF

The scientific challenge of developing advanced data assimilation applications is a daunting task. Independently developed components may have incompatible interfaces or may be written in different computer languages. The high-performance computer (HPC) platforms required by numerically intensive Earth system applications are complex, varied, rapidly evolving and multi-part systems themselves. Since the market for high-end platforms is relatively small, there is little robust middleware available to buffer the modeler from the difficulties of HPC programming. To complicate matters further, the collaborations required to develop large Earth system applications often span initiatives, institutions and agencies, involve geoscience, software engineering, and computer science communities, and cross national borders.The Earth System Modeling Framework (ESMF) project is a concerted response to these challenges. Its goal is to increase software reuse, interoperability, ease of use and performance in Earth system models through the use of a common software framework, developed in an open manner by leaders in the modeling community. The ESMF addresses the technical and to some extent the cultural - aspects of Earth system modeling, laying the groundwork for addressing the more difficult scientific aspects, such as the physical compatibility of components, in the future. In this talk we will discuss the general philosophy and architecture of the ESMF, focussing on those capabilities useful for developing advanced data assimilation applications.

daSilva, A.↗

Comparison of chemical lateral boundary conditions for air quality predictions over the contiguous United States during pollutant intrusion events

The National Air Quality Forecast Capability (NAQFC) operated in U.S.’s National Oceanic and Atmospheric Administration (NOAA) provides the operational forecast guidance for ozone and fine particulate matters with aerodynamic diameters less than 2.5μm (PM2.5) over the contiguous 48 U.S. states (CONUS) using the Community Multi-scale Air Quality (CMAQ) model. The existing NAQFC uses climatological chemical lateral boundary conditions ( CLBCs), which cannot capture pollutant intrusion events originating outside of the model domain. In this study, we developed a model framework to use dynamic CLBCs from the Goddard Earth Observing System Model, version 5 (GEOS) to drive NAQFC. A mapping of the GEOS chemical species to the CMAQ’s Carbon Bond 5 (CB05)-Aero6 species was developed. The utilization of the GEOS dynamic CLBCs in NAQFC showed the best overall performance in simulating the surface observations during the Saharan dust intrusion and Canadian wildfire events in summer 2015.The simulated PM2.5 was improved from 0.18 to 0.37 and the mean bias was reduced from -6.74 μg/m3 to -2.96 μg/m3 over CONUS. Although the effect of CLBCs on the PM2.5 correlation was mainly near the inflow boundary, its impact on the background concentrations reached further inside the domain. The CLBCs could affect background ozone concentrations through the inflows of ozone itself and its precursors, such as CO. It was further found that the aerosol optical thickness (AOT) from satellite retrievals correlated well with the column CO and elemental carbon from GEOS. The satellite-derived AOT CLBCs generally improved the model performance for the wildfire intrusion events during a summer 2018 case study, and demonstrated how satellite observations of atmospheric composition could be used as an alternative method to capture the air quality effects of intrusions when the global model CLBCs, such as GEOS CLBCs, are not available.

Community Multi-scale Air Quality (CMAQ) model↗

Rapid Rebuilding of the Outer Radiation Belt

Recent observations by the radiation monitor (RDM) on the spacecraft Akebono have shown several cases of greater than 2.5 MeV radiation belt electron enhancements occurring on timescales of less than a few hours. Similar enhancements are also seen in detectors on board the NOAA/POES and TWINS 1 satellites. These intervals are shorter than typical radial diffusion or wave-particle interactions can account for. We choose two so-called "rapid rebuilding" events that occur during high speed streams (4 September 2008 and 22 July 2009) and simulated them with the Space Weather Modeling Framework configured with global magnetosphere, radiation belt, ring current, and ionosphere electrodynamics model. Our simulations produce a weaker and delayed dipolarization as compared to observations, but the associated inductive electric field in the simulations is still strong enough to rapidly transport and accelerate MeV electrons resulting in an energetic electron flux enhancement that is somewhat weaker than is observed. Nevertheless, the calculated flux enhancement and dipolarization is found to be qualitatively consistent with the observations. Taken together, the modeling results and observations support the conclusion that storm-time dipolarization events in the magnetospheric magnetic field result in strong radial transport and energization of radiation belt electrons.

Glocer, A.↗

A Storm Surge and Inundation Model of the Back River Watershed at NASA Langley Research Center

This report on a Virginia Institute for Marine Science project demonstrates that the sub-grid modeling technology (now as part of Chesapeake Bay Inundation Prediction System, CIPS) can incorporate high-resolution Lidar measurements provided by NASA Langley Research Center into the sub-grid model framework to resolve detailed topographic features for use as a hydrological transport model for run-off simulations within NASA Langley and Langley Air Force Base. The rainfall over land accumulates in the ditches/channels resolved via the model sub-grid was tested to simulate the run-off induced by heavy precipitation. Possessing both the capabilities for storm surge and run-off simulations, the CIPS model was then applied to simulate real storm events starting with Hurricane Isabel in 2003. It will be shown that the model can generate highly accurate on-land inundation maps as demonstrated by excellent comparison of the Langley tidal gauge time series data (CAPABLE.larc.nasa.gov) and spatial patterns of real storm wrack line measurements with the model results simulated during Hurricanes Isabel (2003), Irene (2011), and a 2009 Nor'easter. With confidence built upon the model's performance, sea level rise scenarios from the ICCP (International Climate Change Partnership) were also included in the model scenario runs to simulate future inundation cases.

Loftis, Jon Derek↗

A Source-Term Based Boundary Layer Bleed/Effusion Model for Passive Shock Control

A modeling framework for boundary layer effusion has been developed based on the use of source (or sink) terms instead of the usual practice of specifying bleed directly as a boundary condition. This framework allows the surface boundary condition (i.e. isothermal wall, adiabatic wall, slip wall, etc.) to remain unaltered in the presence of bleed. This approach also lends itself to easily permit the addition of empirical models for second order effects that are not easily accounted for by simply defining effective transpiration values. Two effusion models formulated for supersonic flows have been implemented into this framework; the Doerffer/Bohning law and the Slater formulation. These models were applied to unit problems that contain key aspects of the flow physics applicable to bleed systems designed for hypersonic air-breathing propulsion systems. The ability of each model to predict bulk bleed properties was assessed, as well as the response of the boundary layer as it passes through and downstream of a porous bleed system. The model assessment was performed with and without the presence of shock waves. Three-dimensional CFD simulations that included the geometric details of the porous plate bleed systems were also carried out to supplement the experimental data, and provide additional insights into the bleed flow physics. Overall, both bleed formulations fared well for the tests performed in this study. However, the sample of test problems considered in this effort was not large enough to permit a comprehensive validation of the models.

Baurle, Robert A.↗

A Multi-scale Modeling System: Developments, Applications and Critical Issues

A multi-scale modeling framework (MMF), which replaces the conventional cloud parameterizations with a cloud-resolving model (CRM) in each grid column of a GCM, constitutes a new and promising approach. The MMF can provide for global coverage and two-way interactions between the CRMs and their parent GCM. The GCM allows global coverage and the CRM allows explicit simulation of cloud processes and their interactions with radiation and surface processes. A new MMF has been developed that is based the Goddard finite volume GCM (fvGCM) and the Goddard Cumulus Ensemble (GCE) model. This Goddard MMF produces many features that are similar to another MMF that was developed at Colorado State University (CSU), such as an improved .surface precipitation pattern, better cloudiness, improved diurnal variability over both oceans and continents, and a stronger, propagating Madden-Julian oscillation (MJO) compared to their parent GCMs using conventional cloud parameterizations. Both MMFs also produce a precipitation bias in the western Pacific during Northern Hemisphere summer. However, there are also notable differences between two MMFs. For example, the CSU MMF simulates less rainfall over land than its parent GCM. This is why the CSU MMF simulated less overall global rainfall than its parent GCM. The Goddard MMF overestimates global rainfall because of its oceanic component. Some critical issues associated with the Goddard MMF are presented in this paper.

Tao, Wei-Kuo↗

The Connected Isotopic Water Cycle in the Community Earth System Model Version 1

Because of the pervasive role of water in the Earth system, the relative abundances of stable isotopologues of water are valuable for understanding atmospheric, oceanic, and biospheric processes, and for interpreting paleoclimate proxy reconstructions. Isotopologues are transported by both largescale and turbulent flows, and the ratio of heavy to light isotopologues changes due to fractionation that can accompany condensation and evaporation processes. Correctly predicting the isotopic distributions requires resolving the relationships between largescale ocean and atmospheric circulation and smallerscale hydrological processes, which can be accomplished within a coupled climate modeling framework. Here we present the water isotopeenabled version of the Community Earth System Model version 1 (iCESM1), which simulates global variations in water isotopic ratios in the atmosphere, land, ocean, and sea ice. In a transient Last Millennium simulation covering the 850-2005 period, iCESM1 correctly captures the latetwentiethcentury structure of δ(exp 18)O and δD over the global oceans, with more limited accuracy over land. The relationship between salinity and seawater δ(exp 18)O is also well represented over the observational period, including interbasin variations. We illustrate the utility of coupled, isotopeenabled simulations using both Last Millennium simulations and freshwater hosing experiments with iCESM1. Closing the isotopic mass balance between all components of the coupled model provides new confidence in the underlying depiction of the water cycle in CESM, while also highlighting areas where the underlying hydrologic balance can be improved. The iCESM1 is poised to be a vital community resource for ongoing model development with both modern and paleoclimate applications.

Brady, E.↗

Effect of Boundary Conditions on Process-Induced Stresses in a Plain Weave Unit Cell

Woven polymer matrix composites (PMCs) are leveraged in aerospace applications for their desirable specific properties, yet they are vulnerable to high residual stresses during manufacturing and their complex geometry makes experimental results difficult to observe. Process modeling is needed to characterize the effects of the curing and predict end stress states. Finite element software can be used to model woven architectures, however accurate representation of processing conditions remains a challenge when it comes to selecting boundary conditions. The effect of BCs on process-induced stress within woven PMCs is studied. The commercial Finite Element Analysis (FEA) software Abaqus is coupled with user-written subroutines in a process modeling framework. A two-dimensionally (2D) woven PMC repeating unit cell (RUC) is modeled with TexGen and Abaqus. Virtual curing is imposed on the bulk matrix. The BC study is conducted with Free, Periodic, Flat, and Flat-Free configurations. Results show that the end stress state is sensitive to the boundary condition assumptions. Flat BC results show great agreement with Periodic BCs. Residual stress results from process modeling are then compared with a linear-elastic thermal cooldown analysis in Abaqus. Cooldown results indicate an overestimation in matrix stresses compared with process modeling.

micromechanics↗

Evaluating Approaches Relating Ecosystem Productivity with Desis Spectral Information

Data from the DLR Earth Sensing Imaging Spectrometer (DESIS), mounted on the International Space Station (ISS), were used to develop and test algorithms for remotely retrieving ecosystem productivity. Twenty DESIS images were used from three widely separated forested study sites representing deciduous and conifer forests. Gross primary production (GPP) values from eddy covariance flux towers at the sites were matched with DESIS spectral reflectances collected on the same days. Multiple algorithms were successful relating spectral reflectance with GPP, including: spectral vegetation indices (SVI) sensitive to chlorophyll content, SVI used in a photosynthetic light-use efficiency model framework, spectral shape characteristics through spectral derivatives and absorption feature analysis, and statistical models leading to multiband hyperspectral indices from partial least squares regression. Successful algorithms were able to achieve R2 better than 0.7 using a diverse set of observations combining data from different sites from multiple years and at multiple times during the year. The demonstrated robustness of the algorithms provides some confidence in using DESIS imagery to map spatial patterns of GPP.

K F Huemmrich↗

Evaluating Approaches Relating Ecosystem Productivity with DESIS Spectral Information

Data from the DLR Earth Sensing Imaging Spectrometer (DESIS), mounted on the International Space Station (ISS), were used to develop and test algorithms for remotely retrieving ecosystem productivity. Twenty DESIS images were used from three widely separated forested study sites representing deciduous and conifer forests. Gross primary production (GPP) values from eddy covariance flux towers at the sites were matched with DESIS spectral reflectances collected on the same days. Multiple algorithms were successful relating spectral reflectance with GPP, including: spectral vegetation indices (SVI) sensitive to chlorophyll content, SVI used in a photosynthetic light-use efficiency model framework, spectral shape characteristics through spectral derivatives and absorption feature analysis, and statistical models leading to multiband hyperspectral indices from partial least squares regression. Successful algorithms were able to achieve R2 better than 0.7 using a diverse set of observations combining data from different sites from multiple years and at multiple times during the year. The demonstrated robustness of the algorithms provides some confidence in using DESIS imagery to map spatial patterns of GPP.

Gross Primary Productivity (GPP)↗

Data-constrained Solar Modeling with GX Simulator

To facilitate the study of solarflares and active regions, we have created a modeling framework, the freelydistributed GX Simulator IDL package, that combines 3D magnetic and plasma structures with thermal andnonthermal models of the chromosphere, transition region, and corona. Its object-based modular architecture,which runs on Windows, Mac, and Unix/Linux platforms, offers the ability to either import 3D density andtemperature distribution models, or to assign numerically defined coronal or chromospheric temperatures anddensities, or their distributions, to each individual voxel. GX Simulator can apply parametric heating modelsinvolving average properties of the magneticfield lines crossing a given voxel, as well as compute and investigatethe spatial and spectral properties of radio,(sub)millimeter, EUV, and X-ray emissions calculated from the model,and quantitatively compare them with observations. The package includes a fully automatic model productionpipeline that, based on minimal users input, downloads the required SDO/HMI vector magneticfield data,performs potential or nonlinear force-freefield extrapolations, populates the magneticfield skeleton withparameterized heated plasma coronal models that assume either steady-state or impulsive plasma heating, andgenerates non-LTE density and temperature distribution models of the chromosphere that are constrained byphotospheric measurements. The standardized models produced by this pipeline may be further customizedthrough specialized IDL scripts, or a set of interactive tools provided by the graphical user interface. Here, wedescribe the GX Simulator framework and its applications.

Solar active regions↗

The Circulation of the Plasmasphere Fluid during the Erosion Event on September 8, 2017

A strong solar wind pressure pulse triggered the magnetic storm on September 7, 2017. Near the end of September 7, the z-component of the interplanetary magnetic field (IMF Bz) dropped from 9 to -10 nT in 30 min. The IMF Bz remained at the level of -10 nT for 2 hours and then had another rapid drop to -31 nT in 30 min. The sudden plunge of IMF Bz and the associated strong convection electric field stirred up the storm main phase with Dst falling from 0 to -122 nT from 2200 UT on September 7 to 0200 UT on the 8th. Severe plasmasphere erosion was observed on September 8 by multiple spacecraft, such as the Van Allen Probes and the Arase satellite. In this study, we examine the fate of the eroded plasmasphere particles during this event by model simulation as well as satellite data analysis. The simulation tool we use is the Space Weather Modeling Framework (SWMF)/Block-Adaptive Tree Solarwind Roe-type Upwind Scheme (BATS-R-US) model coupled with the Comprehensive Inner Magnetosphere-Ionosphere (CIMI) model. One of the distinctive capabilities of the SWMF/BATSRUS-CIMI model is that it treats the cold plasmas in the plasmasphere as a separate fluid in the MHD equations. As a result, the transport and circulation of the plasmasphere fluid in the global magnetosphere can be traced and the impacts of this cold fluid on the global magnetosphere can be evaluated. In this paper, we will show how the drainage plume is formed during the storm and how the plasmasphere fluid is transported to the flank and lobe regions and eventually to the plasma sheet and reenters into the plasmasphere. We will validate our simulation by plasmasphere signatures observed in both the inner and outer magnetosphere.

Mei-Ching Fok↗

Trending in Probability of Collision Measurements via a Bayesian Zero-Inflated Beta Mixed Model

We investigate the performance of a generalized linear mixed model in predicting the Probabilities of Collision (Pc) for conjunction events. Specifically, we apply this model to the log(sub 10) transformation of these probabilities and argue that this transformation yields values that can be considered bounded in practice. Additionally, this bounded random variable, after scaling, is zero-inflated. Consequently, we model these values using the zero-inflated Beta distribution, and utilize the Bayesian paradigm and the mixed model framework to borrow information from past and current events. This provides a natural way to model the data and provides a basis for answering questions of interest, such as what is the likelihood of observing a probability of collision equal to the effective value of zero on a subsequent observation.

Conjuction Assessment↗