Search NASA⌕ Search

SEARCH · Search NASA

Results for “PLANETARY”

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 235 records · Page 13

Estimating the Ultraviolet Emission of M Dwarfs with Exoplanets from Ca ii and Hα

M dwarf stars are excellent candidates around which to search for exoplanets, including temperate, Earth-sized planets. To evaluate the photochemistry of the planetary atmosphere, it is essential to characterize the UV spectral energy distribution of the planet’s host star. This wavelength regime is important because molecules in the planetary atmosphere such as oxygen and ozone have highly wavelength-dependent absorption cross sections that peak in the UV (900–3200 Å). We seek to provide a broadly applicable method of estimating the UV emission of an M dwarf, without direct UV data, by identifying a relationship between noncontemporaneous optical and UV observations. Our work uses the largest sample of M dwarf star far- and near-UV observations yet assembled. We evaluate three commonly observed optical chromospheric activity indices—Hα equivalent widths and log{sub 10} L Hα /L bol , and the Mount Wilson Ca ii H&K S and R HK ' indices—using optical spectra from the HARPS, UVES, and HIRES archives and new HIRES spectra. Archival and new Hubble Space Telescope COS and STIS spectra are used to measure line fluxes for the brightest chromospheric and transition region emission lines between 1200 and 2800 Å. Our results show a correlation between UV emission-line luminosity normalized to the stellar bolometric luminosity and Ca ii R HK ' with standard deviations of 0.31–0.61 dex (factors of ~2–4) about the best-fit lines. We also find correlations between normalized UV line luminosity and Hα log 10 L Hα /L bol and the S index. These relationships allow one to estimate the average UV emission from M0 to M9 dwarfs when UV data are not available.

79 ASTRONOMY AND ASTROPHYSICS↗

Rare Event Sampling Improves Mercury Instability Statistics

Due to the chaotic nature of planetary dynamics, there is a non-zero probability that Mercury's orbit will become unstable in the future. Previous efforts have estimated the probability of this happening between 3 and 5 billion years in the future using a large number of direct numerical simulations with an N-body code, but were not able to obtain accurate estimates before 3 billion years in the future because Mercury instability events are too rare. In this paper we use a new rare-event sampling technique, Quantile Diffusion Monte Carlo (QDMC), to estimate that the probability of a Mercury instability event in the next 2 billion years is approximately 10 –4 in the REBOUND N-body code. We show that QDMC provides unbiased probability estimates at a computational cost of up to 100 times less than direct numerical simulation. QDMC is easy to implement and could be applied to many problems in planetary dynamics in which it is necessary to estimate the probability of a rare event.

79 ASTRONOMY AND ASTROPHYSICS↗

Experimental investigation of elemental and isotopic evaporation processes by laser heating in an aerodynamic levitation furnace

In this report we carried out evaporation experiments on a B-type calcium–aluminium-rich inclusion (CAI) melt in a gas-mixing aerodynamic levitation laser furnace, at 1873 K and an oxygen partial pressure of 10 -9.1 atm, for durations ranging from 60 to 600 s. Evaporation of SiO 2 and MgO follow the same trend as those observed in vacuum furnace experiments at the same temperature and starting composition, showing that their evaporation relative to one another from the melt is independent of pressure, oxygen fugacity, and hydrodynamical regime specific to the furnace. Isotopic ratios of Mg and Si in evaporation residues are used to derive fractionation factors of α 26/24 Mg vap-liq = 0.9906± 0.0004 and α 30/28 Si vap-lip = 0.9943±0.0003, which are both significantly closer to unity than those found for evaporation in a vacuum, which translates to less isotope fractionation. The residues are also less isotopically fractionated than expected for cases in which transport of the gas species away from the melt is diffusion-controlled at 1-atm. By analysing the flow regimes in our furnace, we find that advection by the levitating gas is the primary mode of mass transport away from the melt surface, as opposed to diffusion-limited transport in a vacuum or 1-atm tube furnace. A modified Hertz–Knudsen–Langmuir formulation accounts for this process, and shows that isotopic fractionation of both Si and Mg reflect a saturation factor (ratio of the pressure of the evaporating species to vapour saturation pressure) equal to 0.75. This is in perfect accord with recent measurements of Cu isotopic fractionation using a similar furnace. The fact that three elements (Mg, Si, Cu) with varying equilibrium vapour pressures, activity coefficients in the liquid, and diffusion coefficients in the gas have the same scaling behaviour to saturation pressure is a strong indication that the mechanism controlling evaporation is driven by the hydrodynamical regime imposed in the furnace. Therefore, this class of experiments can be used to constrain processes in which advection dominates over diffusion, such as (but not limited to) planetary ejecta, tektites, giant impacts, nebular condensation in a turbulent flow, or nuclear fallout material. Finally, the possibility to reach high temperatures (in excess of 3500 K) in this furnace allows it to be used to evaluate the activity coefficients of melt components in extreme conditions relevant to molten planetary interiors (i.e., magma oceans), with a specific focus on refractory elements.

