Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Particle representation”

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.

267 records · Page 15
All sourcesNASA NTRSDOE OSTI
NASA NTRS1986

Transfer function analysis of thermospheric perturbations

Applying perturbation theory, a spectral model in terms of vectors spherical harmonics (Legendre polynomials) is used to describe the short term thermospheric perturbations originating in the auroral regions. The source may be Joule heating, particle precipitation or ExB ion drift-momentum coupling. A multiconstituent atmosphere is considered, allowing for the collisional momentum exchange between species including Ar, O2, N2, O, He and H. The coupled equations of energy, mass and momentum conservation are solved simultaneously for the major species N2 and O. Applying homogeneous boundary conditions, the integration is carred out from the Earth's surface up to 700 km. In the analysis, the spherical harmonics are treated as eigenfunctions, assuming that the Earth's rotation (and prevailing circulation) do not significantly affect perturbations with periods which are typically much less than one day. Under these simplifying assumptions, and given a particular source distribution in the vertical, a two dimensional transfer function is constructed to describe the three dimensional response of the atmosphere. In the order of increasing horizontal wave numbers (order of polynomials), this transfer function reveals five components. To compile the transfer function, the numerical computations are very time consuming (about 100 hours on a VAX for one particular vertical source distribution). However, given the transfer function, the atmospheric response in space and time (using Fourier integral representation) can be constructed with a few seconds of a central processing unit. This model is applied in a case study of wind and temperature measurements on the Dynamics Explorer B, which show features characteristic of a ringlike excitation source in the auroral oval. The data can be interpreted as gravity waves which are focused (and amplified) in the polar region and then are reflected to propagate toward lower latitudes.

Mayr, H. G.↗
NASA NTRS2021

Improved representation of the global dust cycle using observational constraints on dust properties and abundance

Even though desert dust is the most abundant aerosol by mass in Earth's atmosphere, atmospheric models struggle to accurately represent its spatial and temporal distribution. These model errors are partially caused by fundamental difficulties in simulating dust emission in coarse-resolution models and in accurately representing dust microphysical properties. Here we mitigate these problems by developing a new methodology that yields an improved representation of the global dust cycle. We present an analytical framework that uses inverse modeling to integrate an ensemble of global model simulations with observational constraints on the dust size distribution, extinction efficiency, and regional dust aerosol optical depth. We then compare the inverse model results against independent measurements of dust surface concentration and deposition flux and find that errors are reduced by approximately a factor of 2 relative to current model simulations of the Northern Hemisphere dust cycle. The inverse model results show smaller improvements in the less dusty Southern Hemisphere, most likely because both the model simulations and the observational constraints used in the inverse model are less accurate. On a global basis, we find that the emission flux of dust with a geometric diameter up to 20 µm (PM20) is approximately 5000 Tg yr−1, which is greater than most models account for. This larger PM20 dust flux is needed to match observational constraints showing a large atmospheric loading of coarse dust. We obtain gridded datasets of dust emission, vertically integrated loading, dust aerosol optical depth, (surface) concentration, and wet and dry deposition fluxes that are resolved by season and particle size. As our results indicate that this dataset is more accurate than current model simulations and the MERRA-2 dust reanalysis product, it can be used to improve quantifications of dust impacts on the Earth system.

Dust cycle↗
NASA NTRS2021

Improved Representation of the Global Dust Cycle Using Observational Constraints on Dust Properties and Abundance

Even though desert dust is the most abundant aerosol by mass in Earth’s atmosphere, atmospheric models 25 struggle to accurately represent its spatial and temporal distribution. These model errors are partially caused byfundamental difficulties in simulating dust emission in coarse-resolution models and in accurately representing dust microphysical properties. Here we mitigate these problems by developing a new methodology that yields an improved representation of the global dust cycle. We present an analytical framework that uses inverse modeling to integrate an ensemble of global model simulations with observational constraints on the dust size distribution, 30 extinction efficiency, and regional dust aerosol optical depth. We then compare the inverse model results against independent measurements of dust surface concentration and deposition flux and find that errors are reduced by approximately a factor of two relative to current model simulations of the Northern Hemisphere dust cycle. The inverse model results show smaller improvements in the less dusty Southern Hemisphere, most likely because both the model simulations and the observational constraints used in the inverse model are less accurate. On a global 35 basis, we find that the emission flux of dust with geometric diameter up to 20 μm (PM20) is approximately 5,000 Tg/year, which is greater than most models account for. This larger PM20 dust flux is needed to match observational constraints showing a large atmospheric loading of coarse dust. We obtain gridded data sets of dust emission, vertically integrated loading, dust aerosol optical depth, (surface) concentration, and wet and dry deposition fluxes that are resolved by season and particle size. As our results indicate that this data set is more accurate than current 40 model simulations and the MERRA-2 dust reanalysis product, it can be used to improve quantifications of dust impacts on the Earth system.

