Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ROUGHNESS”

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 1,315 records · Page 73

Vertical Distribution of Aersols and Water Vapor Using CRISM Limb Observations

Near-infrared spectra taken in a limb-viewing geometry by the Compact Reconnaissance Imaging Spectrometer for Mars (CRISM) on-board the Mars Reconnaissance Orbiter (MRO) provide a useful tool for probing atmospheric structure. Specifically, the observed radiance as a function of wavelength and height above the limb allows the vertical distribution of both dust and ice aerosols to be retrieved. These data serve as an important supplement to the aerosol profiling provided by the MRO/MCS instrument allowing independent validation and giving additional information on particle physical and scattering properties through multi-wavelength studies. A total of at least ten CRISM limb observations have been taken so far covering a full Martian year. Each set of limb observations nominally contains about four dozen scans across the limb giving pole-to-pole coverage for two orbits at roughly 100 and 290 W longitude over the Tharsis and Syrtis/Hellas regions, respectively. At each longitude, limb scans are spaced roughly 10 degrees apart in latitude, with a vertical spatial resolution on the limb of roughly 800 m. Radiative transfer modeling is used to model the observations. We compute synthetic CRISM limb spectra using a discrete-ordinates radiative transfer code that accounts for multiple scattering from aerosols and accounts for spherical geometry of the limb observations by integrating the source functions along curved paths in that coordinate system. Retrieved are 14-point vertical profiles for dust and water ice aerosols with resolution of 0.4 scale heights between one and six scale heights above the surface. After the aerosol retrieval is completed, the abundances of C02 (or surface pressure) and H20 gas are retrieved by matching the depth of absorption bands at 2000 nm for carbon dioxide and at 2600 run for water vapor. In addition to the column abundance of water vapor, limited information on its vertical structure can also be retrieved depending on the signal available from aerosol scattering.

Smith, Michael D.↗

Pneumatically Actuated Miniature Peristaltic Vacuum Pumps

Pneumatically actuated miniature peristaltic vacuum pumps have been proposed for incorporation into advanced miniature versions of scientific instruments that depend on vacuum for proper operation. These pumps are expected to be capable of reaching vacuum-side pressures in the torr to millitorr range (from .133 down to .0.13 Pa). Vacuum pumps that operate in this range are often denoted roughing pumps. In comparison with previously available roughing pumps, these pumps are expected to be an order of magnitude less massive and less power-hungry. In addition, they would be extremely robust, and would operate with little or no maintenance and without need for oil or other lubricants. Portable mass spectrometers are typical examples of instruments that could incorporate the proposed pumps. In addition, the proposed pumps could be used as roughing pumps in general laboratory applications in which low pumping rates could be tolerated. The proposed pumps could be designed and fabricated in conventionally machined and micromachined versions. A typical micromachined version (see figure) would include a rigid glass, metal, or plastic substrate and two layers of silicone rubber. The bottom silicone layer would contain shallow pump channels covered by silicone arches that could be pushed down pneumatically to block the channels. The bottom silicone layer would be covered with a thin layer of material with very low gas permeability, and would be bonded to the substrate everywhere except in the channel areas. The top silicone layer would be attached to the bottom silicone layer and would contain pneumatic- actuation channels that would lie crosswise to the pump channels. This version is said to be micromachined because the two silicone layers containing the channels would be fabricated by casting silicone rubber on micromachined silicon molds. The pneumatic-actuation channels would be alternately connected to a compressed gas and (depending on pump design) either to atmospheric pressure or to a partial vacuum source. The design would be such that the higher pneumatic pressure would be sufficient to push the silicone arches down onto the substrates, blocking the channels. Thus, by connecting pneumatic- actuation channels to the two pneumatic sources in spatial and temporal alternation, waves of opening and closing, equivalent to peristalsis, could be made to move along the pump channels. A pump according to this concept could be manufactured inexpensively. Pneumatic sources (compressors and partial vacuum sources) similar those needed for actuation are commercially available; they typically have masses of .100 g and power demands of the order of several W. In a design-optimization effort, it should be possible to reduce masses and power demands below even these low levels and to integrate pneumatic sources along with the proposed pumps into miniature units with overall dimensions of no more than a few centimeters per side.

Feldman, Sabrina↗

Process for Smoothing an Si Substrate after Etching of SiO2