58 GEOSCIENCES↗

Lattice dynamics, sound velocities, and atomic environments of szomolnokite at high pressure

Complex mixtures of sulfates, silicates, and ice have been observed in a variety of planetary environments on Earth, Mars and the icy satellites of the solar system. Characterizing the properties of the corresponding compositional endmembers is important for understanding the interiors of a range of planetary bodies in which these phases are observed. To measure the electronic and vibrational properties of the pure ferrous iron endmember of the kieserite group, szomolnokite, (FeSO 4 ∙H 2 O), we have performed synchrotron 57 Fe nuclear resonant inelastic and forward scattering experiments in the diamond-anvil cell up to 14.5 GPa. This pressure range covers depths within Earth’s interior relevant to sulfur cycling in subduction zones and the range of pressures expected within icy satellite interiors. We find evidence of crystal lattice softening, changes in elastic properties, and changes in the electric field gradients of iron atoms associated with two structural transitions occurring within the experimental pressure range. Here we apply these findings to icy satellite interiors, including discussion of elastic properties, modeling of ice-sulfate aggregates, and implications for tidal observations.

58 GEOSCIENCES↗

SERENA: Particle Instrument Suite for Determining the Sun-Mercury Interaction from BepiColombo

The ESA-JAXA BepiColombo mission to Mercury will provide simultaneous measurements from two spacecraft, offering an unprecedented opportunity to investigate magnetospheric and exospheric particle dynamics at Mercury as well as their interactions with solar wind, solar radiation, and interplanetary dust. The particle instrument suite SERENA (Search for Exospheric Refilling and Emitted Natural Abundances) is flying in space on-board the BepiColombo Mercury Planetary Orbiter (MPO) and is the only instrument for ion and neutral particle detection aboard the MPO. It comprises four independent sensors: ELENA for neutral particle flow detection, Strofio for neutral gas detection, PICAM for planetary ions observations, and MIPA, mostly for solar wind ion measurements. SERENA is managed by a System Control Unit located inside the ELENA box. In the present paper the scientific goals of this suite are described, and then the four units are detailed, as well as their major features and calibration results. Finally, the SERENA operational activities are shown during the orbital path around Mercury, with also some reference to the activities planned during the long cruise phase.

79 ASTRONOMY AND ASTROPHYSICS↗

Crystal chemistry at high pressure

The chemistry we are taught in school, and we experience in our daily existence occurs at 1 atm. However, pressure spans an astounding 62 orders of magnitude in the Universe in going from the void of interstellar space to the crushing conditions at the center of a neutron star. The way in which pressure affects chemistry is important for Earth and planetary sciences, materials science, in understanding the extreme conditions experienced in nuclear explosions, and it may be key in addressing the future energy needs of our society. In this article we outline how the often neglected pressure variable affects chemistry, beginning from the way in which atomic energy levels are altered. This lays the foundation for understanding the unique crystal and electronic structures that emerge when matter is squeezed, as well as pressure’s effect on chemical reactivity. Finally, we give an overview of the main concepts behind conventional, or phonon-mediated, superconductivity, and describe how the pressure variable may be key in discovering and designing light-element based materials whose superconducting critical temperatures approach room temperature. Here, we discuss some of the main families of superconducting hydrides that have been predicted computationally, and the experimental successes in this exciting and rapidly developing field.

Chemical bonding↗

Atmospheric humidity deficits tell us how soil moisture deficits down-regulate ecosystem evaporation