aerosols, dust↗
NASA NTRS2008

NASA Tech Briefs, August 2008

Customizable Digital Receivers for Radar Two-Camera Acquisition and Tracking of a Flying Target Visual Data Analysis for Satellites A Data Type for Efficient Representation of Other Data Types Hand-Held Ultrasonic Instrument for Reading Matrix Symbols Broadband Microstrip-to-Coplanar Strip Double-Y Balun A Topographical Lidar System for Terrain-Relative Navigation Programmable Low-Voltage Circuit Breaker and Tester Electronic Switch Arrays for Managing Microbattery Arrays Topics covered include: Lower-Dark-Current, Higher-Blue-Response CMOS Imagers; Fabricating Large-Area Sheets of Single-Layer Graphene by CVD; Support for Diagnosis of Custom Computer Hardware; Providing Goal-Based Autonomy for Commanding a Spacecraft; Dynamic Method for Identifying Collected Sample Mass; Optimal Planning and Problem-Solving; Attitude-Control Algorithm for Minimizing Maneuver Execution Errors; Grants Document-Generation System; Heat-Storage Modules Containing LiNO3 3H2O and Graphite Foam; Precipitation-Strengthened, High-Temperature, High-Force Shape Memory Alloys; Improved Relief Valve Would Be Less Susceptible to Failure; Safety Modification of Cam-and-Groove Hose Coupling; Using Composite Materials in a Cryogenic Pump; Using Electronic Noses to Detect Tumors During Neurosurgery; Producing Newborn Synchronous Mammalian Cells; Smaller, Lower-Power Fast-Neutron Scintillation Detectors; Rotationally Vibrating Electric-Field Mill; Estimating Hardness from the USDC Tool-Bit Temperature Rise; Particle-Charge Spectrometer; Automated Production of Movies on a Cluster of Computers; FIDO-Class Development Rover; and Tone-Based Command of Deep Space Probes Using Ground Antennas.

Source record↗
NASA NTRS2006

Various Ambiguities in Re-constructing Laser Pulse Parameters

We think that mode lock laser pulses are generated by the summation process that take place between the monochromatic EM filed frequencies as if they interact with each other as shown in equation 1. In reality, the pulse generation is a collaborative interaction process between EM fields and various material medium. When we carry out the actual mode lock analysis, we do take into account of interpaly between all the temporal dynamics of the cavity gain medium, cavity round trip time and the response time of the intra cavity element (saturable absorber, Kerr medium, etc.). that really enforces the locking of the phase of the cavity spontaneous emissions. On a conceptual level, this simplistic representation of the mode locking by Eq.1 ignores all these critical physical processes. When we try to analyze a pulsed field, again we start by representing it very much like this equation, even though we can only detect the square modulus of this complex field and loose a lot of phase related information to the detectors quantum whims and their time constants. The key parameters for a light pulse are as follows. Foremost is the (i) carrier frequency, which cannot be described or imagined without its state of undulation expressed as its (ii) phase. Next is our imagined time finite (iii) carrier envelope that provides the temporal boundary of the field amplitude strength of the undulating E-field. The final parameter is the (iv) state of polarization or the unique plane along which the strength of the E-field gradient undulates. None of these filed characteristics are made self-evident to us by the fields themselves. We do not see light. Light does not see light. Light beams pass through each other without altering each others energy distribution unless there are interacting material molecules (dipoles) within the physical volume of superposition of the beams. In contrast, we can sense the material particles. Material particles sense each other and they cannot pass through each other without interacting with (scattering from) each other. Thus the interpretation of the superposition phenomenon of multiple fields on detectors should not be lumped under the mysterious "wave-particle duality" philosophy. The phenomenon of superposition can be understood better when we focus on the actual process experienced by the detecting dipoles when allowed by QM rules, they respond to and sum all the induced stimulations due to all the superposed fields followed by the proportionate energy absorption giving rise to the fringes we observe. We will present various experimental results to illustrate our arguments. Our position is that such detector behavior driven interpretations rather than the generally implied field-field interaction driven explanations, will help us better understand the ultimate nature of light and hence invent better and newer devices and instruments.