A reactive-ion etching (RIE) process for smoothing a silicon substrate has been devised. The process is especially useful for smoothing those silicon areas that have been exposed by etching a pattern of holes in a layer of silicon dioxide that covers the substrate. Applications in which one could utilize smooth silicon surfaces like those produced by this process include fabrication of optical waveguides, epitaxial deposition of silicon on selected areas of silicon substrates, and preparation of silicon substrates for deposition of adherent metal layers. During etching away of a layer of SiO2 that covers an Si substrate, a polymer becomes deposited on the substrate, and the substrate surface becomes rough (roughness height approximately equal to 50 nm) as a result of over-etching or of deposition of the polymer. While it is possible to smooth a silicon substrate by wet chemical etching, the undesired consequences of wet chemical etching can include compromising the integrity of the SiO2 sidewalls and undercutting of the adjacent areas of the silicon dioxide that are meant to be left intact. The present RIE process results in anisotropic etching that removes the polymer and reduces height of roughness of the silicon substrate to less than 10 nm while leaving the SiO2 sidewalls intact and vertical. Control over substrate versus sidewall etching (in particular, preferential etching of the substrate) is achieved through selection of process parameters, including gas flow, power, and pressure. Such control is not uniformly and repeatably achievable in wet chemical etching. The recipe for the present RIE process is the following: Etch 1 - A mixture of CF4 and O2 gases flowing at rates of 25 to 75 and 75 to 125 standard cubic centimeters per minute (stdcm3/min), respectively; power between 44 and 55 W; and pressure between 45 and 55 mtorr (between 6.0 and 7.3 Pa). The etch rate lies between approximately equal to 3 and approximately equal to 6 nm/minute. Etch 2 - O2 gas flowing at 75 to 125 stdcm3/min, power between 44 and 55 W, and pressure between 50 and 100 mtorr (between 6.7 and 13.3 Pa).

Turner, Tasha↗

Comparison of PARASOL Observations with Polarized Reflectances Simulated Using Different Ice Habit Mixtures

Insufficient knowledge of the habit distribution and the degree of surface roughness of ice crystals within ice clouds is a source of uncertainty in the forward light scattering and radiative transfer simulations required in downstream applications involving these clouds. The widely used MODerate Resolution Imaging Spectroradiometer (MODIS) Collection 5 ice microphysical model assumes a mixture of various ice crystal shapes with smooth-facets except aggregates of columns for which a moderately rough condition is assumed. When compared with PARASOL (Polarization and Anisotropy of Reflectances for Atmospheric Sciences coupled with Observations from a Lidar) polarized reflection data, simulations of polarized reflectance using smooth particles show a poor fit to the measurements, whereas very rough-faceted particles provide an improved fit to the polarized reflectance. In this study a new microphysical model based on a mixture of 9 different ice crystal habits with severely roughened facets is developed. Simulated polarized reflectance using the new ice habit distribution is calculated using a vector adding-doubling radiative transfer model, and the simulations closely agree with the polarized reflectance observed by PARASOL. The new general habit mixture is also tested using a spherical albedo differences analysis, and surface roughening is found to improve the consistency of multi-angular observations. It is suggested that an ice model incorporating an ensemble of different habits with severely roughened surfaces would potentially be an adequate choice for global ice cloud retrievals.

Cole, Benjamin H.↗

Silicon Micromachined Microlens Array for THz Antennas

5 5 silicon microlens array was developed using a silicon micromachining technique for a silicon-based THz antenna array. The feature of the silicon micromachining technique enables one to microfabricate an unlimited number of microlens arrays at one time with good uniformity on a silicon wafer. This technique will resolve one of the key issues in building a THz camera, which is to integrate antennas in a detector array. The conventional approach of building single-pixel receivers and stacking them to form a multi-pixel receiver is not suited at THz because a single-pixel receiver already has difficulty fitting into mass, volume, and power budgets, especially in space applications. In this proposed technique, one has controllability on both diameter and curvature of a silicon microlens. First of all, the diameter of microlens depends on how thick photoresist one could coat and pattern. So far, the diameter of a 6- mm photoresist microlens with 400 m in height has been successfully microfabricated. Based on current researchers experiences, a diameter larger than 1-cm photoresist microlens array would be feasible. In order to control the curvature of the microlens, the following process variables could be used: 1. Amount of photoresist: It determines the curvature of the photoresist microlens. Since the photoresist lens is transferred onto the silicon substrate, it will directly control the curvature of the silicon microlens. 2. Etching selectivity between photoresist and silicon: The photoresist microlens is formed by thermal reflow. In order to transfer the exact photoresist curvature onto silicon, there needs to be etching selectivity of 1:1 between silicon and photoresist. However, by varying the etching selectivity, one could control the curvature of the silicon microlens. The figure shows the microfabricated silicon microlens 5 x5 array. The diameter of the microlens located in the center is about 2.5 mm. The measured 3-D profile of the microlens surface has a smooth curvature. The measured height of the silicon microlens is about 280 microns. In this case, the original height of the photoresist was 210 microns. The change was due to the etching selectivity of 1.33 between photoresist and silicon. The measured surface roughness of the silicon microlens shows the peak-to-peak surface roughness of less than 0.5 microns, which is adequate in THz frequency. For example, the surface roughness should be less than 7 microns at 600 GHz range. The SEM (scanning electron microscope) image of the microlens confirms the smooth surface. The beam pattern at 550 GHz shows good directivity.