Knowing how actual evaporation is down-regulated from potential evaporation during periods with soil moisture deficits is one of the greatest challenges towards computing evaporation everywhere on a regular basis. We propose the hypothesis that vegetated landscapes transmit information on soil moisture deficits through its down-regulation of evaporation, which in turn affects the humidification and growth of the planetary boundary layer. To test this hypothesis, here, we examined how the evaporative fraction, defined as the ratio between actual and potential evaporation, corresponds with the soil moisture stress index, defined as the relative humidity with a power law dependence on vapor pressure deficit times a parameter, β. We tested the parameterized soil moisture stress index with direct eddy covariance measurements of actual evaporation and computations of potential evaporation based on meteorological conditions, averaged on monthly time steps. The analysis was conducted on a dataset obtained from 144 FLUXNET sites. These sites spanned much of the world's climates and biomes and contained over 5900 months of observations. Observations of the evaporative fraction fit the model of the soil moisture stress index best for semi-arid and arid ecosystems, and least for the humid and wet ecosystems. Consequently, a significant relationship between measurements of the evaporative fraction and soil moisture stress index held for about one-half of the population of sites. Under this condition, the median value of the parameter, β, was 1.41. To investigate the mechanism of this empirical soil moisture stress index, we diagnosed it with a coupled evaporation-planetary boundary layer model. The soil moisture stress index is strongly related to surface resistance, as defined by inverting Penman-Monteith equation. Consequently, this index provides an independent estimate of surface resistance based on easy to measure mean monthly weather conditions like relative humidity and temperature. Thereby, this soil moisture stress index has potential to be applied to weather, climate and biogeochemical models and the interpretation of satellite derived evaporation products, like the one provided by the ECOSTRESS mission.

54 ENVIRONMENTAL SCIENCES↗

Isotopic signatures as tools to reconstruct the primordial architecture of the Solar System

The significant reorganization of the early Solar System due to giant planet migration has hampered our understanding of where planetary bodies formed. Previously employed proxies for reconstructing the primordial planetary architecture, such as water content or oxidation state, are complicated by post-accretionary processes. Furthermore we investigate basaltic achondrites for their nucleosynthetic isotope signatures in the elements neodymium (Nd) and zirconium (Zr) and show that they are—similar to previously investigated chondritic meteorites—characterized by a relative deficit in isotopes produced by the s-process of nucleosynthesis. Importantly, these data are well correlated with nucleosynthetic signatures observed in other elements, demonstrating that s-process matter was heterogeneously distributed throughout the early Solar System. By comparing these isotopic signatures with potential proxies for Solar System reconstruction and computer modeling, we here argue that this isotopic heterogeneity in bulk meteoritic materials is linked to the original heliocentric distance of formation. Such scaling of nucleosynthetic signatures with heliocentric distance could permit reconstruction of the primordial architecture of the Solar System by ‘cosmolocating’ the accretion orbits of meteoritic parent bodies as a function of incorporated s-process matter.

58 GEOSCIENCES↗

A window into atmospheric escape on early Mars provided by argon isotopes in a Martian meteorite

Novel analytical approaches to determining the trapped 38 Ar/ 36 Ar ratio of gases contained within Martian meteorites are presented and applied to the Martian regolith breccia Northwest Africa (NWA) 7034 and paired stone NWA 11220. The resulting data indicate that extensive mass-dependent fractionation of atmospheric Ar may have occurred within 150 million years of planetary formation, ostensibly as a result of diffusion-limited hydrodynamic escape. The inferred fractional loss of Ar and lighter atmospheric constituents exceeds 50%. In conclusion, these data suggest that volatiles derived from planetary outgassing and/or impactors may dominate the present abundance of atmospheric Ar.

07 ISOTOPE AND RADIATION SOURCES↗

Titanium and chromium isotopic compositions of calcium-aluminum-rich inclusions: Implications for the sources of isotopic anomalies and the formation of distinct isotopic reservoirs in the early Solar System

