Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “subsidence”

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 649 records · Page 36

Impacts of Snow Darkening Effects by Light Absorbing Aerosols on Hydroclimate of Eurasia During Boreal Spring and Summer

In this study, we used the NASA GEOS-5 climate model to investigate the impact of snow darkening by light-absorbing aerosols on the hydro-climate if Eurasia during boreal spring and summer. Two sets of 10-member ensemble model integrations with prescribed sea surface temperature were carried out for 10 simulated years (2002-2011); one includes snow-darkening effects (SDE) by light-absorbing aerosols, and one does not (NSDE). Differences between the two experiments in the hydro-climates over Eurasia were evaluated. Results show that SDE warming is most pronounced during the melting season due to strong snow-albedo feedback in the vicinity of the retreating seasonal snowline. SDE spurs a Wet-First-Dry-Later (WFDL) modulation of the surface energy and water balances, characterized by an accelerated (days-to-weeks) snowmelt along with excessive runoff and a warming and wetting of the land during the early part of the melting season. The fast snowmelt is followed by a prolonged (> months) warming and desiccation of the land, stemming from a deficit of late-season snowmelt water and increased surface evaporation. The prolonged warming and drying of the land atmosphere system persist through the summer, sustained by atmospheric conditions favorable for the development of atmospheric blocking, i.e., higher mid-to-upper tropospheric geopotential height, lower relative humidity, reduced cloudiness and enhanced atmospheric subsidence. Overall, SDE by light absorbing aerosols leads to a warmer and drier boreal summer hydro-climate, increasing the frequency of the top 5% and 1% extreme hot days (as defined by NSDE statistics) over western and northern Eurasia by approximately three-fold and tenfold, respectively.

Snow darkening↗

Aerosol-Radiation-Cloud Interactions in the South-East Atlantic: First Results from the ORACLES-2016 Deployment and Plans for Future Activities

Southern Africa produces almost a third of the Earths biomass burning (BB) aerosol particles. Particles lofted into the mid-troposphere are transported westward over the South-East (SE) Atlantic, home to one of the three permanent subtropical stratocumulus (Sc) cloud decks in the world. The SE Atlantic stratocumulus deck interacts with the dense layers of BB aerosols that initially overlay the cloud deck, but later subside and may mix into the clouds. These interactions include adjustments to aerosol-induced solar heating and microphysical effects, and their global representation in climate models remains one of the largest uncertainties in estimates of future climate. Hence, new observations over the SE Atlantic have significant implications for regional and global climate change predictions.Our understanding of aerosol-cloud interactions in the SE Atlantic is severely limited. Most notably, we are missing knowledge on the absorptive and cloud nucleating properties of aerosols, including their vertical distribution relative to clouds, on the locations and degree of aerosol mixing into clouds, on the processes that govern cloud property adjustments, and on the importance of aerosol effects on clouds relative to co-varying synoptic scale meteorology.We describe first results from various synergistic, international research activities aimed at studying aerosol-cloud interactions in the region:NASAs airborne ORACLES (ObseRvations of Aerosols Above Clouds and Their IntEractionS) deployment in AugustSeptember of 2016,the DoEs LASIC (Layered Atlantic Smoke Interactions with Clouds) deployment of the ARM Mobile Facility to Ascension Island (June 2016 October 2017), the ground-based components of CNRS AEROCLO-sA (Aerosols Clouds and Fog over the west coast of southern Africa), and ongoing regional-scale integrative, process-oriented science efforts as part of SEALS-sA (Sea Earth Atmosphere Linkages Study in southern Africa).We expect to describe experimental setups as well as showcase initial aerosol and cloud property distributions. Furthermore, we discuss the implementation of future activities in these programs in coordination with the UK Met Offices CLARIFY (CLoud-Aerosol-Radiation Interactions and Forcing) experiment in 2017.

aerosol↗

Sea-Level Rise: Projections for Maryland 2018