Lee, Choonsup↗

Rewriting the Landform History of One of Africa's Three Largest Basins

The Kalahari Basin in southern Africa - one of the largest basins in Africa, along with the Congo and Chad basins - has attracted attention since David Livingstone traveled through the area in the 1840s. It is a semiarid desert with a large freshwater swampland known as the Okavango Swamp (150 km radius). This prominent megafan (a fan with radii >100 km), with its fingers of dark green forests projecting into the dun colors of the dunes of the Kalahari semi-desert, has been well photographed by astronauts over the years. The study area in the northern Kalahari basin is centered on the Okavango megafan of northwest Botswana, whose swampland has become well known as an African wildlife preserve of importance to biology and tourism alike. The Okavango River is unusual because it has deposited not one but two megafans along its course: the Okavango megafan and the Cubango megafan. The Okavango megafan is one of only three well-known megafans in Africa. Megafans on Earth were once thought to be rare, but recent research has documented 68 in Africa alone. Eleven megafans, plus three more candidates, have been documented in the area immediately surrounding the Okavango feature. These 11 megafans occupy the flattest and smoothest terrains adjacent to the neighboring upland and stand out as the darkest areas in the roughness map of the area. Megafan terrains occupy at least 200,000 sq km of the study area. The roughness map shown is based on an algorithm used first on Mars to quantify topographic roughness. Research of Earth's flattest terrains is just beginning with the aid of such maps, and it appears that these terrains are analogous to the flattest regions of Mars. Implications: 1. The variability in depositional style in each subbasin may apply Africa-wide: rift megafan length is dominated by rift width, whereas Owambo subbasin megafans are probably controlled by upland basin size; Zambezi subbasin megafans appear more like foreland basin types, with the position of the trunk river controlling size. 2. These perspectives were successfully applied to identify the largest megafan in the group (Cubango), a fan that was sufficiently overprinted by dunes and dry lakelets not to be detectable remotely. Such undertsanding can probably be applied on Mars, where Earth experience suggests megafans ought to exist. 3. Sweep angles of rivers on megafans drastically change the hydrology in some subbasins: when the Cubango and Kunene rivers were oriented to the Etosha Pan, it was probably a permanent water body. Now that the rivers are oriented away from the basin, 93 percent of the discharge area from the pan's northerly (main) source area is gone. 4. Biotic contact between major river systems was probably controlled by megafans situated on divides: various fish species that originated in the Congo basin are now found in the Upper Zambezi R., and vice versa, apparently because of river switching behavior on the Cassai megafan that has mediated migrations both to the south and the north.

Wilkinson, Justin↗

Low-Reynolds Number Aerodynamics of an 8.9 Percent Scale Semispan Swept Wing for Assessment of Icing Effects

Aerodynamic assessment of icing effects on swept wings is an important component of a larger effort to improve three-dimensional icing simulation capabilities. An understanding of ice-shape geometric fidelity and Reynolds and Mach number effects on the iced-wing aerodynamics is needed to guide the development and validation of ice-accretion simulation tools. To this end, wind-tunnel testing and computational flow simulations were carried out for an 8.9%-scale semispan wing based upon the Common Research Model airplane configuration. The wind-tunnel testing was conducted at the Wichita State University 7 ft x 10 ft Beech wind tunnel from Reynolds numbers of 0.8×10(exp 6) to 2.4×10(exp 6) and corresponding Mach numbers of 0.09 to 0.27. This paper presents the results of initial studies investigating the model mounting configuration, clean-wing aerodynamics and effects of artificial ice roughness. Four different model mounting configurations were considered and a circular splitter plate combined with a streamlined shroud was selected as the baseline geometry for the remainder of the experiments and computational simulations. A detailed study of the clean-wing aerodynamics and stall characteristics was made. In all cases, the flow over the outboard sections of the wing separated as the wing stalled with the inboard sections near the root maintaining attached flow. Computational flow simulations were carried out with the ONERA elsA software that solves the compressible, three-dimensional RANS equations. The computations were carried out in either fully turbulent mode or with natural transition. Better agreement between the experimental and computational results was obtained when considering computations with free transition compared to turbulent solutions. These results indicate that experimental evolution of the clean wing performance coefficients were due to the effect of three-dimensional transition location and that this must be taken into account for future data analysis. This research also confirmed that artificial ice roughness created with rapid-prototype manufacturing methods can generate aerodynamic performance effects comparable to grit roughness of equivalent size when proper care is exercised in design and installation. The conclusions of this combined experimental and computational study contributed directly to the successful implementation of follow-on test campaigns with numerous artificial ice-shape configurations for this 8.9% scale model.

