pubbl0.bib
@comment{{This file has been generated by bib2bib 1.98}}
@comment{{Command line: /usr/bin/bib2bib --quiet -c 'not journal:"Discussions"' -c 'not title:"Correction to"' -c 'abstract:"turbulent" or abstract:"large-eddy"' -c $type="ARTICLE" -oc pubbl0.txt -ob pubbl0.bib spiga.link.bib}}
@article{2018SSRv..214..109S,
author = {{Spiga}, A. and {Banfield}, D. and {Teanby}, N.~A. and {Forget}, F. and
{Lucas}, A. and {Kenda}, B. and {Rodriguez Manfredi}, J.~A. and
{Widmer-Schnidrig}, R. and {Murdoch}, N. and {Lemmon}, M.~T. and
{Garcia}, R.~F. and {Martire}, L. and {Karatekin}, {\"O}. and
{Le Maistre}, S. and {Van Hove}, B. and {Dehant}, V. and {Lognonné}, P. and
{Mueller}, N. and {Lorenz}, R. and {Mimoun}, D. and {Rodriguez}, S. and
{Beucler}, {\'E}. and {Daubar}, I. and {Golombek}, M.~P. and
{Bertrand}, T. and {Nishikawa}, Y. and {Millour}, E. and {Rolland}, L. and
{Brissaud}, Q. and {Kawamura}, T. and {Mocquet}, A. and {Martin}, R. and
{Clinton}, J. and {Stutzmann}, {\'E}. and {Spohn}, T. and {Smrekar}, S. and
{Banerdt}, W.~B.},
title = {{Atmospheric Science with InSight}},
journal = {\ssr},
keywords = {Mars, InSight, Atmospheric science, Planetary atmospheres},
year = 2018,
volume = 214,
eid = {109},
pages = {109},
abstract = {{In November 2018, for the first time a dedicated geophysical station,
the InSight lander, will be deployed on the surface of Mars. Along with
the two main geophysical packages, the Seismic Experiment for Interior
Structure (SEIS) and the Heat-Flow and Physical Properties Package
(HP$^{3}$), the InSight lander holds a highly sensitive pressure
sensor (PS) and the Temperature and Winds for InSight (TWINS)
instrument, both of which (along with the InSight FluxGate (IFG)
Magnetometer) form the Auxiliary Sensor Payload Suite (APSS). Associated
with the RADiometer (RAD) instrument which will measure the surface
brightness temperature, and the Instrument Deployment Camera (IDC) which
will be used to quantify atmospheric opacity, this will make InSight
capable to act as a meteorological station at the surface of Mars. While
probing the internal structure of Mars is the primary scientific goal of
the mission, atmospheric science remains a key science objective for
InSight. InSight has the potential to provide a more continuous and
higher-frequency record of pressure, air temperature and winds at the
surface of Mars than previous in situ missions. In the paper, key
results from multiscale meteorological modeling, from Global Climate
Models to Large-Eddy Simulations, are described as a reference for
future studies based on the InSight measurements during operations. We
summarize the capabilities of InSight for atmospheric observations, from
profiling during Entry, Descent and Landing to surface measurements
(pressure, temperature, winds, angular momentum), and the plans for how
InSight's sensors will be used during operations, as well as possible
synergies with orbital observations. In a dedicated section, we describe
the seismic impact of atmospheric phenomena (from the point of view of
both ``noise'' to be decorrelated from the seismic signal and ``signal'' to
provide information on atmospheric processes). We discuss in this
framework Planetary Boundary Layer turbulence, with a focus on
convective vortices and dust devils, gravity waves (with idealized
modeling), and large-scale circulations. Our paper also presents
possible new, exploratory, studies with the InSight instrumentation:
surface layer scaling and exploration of the Monin-Obukhov model,
aeolian surface changes and saltation / lifing studies, and monitoring
of secular pressure changes. The InSight mission will be instrumental in
broadening the knowledge of the Martian atmosphere, with a unique set of
measurements from the surface of Mars.
}},
doi = {10.1007/s11214-018-0543-0},
adsurl = {https://ui.adsabs.harvard.edu/abs/2018SSRv..214..109S},
localpdf = {REF/2018SSRv..214..109S.pdf},
adsnote = {Provided by the SAO/NASA Astrophysics Data System}
}
@article{2017SSRv..211..501K,
author = {{Kenda}, B. and {Lognonné}, P. and {Spiga}, A. and {Kawamura}, T. and
{Kedar}, S. and {Banerdt}, W.~B. and {Lorenz}, R. and {Banfield}, D. and
{Golombek}, M.},
title = {{Modeling of Ground Deformation and Shallow Surface Waves Generated by Martian Dust Devils and Perspectives for Near-Surface Structure Inversion}},
journal = {\ssr},
keywords = {Dust devils, Mars, Ground tilt, Subsurface, Large-eddy simulation, Insight},
year = 2017,
volume = 211,
pages = {501-524},
abstract = {{We investigated the possible seismic signatures of dust devils on Mars,
both at long and short period, based on the analysis of Earth data and
on forward modeling for Mars. Seismic and meteorological data collected
in the Mojave Desert, California, recorded the signals generated by dust
devils. In the 10-100 s band, the quasi-static surface deformation
triggered by pressure fluctuations resulted in detectable ground-tilt
effects: these are in good agreement with our modeling based on
Sorrells' theory. In addition, high-frequency records also exhibit a
significant excitation in correspondence to dust devil episodes. Besides
wind noise, this signal includes shallow surface waves due to the
atmosphere-surface coupling and is used for a preliminary inversion of
the near-surface S-wave profile down to 50 m depth. In the case of Mars,
we modeled the long-period signals generated by the pressure field
resulting from turbulence-resolving Large-Eddy Simulations. For typical
dust-devil-like vortices with pressure drops of a couple Pascals, the
corresponding horizontal acceleration is of a few nm/s$^{2}$ for
rocky subsurface models and reaches 10-20 nm/s$^{2}$ for weak
regolith models. In both cases, this signal can be detected by the
Very-Broad Band seismometers of the InSight/SEIS experiment up to a
distance of a few hundred meters from the vortex, the amplitude of the
signal decreasing as the inverse of the distance. Atmospheric vortices
are thus expected to be detected at the InSight landing site; the
analysis of their seismic and atmospheric signals could lead to
additional constraints on the near-surface structure, more precisely on
the ground compliance and possibly on the seismic velocities.
}},
doi = {10.1007/s11214-017-0378-0},
adsurl = {https://ui.adsabs.harvard.edu/abs/2017SSRv..211..501K},
localpdf = {REF/2017SSRv..211..501K.pdf},
adsnote = {Provided by the SAO/NASA Astrophysics Data System}
}
@article{2017Icar..291...82P,
author = {{Pottier}, A. and {Forget}, F. and {Montmessin}, F. and {Navarro}, T. and
{Spiga}, A. and {Millour}, E. and {Szantai}, A. and {Madeleine}, J.-B.
},
title = {{Unraveling the martian water cycle with high-resolution global climate simulations}},
journal = {\icarus},
keywords = {Mars atmosphere, Atmospheres dynamics, Mars, Climate, Meteorology},
year = 2017,
volume = 291,
pages = {82-106},
abstract = {{Global climate modeling of the Mars water cycle is usually performed at
relatively coarse resolution (200 - 300km), which may not be sufficient
to properly represent the impact of waves, fronts, topography effects on
the detailed structure of clouds and surface ice deposits. Here, we
present new numerical simulations of the annual water cycle performed at
a resolution of 1{\deg} {\times} 1{\deg} ({\sim} 60 km in latitude). The
model includes the radiative effects of clouds, whose influence on the
thermal structure and atmospheric dynamics is significant, thus we also
examine simulations with inactive clouds to distinguish the direct
impact of resolution on circulation and winds from the indirect impact
of resolution via water ice clouds. To first order, we find that the
high resolution does not dramatically change the behavior of the system,
and that simulations performed at {\sim} 200 km resolution capture well
the behavior of the simulated water cycle and Mars climate.
Nevertheless, a detailed comparison between high and low resolution
simulations, with reference to observations, reveal several significant
changes that impact our understanding of the water cycle active today on
Mars. The key northern cap edge dynamics are affected by an increase in
baroclinic wave strength, with a complication of northern summer
dynamics. South polar frost deposition is modified, with a westward
longitudinal shift, since southern dynamics are also influenced.
Baroclinic wave mode transitions are observed. New transient phenomena
appear, like spiral and streak clouds, already documented in the
observations. Atmospheric circulation cells in the polar region exhibit
a large variability and are fine structured, with slope winds. Most
modeled phenomena affected by high resolution give a picture of a more
turbulent planet, inducing further variability. This is challenging for
long-period climate studies.
}},
doi = {10.1016/j.icarus.2017.02.016},
adsurl = {https://ui.adsabs.harvard.edu/abs/2017Icar..291...82P},
localpdf = {REF/2017Icar..291...82P.pdf},
adsnote = {Provided by the SAO/NASA Astrophysics Data System}
}
@article{2017JGRE..122..134L,
author = {{Lefèvre}, M. and {Spiga}, A. and {Lebonnois}, S.},
title = {{Three-dimensional turbulence-resolving modeling of the Venusian cloud layer and induced gravity waves}},
journal = {Journal of Geophysical Research (Planets)},
keywords = {3-D mesoscale modeling, Venus, convective cloud layer, gravity waves},
year = 2017,
volume = 122,
pages = {134-149},
abstract = {{The impact of the cloud convective layer of the atmosphere of Venus on
the global circulation remains unclear. The recent observations of
gravity waves at the top of the cloud by the Venus Express mission
provided some answers. These waves are not resolved at the scale of
global circulation models (GCM); therefore, we developed an
unprecedented 3-D turbulence-resolving large-eddy simulations (LES)
Venusian model using the Weather Research and Forecast terrestrial
model. The forcing consists of three different heating rates: two
radiative ones for solar and infrared and one associated with the
adiabatic cooling/warming of the global circulation. The rates are
extracted from the Laboratoire de Météorlogie Dynamique
Venus GCM using two different cloud models. Thus, we are able to
characterize the convection and associated gravity waves in function of
latitude and local time. To assess the impact of the global circulation
on the convective layer, we used rates from a 1-D radiative-convective
model. The resolved layer, taking place between 1.0 {\times}
10$^{5}$ and 3.8 {\times} 10$^{4}$ Pa (48-53 km), is
organized as polygonal closed cells of about 10 km wide with vertical
wind of several meters per second. The convection emits gravity waves
both above and below the convective layer leading to temperature
perturbations of several tenths of kelvin with vertical wavelength
between 1 and 3 km and horizontal wavelength from 1 to 10 km. The
thickness of the convective layer and the amplitudes of waves are
consistent with observations, though slightly underestimated. The global
dynamics heating greatly modify the convective layer.
}},
doi = {10.1002/2016JE005146},
adsurl = {https://ui.adsabs.harvard.edu/abs/2017JGRE..122..134L},
localpdf = {REF/2017JGRE..122..134L.pdf},
adsnote = {Provided by the SAO/NASA Astrophysics Data System}
}
@article{2016SSRv..203..245S,
author = {{Spiga}, A. and {Barth}, E. and {Gu}, Z. and {Hoffmann}, F. and
{Ito}, J. and {Jemmett-Smith}, B. and {Klose}, M. and {Nishizawa}, S. and
{Raasch}, S. and {Rafkin}, S. and {Takemi}, T. and {Tyler}, D. and
{Wei}, W.},
title = {{Large-Eddy Simulations of Dust Devils and Convective Vortices}},
journal = {\ssr},
keywords = {Dust devils, Large-Eddy Simulations, Convective vortices, Convective boundary layer},
year = 2016,
volume = 203,
pages = {245-275},
abstract = {{In this review, we address the use of numerical computations called
Large-Eddy Simulations (LES) to study dust devils, and the more general
class of atmospheric phenomena they belong to (convective vortices). We
describe the main elements of the LES methodology. We review the
properties, statistics, and variability of dust devils and convective
vortices resolved by LES in both terrestrial and Martian environments.
The current challenges faced by modelers using LES for dust devils are
also discussed in detail.
}},
doi = {10.1007/s11214-016-0284-x},
adsurl = {https://ui.adsabs.harvard.edu/abs/2016SSRv..203..245S},
localpdf = {REF/2016SSRv..203..245S.pdf},
adsnote = {Provided by the SAO/NASA Astrophysics Data System}
}
@article{2013JGRE..118.1468C,
author = {{Cola{\"i}tis}, A. and {Spiga}, A. and {Hourdin}, F. and {Rio}, C. and
{Forget}, F. and {Millour}, E.},
title = {{A thermal plume model for the Martian convective boundary layer}},
journal = {Journal of Geophysical Research (Planets)},
archiveprefix = {arXiv},
eprint = {1306.6215},
primaryclass = {physics.ao-ph},
keywords = {Mars, atmosphere, convection, boundary layer, large-eddy simulations, PBL parameterization},
year = 2013,
volume = 118,
pages = {1468-1487},
abstract = {{The Martian planetary boundary layer (PBL) is a crucial component of the
Martian climate system. Global climate models (GCMs) and mesoscale
models (MMs) lack the resolution to predict PBL mixing which is
therefore parameterized. Here we propose to adapt the ``thermal plume''
model, recently developed for Earth climate modeling, to Martian GCMs,
MMs, and single-column models. The aim of this physically based
parameterization is to represent the effect of organized turbulent
structures (updrafts and downdrafts) on the daytime PBL transport, as it
is resolved in large-eddy simulations (LESs). We find that the
terrestrial thermal plume model needs to be modified to satisfyingly
account for deep turbulent plumes found in the Martian convective PBL.
Our Martian thermal plume model qualitatively and quantitatively
reproduces the thermal structure of the daytime PBL on Mars:
superadiabatic near-surface layer, mixing layer, and overshoot region at
PBL top. This model is coupled to surface layer parameterizations taking
into account stability and turbulent gustiness to calculate
surface-atmosphere fluxes. Those new parameterizations for the surface
and mixed layers are validated against near-surface lander measurements.
Using a thermal plume model moreover enables a first-order estimation of
key turbulent quantities (e.g., PBL height and convective plume
velocity) in Martian GCMs and MMs without having to run costly LESs.
}},
doi = {10.1002/jgre.20104},
adsurl = {https://ui.adsabs.harvard.edu/abs/2013JGRE..118.1468C},
localpdf = {REF/2013JGRE..118.1468C.pdf},
adsnote = {Provided by the SAO/NASA Astrophysics Data System}
}
@article{2012GI......1..151S,
author = {{Spiga}, A.},
title = {{Comment on ''Observing desert dust devils with a pressure logger`` by Lorenz (2012) - insights on measured pressure fluctuations from large-eddy simulations}},
journal = {Geoscientific Instrumentation, Methods and Data Systems},
year = 2012,
volume = 1,
pages = {151-154},
abstract = {{Lorenz et al. (2012) proposes to use pressure loggers for long-term
field measurements in terrestrial deserts. The dataset obtained through
this method features both pressure drops (reminiscent of dust devils)
and periodic convective signatures. Here we use large-eddy simulations
to provide an explanation for those periodic convective signatures and
to argue that pressure measurements in deserts have broader applications
than monitoring dust devils.
}},
doi = {10.5194/gi-1-151-2012},
adsurl = {https://ui.adsabs.harvard.edu/abs/2012GI......1..151S},
localpdf = {REF/2012GI......1..151S.pdf},
adsnote = {Provided by the SAO/NASA Astrophysics Data System}
}
@article{2011P&SS...59..915S,
author = {{Spiga}, A.},
title = {{Elements of comparison between Martian and terrestrial mesoscale meteorological phenomena: Katabatic winds and boundary layer convection}},
journal = {\planss},
year = 2011,
volume = 59,
pages = {915-922},
abstract = {{Terrestrial and Martian atmospheres are both characterised by a large
variety of mesoscale meteorological events, occurring at horizontal
scales of hundreds of kilometres and below. Available measurements from
space exploration and recently developed high-resolution numerical tools
have given insights into Martian mesoscale phenomena, as well as
similarities and differences with their terrestrial counterparts. The
remarkable intensity of Martian mesoscale events compared to terrestrial
phenomena mainly results from low density and strong radiative control.
This is exemplified in the present paper by discussing two mesoscale
phenomena encountered in the lowest atmospheric levels of both planets
with notable differences: nighttime katabatic winds (drainage flow down
sloping terrains) and daytime boundary layer convection (vertical growth
of mixed layer over heated surfaces). While observations of katabatic
events are difficult on Earth, except over vast ice sheets, intense
clear-cut downslope circulations are expected to be widespread on Mars.
Convective motions in the daytime Martian boundary layer are primarily
driven by radiative contributions, usually negligible on Earth where
sensible heat flux dominates, and exhibit turbulent variances one order
of magnitude larger. Martian maximum heat fluxes are not attained close
to the surface as on Earth but a few hundreds of metres above, which
implies generalised definitions for mixing layer scales such as vertical
velocity w$_{*}$. Measurements on Mars of winds in uneven
topographical areas and of heat fluxes over flat terrains could be
useful to assess general principles of mesoscale meteorology applicable
to both terrestrial and Martian environments.
}},
doi = {10.1016/j.pss.2010.04.025},
adsurl = {https://ui.adsabs.harvard.edu/abs/2011P%26SS...59..915S},
localpdf = {REF/2011P_26SS...59..915S.pdf},
adsnote = {Provided by the SAO/NASA Astrophysics Data System}
}
@article{2010QJRMS.136..414S,
author = {{Spiga}, A. and {Forget}, F. and {Lewis}, S.~R. and {Hinson}, D.~P.
},
title = {{Structure and dynamics of the convective boundary layer on Mars as inferred from large-eddy simulations and remote-sensing measurements}},
journal = {Quarterly Journal of the Royal Meteorological Society},
year = 2010,
volume = 136,
pages = {414-428},
abstract = {{The structure of the Martian convective boundary layer (BL) is decribed
by means of a novel approach involving both modelling and data analysis.
Mars Express radio-occultation (RO) temperature profiles are compared to
large-eddy simulations (LESs) performed with the Martian mesoscale
model. The model combines the Martian radiative transfer, soil and
surface layer schemes designed at Laboratoire de
M{\~A}{\copy}t{\~A}{\copy}orologie Dynamique (LMD) with the most
recent version of the Weather Research and Forecast (WRF) fully
compressible non-hydrostatic dynamical core. The key roles of the
vertical resolution and, to lesser extent, of the domain horizontal
extent have been investigated to ensure the robustness of the LES
results. The dramatic regional variations of the BL depth are
quantitatively reproduced by the Martian LES. Intense BL dynamics are
found to underlie the measured depths (up to 9 km): vertical speed up to
20 m s-1, heat flux up to 2.7 K m s-1 and turbulent kinetic energy up to 26
m2 s-2. Under specific conditions, both the model and the
measurements show a distinctive positive correlation between surface
topography and BL depth. Our interpretation is that, in the tenuous CO2
Martian near-surface environment, the daytime BL is to first order
controlled by the infrared radiative heating, fairly independent of
elevation, which implies a simple correlation between the BL potential
temperature and the inverse pressure (``pressure effect''). No prominent
``pressure effect'' is in action on Earth where sensible heat flux
dominates the BL energy budget. Both RO observations and numerical
simulations confirm the terrain-following behaviour of near-surface
temperature on Mars induced by the dominant radiative influence. The
contribution of the Martian sensible heat flux is not negligible and
results in a given isotherm in the BL being comparatively closer to the
ground at higher surface elevation. The strong radiative control of the
Martian convective BL implies a generalised formulation for the BL
dimensionless quantities. Based on this formulation and the variety of
simulated BL depths by the LES, new similarity relationships for the
Martian convective BL in quasi-steady midday conditions are derived.
Rigorous comparisons between the Martian and terrestrial BL and fast
computations of the mean Martian BL turbulent statistics are now made
possible by such similarity laws.
}},
adsurl = {https://ui.adsabs.harvard.edu/abs/2010QJRMS.136..414S},
localpdf = {REF/2010QJRMS.136..414S.pdf},
adsnote = {Provided by the SAO/NASA Astrophysics Data System}
}
@article{2010IJMSE...5..146S,
author = {{Spiga}, A. and {Lewis}, S.~R.},
title = {{Martian mesoscale and microscale wind variability of relevance for dust lifting}},
journal = {International Journal of Mars Science and Exploration},
year = 2010,
volume = 5,
pages = {146-158},
abstract = {{Background: Mars is both a windy and dusty environment. Ariborne dust is
a crucial climate component on Mars. It impacts atmospheric circulations
at large-, meso- and micro-scales, which in turn control dust lifting
from the surface and transport in the atmosphere. Dust lifting processes
and feedbacks on atmospheric circulations are currently not well
understood. Method: Our purpose is to show how mesoscale models and
large-eddy simulations help to explore small-scale circulation patterns
which are potentially important for lifting dust into the atmosphere but
which are unresolved by global climate models. We focus on variations of
friction velocity, u*, relevant for dust lifting, in particular
investigating maximum values and the spatial and temporal variability of
u*. Conclusion: Meteorological scales between 100 km and 10 km can be
studied by high-resolution global circulation and limited-area mesoscale
models, which both show strong topographic control of the daytime and
nighttime near-surface winds. Scales below 10 km and 1 km are dominated
by turbulent gusts and dust devils, two distinct convective boundary
layer processes likely to lift dust from the surface. In low-latitude
regions, boundary layer depth and friction velocity u* are correlated
with surface altimetry. Further studies will be carried out to
parameterize lifting by boundary layer processes and dust radiative
effects once transported in the atmosphere.
}},
doi = {10.1555/mars.2010.0006},
adsurl = {https://ui.adsabs.harvard.edu/abs/2010IJMSE...5..146S},
localpdf = {REF/2010IJMSE...5..146S.pdf},
adsnote = {Provided by the SAO/NASA Astrophysics Data System}
}
@article{2009JGRE..114.2009S,
author = {{Spiga}, A. and {Forget}, F.},
title = {{A new model to simulate the Martian mesoscale and microscale atmospheric circulation: Validation and first results}},
journal = {Journal of Geophysical Research (Planets)},
keywords = {Planetary Sciences: Solar System Objects: Mars, Atmospheric Processes: Mesoscale meteorology, Atmospheric Processes: Boundary layer processes, Atmospheric Processes: Regional modeling, Planetary Sciences: Fluid Planets: Meteorology (3346)},
year = 2009,
volume = 114,
eid = {E02009},
pages = {E02009},
abstract = {{The Laboratoire de Météorologie Dynamique (LMD) Mesoscale
Model is a new versatile simulator of the Martian atmosphere and
environment at horizontal scales ranging from hundreds of kilometers to
tens of meters. The model combines the National Centers for
Environmental Prediction(NCEP)-National Center for Atmospheric Research
(NCAR) fully compressible nonhydrostatic Advanced Research Weather
Research and Forecasting (ARW-WRF) dynamical core, adapted to Mars, with
the LMD-general circulation model (GCM) comprehensive set of physical
parameterizations for the Martian dust, CO$_{2}$, water, and
photochemistry cycles. Since LMD-GCM large-scale simulations are also
used to drive the mesoscale model at the boundaries of the chosen domain
of interest, a high level of downscaling consistency is reached. To
define the initial state and the atmosphere at the domain boundaries, a
specific ``hybrid'' vertical interpolation from the coarse-resolution
GCM fields to the high-resolution mesoscale domain is used to ensure the
stability and the physical relevancy of the simulations. Used in
synoptic-scale mode with a cyclic domain wrapped around the planet, the
mesoscale model correctly replicates the main large-scale thermal
structure and the zonally propagating waves. The model diagnostics of
the near-surface pressure, wind, and temperature daily cycles in Chryse
Planitia are in accordance with the Viking and Pathfinder measurements.
Afternoon gustiness at the respective landing sites is adequately
accounted for on the condition that convective adjustment is turned off
in the mesoscale simulations. On the rims of Valles Marineris, intense
daytime anabatic (\~{}30 m s$^{-1}$) and nighttime katabatic (\~{}40 m
s$^{-1}$) winds are predicted. Within the canyon corridors,
topographical channeling can amplify the wind a few kilometers above the
ground, especially during the night. Through large-eddy simulations in
Gusev Crater, the model describes the mixing layer growth during the
afternoon, and the associated dynamics: convective motions, overlying
gravity waves, and dust devil-like vortices. Modeled temperature
profiles are in satisfactory agreement with the Miniature Thermal
Emission Spectrometer (Mini-TES) measurements. The ability of the model
to transport tracers at regional scales is exemplified by the model's
prediction for the altitude of the Tharsis topographical water ice
clouds in the afternoon. Finally, a nighttime ``warm ring'' at the base
of Olympus Mons is identified in the simulations, resulting from
adiabatic warming by the intense downslope winds along the flanks of the
volcano. The surface temperature enhancement reaches +20 K throughout
the night. Such a phenomenon may have adversely affected the thermal
inertia derivations in the region.
}},
doi = {10.1029/2008JE003242},
adsurl = {https://ui.adsabs.harvard.edu/abs/2009JGRE..114.2009S},
localpdf = {REF/2009JGRE..114.2009S.pdf},
adsnote = {Provided by the SAO/NASA Astrophysics Data System}
}