In fulfillment of requirements of the Maryland Commission on Climate Change Act of 2015, this report provides updated projections of the amount of sea-level rise relative to Maryland coastal lands that is expected into the next century. These projections represent the consensus of an Expert Group drawn from the Mid-Atlantic region. The framework for these projections is explicitly tied to the projections of global sea-level rise included in the Intergovernmental Panel on Climate Change Fifth Assessment (2014) and incorporates regional factors such as subsidence, distance from melting glaciers and polar ice sheets, and ocean currents. The probability distribution of estimates of relative sea-level rise from the baseline year of 2000 are provided over time and, after 2050, for three different greenhouse gas emissions pathways: Growing Emissions (RCP8.5), Stabilized Emissions (RCP4.5), and meeting the Paris Agreement (RCP2.6). This framework has been recently used in developing relative sea-level rise projections for California, Oregon, Washington, New Jersey, and Delaware as well as several metropolitan areas. The Likely range (66% probability) of the relative rise of mean sea level expected in Maryland between 2000 and 2050 is 0.8 to 1.6 feet, with about a one-in-twenty chance it could exceed 2.0 feet and about a one-in-one hundred chance it could exceed 2.3 feet. Later this century, rates of sea-level rise increasingly depend on the future pathway of global emissions of greenhouse gases during the next sixty years. If emissions continue to grow well into the second half of the 21st century, the Likely range of sea-level rise experienced in Maryland is 2.0 to 4.2 feet over this century, two to four times the sea-level rise experienced during the 20th century. Moreover, there is a one-in-twenty chance that it could exceed 5.2 feet. If, on the other hand, global society were able to bring net greenhouse gas emissions to zero in time to meet the goals of the Paris Climate Agreement and reduce emissions sufficient to limit the increase in global mean temperature to less than 2°Celsius over pre-industrial levels, the Likely range for 2100 is 1.2 to 3.0 feet, with a one-in-twenty chance that it would exceed 3.7 feet. The difference in sea-level rise between these contrasting scenarios would diverge even more during the next century, with the failure to reduce emissions in the near term resulting in much greater sea-level rise 100 years from now. Moreover, recent research suggests that, without imminent and substantial reductions in greenhouse gas emissions, the loss of polar ice sheets-and thus the rate of sea-level rise-may be more rapid than assumed in these projections, particularly under the Growing Emissions scenario. These probabilistic sea-level rise projections can and should be used in planning and regulation, infrastructure siting and design, estimation of changes in tidal range and storm surge, developing inundation mapping tools, and adaptation strategies for high-tide flooding and saltwater intrusion.

Boesch, D. F.↗

Retrieval of Daily Maps of PM2.5 Aerosol in the Problematic California Valleys: Bright, Speckled Reflectances, Thin AOT, but High Pollution