Aerodynamics↗

Low-Reynolds Number Aerodynamics of an 8.9 Percent Scale Semispan Swept Wing for Assessment of Icing Effects

Aerodynamic assessment of icing effects on swept wings is an important component of a larger effort to improve three-dimensional icing simulation capabilities. An understanding of ice-shape geometric fidelity and Reynolds and Mach number effects on the iced-wing aerodynamics is needed to guide the development and validation of ice-accretion simulation tools. To this end, wind-tunnel testing and computational flow simulations were carried out for an 8.9 percent-scale semispan wing based upon the Common Research Model airplane configuration. The wind-tunnel testing was conducted at the Wichita State University 7 by 10 ft Beech wind tunnel from Reynolds numbers of 0.8×10(exp 6) to 2.4×10(exp 6) and corresponding Mach numbers of 0.09 to 0.27. This paper presents the results of initial studies investigating the model mounting configuration, clean-wing aerodynamics and effects of artificial ice roughness. Four different model mounting configurations were considered and a circular splitter plate combined with a streamlined shroud was selected as the baseline geometry for the remainder of the experiments and computational simulations. A detailed study of the clean-wing aerodynamics and stall characteristics was made. In all cases, the flow over the outboard sections of the wing separated as the wing stalled with the inboard sections near the root maintaining attached flow. Computational flow simulations were carried out with the ONERA elsA software that solves the compressible, threedimensional RANS equations. The computations were carried out in either fully turbulent mode or with natural transition. Better agreement between the experimental and computational results was obtained when considering computations with free transition compared to turbulent solutions. These results indicate that experimental evolution of the clean wing performance coefficients were due to the effect of three-dimensional transition location and that this must be taken into account for future data analysis. This research also confirmed that artificial ice roughness created with rapid-prototype manufacturing methods can generate aerodynamic performance effects comparable to grit roughness of equivalent size when proper care is exercised in design and installation. The conclusions of this combined experimental and computational study contributed directly to the successful implementation of follow-on test campaigns with numerous artificial ice-shape configurations for this 8.9 percent scale model.

Aerodynamics↗

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↗

Design and Development of a Laboratory-Scale Ice Adhesion Testing Device

When an aircraft traverses through clouds containing supercooled water droplets, in-flight icing can occur that negatively affects vehicle performance by increasing weight and drag leading to loss of lift. Super-cooled water droplets present in clouds that impact vehicle surfaces can lead to inflight icing any time during the year.1 Most events occur at temperatures ranging from 0 to -20degC. Ice generated on the aircraft can vary between clear/glaze, rime, and mixed (Fig. 1) depending on air temperature (-5 to -20degC), liquid water content (0.3-0.6 g/m3), and droplet size (median volumetric diameter of 15-40 μm). Current strategies to remove ice are based on active technologies such as pneumatic boots, heated surfaces, and deicing agents (i.e., ethylene- and propylene-based glycols). The latter have potential environmental concerns. A passive approach to mitigate accreting ice that is actively being investigated are protective coatings. An ice mitigating coating could potentially be used as a stand-alone material, but more likely in combination with an active approach. In the latter scenario, potential reduction in power consumption by the active approach may be realized. To determine the ice adhesion strength of impact ice that is representative of the aircraft environment is not a trivial matter. Test methods utilizing slowly formed ice (i.e., freezer ice) do not accurately simulate this environment. Likewise, some testing methodologies involve sample relocation from the icing environment to the test chamber that can result in thermal shock to the sample, thus affecting the results. The Adverse Environment Rotor Test Stand (AERTS) located at Pennsylvania State University (PSU) has been demonstrated to simulate impact icing conditions within the icing envelope for the determination of ice adhesion shear strength (IASS) without removal/relocation of the sample.2 Due to the confidence in results obtained from AERTS, this instrument is in high demand and requires a significant amount of lead time and capital investment to obtain IASS results. As a solution for quickly and economically screening coatings in a controlled manner under impact icing conditions, a laboratory-scale ice adhesion test and dead blades were then removed from the rotor/blade assembly to obtain the final mass. The IASS of the live blade was determined from the difference in mass (before and after testing) of the live and dead blades, the ice shed area, and the rpm of the shed event. The same live blade sample was tested in triplicate at all three test temperatures. Surface roughness was determined using a Bruker Dektak XT Stylus Profilometer. Measurements were conducted using a 12.5 μm tip at a vertical range of 65.5 μm with an applied force of 3 mg. Data were collected over a 1.0 mm length at a resolution of 0.056 μm/point. Five single line scans at different locations were collected and processed using a two-point leveling subtraction. The resultant Ra (arithmetic roughness) and Rq (root mean square roughness) average values were calculated.