As the earliest-dated solids in our Solar System, calcium-aluminum-rich inclusions (CAIs) provide a record of their formation environment near the young Sun and hold clues to the formation of planetary-scale isotopic reservoirs in the solar protoplanetary disk. Although CAIs from several CV, CK, CM, CO, and ordinary chondrites have been analyzed previously for their Ti isotopic compositions, CAIs from just three CV chondrites have been analyzed for their Cr isotopic compositions, and only a handful of CAIs have been measured for both their Ti and Cr isotopic compositions. We report mass-independent Ti and Cr isotopic anomalies in several CAIs from CV and CK chondrites; this is the first report of the Cr isotopic composition of a CAI from a CK chondrite. With this data, we aim to better constrain the compositional range of CAIs in ε 50 Ti versus ε 54 Cr space, thereby facilitating the isotopic characterization of the material inherited by the solar protoplanetary disk and the role of CAIs in the formation of distinct planetary-scale isotopic reservoirs in our early Solar System. The narrow range in isotopic anomalies in CAIs when compared to other early-formed refractory inclusions such as platy hibonite crystals (PLACs) and spinel-hibonite inclusions (SHIBs) suggests that CAIs record the mixing of these precursor materials and the averaging of their larger isotopic anomalies. The isotopic composition of CAIs is therefore likely the result of a combination of factors, including mixing of material inherited from their formation region, heterogeneous carrier phase distribution, and thermal processing in the disk. The ε 50 Ti and ε 54 Cr isotopic compositions of CAIs are not correlated, further demonstrating that these isotopic anomalies have different carrier phases. The Ti and Cr isotopic compositions of CAIs additionally show that CAIs alone cannot be responsible for the compositional difference between the non-carbonaceous chondritic (NC) and carbonaceous chondritic (CC) isotopic reservoirs but nevertheless do play a role in the formation of these large-scale isotopic reservoirs in the early Solar System.

58 GEOSCIENCES↗

Development towards stable chlorine isotope measurements of astromaterials using the modified Middleton source of an accelerator mass spectrometer

The volatile element Cl can be lost during the formation and evolution of planetary bodies, leading to fractionation of its two stable isotopes 35 Cl and 37 Cl. Chlorine isotope variations (reported as δ 37 Cl in parts per thousand (‰) relative to Standard Mean Ocean Chloride, SMOC) are documented to exceed 80‰ between different lunar rock samples and have been variably interpreted as the fingerprint of degassing during accretion, magma ocean, or volcanic portions of lunar history. The large intersample and intrasample variations observed by both bulk isotope ratio mass spectrometry (IRMS) and in-situ secondary ion mass spectrometry (SIMS) methods are difficult to interpret in part because of a paucity of bulk Cl isotope measurements. This lack of high-precision bulk data is due to the relative rarity of IRMS laboratories capable of making these high precision measurements on small samples of precious planetary materials such as those returned by human or robotic exploration. Here we present a new method for performing high precision δ 37 Cl measurements using the modified Middleton ion source of an existing accelerator mass spectrometer. For samples with as little as 1 μg Cl–the equivalent of 2–4 mm 3 of a typical lunar rock sample–the average cathode accuracy is ~1‰. Finally, cathode reproducibility is typically ~1‰ (2σ) for samples with at least 10 μg of Cl, increasing to ~3–6‰ for aliquots with ~1–2 μg Cl, similar to published SIMS results and sufficient to study astromaterials from the Moon, Mars, or 4 Vesta, which have tens of ‰ observed variations.

07 ISOTOPE AND RADIATION SOURCES↗

Geostrophic to magnetically-damped turbulence in liquid metal rotating magnetoconvection

Understanding planetary core convection dynamics requires the study of convective flows in which the Coriolis and Lorentz forces attain a leading-order, so-called magnetostrophic balance. Experimental investigations of rotating magnetoconvection (RMC) in the magnetostrophic regime are therefore essential to broadly characterize the properties of local-scale planetary core flow. Toward this end, we present here the first thermovelocimetric measurements of magnetostrophic, liquid metal convection, which are made using liquid gallium as the working fluid, at moderate rotation rates (Ekman numbers 10 -4 ≤ Ek ≤ 10 -5 ) and in the presence of dynamically strong magnetic fields (Elsasser number Λ = 1). Complementary rotating convection (RC) experiments are performed at the same rotation rates to serve as reference cases. Our RMC velocity measurements adequately follow a geostrophic turbulent scaling for cases in which local-scale convective inertial forces exceed the Lorentz forces in the fluid bulk. In cases where Lorentz forces exceed local-scale inertia (local or convective interaction parameter N ℓ or N C ≳ 3), the root-mean-square RMC velocities are magnetically damped, yielding values below the geostrophic turbulent RC scaling prediction. An enhancement in heat transfer is observed, which we attribute to the increased coherence of vertically aligned magnetostrophic convective flow. Finally, our results and comparisons to dynamo simulations suggest that the system is in the strong-field regime, where convection-scale flows are magnetically damped with N C ≳ 3).