The San Joaquin Valley suffers from severe episodes of respirable aerosol (PM2.5) in wintertime.We provide maps of aerosol episodes using daily snapshots of PM2.5 and its changing features despite numerous difficulties inherent to sampling the region. Linear relationships relating aerosol optical thickness (AOT) to PM2.5 give an explained variance of approximately 3.The GEO-CAPE mission has as a goal the provision of relevant measures of respirable aerosol to the community,but has not formulated a science goal beyond the limited goal of retrieval of AOT, bringing the usefulness of GEO-CAPE into doubt.Our special focus was on the DISCOVER-AQ period, Jan-Feb 2013, which had many supporting measurements.Both high pollution and retrieval difficulties tend to occur in many Mediterranean agricultural regions like the San Joauin. One difficulty is the relatively bright surfaces with considerable exposed soil. NASAs MAIAC and MODIS Deep Blue retrieval techniques are shown to have considerable skill even at low aerosol optical thickness (AOT) values, as evaluated by concurrent AERONET sunphotometer measurements.More significantly, these AOT values can correspond to high daytime PM2.5 since aerosol mixed layer depth is thin and variable, 200m 600 m. The thin layers derive from typical subsidence of dry air between more stormy periods. This situation provides an advantage: water vapor column is also almost completely limited to a similar mixed layer depth, and can thus serve as a measure of aerosol dilution.Using the MAIAC Water Vapor Column:In order to make the maps below, we used the MAIAC data but subtracted partial water-vapor columns estimated from MERRA Reanalysis Data availabe from the GSFC GMAO using kriging. We did not use the mixed-layer estimates from MERRA, since such analyses were found problematic during our forecasting exercises for DISCOVER-AQ. Observations from the aircraft soundings suggested that this overlying moisture was mostly due to larger scale flows, not ML venting.However, the specific humidity at the surface and a nearly well-mixed ML was analyzed by kriging from the surface network (MesoWest, University of Utah). These were thought to be truer, uninfluenced by physical process modeling that combines with data observations. (TBD: How different are they?) Procedure: Subtract overlying partial water columns from MAIAC column water and divide this by a surface value of water vapor. (MAIAC column is expressed in cm of water, i.e., water vapor at surface conditions.This method appears to bring out useful details in the distribution of submicron particles in the very problematic Wintertime San Joaquin Valley, and allow analysis of pollution episodes throughout the valley, rather than long-term averages.

san joaquin valley↗

An Automated Detection Methodology for Dry Well-Mixed Layers

The intense surface heating over arid land surfaces produces dry well-mixed layers (WML) via dry convection. These layers are characterized by nearly constant potential temperature and low, nearly constant water vapor mixing ratio. To further the study of dry WMLs, we created a detection methodology and supporting software to automate the identification and characterization of dry WMLs from multiple data sources including rawinsondes, remote sensing platforms, and model products. The software is a modular code written in Python, an open source language. Radiosondes from a network of synoptic stations in North Africa were used to develop and test the WML detection process. The detection involves an iterative decision tree that ingests a vertical profile from an input data file, performs a quality check for sufficient data density, and then searches upward through the column for successive points where the simultaneous changes in water vapor mixing ratio and potential temperature are less than the specified maxima. If points in the vertical profile meet the dry WML identification criteria, statistics are generated detailing the characteristics of each layer in the profile. At the end of the vertical profile analysis, there is an option to plot analyzed profiles in a variety of file formats. Initial results show that the detection methodology can be successfully applied across a wide variety of input data and North African environments and for all seasons. It is sensitive enough to identify dry WMLs from other types of isentropic phenomena such as subsidence layers and distinguish the current day’s dry WML from previous days.

Stephen D. Nicholls↗

P43K-3884 Characterizing Martian Crater Circulations with the NASA Ames Mars GCM

Observations made in Gale Crater by instruments on the MSL Curiosity Rover show that the diurnal amplitude of the surface pressure is increased and the depth of the Convective Boundary Layer (CBL) is decreased relative to other lander locations on flatter regions of Mars (Haberle et al., 2014; Moores et al., 2015). Mesoscale modeling studies of Gale Crater suggest that crater circulations produce these effects. Tyler & Barnes (2013) show that local upslope/downslope flows along the crater rim and Mt. Sharp amplify the diurnal pressure cycle. These same flows are thought to be at least partly responsible for the suppression of the CBL because upward air flow at the rim and in the center (due to Mt. Sharp) forces subsidence over the lowest regions of the crater during the day. Regional flows, largely due to the location of Gale near the dichotomy boundary, may also play a role in shaping the circulation internal to the crater. Whether the behavior of the CBL and the amplified diurnal pressure cycle are phenomena observed in craters morphologically different from Gale (i.e. bowl-shaped, irregular, degraded) is not yet understood. We will explore these questions by characterizing the behavior of these processes as they are shaped by the morphology of craters greater than 100 km in diameter. We use the NASA Ames Mars Global Circulation Model (GCM) that now utilizes the NOAA/GFDL cubed-sphere finite-volume dynamical core to examine ~100 craters of varying size and shape from a database of known Martian craters (Robbins & Hynek, 2014). Run at 7.5 km resolution, the GCM is capable of resolving surface winds, temperature, and pressure inside craters of this size allowing for the analysis of dozens of craters simulated at various seasons and within the context of synoptic and global-scale phenomena.

Batterson, Courtney MaryLou↗

Abundance Measurements of Titan's Stratospheric HCN, HC3N, C3H4, and CH3CN from ALMA Observations

Previous investigations have employed more than 100 close observations of Titan by the Cassini orbiter to elucidate connections between the production and distribution of Titan's vast, organic-rich chemical inventory and its atmospheric dynamics. However, as Titan transitions into northern summer, the lack of incoming data from the Cassini orbiter presents a potential barrier to the continued study of seasonal changes in Titan's atmosphere. In our previous work, we demonstrated that the Atacama Large Millimeter/submillimeter Array (ALMA) is well suited for measurements of Titan's atmosphere in the stratosphere and lower mesosphere (~100 - 500 km) through the use of spatially resolved (beam sizes <1") flux calibration observations of Titan. Here, we derive vertical abundance profiles of four of Titan's trace atmospheric species from the same 3 independent spatial regions across Titan's disk during the same epoch (2012-2015): HCN, HC3N, C3H4, and CH3CN. We find that Titan's minor constituents exhibit large latitudinal variations, with enhanced abundances at high latitudes compared to equatorial measurements; this includes CH3CN, which eluded previous detection by Cassini in the stratosphere, and thus spatially resolved abundance measurements were unattainable. Even over the short 3-year period, vertical profiles and integrated emission maps of these molecules allow us to observe temporal changes in Titan's atmospheric circulation during northern spring. Our derived abundance profiles are comparable to contemporary measurements from Cassini infrared observations, and we find additional evidence for subsidence of enriched air onto Titan's south pole during this time period. Continued observations of Titan with ALMA beyond the summer solstice will enable further study of how Titan's atmospheric composition and dynamics respond to seasonal changes.

Thelen, Alexander E.↗

Evaluation of Modeled Precipitation in Oceanic Extratropical Cyclones Using IMERG

Using the high spatial and temporal resolution precipitation dataset Integrated Multi-satellitE Retrievals for GPM (IMERG), extratropical cyclone precipitation is evaluated in two reanalyses and two climate models. Based on cyclone-centered composites, all four models overestimate precipitation in the western subsiding and dry side of the cyclones, and underestimate the precipitation in the eastern ascending and moist side. By decomposing the composites into frequency of occurrence and intensity (mean precipitation rate when precipitating), the analysis reveals a tendency for all four models to overestimate frequency and underestimate intensity, with the former issue dominating in the western half and the latter in the eastern half of the cyclones. Differences in frequency are strongly dependent on cyclone environmental moisture, while the differences in intensity are strongly impacted by the strength of ascent within the cyclone. There are some uncertainties associated with the observations: IMERG might under-report frozen precipitation and possibly exaggerate rates in vigorously ascending regions. Nevertheless, the analysis suggests that all models produce extratropical cyclone precipitation too often and too lightly. These biases have consequences when evaluating the changes in precipitation characteristics with changes in cyclone properties: the models disagree on the magnitude of the change in precipitation intensity with a change in environmental moisture and in precipitation frequency with a change in cyclone strength. This complicates accurate predictions of precipitation changes in a changing climate.

Catherine M Naud↗

Changes in Tropical Clouds and Atmospheric Circulation Associated with Rapid Adjustment Induced by Increased Atmospheric CO2 – A Multiscale Modeling Framework Study

The radiative heating increase due to increased CO2 concentration is the primary source for the rapid adjustment of atmospheric circulation and clouds. In this study, we investigate the rapid adjustment resulting from doubling of CO2 and its physical mechanism using a multiscale modeling framework (MMF). The MMF includes an advanced higher-order turbulence closure in its cloud-resolving model component and simulates realistic shallow and deep cloud climatology and boundary layer turbulence. The rapid adjustment over the tropics is characterized by 1) reduced ascent and descent strengths over the ocean, 2) increased lower tropospheric stability (LTS) over the subsidence region, 3) shoaling of planetary boundary layers over the ocean, 4) increased deep convection over lands and shift of cloud coverage from the ocean to lands, and 5) reduced sensible (SH) and latent heat (LH) fluxes over the oceanic deep convective regions. Unlike conventional general circulation models and another MMF, a reduction in the global-mean shortwave cloud radiative cooling is not simulated, due to the increase in low clouds at lower altitudes over the ocean, resulting from reduced cloud-top entrainment due to strengthened inversion. Changes in regional circulation play a key role in cloud changes and shift of cloud coverage to lands. Weaker energy transport resulting from water vapor and cloud CO2 masking effects reduces the upward motion and convective clouds in the oceanic regions. The ocean-land transports are linked to the partitioning of surface SH and LH fluxes that increases humidity over lands and enhances deep convection over the tropical lands.