Smith, Joseph G., Jr.↗

Swarm Flyby Gravimetry

This study describes a new technology for discerning the gravity fields and mass distribution of a solar system small body, without requiring dedicated orbiters or landers. Instead of a lander, a spacecraft releases a collection of small, simple probes during a flyby past an asteroid or comet. By tracking those probes from the host spacecraft, one can estimate the asteroid's gravity field and infer its underlying composition and structure. This approach offers a diverse measurement set,equivalent to planning and executing many independent and unique flyby encounters of a single spacecraft. This report assesses a feasible hardware implementation, derives the underlying models,and analyzes the performance of this concept via simulation.In terms of hardware, a small, low mass, low cost implementation is presented, which consists of a dispenser and probes. The dispenser contains roughly 12 probes in a tube and has a total size commensurate with a 6U P-Pod. The probes are housed in disc shaped sabots. When commanded,the dispenser ejects the top-most probe using a linear motor. The ejected probe separates from its sabots and unfolds using internal springs. There are two types of probes, each designed for a particular tracking modality. The reflective probe type, tracked by a telescope, unfolds to forma diffusely reflective sphere. The retroreflector probe type, tracked by a lidar, unfolds to form a corner-cube retroreflector assembly. Both types are designed to spherical so that their attitude doesn't affect the spacecraft's tracking performance.This analysis indicates that the point-mass term of small bodies larger than roughly 500 m in diameter can be observed from a host spacecraft that tracks locally deployed probes throughout a flyby to an uncertainty of better than 5%. The conditions by which this measurement is possible depends on the characteristics of the asteroid (size, type), the flyby velocity, and the type of tracking available (angles-only or angles+ranging). For most encounters, a few (1-3) well placed probes can be very effective, with marginal improvement for additional probes. Given realistic deployment errors, an encounter may require roughly 10-12 probes to ensure that 1-3 achieve their target. Long duration tracking of probes flying by large asteroids (>5 km diameter) can sometimes provide observability of the gravity field's first spherical harmonic, J( sub 2). In summary, this method offers a feasible, affordable approach to enabling or augmenting flyby science.

Atchison, Justin↗

Use of Polarimetric Radar Measurements to Constrain Simulated Convective Cell Evolution: A Pilot Study with Lagrangian Tracking

To probe the potential value of a radar-driven field campaign to constrain simulation of isolated convection subject to a strong aerosol perturbation, convective cells observed by the operational KHGX weather radar in the vicinity of Houston, Texas, are examined individually and statistically. Cells observed in a single case study of onshore flow conditions during July 2013 are first examined and compared with cells in a regional model simulation. Observed and simulated cells are objectively identified and tracked from observed or calculated positive specific differential phase (K(sub DP)) above the melting level, which is related to the presence of supercooled liquid water. Several observed and simulated cells are subjectively selected for further examination. Below the melting level, we compare sequential cross sections of retrieved and simulated raindrop size distribution parameters. Above the melting level, we examine time series of KDP and radar differential reflectivity (Z(sub DR)) statistics from observations and calculated from simulated supercooled rain properties, alongside simulated vertical wind and supercooled rain mixing ratio statistics. Results indicate that the operational weather radar measurements offer multiple constraints on the properties of simulated convective cells, with substantial value added from derived K(sub DP) and retrieved rain properties. The value of collocated three-dimensional lightning mapping array measurements, which are relatively rare in the continental US, supports the choice of Houston as a suitable location for future field studies to improve the simulation and understanding of convective updraft physics. However, rapid evolution of cells between routine volume scans motivates consideration of adaptive scan strategies or radar imaging technologies to amend operational weather radar capabilities. A 3-year climatology of isolated cell tracks, prepared using a more efficient algorithm, yields additional relevant information. Isolated cells are found within the KHGX domain on roughly 40 % of days year-round, with greatest concentration in the northwest quadrant, but roughly 5-fold more cells occur during June through September. During this enhanced occurrence period, the cells initiate following a strong diurnal cycle that peaks in the early afternoon, typically follow a south-to-north flow, and dissipate within 1 h, consistent with the case study examples. Statistics indicate that ∼ 150 isolated cells initiate and dissipate within 70 km of the KHGX radar during the enhanced occurrence period annually, and roughly 10 times as many within 200 km, suitable for multi-instrument Lagrangian observation strategies. In addition to ancillary meteorological and aerosol measurements, robust vertical wind speed retrievals would add substantial value to a radar-driven field campaign.