Convective turbulence↗

Advancing Aerosol Chemical Characterization and Vertical Profiling over the Southern Great Plains Using Uncrewed Aerial Sampling and Offline Aerosol Mass Spectrometry

Recent advancements in uncrewed aerial systems (UASs) and particulate matter (PM) analytical techniques have provided opportunities for atmospheric research. In this study, we deployed the Department of Energy’s fixed-wing ArcticShark UAS to examine PM 2.5 composition at varying altitudes─within and above the planetary boundary layer (PBL)─over the Southern Great Plains atmospheric observatory (SGP). A total of 22 flights were conducted across March, June, and August 2023. Composite filter samples were collected during each flight and analyzed with offline aerosol mass spectrometry (AMS), complemented by on-board real-time sensors and ground-based instrumentation, to provide a comprehensive view of regional aerosol characteristics. Results show clear vertical and seasonal differences in the aerosol composition. Relative to ground-level measurements, aloft samples exhibited shifts in the distribution of organic and inorganic PM, with the organic composition varying distinctly across seasons. Particulate organic nitrogen (ON) was elevated, with bulk compositions similar in March and June but strongly altered in August, likely driven by biomass burning and enhanced photochemical activity. Combined AMS and chemical ionization mass spectrometry analyses detected amines, amides, and amino acids. PM above the planetary boundary layer was enriched in oxidized organic aerosols, while ground-level PM contained higher nitrate and sulfate. Seasonal differences in aqueous-phase processing were also observed, which were strongest in March during persistent cloud cover and weaker in the drier August period, suggesting a shift from aqueous- to gas-phase SOA formation. In conclusion, these findings highlight the value of UAS in advancing PM measurements and vertical profiling of aerosol composition.

54 ENVIRONMENTAL SCIENCES↗

Liquid-Vapor Coexistence and Critical Point of Mg 2 SiO 4 From Ab Initio Simulations

Hypervelocity impact-driven vaporization is characteristic of late-stage planet formation. Yet the behavior and properties of liquid-vapor mixtures of planetary materials of interest are typically unknown. Multiphase equations of state used in hydrodynamic simulations of planet impacts therefore lack reliable data for this important phenomenon. Here, we present the first constraints on the liquid-vapor critical point and coexistence phase boundary of Mg2SiO4 computed from ab initio molecular dynamics simulations. We found that the vapor is depleted in magnesium and enriched in silica and oxygen, while the coexisting liquid is enriched in magnesium and depleted in oxygen, from which we infer vaporization is incongruent. The critical point was estimated from an equation of state fit to the data. The results are in line with recent calculations of MgSiO3 and together confirm that extant multiphase equation of state (EOS) models used in planetary accretion modeling significantly underestimate the amount of supercritical material postimpact.

58 GEOSCIENCES↗

Temperature and Density on the Forsterite Liquid‐Vapor Phase Boundary

Abstract The physical processes during planet formation span a large range of pressures and temperatures. Giant impacts, such as the one that formed the Moon, achieve peak pressures of 100s of GPa. The peak shock states generate sufficient entropy such that subsequent decompression to low pressures intersects the liquid‐vapor phase boundary. The entire shock‐and‐release thermodynamic path must be calculated accurately in order to predict the post‐impact structures of planetary bodies. Forsterite (Mg 2 SiO 4 ) is a commonly used mineral to represent the mantles of differentiated bodies in hydrocode models of planetary collisions. Here, we performed shock experiments on the Sandia Z Machine to obtain the density and temperature of the liquid branch of the liquid‐vapor phase boundary of forsterite. This work is combined with previous work constraining pressure, density, temperature, and entropy of the forsterite principal Hugoniot. We find that the vapor curves in previous forsterite equation of state models used in giant impacts vary substantially from our experimental results, and we compare our results to a recently updated equation of state. We have also found that due to under‐predicted entropy production on the principal Hugoniot and elevated temperatures of the liquid vapor phase boundary of these past models, past impact studies may have underestimated vapor production. Furthermore, our results provide experimental support to the idea that giant impacts can transform much of the mantles of rocky planets into supercritical fluids.