Xu, Kuan-Man↗

Modal Analysis of a Two-Parachute System

The Orion capsule is designed to land under a nominal configuration of three main parachutes; however, the system is required to be fault tolerant and land successfully if one of the main parachutes fails to open. The Capsule Parachute Assembly System (CPAS) Team performed a series of drop tests in order to characterize the performance of the system with two main parachutes. During the series of drop tests, several distinct dynamical modes were observed. The most consequential of these is the pendulum mode. Three other modes are benign: flyout (scissors), maypole, and breathing. The actual multi-body system is nonlinear, flexible, and possesses significant cross-coupling. Rather than perform analysis of this highly complex system directly, we conduct analysis of each dynamical mode observed during flight, based on first principles. This approach is analogous to traditional aircraft flight dynamics analysis in which the full nonlinear behavior of the airframe is decomposed into longitudinal dynamics (phugoid and short-period modes) and lateral dynamics (spiral, roll-subsidence, and dutch-roll modes). This analysis is intended to supplement multi-body nonlinear simulations in order to provide further insight into the system.

Pei, Jing↗

Modal Analysis of the Orion Capsule Two Parachute System

As discussed in Ref [1], it is apparent from flight tests that the system made up of two main parachutes and a capsule can undergo several distinct dynamical behaviors. The most significant and problematic of these is the pendulum mode in which the system develops a pronounced swinging motion with an amplitude of up to 24 deg. Large excursions away from vertical by the capsule could cause it to strike the ground at a large horizontal or vertical speed and jeopardize the safety of the astronauts during a crewed mission. In reference [1], Ali et al. summarized a series of efforts taken by the Capsule Parachute Assembly System (CPAS) Program to understand and mitigate the pendulum issue. The period of oscillation and location of the system's pivot point are determined from post-flight analysis. Other noticeable but benign modes include: 1) flyout (scissors) mode, where the parachutes move back and forth symmetrically with respect to the vertical axis similar to the motion of a pair of scissors; 2) maypole mode, where the two parachutes circle around the vertical axis at a nearly constant radius and period; and 3) breathing mode, in which deformation of the non-rigid canopies affects the axial acceleration of the system in an oscillatory manner. Because these modes are relatively harm- less, little effort has been devoted to analyzing them in comparison with the pendulum motion. Motions of the actual system made up of two parachutes and a capsule are extremely complicated due to nonlinearities and flexibility effects. Often it is difficult to obtain insight into the fundamental dynamics of the system by examining results from a multi-body simulation based on nonlinear equations of motion (EOMs). As a part of this study, the dynamics of each mode observed during flight is derived from first principles on an individual basis by making numerous simplifications along the way. The intent is to gain a better understanding into the behavior of the complex multi-body system by studying the reduced set of differential equations associated with each mode. This approach is analogous to the traditional modal analysis technique used to study airplane flight dynamics, in which the full nonlinear behavior of the airframe is decomposed into the phugoid and short period modes for the longitudinal dynamics and the spiral, roll-subsidence, and dutch-roll modes for the lateral dynamics. It is important to note that the study does not address the mechanisms that cause the system to transition from one mode to another, nor does it discuss motions during which two or more modes occur simultaneously.