Fridlind, Ann M.↗

A Fresh Look at the Stratigraphy of Northern Australe

The roughly circular collection of mare deposits centered at ~38.9°S, 93°E is often re- ferred to as Mare Australe. It is located outside of the Procellarum KREEP Terrain. The circular arrangement of Australe’s mare patches has suggested an ancient, heavily degraded or relaxed impact basin roughly 900 km in diameter. The mare deposits are generally thought to have erupted into smaller post- basin craters. The type, volume, and distribution of mare eruptions potentially resembles the early stages of basin-filling mare events, but which are preserved in Australe and some farside locations. Gravity data suggest that if there was a basin, it is much smaller than originally proposed (now ~600 km) and located in the northern part of Mare Australe, between Humboldt, Milne, and Jenner craters. As a whole, Mare Australe lacks the topography typically associated with a basin; however, northern Australe has a slight topographic depression that roughly corresponds to the basin-like Bouguer gravity signature in the same area. The compositions exposed in Humboldt crater suggest that a preexisting basin might have excavated deeper crustal material. However, the underlying cause of the circularity of Mare Australe’s deposits, particularly those extending outside of the potential impact basin setting, is not yet understood. Thus, Australe may preserve fundamental information about mare volcanism potentially uncoupled from basin formation and structure. The objectives of this study are to use new high- resolution data (images, gravity, topography, and com- position) to reassess Australe’s mare deposits, deter- mine the timing and style of volcanism, identify discrete basalt deposits, and to further characterize the evolution of magmatism and subsurface structure in this area. Here, we focus on the northern Australe deposits (between Humboldt, Jenner, and Milne). As originally noted by Whitford-Stark (1979), Humboldt crater and its ejecta make an excellent stratigraphic marker that can be traced across much of the Australe region. The ejecta serves as a stratigraphic constraint for absolute model ages (AMAs) derived from crater size-frequency distributions (CSFDs).

Stopar, J. D.↗

Preface to the Special Collection: Recollections in Space Physics

Space Physics is a comparatively young scientific discipline, tightly linked to the era of satellite based investigations and the discoveries that came with it. As such, we as a community are fortunate to have met, been taught and mentored by, and even become friends with, many of those who were around to witness the birth of, and in many cases establish, our field. Each of these "Pioneers" have remarkable stories to tell. These accounts from the dawn of the space age provide a glimpse into that era of scientific discovery and are an important part of our collective history that deserve to be shared and commemorated. These are stories of perseverance, ingenuity, luck, and sometimes failure. We are fortunate to live in an era in which these distinguished scientists are still with us and able to share their experiences.To commemorate the 100th Anniversary of the American Geophysical Union (AGU), we have solicited a special collection of recollection papers to memorialize these experiences. This is not the first such effort to capture these stories. In 1994 (Vol. 99, No. A10, pp. 19,099-19,212) and again in 1996 (Vol. 101, No. A5, pp. 10,477-10,585) JGRSpace Physics published special sections entitled, "Pioneers of Space Physics," in celebration of AGU's 75th anniversary. Later, in 1997, Gilmore and Spreiter (1997) added a set of recollection papers on the discovery of the magnetosphere. Together, these papers containing personal accounts from the pioneers of the space age covered the period of roughly 1958-1967.It is our great honor to present retrospective papers from several of our distinguished colleagues who contributed significantly to the explosion in understanding of space physics and aeronomy that occurred from roughly 1967-1980. We chose to start at the end of the previous set of recollections and cover the decade of the 1970s, during which our field greatly expanded. It is also a time in which the second generation of space scientists entered the field and made foundational discoveries that continue to define our discipline.As with the previous papers, we asked the authors to "recount some of the events leading to the emergence of space physics and to put the events into a professional as well as personal perspective." (Gombosi et al.,1994). The authors of this special collection of recollections were selected following the same guidelines as the previous efforts. We started from a long list of senior colleagues and narrowed the list to roughly two dozen distinguished scientists based on discipline and geographic balance. Some of our colleagues, when asked, felt they would be unable to devote the time or energy to such a recollection and declined. Sadly,in the intervening 25 years since the original effort, some of our colleagues who were most active during the early years of the space age have passed. Therefore, this new special collection, 25 years after the first one, is timely with the AGU centennial celebrations, but late in fully capturing the stories of the pioneers of space physics. This is unfortunate and we encourage future editors of this journal to commission specialsections of legacy perspectives with a faster cadence than a quarter of a century.

