
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
projectedgradient algorithm  faithfulness to the recorded traces is
constrained by the measured noise level. Wavelet reestimation 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 multiplepoint statistics framework is often
computationally prohibitive for highdimensional problems. To partly address
this, we introduce and evaluate a new trainingimage 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) lowdimensional parameterization, thereby allowing for efficient
probabilistic inversion using stateoftheart 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 steadystate flow
and 3D transient hydraulic tomography with and without direct conditioning data
are used to illustrate the effectiveness of our proposed SGANbased 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 leakoff, previously demonstrated for the PKN and KGD models, is
extended to obtain solutions for a pennyshaped crack. The numerical scheme is
capable of dealing with both the viscosity and toughness dominated regimes,
with the fracture being driven by a powerlaw 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,
semianalytical approximations of the solution are provided for the zero
leakoff case. A comparison with other results available in the literature is
undertaken.

A multistep steadystate (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
steadystate 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. Multistep 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 particlescale 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 normalbed pressure) evaluated at the base of the transport layer
(i.e., the effective elevation of energetic particlebed 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 solidliquid
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 nonstationarity 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 nonstationarity
from effects due to nonlinearity (and the nonGaussianity it causes), a new
doubly stochastic advectiondiffusiondecay 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 nonstationarity 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,
timesmoothed, and spacesmoothed covariances. It is shown that the stronger
the nonstationarity, the less useful the static covariance matrix becomes and
the more beneficial the timesmoothed covariances are. Timesmoothing of
backgrounderror covariances proved to be systematically more useful than their
spacesmoothing. Under nonstationarity, a filter that extends the (previously
proposed by the authors) Hierarchical Bayes Ensemble Filter and accommodates
the three covarianceblending 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 electronelectron scattering, usually neglected at high
temperatures in transition metals. Our calculations show a quasilinear
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 electronelectron
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 selfsimilarity in random binary rooted trees. In a wellunderstood
case of GaltonWatson trees, a distribution on a space of trees is said to be
selfsimilar if it is invariant with respect to the operation of pruning, which
cuts the tree leaves. This only happens for the critical GaltonWatson tree (a
constant process progeny), which also exhibits other special symmetries. We
extend the pruneinvariance setup to arbitrary binary trees with edge lengths.
In this general case the class of selfsimilar 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
selfsimilarity definitions (including mean selfsimilarity and selfsimilarity
in edgelengths) 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 GaltonWatson 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 extremeinvariant processes,
constructed as the Harris paths of a selfsimilar hierarchical branching
process, whose local minima has the same (linearly scaled) distribution as the
original process.

Using particlescale simulations of nonsuspended 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 Shieldslike
threshold diagram in which a Stokeslike 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 particlebed 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
presentday 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
solidstate 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
4001000 meters. The resulting massflux formulation for when the nitrogen ice
flows as a laminar dry glacier is characterized by an ArrheniusGlen 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 Xray
diffraction results of samples with both crystalline and amorphous components.
AMORPH fits Xray 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 intrauser 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 FamilyVicsek 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 (driftwave quanta) is studied within the
generalized HasegawaMima collisionlessplasma 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 geometricaloptics (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 RayleighKuo
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
wavekinetic simulations. On the other hand, the GO theory predicts that zonal
flows above the RayleighKuo threshold can only grow; hence, the deterioration
of intense zonal flows cannot be captured within a GO model. In particular,
this means that the socalled 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 coseismic 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 offfault 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 stateof theart
observation of surface deformation, provided by optical image correlation, and
first principle physicsbased numerical modeling to determine the most likely
rupture path. We quantify in detail the observed horizontal coseismic
deformation and identify specific offfault 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 offfault 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 offfault deformation fields. Our results provide a unique insight into
the rupture pathway, by observing, and modeling, both on and offfault
deformation. We propose this combined approach here to narrow down the possible
rupture scenarios for large continental earthquakes accompanied by coseismic
offfault damage.

Please refer to the Abstract part in the paper.

We study the relationship between the El Ni\~noSouthern Oscillation (ENSO)
and the Indian summer monsoon in ensemble simulations from stateoftheart
climate models, the Max Planck Institute Earth System Model (MPIESM) and the
Community Earth System Model (CESM). We consider two simple variables: the
TahitiDarwin sealevel 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 MPIESM under historical forcing between 1890 and 2005,
which is in contrast with scientific consensus. In the MPIESM no similar
increase is present between 2006 and 2099 under the Representative
Concentration Pathway 8.5 (RCP8.5), and in a 110yearlong 1percent 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 MPIESM 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 porositydependent 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
spacetime trade off methods at field scale. This study proved the feasibility
and reliability of soil water content monitoring through a sevenmonth
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 soilcrop 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 datataking
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 HasegawaMima 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 TIunstable under the RayleighKuo
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 nonsinusoidal
ZFs, the results are qualitatively consistent. As a corollary, there is no TI
in the geometricaloptics 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 geometricaloptics assumption cannot adequately describe the ZF
stability.

Ultralowfrequency (ULF) oscillations of the magnetospheric cusps are
observed permanently in the nearmid sector of auroral oval in the form of
socalled 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 northsouth 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, ultralowfrequency
oscillations, northsouth 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 (timeevolving) melting rate,
which is tightly connected to the turbulent convective dynamics in the bulk of
the melt fluid and the heat transfer at the liquidsolid interface. In this
work, we consider a convectivemelting model, constructed as a generalization
of the RayleighB\'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 LatticeBoltzmann numerical simulations employing an enthalpy
formulation of the governing equations, we explore the model dynamics in two
and threedimensional 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
convectivemelting system behavior shares strong resemblances with that of the
RayleighB\'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 lowStefan setups from more convenient
highStefan ones.

To optimally monitor earthquakegenerating 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 handlabeled 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 TanDEMX 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 TanDEMX compared to SRTM and ASTER GDEM. In areas of
open vegetation, TanDEMX 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 TanDEMX 12m,
TanDEMX 30m and SRTM was observed in the steep slopes of the coastal ranges,
related to the movingwindow 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 TanDEMX (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.