Pei, J.↗

Untangling causality in midlatitude aerosol–cloud adjustments

Aerosol–cloud interactions represent the leading uncertainty in our ability to infer climate sensitivity from the observational record. The forcing from changes in cloud albedo driven by increases in cloud droplet number (Nd) (the first indirect effect) is confidently negative and has narrowed its probable range in the last decade, but the sign and strength of forcing associated with changes in cloud macrophysics in response to aerosol (aerosol–cloud adjustments) remain uncertain. This uncertainty reflects our inability to accurately quantify variability not associated with a causal link flowing from the cloud microphysical state to the cloud macrophysical state. Once variability associated with meteorology has been removed, covariance between the liquid water path (LWP) averaged across cloudy and clear regions (here characterizing the macrophysical state) and Nd (characterizing the microphysical) is the sum of two causal pathways linking Nd to LWP: Nd altering LWP (adjustments) and precipitation scavenging aerosol and thus depleting Nd. Only the former term is relevant to constraining adjustments, but disentangling these terms in observations is challenging. We hypothesize that the diversity of constraints on aerosol–cloud adjustments in the literature may be partly due to not explicitly characterizing covariance flowing from cloud to aerosol and aerosol to cloud. Here, we restrict our analysis to the regime of extratropical clouds outside of low-pressure centers associated with cyclonic activity. Observations from MAC-LWP (Multisensor Advanced Climatology of Liquid Water Path) and MODIS are compared to simulations in the Met Office Unified Model (UM) GA7.1 (the atmosphere model of HadGEM3-GC3.1 and UKESM1). The meteorological predictors of LWP are found to be similar between the model and observations. There is also agreement with previous literature on cloud-controlling factors finding that increasing stability, moisture, and sensible heat flux enhance LWP, while increasing subsidence and sea surface temperature decrease it. A simulation where cloud microphysics are insensitive to changes in Nd is used to characterize covariance between Nd and LWP that is induced by factors other than aerosol–cloud adjustments. By removing variability associated with meteorology and scavenging, we infer the sensitivity of LWP to changes in Nd. Application of this technique to UM GA7.1 simulations reproduces the true model adjustment strength. Observational constraints developed using simulated covariability not induced by adjustments and observed covariability between Nd and LWP predict a 25 %–30 % overestimate by the UM GA7.1 in LWP change and a 30 %–35 % overestimate in associated radiative forcing.