Kepko, Emil L.↗

The benefit of brightness temperature assimilation for the SMAP Level-4 surface and root-zone soil moisture analysis

The Soil Moisture Active Passive (SMAP) Level-4 (L4) product provides global estimates of surface soil moisture (SSM) and root-zone soil moisture (RZSM) via the assimilation of SMAP brightness temperature (Tb) observations into the NASA Catchment Land Surface Model (CLSM). Here, using in situ measurements from 2474 sites in China, we evaluate the performance of soil moisture estimates from the L4 data assimilation (DA) system and from a baseline “open-loop” (OL) simulation of CLSM without Tb assimilation. Using random forest regression, the efficiency of the L4 DA system (i.e., the performance improvement in DA relative to OL) is attributed to eight control factors related to the CLSM as well as τ–ω radiative transfer model (RTM) components of the L4 system. Results show that the Spearman rank correlation (R) for L4 SSM with in situ measurements increases for 77 % of the in situ measurement locations (relative to that of OL), with an average R increase of approximately 14 % (ΔR=0.056). RZSM skill is improved for about 74 % of the in situ measurement locations, but the average R increase for RZSM is only 7 % (ΔR=0.034). Results further show that the SSM DA skill improvement is most strongly related to the difference between the RTM-simulated Tb and the SMAP Tb observation, followed by the error in precipitation forcing data and estimated microwave soil roughness parameter h. For the RZSM DA skill improvement, these three dominant control factors remain the same, although the importance of soil roughness exceeds that of the Tb simulation error, as the soil roughness strongly affects the ingestion of DA increments and further propagation to the subsurface. For the skill of the L4 and OL estimates themselves, the top two control factors are the precipitation error and the SSM–RZSM coupling strength error, both of which are related to the CLSM component of the L4 system. Finally, we find that the L4 system can effectively filter out errors in precipitation. Therefore, future development of the L4 system should focus on improving the characterization of the SSM–RZSM coupling strength.

Jianxiu Qiu↗

Estimating Everglades Peat Reserves and Losses on a Landscape-Scale

The Florida Everglades formed from sawgrass and other aquatic plant remains that have accumulated over thousands of years and is a patterned peatland. The region was initially drained for agricultural and urban development in the late 1800s by the lowering of Lake Okeechobee water levels through canal drainage. Over the course of the past century, the Everglades has further been altered by the construction of canals and levees to provide control of Lake Okeechobee and Everglades water levels. The drainage has caused the loss of much of the peat of the Everglades by peat fires and microbial oxidation, much of it occurring between the 1880s and the 1950s. Plans for Everglades restoration are designed to provide additional water to flow through the remaining Everglades. Additional water into the Everglades would prevent further peat loss by reducing oxidation and encouraging the accretion of peat that has been lost. This potentially could result in the sequestration of significant quantities of atmospheric carbon. We used historic records to determine how much peat was there originally and how much has been lost. Over the past century and a half, approximately 50% of the original Everglades has been lost to agricultural and urban development. In addition, much of the peat in the remaining Everglades Protection Area has been lost. We used Geographic Information Systems technologies to determine the volume and mass of peat, the mass of carbon for each of the predrainage landscapes as well as the current Everglades regions. To do this, historical and current surveys within the Everglades were used to create elevation maps. With these maps, along with a bedrock contour map of south Florida, and soils data from the USEPA, we estimated Everglades peat characteristics. These calculations provide rough, quantitative estimates of the historical Everglades, its current condition, and the changes in the peat soil that have taken place since the 19th century. Given the uncertainties in the hindcasting of some of these data sets, this analysis provides rough estimates for the peat volumes and related constituents. Our calculations indicate that the current Everglades roughly contains less than 24% of theoriginal peat volume, 17% of its mass and 19% of its carbon. The work described was fully supported by the South Florida Water Management District, West Palm Beach, FL

Thomas Walter Dreschel↗

Mars Entry, Descent, Landing, and Ascent Systems Sensitivities to Landing Site and Atmospheric Dust