Roychoudhuri, Chandrasekhar↗
NASA NTRS2007

Multi-Scale Modeling of Magnetospheric Reconnection

One of the major challenges in modeling the magnetospheric magnetic reconnection is to quantify the interaction between large-scale global magnetospheric dynamics and microphysical processes in diffusion regions near reconnection sites. There is still considerable debate as to what degree microphysical processes on kinetic scales affect the global evolution and how important it is to substitute numerical dissipation and/or ad hoc anomalous resistivity by a physically motivated model of dissipation. Comparative studies of magnetic reconnection in small scale geometries demonstrated that MHD simulations that included non-ideal processes in terms of a resistive term $\eta J$ did not produce the fast reconnection rates observed in kinetic simulations. For a broad range of physical parameters in collisionless magnetospheric plasma, the primary mechanism controlling the dissipation in the vicinity of the reconnection site is non-gyrotropic effects with spatial scales comparable with the particle Larmor radius. We utilize the global MHD code BATSRUS and incorporate nongyrotropic effects in diffusion regions in terms of corrections to the induction equation. We developed an algorithm to search for magnetotail reconnection sites, specifically where the magnetic field components perpendicular to the local current direction approaches zero and form an X-type configuration. Spatial scales of the diffusion region and magnitude of the reconnection electric field are calculated selfconsistently using MHD plasma and field parameters in the vicinity of the reconnection site. The location of the reconnection sites is updated during the simulations. To clarify the role of nongyrotropic effects in diffusion region on the global magnetospheric dynamic we perform simulations with steady southward IMF driving of the magnetosphere. Ideal MHD simulations with magnetic reconnection supported by numerical resistivity produce steady configuration with almost stationary near-earth neutral line (NENL). Simulations with non-gyrotropic corrections demonstrate dynamic quasi-periodic response to the steady driving condition. The loading/unloading cycle in non-gyrotropic MHD results has a non-stationary reconnection site in the magnetotail, with the retreating during the stretching phase and then a new NENL forming in the resulting thin plasma sheet. We expect that this model will lead to improved representations of space weather event in the magnetosphere.

Kuznetsova, M. M.↗
NASA NTRS2007

New, Improved Bulk-microphysical Schemes for Studying Precipitation Processes in WRF: Comparisons with Other Schemes - Part 1