Daniel T. McCoy↗

Evaluation of CFD as a Surrogate for Wind-Tunnel Testing for Mach 2.4 to 4.6 - Project Overview

The debate over when wind-tunnel testing (WTT) will be replaced by Computational Fluid Dynamics (CFD) comes and goes. More recently the debate has subsided with a more collaborative spirit between practitioners of these two disciplines resulting in significant improvements in the outcomes of both. There may come a time, however, when CFD has sufficient accuracy to supplant WTT as the dominant or perhaps only tool for aerodynamic simulation. If and/or when that happens, financial pressure result in efforts to close or severely limit the operations of wind tunnels. Presumably additional resources will go toward CFD in order to generate aerodynamic databases, load environments, and new aero/fluid-dynamic knowledge. It is therefore important to develop appropriate processes by which wind-tunnel closure decisions are made to ensure that facilities critical to industry and government research and development aren’t closed prematurely without ensuring that the existing CFD tools have sufficient accuracy and low-enough cost (and enough experts and computational facilities) to take on the traditional role of wind tunnels. This paper will describe a project intended to answer the specific question of whether CFD can replace WTT for the limited Mach-number range 2.4 to 4.6. The wind tunnel being examined in this context is the high-speed leg of the Unitary Plan Wind Tunnel at NASA’s Langley Research Center (LaRC UPWT).

James C. Ross↗

Mapping the Zonal Structure of Titan’s Northern Polar Vortex