Plans for human missions to Mars continue to go through several architectural changes, dating all the way back to the 1950s [1]. The continuous study, reformulation, and refinement of Mars architectures and system concepts is necessary in order to incorporate evolving mission objectives, technology advancements, and growth in the body of knowledge regarding human factors and the various environments of human space travel. This evolution has continued into the 21st century, with architectures concepts such as NASA’s Design Reference Architecture 5.0 in 2009 [2], The Evolvable Mars Campaign in 2016 [3], and as recently as 2020, an architecture focused on a crewed mission as early as the 2030s that aims to minimizing crewed duration and infrastructure investment for the first mission [4]. Within these architectures, numerous studies around the current concept designs for Mars entry, descent, landing, and ascent (EDLA) systems have been performed over the last half a decade [[5], [6], [7], [8], [9], [10]]. Despite the breadth and depth of these studies, landing site, a key design parameter relevant to the design of EDLA systems, has remained nebulous over the years, largely due to the ever evolving mission objectives and architecture concept over the decades. However, the specific landing site has direct impact on the altitude and atmospheric conditions, which subsequently impact the design of EDLA systems. To accommodate the lack of a specific landing site, a baseline reference altitude of 0 km relative to the Mars Orbiter Laser Altimeter (MOLA), which is similar to an Earth sea level reference, has typically been selected and fixed for these past studies. Similarly, a baseline reference atmosphere has typically been utilized in these studies, either the Mars Global Reference Atmospheric Model [11], or a general 1982 standard warm, high pressure atmosphere model derived from Viking lander data [12]. Fig. 1 shows the range of elevations across the surface of Mars. Current architectures are focused on latitudes greater than 30 degrees north in an effort to ensure access to frozen water ice. From the figure, it is quickly evident that elevations over the range of longitudes at or above this latitude are significantly varied between roughly -4 to +4 km MOLA. Based on these observations, initial qualitative assessments of the impacts of landing site elevation on EDLA systems were performed. Fig. 2 depicts the expected trends in EDLA system mass for variations in both land site latitude and elevation. The background coloring is a qualitative evaluation of the current knowledge on ice water availability at the latitudes. Further modeling and simulation was performed to obtain numerical predictions of sensitivities to these parameters. Results indicate up to +2% to -4% mass variation from the current baseline Mars Ascent Vehicle concept, with the potential for greater than 6% mass variation at latitudes greater than 70 degrees north. However, variations in the descent system due to landing site were much more significant with mass variations in the range of -15% to +30% around the current baseline Mars Descent System concept. Additionally, Fig. 3 provides a depiction of potential Mars atmospheric density variations with both dust and Martian season. Though the Martian atmosphere is relatively thin compared to Earth’s atmosphere, it still poses significant impact on the design of EDLA systems. Variations in atmospheric density indicated by this slice of data would have significant impact on the design of EDLA systems. However, due to the highly complex nature of atmospheric flight, a qualitative assessment could not be performed. Rather, sensitivity results relied on modeling and simulation to provide numerical results for sensitivities on the current EDLA design concepts under consideration. Understanding these sensitivities is vital to the overall systems design of a Mars architecture. The resulting mass impacts on EDLA systems, driven by landing site elevation and dust level variations, has rippling impacts throughout the architecture that, ultimately, impact the viability of the architecture. Results showed fairly minor mass impacts to the current Mars Ascent Vehicle baseline configuration, roughly -0.6% to + 0.4% mass variations, while the Mars Descent System say greater variations due to dust, roughly -4% to +6% around the baseline concept.

EDLA↗

The Harmonic Linearized Navier-Stokes Equations for Transition Prediction in Three-Dimensional Flows

The conventional method to predict the onset of laminar-turbulent transition in convectively unstable boundary-layer flows is based on the logarithmic amplification ratio, the so-called N-factor, of the linear instability waves. To calculate the N-factor, the flow variables are decomposed into a laminar basic state solution and the linear disturbances, which are assumed to be harmonic in time. The most commonly used linear stability analysis approaches include the locally parallel linear stability theory (LST) and the non-local, weakly nonparallel parabolized stability equations (PSE). However, these methods do not account for strong streamwise gradients that are encountered in several configurations of interest, as roughness elements, steps, gaps, or corners. To solve the linear evolution of disturbances along such strongly nonparallel regions, the harmonic linearized Navier-Stokes equations (HLNSE) need to be solved. The discretization of the HLNSE for spanwise/azimuthally inhomogeneous laminar basic states yields a linear system of complex arithmetic with a leading dimension of the order of 107 to 108. A combined multithread and multiprocessor algorithm is implemented for the direct solution of such linear system. Results for a supersonic boundary layer over a three-dimensional roughness patch show good agreement with experimental measurements when the evolution of the instability waves over the roughness patch is included via the HLNSE.

Boundary Layer Stability↗