Advances in computing power allow atmospheric prediction models to be mn at progressively finer scales of resolution, using increasingly more sophisticated physical parameterizations and numerical methods. The representation of cloud microphysical processes is a key component of these models, over the past decade both research and operational numerical weather prediction models have started using more complex microphysical schemes that were originally developed for high-resolution cloud-resolving models (CRMs). A recent report to the United States Weather Research Program (USWRP) Science Steering Committee specifically calls for the replacement of implicit cumulus parameterization schemes with explicit bulk schemes in numerical weather prediction (NWP) as part of a community effort to improve quantitative precipitation forecasts (QPF). An improved Goddard bulk microphysical parameterization is implemented into a state-of the-art of next generation of Weather Research and Forecasting (WRF) model. High-resolution model simulations are conducted to examine the impact of microphysical schemes on two different weather events (a midlatitude linear convective system and an Atllan"ic hurricane). The results suggest that microphysics has a major impact on the organization and precipitation processes associated with a summer midlatitude convective line system. The 31CE scheme with a cloud ice-snow-hail configuration led to a better agreement with observation in terms of simulated narrow convective line and rainfall intensity. This is because the 3ICE-hail scheme includes dense ice precipitating (hail) particle with very fast fall speed (over 10 m/s). For an Atlantic hurricane case, varying the microphysical schemes had no significant impact on the track forecast but did affect the intensity (important for air-sea interaction)

Tao, W.-K.↗
NASA NTRS2013

Jettison Engineering Trajectory Tool

The Jettison Engineering Trajectory Tool (JETT) performs the jettison analysis function for any orbiting asset. It provides a method to compute the relative trajectories between an orbiting asset and any jettisoned item (intentional or unintentional) or sublimating particles generated by fluid dumps to assess whether an object is safe to jettison, or if there is a risk with an item that was inadvertently lost overboard. The main concern is the interaction and possible recontact of the jettisoned object with an asset. This supports the analysis that jettisoned items will safely clear the vehicle, ensuring no collisions. The software will reduce the jettison analysis task from one that could take days to complete to one that can be completed in hours, with an analysis that is more comprehensive than the previous method. It provides the ability to define the jettison operation relative to International Space Station (ISS) structure, and provides 2D and 3D plotting capability to allow an analyst to perform a subjective clearance assessment with ISS structures. The developers followed the SMP to create the code and all supporting documentation. The code makes extensive use of the object-oriented format of Java and, in addition, the Model-View-Controller architecture was used in the organization of the code, allowing each piece to be independent of updates to the other pieces. The model category is for maintaining data entered by the user and generated by the analysis. The view category provides capabilities for data entry and displaying all or a portion of the analysis data in tabular, 2D, and 3D representation. The controller category allows for handling events that affect the model or view(s). The JETT utilizes orbital mechanics with complex algorithms. Since JETT is written in JAVA, it is essentially platform-independent.

Zaczek, Mariusz↗
NASA NTRS2020

Analysis of Sample Acquisition Dynamics Using Discrete Element Method

The analysis presented in this paper is conducted in the framework of the Ocean Worlds Autonomy Testbed for Exploration Research and Simulation (OceanWATERS) project, currently under development at NASA Ames Research Center. OceanWATERS aims at designing a simulation environment which allows for testing autonomy of scientific lander missions to the icy moons of our solar system. Mainly focused on reproducing the end effector interaction with the inherent terrain, this paper introduces a novel discrete element method (DEM)-based approach to determine forces and torques acting on the lander’s scoop during the sample acquisition process. An accurate force feedback from the terrain on the scoop is required by fault-detection and autonomous decision-making algorithms to identify when the requested torque on the robotic arm’s joints exceeds the maximum available torque. Knowledge of the terrain force feedback significantly helps evaluating the arm’s links structural properties and properly selecting actuators for the joints. Models available in literature constitute a partial representation of the dynamics of the interaction. As an example, Balovnev derived an analytical expression of the vertical and horizontal force acting on a bucket while collecting a sample as a function of its geometry and velocity, soil parameters and reached depth. Although the model represents an adequate approximation of the two force components, it ignores the direction orthogonal to the scoop motion and neglects the torque. This work relies on DEM analysis to compensate for analytical models’ deficiencies and inaccuracies, i. e. provide force and torque 3D vectors, defined in the moving reference (body) frame attached to the scoop, at each instant of the sample collection process. Results from the first presented analysis relate to the specific OceanWATERS sampling strategy, which consists of collecting the sample through five consecutive passes with increasing depth, each pass following the same circularlinear- circular trajectory. Data is collected given a specific scoop design interacting with two types of bulk materials, which may characterize the surface of icy planetary bodies: snow and ice. Although specifically concerned with the OceanWATERS design, this first analysis provides the expected force trends for similar sampling strategies and allows to deduce phenomenological information about the general scooping process. In order to further instruct the community on the use of DEM tools as a solution to the sampling collection problem, two more analyses have been carried out, mainly focused on reducing the DEM computation time, which increases with a decrease in particle size. After running a set of identical simulations, where the only changing parameter is the size of the spherical particle, it is observed that the resulting force trajectories, starting from a given particle size, converge to the true trend. It is deducible that a further decrease in size yields negligible improvements in the accuracy, while it sensibly increases computation time. A final analysis aims at discussing limitations of approximating bulk material particles having a complex shape, e. g. ice fragments, with spheres, by comparing force trends resulting in the two cases for the same simulation scenario.

Catanoso, Damiana↗
NASA NTRS2003

Mathematical model of diffusion-limited evolution of multiple gas bubbles in tissue

Models of gas bubble dynamics employed in probabilistic analyses of decompression sickness incidence in man must be theoretically consistent and simple, if they are to yield useful results without requiring excessive computations. They are generally formulated in terms of ordinary differential equations that describe diffusion-limited gas exchange between a gas bubble and the extravascular tissue surrounding it. In our previous model (Ann. Biomed. Eng. 30: 232-246, 2002), we showed that with appropriate representation of sink pressures to account for gas loss or gain due to heterogeneous blood perfusion in the unstirred diffusion region around the bubble, diffusion-limited bubble growth in a tissue of finite volume can be simulated without postulating a boundary layer across which gas flux is discontinuous. However, interactions between two or more bubbles caused by competition for available gas cannot be considered in this model, because the diffusion region has a fixed volume with zero gas flux at its outer boundary. The present work extends the previous model to accommodate interactions among multiple bubbles by allowing the diffusion region volume of each bubble to vary during bubble evolution. For given decompression and tissue volume, bubble growth is sustained only if the bubble number density is below a certain maximum.

Evaluation Studies↗
NASA NTRS2015

Thermodynamic Derivation of the Activation Energy for Ice Nucleation

Cirrus clouds play a key role in the radiative and hydrological balance of the upper troposphere. Their correct representation in atmospheric models requires an understanding of the microscopic processes leading to ice nucleation. A key parameter in the theoretical description of ice nucleation is the activation energy, which controls the flux of water molecules from the bulk of the liquid to the solid during the early stages of ice formation. In most studies it is estimated by direct association with the bulk properties of water, typically viscosity and self-diffusivity. As the environment in the ice-liquid interface may differ from that of the bulk, this approach may introduce bias in calculated nucleation rates. In this work a theoretical model is proposed to describe the transfer of water molecules across the ice-liquid interface. Within this framework the activation energy naturally emerges from the combination of the energy required to break hydrogen bonds in the liquid, i.e., the bulk diffusion process, and the work dissipated from the molecular rearrangement of water molecules within the ice-liquid interface. The new expression is introduced into a generalized form of classical nucleation theory. Even though no nucleation rate measurements are used to fit any of the parameters of the theory the predicted nucleation rate is in good agreement with experimental results, even at temperature as low as 190 K, where it tends to be underestimated by most models. It is shown that the activation energy has a strong dependency on temperature and a weak dependency on water activity. Such dependencies are masked by thermodynamic effects at temperatures typical of homogeneous freezing of cloud droplets; however, they may affect the formation of ice in haze aerosol particles. The new model provides an independent estimation of the activation energy and the homogeneous ice nucleation rate, and it may help to improve the interpretation of experimental results and the development of parameterizations for cloud formation.

Cirrus clouds↗
NASA NTRSDate not supplied

An Ensemble Investigation of the Causes for Regional Air-Quality Model Critical Load Exceedances Prediction Variability in European and North American Domains Using Diagnostics From Phase 4 of the Air Quality Model Evaluation International Initiative

We summarize tentative findings from multi air quality model ensembles for the years 2009 and 2010 in Europe (EU), and 2010 and 2016 in North America (NA), under AQMEII-4. The model predictions of sulphur and nitrogen deposition were used to estimate exceedances of critical loads for acidification and eutrophication, to show the extent to which the ensemble members agree in the magnitude and the trend of ecologically meaningful impacts. Model exceedance variability was analyzed using AQMEII-4 diagnostics. Evaluation against concentration and wet deposition observations, coupled with these diagnostics, identified specific process representations as the causes for variability between model predictions and for reduced model performance. All models predicted reductions in ecosystem acidification impacts in North America between the years 2010 and 2016, in accord with SO2 emissions reduction legislation which started in 2010 (SO2 SIP) However, all models in EU and NA domains had net negative biases for wet deposition of sulphur and nitrogen relative to observations. The wet S deposition average mean bias for the NA ensemble was -0.17 eq ha-1 d-1, and for the EU ensemble -1.15 eq ha-1 d-1. The NA daily wet deposition average mean bias for NH4+ was -0.37 eq ha-1d-1; EU -1.19 eq ha-1 d-1. The daily NA wet NO3- deposition average mean bias was -0.24 eq ha-1d-1; EU -0.69 eq ha-1 d-1. The members of the ensemble diverged (factor of 10) in their North American predictions for Ndep and consequently their eutrophication exceedances. The models with the highest eutrophication predictions also predicted the highest levels of gas-phase ammonia dry deposition (standard deviation of ammonia dry deposition flux across ensemble members was larger than the ensemble average). These models also had negative biases of predicted ammonia concentrations; average mean biases of -0.63 (satellite NH3) and -0.85 ppbv (surface NH3) compared to ensemble averages of -0.30 and -0.34 ppbv. Diagnostics showed that these differences resulted from the manner in which bidirectional ammonia fluxes were parameterized within these models. The second largest source of NA eutrophication prediction variability were models with positive biases in particulate ammonium and nitrate concentrations, and higher particle nitrogen deposition levels ( particle ammonium concentration bias +0.35 ug m-3; ensemble bias +0.15 ug m-3). We believe two factors may have led to these latter overestimates: higher levels of fine mode particle nitrate formation compared to other models (due to the use of an inorganic heterogeneous chemistry algorithm which did not take base cation chemistry into account), and updates to particle dry deposition velocities carried out in the absence of concurrent updates to wet scavenging algorithms. The relative importance of dry gas, dry particulate, and wet deposition towards total sulphur and nitrogen deposition totals differed between EU and North American domains, though all models had negative biases in wet deposition as noted above. Parallel and subsequent work suggests that multiphase hydrometeor scavenging may improve model wet deposition performance. An increased research focus is recommended for four model processes: multiphase hydrometeor scavenging, ammonia bidirectional fluxes, base cation chemistry and emissions, and particle dry deposition.

regional air-quality model↗
NASA NTRS2022

Satellite-Based Evaluation of AeroCom Model Bias in Biomass Burning Regions

Global models are widely used to simulate biomass burning aerosol (BBA). Exhaustive evaluations on model representation of aerosol distributions and properties are fundamental to assess health and climate impacts of BBA. Here we conducted a comprehensive comparison of Aerosol Comparisons between Observations and Models (AeroCom) project model simulations with satellite observations. A total of 59 runs by 18 models from three AeroCom Phase-III experiments (i.e., biomass burning emissions, CTRL16, and CTRL19) and 14 satellite products of aerosols were used in the study. Aerosol optical depth (AOD) at 550 nm was investigated during the fire season over three key fire regions reflecting different fire dynamics (i.e., deforestation-dominated Amazon, Southern Hemisphere Africa where savannas are the key source of emissions, and boreal forest burning in boreal North America). The 14 satellite products were first evaluated against AErosol RObotic NETwork (AERONET) observations, with large uncertainties found. But these uncertainties had small impacts on the model evaluation that was dominated by modeling bias. Through a comparison with Polarization and Directionality of the Earth’s Reflectances measurements with the Generalized Retrieval of Aerosol and Surface Properties algorithm (POLDER-GRASP), we found that the modeled AOD values were biased by −93 % to 152 %, with most models showing significant underestimations even for the state-of-the-art aerosol modeling techniques (i.e., CTRL19). By scaling up BBA emissions, the negative biases in modeled AOD were significantly mitigated, although it yielded only negligible improvements in the correlation between models and observations, and the spatial and temporal variations in AOD biases did not change much. For models in CTRL16 and CTRL19, the large diversity in modeled AOD was in almost equal measures caused by diversity in emissions, lifetime, and the mass extinction coefficient (MEC). We found that in the AeroCom ensemble, BBA lifetime correlated significantly with particle deposition (as expected) and in turn correlated strongly with precipitation. Additional analysis based on Cloud-Aerosol LIdar with Orthogonal Polarization (CALIOP) aerosol profiles suggested that the altitude of the aerosol layer in the current models was generally too low, which also contributed to the bias in modeled lifetime. Modeled MECs exhibited significant correlations with the Ångström exponent (AE, an indicator of particle size). Comparisons with the POLDER-GRASP-observed AE suggested that the models tended to overestimate the AE (underestimated particle size), indicating a possible underestimation of MECs in models. The hygroscopic growth in most models generally agreed with observations and might not explain the overall underestimation of modeled AOD. Our results imply that current global models contain biases in important aerosol processes for BBA (e.g., emissions, removal, and optical properties) that remain to be addressed in future research.

biomass burning aerosols↗
NASA NTRS2017

Remote Sensing of Crystal Shapes in Ice Clouds

Ice crystals in clouds exist in a virtually limitless variation of geometries. The most basic shapes of ice crystals are columnar or plate-like hexagonal prisms with aspect ratios determined by relative humidity and temperature. However, crystals in ice clouds generally display more complex structures owing to aggregation, riming and growth histories through varying temperature and humidity regimes. Crystal shape is relevant for cloud evolution as it affects microphysical properties such as fall speeds and aggregation efficiency. Furthermore, the scattering properties of ice crystals are affected by their general shape, as well as by microscopic features such as surface roughness, impurities and internal structure. To improve the representation of ice clouds in climate models, increased understanding of the global variation of crystal shape and how it relates to, e.g., location, cloud temperature and atmospheric state is crucial. Here, the remote sensing of ice crystal macroscale and microscale structure from airborne and space-based lidar depolarization observations and multi-directional measurements of total and polarized reflectances is reviewed. In addition, a brief overview is given of in situ and laboratory observations of ice crystal shape as well as the optical properties of ice crystals that serve as foundations for the remote sensing approaches. Lidar depolarization is generally found to increase with increasing cloud height and to vary with latitude. Although this variation is generally linked to the variation of ice crystal shape, the interpretation of the depolarization remains largely qualitative and more research is needed before quantitative conclusions about ice shape can be deduced. The angular variation of total and polarized reflectances of ice clouds has been analyzed by numerous studies in order to infer information about ice crystal shapes from them. From these studies it is apparent that pristine crystals with smooth surfaces are generally inconsistent with the data and thus crystal impurity, distortion or surface roughness is prevalent. However, conclusions about the dominating ice shapes are often inconclusive and contradictory and are highly dependent on the limited selection of shapes included in the investigations. Since ice crystal optical properties are mostly determined by the aspect ratios of the crystal components and their microscale structure, it is advised that remote sensing applications focus on the variation of these ice shape characteristics, rather than on the macroscale shape or habit. Recent studies use databases with nearly continuous ranges of crystal component aspect ratio and-or roughness levels to infer the variation of ice crystal shape from satellite and airborne remote sensing measurements. Here, the rationale and results of varying strategies for the remote sensing of ice crystal shape are reviewed. Observed systematic variations of ice crystal geometry with location, cloud height and atmospheric state suggested by the data are discussed. Finally, a prospective is given on the future of the remote sensing of ice cloud particle shapes.

clouds↗
NASA NTRS2022

Newly Identified Climatically and Environmentally Significant High-Latitude Dust Sources

Dust particles from high latitudes have a potentially large local, regional, and global significance to climate and the environment as short-lived climate forcers, air pollutants, and nutrient sources. Identifying the locations of local dust sources and their emission, transport, and deposition processes is important for understanding the multiple impacts of high-latitude dust (HLD) on the Earth’s systems. Here, we identify, describe, and quantify the source intensity (SI) values, which show the potential of soil surfaces for dust emission scaled to values 0 to 1 concerning globally best productive sources, using the Global Sand and Dust Storms Source Base Map (G-SDS-SBM). This includes 64 HLD sources in our collection for the northern (Alaska, Canada, Denmark, Greenland, Iceland, Svalbard, Sweden, and Russia) and southern (Antarctica and Patagonia) high latitudes. Activity from most of these HLD sources shows seasonal character. It is estimated that high-latitude land areas with higher (SI ≥ 0.5), very high (SI ≥ 0.7), and the highest potential (SI ≥ 0.9) for dust emission cover > 1 670 000 km 2 , > 560 000 km 2 , and > 240 000 km 2 , respectively. In the Arctic HLD region (≥ 60◦ N), land area with SI ≥ 0.5 is 5.5 % (1 035 059 km 2 ), area with SI ≥ 0.7 is 2.3 % (440 804 km 2 ), and area with SI ≥ 0.9 is 1.1 % (208 701 km 2 ). Minimum SI values in the northern HLD region are about 3 orders of magnitude smaller, indicating that the dust sources of this region greatly depend on weather conditions. Our spatial dust source distribution analysis modeling results showed evidence supporting a northern HLD belt, defined as the area north of 50◦ N, with a “transitional HLD-source area” extending at latitudes 50–58◦ N in Eurasia and 50–55◦ N in Canada and a “cold HLD-source area” including areas north of 60◦ N in Eurasia and north of 58◦ N in Canada, with currently “no dust source” area between the HLD and low-latitude dust (LLD) dust belt, except for British Columbia. Using the global atmospheric transport model SILAM, we estimated that 1.0 % of the global dust emission originated from the high-latitude regions. About 57 % of the dust deposition in snow- and ice-covered Arctic regions was from HLD sources. In the southern HLD region, soil surface conditions are favorable for dust emission during the whole year. Climate change can cause a decrease in the duration of snow cover, retreat of glaciers, and an increase in drought, heatwave intensity, and frequency, leading to the increasing frequency of topsoil conditions favorable for dust emission, which increases the probability of dust storms. Our study provides a step forward to improve the representation of HLD in models and to monitor, quantify, and assess the environmental and climate significance of HLD.

high latitude dust↗