Saturn exhibits an obliquity of 26.7° such that the largest moon, Titan, experiences seasonal variations including the formation of a polar vortex in the winter hemisphere. Titan’s polar vortex is characterised by cold stratospheric temperatures due to the lack of insolation over the winter pole, and an increase in trace gas abundance as a result of complex organic chemistry in the upper atmosphere combined with polar subsidence. Meridional variations in temperature and gas abundance across the vortex have previously been investigated, but there has not yet been any in-depth study of the zonal variations in the temperature or composition of the northern vortex. Here we present the first comprehensive two-dimensional seasonal mapping of Titan’s northern winter vortex. Using 18 nadir mapping sequences observed by the Composite InfraRed Spectrometer (CIRS) instrument on-board Cassini, we investigate the evolution of the vortex over almost half a Titan year, from late winter through to mid summer (L(sub s) = 326 - 86°, 2007-2017). We find the stratospheric symmetry axis to be tilted from the solid body rotation axis by around 3.5°, although our results for the azimuthal orientation of the tilt are inconclusive. We find that the northern vortex appears to remain zonally uniform in both temperature and composition at all times. A comparison with vortices observed on Earth, Mars, and Venus shows that large-scale wave mechanisms that are important on other terrestrial planets are not as significant in Titan’s atmosphere. This allows the northern vortex to be more symmetrical and persist longer throughout the annual cycle compared to other terrestrial planets.

Jason Sharkey↗

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–1 K 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).

Marc Mallet↗

Martian subsurface cryosalt expansion and collapse as trigger for landslides

On Mars, seasonal martian flow features known as recurring slope lineae (RSL) are prevalent on sun-facing slopes and are associated with salts. On Earth, subsurface interactions of gypsum with chlorides and oxychlorine salts wreak havoc: instigating sinkholes, cave collapse, debris flows, and upheave. Here, we illustrate (i) the disruptive potential of sulfate-chloride reactions in laboratory soil crust experiments, (ii) the formation of thin films of mixed ice-liquid water “slush” at −40° to −20°C on salty Mars analog grains, (iii) how mixtures of sulfates and chlorine salts affect their solubilities in low-temperature environments, and (iv) how these salt brines could be contributing to RSL formation on Mars. Our results demonstrate that interactions of sulfates and chlorine salts in fine-grained soils on Mars could absorb water, expand, deliquesce, cause subsidence, form crusts, disrupt surfaces, and ultimately produce landslides after dust loading on these unstable surfaces.

J. L. Bishop↗

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↗

Potential vorticity structure of Titan's polar vortices from Cassini CIRS observations

The Cassini mission has provided the best opportunity to date to extensively study the seasonal variation in Titan’s atmosphere, with observations spanning almost half a Titan year (𝐿𝑠 = 293 − 93°). An important feature in the Titan middle-atmosphere is the formation of a polar vortex. Observations have shown that an initially well-developed northern vortex enriched with trace gas species gradually breaks down after spring equinox as a new vortex emerges in southern winter. Here we use Cassini CIRS observations to derive the temperature and composition of the middle-atmosphere. We use the gradient wind equation to first estimate the mean zonal winds, and then the Potential Vorticity (PV) throughout Titan’s atmosphere over the timespan of the Cassini mission. PV is a useful diagnostic quantity for studying the dynamics of polar vortices because it is materially conserved for adiabatic and frictionless flows, and can be inverted to find all other dynamical fields. Our results show the formation of a strong zonal jet in the winter hemisphere, with wind velocities reaching 220 m/s, which is consistent with previous studies. An annular PV structure is also observed over the winter poles, whereby a ring of PV encircles a local minima over the pole. Such distributions are often found to be unstable without a restoring force, yet they are seen here in numerous observations in both the northern and southern hemispheres. A comparison with the annular Martian vortices shows that latent heat release from condensation or subsidence-induced adiabatic heating may explain the origin and stability of the annulus. Finally, we investigate the evolution of the size of the vortices and the role of strong PV gradients as a dynamical mixing barrier for trace gas species across the vortex edge. We find that longer lived gases are less confined to the vortex than those with shorter photochemical lifetimes.

Conor Andrew Nixon↗