58 GEOSCIENCES↗

The Lithophile Element Budget of Earth's Core

Abstract The relative composition of Earth's core and mantle were set during core formation. By determining how elements partition between metal and silicate at high pressures and temperatures, measurements of the mantle composition and geophysical observations of the core can be used to understand the mechanisms by which Earth formed. Here we present the results of metal‐silicate partitioning experiments for a range of nominally lithophile elements (Al, Ca, K, Mg, O, Si, Th, and U) and S to 85 GPa and up to 5400 K. With our results and a compilation of literature data, we developed a parameterization for partitioning that accounts for compositional dependencies in both the metal and silicate phases. Using this parameterization in a range of planetary growth models, we find that, in general, lithophile element partitioning into the metallic phase is enhanced at high temperatures. The relative abundances of FeO, SiO 2 , and MgO in the mantle vary significantly between planetary growth models, and the mantle abundances of these elements can be used to provide important constraints on Earth's accretion. To match Earth's core mass and mantle composition, Earth's building blocks must have been enriched in Fe and depleted in Si compared with CI chondrites. Finally, too little Mg, Si, and O are partitioned into the core for precipitation of oxides to be a major source of energy for the geodynamo. In contrast, several ppb of U can be partitioned into the core at high temperatures, and this energy source must be accounted for in thermal evolution models.

Chidester, B. A.↗

The Effect of Surface Heating Heterogeneity on Boundary Layer Height and Its Dependence on Background Wind Speed

The planetary boundary layer height (PBLH) is a fundamental variable of the planetary boundary layer (PBL). The effect of surface heterogeneity is challenging in the PBL parameterization for numerical weather and climate models. Here we use large eddy simulation data to investigate the impact of surface heating heterogeneity on the PBLH spatial variations and mean values over a domain of typical mesoscale model grid size. It is found that this impact depends on background wind condition. Variance decomposition and Fourier spectra both reveal that for the fully developed sheared convective PBL, surface heterogeneity contributes little to the spatial variation of PBLH, except for the case with heterogeneity scale of 14.4 km. This suggests that surface heterogeneity with length scales less than ~10 km, which is typical for mesoscale modeling, does not induce significant subgrid spatial variation of PBLH. The domain-averaged PBLH is found to generally increase as surface heterogeneity scale decreases, until reaching the scale on the order of PBLH. The findings can be used to inform the parameterization development of PBL processes over heterogeneous surface.

54 ENVIRONMENTAL SCIENCES↗

A Scale‐Dependent Analysis of the Barotropic Vorticity Budget in a Global Ocean Simulation

Abstract The climatological mean barotropic vorticity budget is analyzed to investigate the relative importance of surface wind stress, topography, planetary vorticity advection, and nonlinear advection in dynamical balances in a global ocean simulation. In addition to a pronounced regional variability in vorticity balances, the relative magnitudes of vorticity budget terms strongly depend on the length‐scale of interest. To carry out a length‐scale dependent vorticity analysis in different ocean basins, vorticity budget terms are spatially coarse‐grained. At length‐scales greater than 1,000 km, the dynamics closely follow the Topographic‐Sverdrup balance in which bottom pressure torque, surface wind stress curl and planetary vorticity advection terms are in balance. In contrast, when including all length‐scales resolved by the model, bottom pressure torque and nonlinear advection terms dominate the vorticity budget (Topographic‐Nonlinear balance), which suggests a prominent role of oceanic eddies, which are of km in size, and the associated bottom pressure anomalies in local vorticity balances at length‐scales smaller than 1,000 km. Overall, there is a transition from the Topographic‐Nonlinear regime at scales smaller than 1,000 km to the Topographic‐Sverdrup regime at length‐scales greater than 1,000 km. These dynamical balances hold across all ocean basins; however, interpretations of the dominant vorticity balances depend on the level of spatial filtering or the effective model resolution. On the other hand, the contribution of bottom and lateral friction terms in the barotropic vorticity budget remains small and is significant only near sea‐land boundaries, where bottom stress and horizontal viscous friction generally peak.

Khatri, Hemant↗