• In this work, an efficient numerical scheme is presented for seismic blind deconvolution in a multichannel scenario. The proposed method iterate with wo steps: first, wavelet estimation across all channels and second, refinement of the reflectivity estimate simultaneously in all channels using sparse deconvolution. The reflectivity update step is formulated as a basis pursuit denoising problem and a sparse solution is obtained with the spectral projected-gradient algorithm - faithfulness to the recorded traces is constrained by the measured noise level. Wavelet re-estimation has a closed form solution when performed in the frequency domain by finding the minimum energy wavelet common to all channels. Nothing is assumed known about the wavelet apart from its time duration. In tests with both synthetic and real data, the method yields sparse reflectivity series and stable wavelet estimates results compared to existing methods with significantly less computational effort.
  • Probabilistic inversion within a multiple-point statistics framework is often computationally prohibitive for high-dimensional problems. To partly address this, we introduce and evaluate a new training-image based inversion approach for complex geologic media. Our approach relies on a deep neural network of the generative adversarial network (GAN) type. After training using a training image (TI), our proposed spatial GAN (SGAN) can quickly generate 2D and 3D unconditional realizations. A key characteristic of our SGAN is that it defines a (very) low-dimensional parameterization, thereby allowing for efficient probabilistic inversion using state-of-the-art Markov chain Monte Carlo (MCMC) methods. In addition, available direct conditioning data can be incorporated within the inversion. Several 2D and 3D categorical TIs are first used to analyze the performance of our SGAN for unconditional geostatistical simulation. Training our deep network can take several hours. After training, realizations containing a few millions of pixels/voxels can be produced in a matter of seconds. This makes it especially useful for simulating many thousands of realizations (e.g., for MCMC inversion) as the relative cost of the training per realization diminishes with the considered number of realizations. Synthetic inversion case studies involving 2D steady-state flow and 3D transient hydraulic tomography with and without direct conditioning data are used to illustrate the effectiveness of our proposed SGAN-based inversion. For the 2D case, the inversion rapidly explores the posterior model distribution. For the 3D case, the inversion recovers model realizations that fit the data close to the target level and visually resemble the true model well.
  • In this paper, we examine the applicability of the approximation, $\overline{f\,g}\approx \overline f\,\overline g\,$, within Backus (1962) averaging. This approximation is a crucial step in the method proposed by Backus (1962), which is widely used in studying wave propagation in layered Hookean solids. According to this approximation, the average of the product of a rapidly varying function and a slowly varying function is approximately equal to the product of the averages of those two functions. Considering that the rapidly varying function represents the mechanical properties of layers, we express it as a step function. The slowly varying function is continuous, since it represents the components of the stress or strain tensors. In this paper, beyond the upper bound of the error for that approximation, which is formulated by Bos et al. (2017), we provide a statistical analysis of the approximation by allowing the function values to be sampled from general distributions. Even though, according to the upper bound, Backus (1962) averaging might not appear as a viable approach, we show that$-$for cases representative of physical scenarios modelled by such an averaging$-$the approximation is typically quite good. We identify the cases for which there can be a deterioration in its efficacy. In particular, we examine a special case for which the approximation results in spurious values. However, such a case$-$though physically realizable$-$is not likely to appear in seismology, where Backus (1962) averaging is commonly used. Yet, such values might occur in material sciences, in general, for which Backus (1962) averaging is also considered.
  • A universal particle velocity based algorithm for simulating hydraulic fracture with leak-off, previously demonstrated for the PKN and KGD models, is extended to obtain solutions for a penny-shaped crack. The numerical scheme is capable of dealing with both the viscosity and toughness dominated regimes, with the fracture being driven by a power-law fluid. The computational approach utilizes two dependent variables; the fracture aperture and the reduced particle velocity. The latter allows for the application of a local condition of the Stefan type (the speed equation) to trace the fracture front. The obtained numerical solutions are carefully tested using various methods, and are shown to achieve a high level of accuracy. Simple, accurate, semi-analytical approximations of the solution are provided for the zero leak-off case. A comparison with other results available in the literature is undertaken.
  • A multi-step steady-state (MSSS) method is proposed here for the measurement of low permeability. This new method can accurately and easily measure very low permeabilities of rock samples using a new setup, where the targeted rock sample and ordinary apparatus components are connected with a reference rock sample to form a series circuit. Any conventional rock sample with high permeability could be used as a reference rock sample such that the traditional steady-state measurement is feasible to accurately determine its permeability as a reference value in the MSSS method. The challenging measurement of tiny mass flux rate by advanced pump system is avoided and the permeability of targeted rock sample can be directly computed using the pressure drops, sectional areas and lengths of the two connected rock samples, and the known permeability of the reference rock sample, based on the mass conservation principle in a series circuit at steady state. Multi-step measurements using additional reference rock samples will be needed if the pressure drop across the first reference rock sample is too small to be accurately measured due to high permeability ratio when it is connected with the targeted rock sample to form a series circuit. The relative pressure drops can be small since the measurement of flow speed is unnecessary, which improves the accuracy in studying the dependence of gas permeability on the pore pressure. Consequently, the advantages of the MSSS method include low expense, simplicity, high accuracy and efficiency.
  • A key interest in geomorphology is to predict how the shear stress $\tau$ exerted by a turbulent flow of air or liquid onto an erodible sediment bed affects the transport load $M\tilde g$ (i.e., the submerged weight of transported nonsuspended sediment per unit area) and its average velocity when exceeding the sediment transport threshold $\tau_t$. Most transport rate predictions in the literature are based on the scaling $M\tilde g\propto\tau-\tau_t$, the physical origin of which, however, has remained controversial. Here we test the universality and study the origin of this scaling law using particle-scale simulations of nonsuspended sediment transport driven by a large range of Newtonian fluids. We find that the scaling coefficient is a universal approximate constant and can be understood as an inverse granular friction coefficient (i.e., the ratio between granular shear stress and normal-bed pressure) evaluated at the base of the transport layer (i.e., the effective elevation of energetic particle-bed rebounds). Usually, the granular flow at this base is gaslike and rapidly turns into the solidlike granular bed underneath: a liquidlike regime does not necessarily exist, which is accentuated by a nonlocal granular flow rheology in both the transport layer and bed. Hence, this transition fundamentally differs from the solid-liquid transition (i.e., yielding) in dense granular flows even though both transitions are described by a friction law. Combining this result with recent insights into the nature of $\tau_t$, we conclude that the transport load scaling is a signature of a steady rebound state and unrelated to entrainment of bed sediment.
  • Effects of non-stationarity on the performance of hybrid ensemble filters are studied (by hybrid filters we mean those which blend ensemble covariances with some other regularizing covariances). To isolate effects of non-stationarity from effects due to nonlinearity (and the non-Gaussianity it causes), a new doubly stochastic advection-diffusion-decay model (DSADM) is proposed. The model is hierarchical: it is a linear stochastic partial differential equation whose coefficients are random fields defined through their own stochastic partial differential equations. DSADM generates conditionally Gaussian spatiotemporal random fields with a tunable degree of non-stationarity in space and time. DSADM allows the use of the exact Kalman filter as a baseline benchmark. In numerical experiments with DSADM as the "model of truth", the relative importance of the three kinds of covariance blending is studied: with static, time-smoothed, and space-smoothed covariances. It is shown that the stronger the non-stationarity, the less useful the static covariance matrix becomes and the more beneficial the time-smoothed covariances are. Time-smoothing of background-error covariances proved to be systematically more useful than their space-smoothing. Under non-stationarity, a filter that extends the (previously proposed by the authors) Hierarchical Bayes Ensemble Filter and accommodates the three covariance-blending techniques is shown to outperform all other configurations of the filters tested. The R code of the model and the filters is available from github.com/cyrulnic/NoStRa.
  • We compute the thermal conductivity and electrical resistivity of solid hcp Fe to pressures and temperatures of Earth's core. We find significant contributions from electron-electron scattering, usually neglected at high temperatures in transition metals. Our calculations show a quasi-linear relation between electrical resistivity and temperature for hcp Fe at extreme high pressures. We obtain thermal and electrical conductivities that are consistent with experiments considering reasonable error. The predicted thermal conductivity is reduced from previous estimates that neglect electron-electron scattering. Our estimated thermal conductivity for the outer core is 77$\pm$10 W/m/K, and is consistent with a geodynamo driven by thermal convection.
  • We study self-similarity in random binary rooted trees. In a well-understood case of Galton-Watson trees, a distribution on a space of trees is said to be self-similar if it is invariant with respect to the operation of pruning, which cuts the tree leaves. This only happens for the critical Galton-Watson tree (a constant process progeny), which also exhibits other special symmetries. We extend the prune-invariance setup to arbitrary binary trees with edge lengths. In this general case the class of self-similar processes becomes much richer and covers a variety of practically important situations. The main result is construction of the hierarchical branching processes that satisfy various self-similarity definitions (including mean self-similarity and self-similarity in edge-lengths) depending on the process parameters. Taking the limit of averaged stochastic dynamics, as the number of trajectories increases, we obtain a deterministic system of differential equations that describes the process evolution. This system is used to establish a phase transition that separates fading and explosive behavior of the average process progeny. We describe a class of critical Tokunaga processes that happen at the phase transition boundary. They enjoy multiple additional symmetries and include the celebrated critical binary Galton-Watson tree with independent exponential edge length as a special case. Finally, we discuss a duality between trees and continuous functions, and introduce a class of extreme-invariant processes, constructed as the Harris paths of a self-similar hierarchical branching process, whose local minima has the same (linearly scaled) distribution as the original process.
  • Using particle-scale simulations of non-suspended sediment transport for a large range of Newtonian fluids driving transport, including air and water, we determine the bulk transport cessation threshold $\Theta^r_t$ by extrapolating the transport load as a function of the dimensionless fluid shear stress (`Shields number') $\Theta$ to the vanishing transport limit. In this limit, the simulated steady states of continuous transport can be described by simple analytical model equations relating the average transport layer properties to the law of the wall flow velocity profile. We use this model to calculate $\Theta^r_t$ for arbitrary environments and derive a general Shields-like threshold diagram in which a Stokes-like number replaces the particle Reynolds number. Despite the simplicity of our hydrodynamic description, the predicted cessation threshold, both from the simulations and analytical model, quantitatively agrees with measurements for transport in air and viscous and turbulent liquids despite not being fitted to these measurements. We interpret the analytical model as a description of a continuous rebound motion of transported particles and thus $\Theta^r_t$ as the minimal fluid shear stress needed to compensate the average energy loss of transported particles during an average rebound at the bed surface. This interpretation, supported by simulations near $\Theta^r_t$, implies that entrainment mechanisms are needed to sustain transport above $\Theta^r_t$. While entrainment by turbulent events sustains intermittent transport, entrainment by particle-bed impacts sustains continuous transport. Combining our interpretations with the critical energy criterion for incipient motion by Valyrakis and coworkers, we put forward a new conceptual picture of sediment transport intermittency.
  • Observations of Pluto's surface made by the New Horizons spacecraft indicates present-day nitrogen ice glaciation in and around the basin known as Sputnik Planum. Motivated by these observations, we have developed an evolutionary glacial flow model of solid nitrogen ice taking into account its published thermophysical and rheologies properties. This model assumes that glacial ice layers flow laminarly and have low aspect ratios which permits a vertically integrated mathematical formulation. We assess the conditions for the validity of laminar nitrogen ice motion by revisiting the problem of the onset of solid-state buoyant convection of nitrogen ice for a variety of bottom thermal boundary conditions. Subject to uncertainties in nitrogen ice rheology, nitrogen ice layers are estimated to flow laminarly for thicknesses less than 400-1000 meters. The resulting mass-flux formulation for when the nitrogen ice flows as a laminar dry glacier is characterized by an Arrhenius-Glen functional form. The flow model developed is used here to qualitatively answer some questions motivated by observed glacial flow features found on Sputnik Planum. We find that the wavy transverse dark features found along the northern shoreline of Sputnik Planum may be a transitory imprint of shallow topography just beneath the ice surface suggesting the possibility that a major shoreward flow event happened relatively recently within the last few hundred years. Model results also support the interpretation that the prominent darkened features resembling flow lobes observed along the eastern shoreline of the Sputnik Planum basin may be a result of wet nitrogen glacial ice flowing into the basin from the pitted highlands of eastern Tombaugh Regio.
  • AMORPH utilizes a new Bayesian statistical approach to interpreting X-ray diffraction results of samples with both crystalline and amorphous components. AMORPH fits X-ray diffraction patterns with a mixture of narrow and wide components, simultaneously inferring all of the model parameters and quantifying their uncertainties. The program simulates background patterns previously applied manually, providing reproducible results, and significantly reducing inter- and intra-user biases. This approach allows for the quantification of amorphous and crystalline materials and for the characterization of the amorphous component, including properties such as the centre of mass, width, skewness, and nongaussianity of the amorphous component. Results demonstrate the applicability of this program for calculating amorphous contents of volcanic materials and independently modeling their properties in compositionally variable materials.
  • Here, we analyze in natural time $\chi$, the slow propagation of a macroscopic crack in brittle heterogeneous materials through sudden jumps and energy release events which are power law distributed with universal exponents. This macroscopic crack growth is currently believed to exhibit similar characteristics with the seismicity associated with earthquakes. Considering that the crack front is self affine and exhibits Family-Vicsek universal scaling, we show that the variance $\kappa_1 (\equiv\langle \chi^2\rangle -\langle \chi \rangle^2)$ of natural time is equal to 0.0686, which almost coincides with the value $\kappa_1\approx 0.07$ obtained from the seismicity preceding major earthquakes. This sheds light on the determination of the occurrence time of an impending mainshock.
  • The phase space of driftons (drift-wave quanta) is studied within the generalized Hasegawa--Mima collisionless-plasma model in the presence of zonal flows. This phase space is made intricate by the corrections to the drifton ray equations that were recently proposed by Parker [J. Plasma Phys. $\textbf{82}$, 95820602 (2016)] and Ruiz $\textit{et al.}$ [Phys. Plasmas $\textbf{23}$, 122304 (2016)]. Contrary to the traditional geometrical-optics (GO) model of the drifton dynamics, it is found that driftons can be not only trapped or passing, but they can also accumulate spatially while experiencing indefinite growth of their momenta. In particular, it is found that the Rayleigh--Kuo threshold known from geophysics corresponds to the regime when such "runaway" trajectories are the only ones possible. On one hand, this analysis helps visualize the development of the zonostrophic instability, particularly its nonlinear stage, which is studied here both analytically and through wave-kinetic simulations. On the other hand, the GO theory predicts that zonal flows above the Rayleigh--Kuo threshold can only grow; hence, the deterioration of intense zonal flows cannot be captured within a GO model. In particular, this means that the so-called tertiary instability of intense zonal flows cannot be adequately described within the quasilinear wave kinetic equation, contrary to some previous studies.
  • Large continental earthquakes activate multiple faults in a complex fault system, dynamically inducing co-seismic damage around them. The 2016 Mw 7.8 Kaikoura earthquake in the northern South Island of New Zealand has been reported as one of the most complex continental earthquakes ever documented1, which resulted in a distinctive on and off-fault deformation pattern. Previous geophysical studies confirm that the rupture globally propagated northward from epicenter. However, the exact rupture- propagation path is still not well understood because of the geometrical complexity, partly at sea, and the possibility of a blind thrust. Here we use a combination of state-of- the-art observation of surface deformation, provided by optical image correlation, and first principle physics-based numerical modeling to determine the most likely rupture path. We quantify in detail the observed horizontal co-seismic deformation and identify specific off-fault damage zones in the area of the triple junction between the Jordan, the Kekerengu and the Papatea fault segments. We also model dynamic rupture propagation, including the activation of off-fault damage, for two alternative rupture scenarios through the fault triple junction. Comparing our observations with the results from the above two modeled scenarios we show that only one of the scenarios best explains both the on and off-fault deformation fields. Our results provide a unique insight into the rupture pathway, by observing, and modeling, both on and off-fault deformation. We propose this combined approach here to narrow down the possible rupture scenarios for large continental earthquakes accompanied by co-seismic off-fault damage.
  • Please refer to the Abstract part in the paper.
  • We study the relationship between the El Ni\~no--Southern Oscillation (ENSO) and the Indian summer monsoon in ensemble simulations from state-of-the-art climate models, the Max Planck Institute Earth System Model (MPI-ESM) and the Community Earth System Model (CESM). We consider two simple variables: the Tahiti--Darwin sea-level pressure difference and the Northern Indian precipitation. We utilize ensembles converged to the snapshot attractor of the system for analyzing possible changes (i) in the teleconnection between the fluctuations of the two variables, and (ii) in their climatic means. (i) With very high confidence, we detect an increase in the strength of the teleconnection in the MPI-ESM under historical forcing between 1890 and 2005, which is in contrast with scientific consensus. In the MPI-ESM no similar increase is present between 2006 and 2099 under the Representative Concentration Pathway 8.5 (RCP8.5), and in a 110-year-long 1-percent pure CO2 scenario; neither is in the CESM between 1960 and 2100 with historical forcing and RCP8.5. The static susceptibility of the strength of the teleconnection with respect to radiative forcing (assuming an instantaneous and linear response) is at least three times larger in the historical MPI-ESM ensemble than in the others. (ii) In the other ensembles, the climatic mean is strongly displaced in the phase space projection spanned by the two variables. This displacement is nevertheless linear. However, the slope exhibits a strong seasonality, falsifying a hypothesis of a universal relation, an emergent constraint, that would be valid at all time scales between these two climatic means.
  • Certain geological features have been interpreted as evidence of channelized magma flow in the mantle, which is a compacting porous medium. Aharonov et al. (1995) developed a simple model of reactive porous flow and numerically analysed its instability to channels. The instability relies on magma advection against a chemical solubility gradient and the porosity-dependent permeability of the porous host rock. We extend the previous analysis by systematically mapping out the parameter space. Crucially, we augment numerical solutions with asymptotic analysis to better understand the physical controls on the instability. We derive scalings for critical conditions of the instability and analyse the associated bifurcation structure. We also determine scalings for the wavelength and growth rate of the channel structures that emerge. We obtain quantitative theories for and a physical understanding of: first, how advection or diffusion over the reactive time scale set the horizontal length scale of channels; second, the role of viscous compaction of the host rock, which also affects the vertical extent of channelized flow. These scalings allow us to derive estimates of the dimensions of emergent channels that are consistent with the geologic record.
  • Proximal soil sensors are taking hold in the understanding of soil hydrogeological processes involved in precision agriculture. In this context, permanently installed gamma ray spectroscopy stations represent one of the best space-time trade off methods at field scale. This study proved the feasibility and reliability of soil water content monitoring through a seven-month continuous acquisition of terrestrial gamma radiation in a tomato test field. By employing a 1 L sodium iodide detector placed at a height of 2.25 m, we investigated the gamma signal coming from an area having a ~25 m radius and from a depth of approximately 30 cm. Experimental values, inferred after a calibration measurement and corrected for the presence of biomass, were corroborated with gravimetric data acquired under different soil moisture conditions, giving an average absolute discrepancy of about 2%. A quantitative comparison was carried out with data simulated by AquaCrop, CRITeRIA, and IRRINET soil-crop system models. The different goodness of fit obtained in bare soil condition and during the vegetated period highlighted that CRITeRIA showed the best agreement with the experimental data over the entire data-taking period while, in presence of the tomato crop, IRRINET provided the best results.
  • This paper reports the stability conditions for intense zonal flows (ZFs) and the growth rate $\gamma_{\rm TI}$ of the corresponding "tertiary" instability (TI) within the generalized Hasegawa--Mima plasma model. The analytic calculation extends and revises Kuo's analysis of the mathematically similar barotropic vorticity equation for incompressible neutral fluids on a rotating sphere [H.-L. Kuo, J. Meteor. $\textbf{6}$, 105 (1949)]; then, the results are applied to the plasma case. An error in Kuo's original result is pointed out. An explicit analytic formula for TI is derived and compared with numerical calculations. It is shown that a ZF is TI-unstable under the Rayleigh--Kuo criterion known from geophysics plus the condition that the ZF wave number must exceed the inverse ion sound radius. For a sinusoidal ZF, these two necessary conditions combined together are also sufficient for the TI. For non-sinusoidal ZFs, the results are qualitatively consistent. As a corollary, there is no TI in the geometrical-optics limit, i.e., when the DW wavelength is small compared to the ZF scale. This also means that the traditional wave kinetic equation under the geometrical-optics assumption cannot adequately describe the ZF stability.
  • Ultra-low-frequency (ULF) oscillations of the magnetospheric cusps are observed permanently in the near-mid sector of auroral oval in the form of so-called IPCL. In this paper, we posed the question: do electromagnetic waves incident on cusp from the foreshock (it is a special region of cosmic plasma existing before the front of the magnetosphere) affect the IPCL regime? We have proposed a method of experimental investigation, which is based on the idea of the foreshock position dependence on the interplanetary magnetic field (IMF) orientation. A hypothesis has been put forward on the existence of a specific effect of the IPCL north-south asymmetry. We tested our hypothesis using IPCL observation data from the Hornsund (Spitsbergen) and Davis (Antarctica) observatories. We found confirmation of the hypothesis at a high confidence level of statistical significance. Keywords: magnetosphere, ultra-low-frequency oscillations, north-south asymmetry.
  • Melting and, conversely, solidification processes in the presence of convection are key to many geophysical problems. An essential question related to these phenomena concerns the estimation of the (time-evolving) melting rate, which is tightly connected to the turbulent convective dynamics in the bulk of the melt fluid and the heat transfer at the liquid-solid interface. In this work, we consider a convective-melting model, constructed as a generalization of the Rayleigh-B\'enard system, accounting for the basal melting of a solid. As the change of phase proceeds, a fluid layer grows at the heated bottom of the system and eventually reaches a turbulent convection state. By means of extensive Lattice-Boltzmann numerical simulations employing an enthalpy formulation of the governing equations, we explore the model dynamics in two and three-dimensional configurations. The focus of the analysis is on the scaling of global quantities like the heat flux and the kinetic energy with the Rayleigh number, as well as on the interface morphology and the effects of space dimensionality. Independently of dimensionality, we find that the convective-melting system behavior shares strong resemblances with that of the Rayleigh-B\'enard one, and that the heat flux is only weakly enhanced with respect to that case. Such similarities are understood, at least to some extent, considering the resulting slow motion of the melting front (with respect to the turbulent fluid velocity fluctuations) and its generally little roughness (compared to the height of the fluid layer). Varying the Stefan number, accounting for the thermodynamical properties of the material, also seems to have only a mild effect, which implies the possibility to extrapolate results in numerically delicate low-Stefan setups from more convenient high-Stefan ones.
  • To optimally monitor earthquake-generating processes, seismologists have sought to lower detection sensitivities ever since instrumental seismic networks were started about a century ago. Recently, it has become possible to search continuous waveform archives for replicas of previously recorded events (template matching), which has led to at least an order of magnitude increase in the number of detected earthquakes and greatly sharpened our view of geological structures. Earthquake catalogs produced in this fashion, however, are heavily biased in that they are completely blind to events for which no templates are available, such as in previously quiet regions or for very large magnitude events. Here we show that with deep learning we can overcome such biases without sacrificing detection sensitivity. We trained a convolutional neural network (ConvNet) on the vast hand-labeled data archives of the Southern California Seismic Network to detect seismic body wave phases. We show that the ConvNet is extremely sensitive and robust in detecting phases, even when masked by high background noise, and when the ConvNet is applied to new data that is not represented in the training set (in particular, very large magnitude events). This generalized phase detection (GPD) framework will significantly improve earthquake monitoring and catalogs, which form the underlying basis for a wide range of basic and applied seismological research.
  • A first assessment of the TanDEM-X DEMs over Brazilian territory is presented through a comparison with SRTM, ASTER GDEM and ALOS AW3D30 DEMs in seven study areas with distinct geomorphological contexts, vegetation coverage and land use. Visual analysis and elevation histograms point to a finer effective spatial resolution of TanDEM-X compared to SRTM and ASTER GDEM. In areas of open vegetation, TanDEM-X lower elevations indicate a better penetration of the radar signal. DEMs of differences (DoDs) allowed the identification of issues inherent to the production methods of the analyzed DEMs, such as mast oscillations in SRTM data and mismatch between adjacent scenes in ASTER GDEM and ALOS AW3D30. A systematic difference in elevations between TanDEM-X 12m, TanDEM-X 30m and SRTM was observed in the steep slopes of the coastal ranges, related to the moving-window process used to resample the 12m data to a 30m pixel size. Due its simplicity, it is strongly recommended to produce a DoD with SRTM before using ASTER GDEM or ALOS AW3D30 in any analysis, to evaluate if the area of interest is affected by these problems. The DoDs also highlighted changes in land use in the time span between the acquisition of SRTM (2000) and TanDEM-X (2013) data, whether by natural causes or by human interference in the environment.
  • A careful reading of old articles puts Olivier Pauluis' criticisms concerning the definition of isentropic processes in terms of a potential temperature closely associated with the entropy of moist air, together with the third principle of thermodynamics, into perspective.