Assessing Tracer Transport Algorithms and the Impact of Vertical Resolution in a Finite-Volume Dynamical Core

2012 ◽  
Vol 140 (5) ◽  
pp. 1620-1638 ◽  
Author(s):  
James Kent ◽  
Christiane Jablonowski ◽  
Jared P. Whitehead ◽  
Richard B. Rood
Keyword(s):  
Finite Volume ◽  
Vertical Component ◽  
Vertical Motion ◽  
Atmospheric Model ◽  
Horizontal Resolution ◽  
Vertical Resolution ◽  
Tracer Transport ◽  
Dynamical Core ◽  
Lagrangian Coordinate ◽  
The Impact

Abstract Modeling the transport of trace gases is an essential part of any atmospheric model. The tracer transport scheme in the Community Atmosphere Model finite-volume dynamical core (CAM-FV), which is part of the National Center for Atmospheric Research’s (NCAR’s) Community Earth System Model (CESM1), is investigated using multidimensional idealized advection tests. CAM-FV’s tracer transport algorithm makes use of one-dimensional monotonic limiters. The Colella–Sekora limiter, which is applied to increase accuracy where the data are smooth, is implemented into the CAM-FV framework, and compared with the more traditional monotonic limiter of the piecewise parabolic method (the default limiter). For 2D flow, CAM-FV splits dimensions, allowing overshoots and undershoots, with the Colella–Sekora limiter producing larger overshoots than the default limiter. The impact of vertical resolution is also explored. A vertical Lagrangian coordinate is used in CAM-FV, and is periodically remapped back to a fixed Eulerian grid. For purely vertical motion, it is found that less-frequent remapping of the Lagrangian coordinate in CAM-FV improves results. For full 3D tests, the vertical component of the tracer transport dominates the error and limits the overall accuracy. If the vertical resolution is inadequate, increasing the horizontal resolution has almost no effect on accuracy. This is because the vertical resolution currently used in CAM version 5 may not be sufficiently fine enough to resolve some atmospheric tracers and provide accurate vertical advection. Idealized tests using tracers in a gravity wave agree with these results.

scholarly journals A Case Study of Observed and Modeled Barrier Flow in the Denmark Strait in May 2015

2017 ◽  
Vol 145 (6) ◽  
pp. 2385-2404 ◽  
Author(s):  
Alice K. DuVivier ◽  
John J. Cassano ◽  
Steven Greco ◽  
G. David Emmitt
Keyword(s):  
Wrf Model ◽  
Horizontal Resolution ◽  
Surface Energy Budget ◽  
Vertical Resolution ◽  
Structure Comparison ◽  
Denmark Strait ◽  
Pbl Scheme ◽  
Barrier Jets ◽  
The Impact ◽  
Jet Structure

Abstract Mesoscale barrier jets in the Denmark Strait are common in winter months and have the capability to influence open ocean convection. This paper presents the first detailed observational study of a summertime (21 May 2015) barrier wind event in the Denmark Strait using dropsondes and observations from an airborne Doppler wind lidar (DWL). The DWL profiles agree well with dropsonde observations and show a vertically narrow (~250–400 m) barrier jet of 23–28 m s−1 near the Greenland coast that broadens (~300–1000 m) and strengthens farther off coast. In addition, otherwise identical regional high-resolution Weather Research and Forecasting (WRF) Model simulations of the event are analyzed at four horizontal grid spacings (5, 10, 25, and 50 km), two vertical resolutions (40 and 60 levels), and two planetary boundary layer (PBL) parameterizations [Mellor–Yamada–Nakanishi–Niino, version 2.5 (MYNN2.5) and University of Washington (UW)] to determine what model configurations best simulate the observed jet structure. Comparison of the WRF simulations with wind observations from satellites, dropsondes, and the airborne DWL scans indicate that the combination of both high horizontal resolution (5 km) and vertical resolution (60 levels) best captures observed barrier jet structure and speeds as well as the observed cloud field, including some convective clouds. Both WRF PBL schemes produced reasonable barrier jets with the UW scheme slightly outperforming the MYNN2.5 scheme. However, further investigation at high horizontal and vertical resolution is needed to determine the impact of the WRF PBL scheme on surface energy budget terms, particularly in the high-latitude maritime environment around Greenland.

scholarly journals Extending the Modular Earth Submodel System (MESSy v2.54) model hierarchy: the ECHAM/MESSy IdeaLized (EMIL) model setup

2020 ◽  
Vol 13 (11) ◽  
pp. 5229-5257
Author(s):  
Hella Garny ◽  
Roland Walz ◽  
Matthias Nützel ◽  
Thomas Birner
Keyword(s):  
Climate Model ◽  
Polar Vortex ◽  
Equilibrium Temperature ◽  
Earth System ◽  
Tracer Transport ◽  
Dynamical Core ◽  
Vortex Strength ◽  
Model Hierarchy ◽  
Model Setup ◽  
The Impact

Abstract. As models of the Earth system grow in complexity, a need emerges to connect them with simplified systems through model hierarchies in order to improve process understanding. The Modular Earth Submodel System (MESSy) was developed to incorporate chemical processes into an Earth System model. It provides an environment to allow for model configurations and setups of varying complexity, and as of now the hierarchy ranges from a chemical box model to a fully coupled chemistry–climate model. Here, we present a newly implemented dry dynamical core model setup within the MESSy framework, denoted as ECHAM/MESSy IdeaLized (EMIL) model setup. EMIL is developed with the aim to provide an easily accessible idealized model setup that is consistently integrated in the MESSy model hierarchy. The implementation in MESSy further enables the utilization of diagnostic chemical tracers. The setup is achieved by the implementation of a new submodel for relaxation of temperature and horizontal winds to given background values, which replaces all other “physics” submodels in the EMIL setup. The submodel incorporates options to set the needed parameters (e.g., equilibrium temperature, relaxation time and damping coefficient) to functions used frequently in the past. This study consists of three parts. In the first part, test simulations with the EMIL model setup are shown to reproduce benchmarks provided by earlier dry dynamical core studies. In the second part, the sensitivity of the coupled troposphere–stratosphere dynamics to various modifications of the setup is studied. We find a non-linear response of the polar vortex strength to the prescribed meridional temperature gradient in the extratropical stratosphere that is indicative of a regime transition. In agreement with earlier studies, we find that the tropospheric jet moves poleward in response to the increase in the polar vortex strength but at a rate that strongly depends on the specifics of the setup. When replacing the idealized topography to generate planetary waves by mid-tropospheric wave-like heating, the response of the tropospheric jet to changes in the polar vortex is strongly damped in the free troposphere. However, near the surface, the jet shifts poleward at a higher rate than in the topographically forced simulations. Those results indicate that the wave-like heating might have to be used with care when studying troposphere–stratosphere coupling. In the third part, examples for possible applications of the model system are presented. The first example involves simulations with simplified chemistry to study the impact of dynamical variability and idealized changes on tracer transport, and the second example involves simulations of idealized monsoon circulations forced by localized heating. The ability to incorporate passive and chemically active tracers in the EMIL setup demonstrates the potential for future studies of tracer transport in the idealized dynamical model.

scholarly journals The Brazilian Global Atmospheric Model (BAM): Performance for Tropical Rainfall Forecasting and Sensitivity to Convective Scheme and Horizontal Resolution

2016 ◽  
Vol 31 (5) ◽  
pp. 1547-1572 ◽  
Author(s):  
Silvio N. Figueroa ◽  
José P. Bonatti ◽  
Paulo Y. Kubota ◽  
Georg A. Grell ◽  
Hugh Morrison ◽  
...  
Keyword(s):  
South America ◽  
Weather Prediction ◽  
Atmospheric Model ◽  
Horizontal Resolution ◽  
Southeastern Brazil ◽  
Rainfall Forecasting ◽  
Dynamical Core ◽  
Tropical Rainfall ◽  
Global Atmospheric Model ◽  
Convective Scheme

Abstract This article describes the main features of the Brazilian Global Atmospheric Model (BAM), analyses of its performance for tropical rainfall forecasting, and its sensitivity to convective scheme and horizontal resolution. BAM is the new global atmospheric model of the Center for Weather Forecasting and Climate Research [Centro de Previsão de Tempo e Estudos Climáticos (CPTEC)], which includes a new dynamical core and state-of-the-art parameterization schemes. BAM’s dynamical core incorporates a monotonic two-time-level semi-Lagrangian scheme, which is carried out completely on the model grid for the tridimensional transport of moisture, microphysical prognostic variables, and tracers. The performance of the quantitative precipitation forecasts (QPFs) from two convective schemes, the Grell–Dévényi (GD) scheme and its modified version (GDM), and two different horizontal resolutions are evaluated against the daily TRMM Multisatellite Precipitation Analysis over different tropical regions. Three main results are 1) the QPF skill was improved substantially with GDM in comparison to GD; 2) the increase in the horizontal resolution without any ad hoc tuning improves the variance of precipitation over continents with complex orography, such as Africa and South America, whereas over oceans there are no significant differences; and 3) the systematic errors (dry or wet biases) remain virtually unchanged for 5-day forecasts. Despite improvements in the tropical precipitation forecasts, especially over southeastern Brazil, dry biases over the Amazon and La Plata remain in BAM. Improving the precipitation forecasts over these regions remains a challenge for the future development of the model to be used not only for numerical weather prediction over South America but also for global climate simulations.

scholarly journals Impact of optimized mixing heights on simulated regional atmospheric transport of CO<sub>2</sub>

2014 ◽  
Vol 14 (14) ◽  
pp. 7149-7172 ◽  
Author(s):  
R. Kretschmer ◽  
C. Gerbig ◽  
U. Karstens ◽  
G. Biavati ◽  
A. Vermeulen ◽  
...  
Keyword(s):  
Dispersion Model ◽  
Synthetic Data ◽  
Wrf Model ◽  
Horizontal Resolution ◽  
Particle Dispersion ◽  
Random Errors ◽  
Tracer Transport ◽  
European Continent ◽  
Co2 Transport ◽  
The Impact

Abstract. The mixing height (MH) is a crucial parameter in commonly used transport models that proportionally affects air concentrations of trace gases with sources/sinks near the ground and on diurnal scales. Past synthetic data experiments indicated the possibility to improve tracer transport by minimizing errors of simulated MHs. In this paper we evaluate a method to constrain the Lagrangian particle dispersion model STILT (Stochastic Time-Inverted Lagrangian Transport) with MH diagnosed from radiosonde profiles using a bulk Richardson method. The same method was used to obtain hourly MHs for the period September/October 2009 from the Weather Research and Forecasting (WRF) model, which covers the European continent at 10 km horizontal resolution. Kriging with external drift (KED) was applied to estimate optimized MHs from observed and modelled MHs, which were used as input for STILT to assess the impact on CO2 transport. Special care has been taken to account for uncertainty in MH retrieval in this estimation process. MHs and CO2 concentrations were compared to vertical profiles from aircraft in situ data. We put an emphasis on testing the consistency of estimated MHs to observed vertical mixing of CO2. Modelled CO2 was also compared with continuous measurements made at Cabauw and Heidelberg stations. WRF MHs were significantly biased by ~10–20% during day and ~40–60% during night. Optimized MHs reduced this bias to ~5% with additional slight improvements in random errors. The KED MHs were generally more consistent with observed CO2 mixing. The use of optimized MHs had in general a favourable impact on CO2 transport, with bias reductions of 5–45% (day) and 60–90% (night). This indicates that a large part of the found CO2 model–data mismatch was indeed due to MH errors. Other causes for CO2 mismatch are discussed. Applicability of our method is discussed in the context of CO2 inversions at regional scales.

Assessing numerical impacts on stratospheric dynamics and transport using the age of air and the leaky pipe theoretical model

2020 ◽  
Author(s):  
Aman Gupta ◽  
Edwin Gerber ◽  
Peter Lauritzen
Keyword(s):  
Finite Volume ◽  
Climate Models ◽  
Trace Gases ◽  
Atmospheric Modeling ◽  
Numerical Representation ◽  
Vertical Resolution ◽  
Tracer Transport ◽  
Stratospheric Dynamics ◽  
High Vertical Resolution ◽  
Diabatic Forcing

&lt;p&gt;Accurate representation of tracer transport --- the movement of trace constituents by the atmospheric flow --- continues to be a challenge for climate models. Differences in the resolved circulation, biases due to physical parameterizations, and differences in the numerical representation of trace gases result in large variations in transport, even among state-of-the-art climate models. These differences result in disagreement among model projections of the evolution of stratospheric ozone throughout the 21st century particularly in the recovery of the Antarctic ozone hole. In addition to transport, the delicate momentum balance in the upper-troposphere and lower-stratosphere (UTLS) also presents a stiff challenge for model numerics, exposing the impacts of numerical dissipation, the resolution of waves, and the consequences of imperfect momentum conservation. Biases in this region impact the global circulation, e.g., influencing the extratropics jets and stratospheric polar vortices, and alter the transport and exchange of trace gases between and through the troposphere and stratosphere.&lt;/p&gt;&lt;p&gt;&amp;#160;&lt;/p&gt;&lt;p&gt;In this study, we compare 2 modern dynamical cores (dycores) that employ very different numerics: the cubed sphere finite volume (CSFV) core from GFDL and the spectral element (SE) core from NCAR-CAM5. We force these dycores using identical Held-Suarez diabatic forcing in the troposphere and Polvani-Kushner diabatic forcing in the stratosphere, varying the horizontal and vertical resolution. We observe significant differences in circulation, between the two models at high vertical resolution in the lower and middle tropical stratosphere. While the finite volume core is relatively insensitive to any changes in vertical resolution, the PS and SE dycores resolve considerably different tropical stratospheric dynamics at high vertical resolution (80 levels). These models develop QBO-like westerly winds in the tropics and induce a secondary meridional circulation in the tropical stratosphere, which sets of transport between the models. Using the theoretical leaky pipe transport model we analyze and separate out the transport differences due to differences is diabatic circulation and isentropic mixing and infer that this secondary circulation strikingly modulates stratospheric tracer transport (age of air) by altering the tropical-extratropical mixing, and impacts the extratropical circulation through the subtropical jets. Implications for comprehensive atmospheric modeling are discussed.&lt;/p&gt;

scholarly journals Climate Simulations with an Isentropic Finite-Volume Dynamical Core

2012 ◽  
Vol 25 (8) ◽  
pp. 2843-2861 ◽  
Author(s):  
Chih-Chieh Chen ◽  
Philip J. Rasch
Keyword(s):  
Finite Volume ◽  
Climate Models ◽  
Transport Processes ◽  
Climate Simulation ◽  
Modeling Framework ◽  
Dynamical Core ◽  
Vertical Coordinate ◽  
Upper Troposphere ◽  
Isentropic Model ◽  
The Impact

Abstract This paper discusses the impact of changing the vertical coordinate from a hybrid pressure to a hybrid-isentropic coordinate within the finite-volume (FV) dynamical core of the Community Atmosphere Model (CAM). Results from a 20-yr climate simulation using the new model coordinate configuration are compared to control simulations produced by the Eulerian spectral and FV dynamical cores of CAM, which both use a pressure-based (σ − P) coordinate. The same physical parameterization package is employed in all three dynamical cores. The isentropic modeling framework significantly alters the simulated climatology and has several desirable features. The revised model produces a better representation of heat transport processes in the atmosphere leading to much improved atmospheric temperatures. The authors show that the isentropic model is very effective in reducing the long-standing cold temperature bias in the upper troposphere and lower stratosphere, a deficiency shared among most climate models. The warmer upper troposphere and stratosphere seen in the isentropic model reduces the global coverage of high clouds, which is in better agreement with observations. The isentropic model also shows improvements in the simulated wintertime mean sea level pressure field in the Northern Hemisphere.

scholarly journals Sensitivity of Simulated Climate to Horizontal and Vertical Resolution in the ECHAM5 Atmosphere Model

2006 ◽  
Vol 19 (16) ◽  
pp. 3771-3791 ◽  
Author(s):  
E. Roeckner ◽  
R. Brokopf ◽  
M. Esch ◽  
M. Giorgetta ◽  
S. Hagemann ◽  
...  
Keyword(s):  
General Circulation ◽  
Circulation Model ◽  
Horizontal Resolution ◽  
Boreal Winter ◽  
Cold Bias ◽  
Vertical Resolution ◽  
Max Planck ◽  
Mean Climate ◽  
The Impact ◽  
Mean Square Errors

Abstract The most recent version of the Max Planck Institute for Meteorology atmospheric general circulation model, ECHAM5, is used to study the impact of changes in horizontal and vertical resolution on seasonal mean climate. In a series of Atmospheric Model Intercomparison Project (AMIP)-style experiments with resolutions ranging between T21L19 and T159L31, the systematic errors and convergence properties are assessed for two vertical resolutions. At low vertical resolution (L19) there is no evidence for convergence to a more realistic climate state for horizontal resolutions higher than T42. At higher vertical resolution (L31), on the other hand, the root-mean-square errors decrease monotonically with increasing horizontal resolution. Furthermore, except for T42, the L31 versions are superior to their L19 counterparts, and the improvements become more evident at increasingly higher horizontal resolutions. This applies, in particular, to the zonal mean climate state and to the stationary wave patterns in boreal winter. As in previous studies, increasing horizontal resolution leads to a warming of the troposphere, most prominently at midlatitudes, and to a poleward shift and intensification of the midlatitude westerlies. Increasing the vertical resolution has the opposite effect, almost independent of horizontal resolution. Whereas the atmosphere is colder at low and middle latitudes, it is warmer at high latitudes and close to the surface. In addition, increased vertical resolution results in a pronounced warming in the polar upper troposphere and lower stratosphere, where the cold bias is reduced by up to 50% compared to L19 simulations. Consistent with these temperature changes is a decrease and equatorward shift of the midlatitude westerlies. The substantial benefits in refining both horizontal and vertical resolution give some support to scaling arguments deduced from quasigeostrophic theory implying that horizontal and vertical resolution ought to be chosen consistently.

scholarly journals Ice injected into the tropopause by deep convection – Part 2: Over the Maritime Continent

2021 ◽  
Vol 21 (3) ◽  
pp. 2191-2210
Author(s):  
Iris-Amata Dion ◽  
Cyrille Dallet ◽  
Philippe Ricaud ◽  
Fabien Carminati ◽  
Thibaut Dauhut ◽  
...  
Keyword(s):  
Water Content ◽  
Deep Convection ◽  
Companion Paper ◽  
Horizontal Resolution ◽  
High Temporal Resolution ◽  
Vertical Resolution ◽  
Maritime Continent ◽  
Ice Water Content ◽  
Ice Water ◽  
The Impact

Abstract. The amount of ice injected into the tropical tropopause layer has a strong radiative impact on climate. A companion paper (Part 1) used the amplitude of the diurnal cycle of ice water content (IWC) as an estimate of ice injection by deep convection, showed that the Maritime Continent (MariCont) region provides the largest injection to the upper troposphere (UT; 146 hPa) and to the tropopause level (TL; 100 hPa). This study focuses on the MariCont region and extends that approach to assess the processes, the areas and the diurnal amount and duration of ice injected over islands and over seas during the austral convective season. The model presented in the companion paper is again used to estimate the amount of ice injected (ΔIWC) by combining ice water content (IWC) measured twice a day by the Microwave Limb Sounder (MLS; Version 4.2) from 2004 to 2017 and precipitation (Prec) measurements from the Tropical Rainfall Measurement Mission (TRMM; Version 007) binned at high temporal resolution (1 h). The horizontal distribution of ΔIWC estimated from Prec (ΔIWCPrec) is presented at 2∘×2∘ horizontal resolution over the MariCont. ΔIWC is also evaluated by using the number of lightning events (Flash) from the TRMM-LIS instrument (Lightning Imaging Sensor, from 2004 to 2015 at 1 h and 0.25∘ × 0.25∘ resolution). ΔIWCPrec and ΔIWC estimated from Flash (ΔIWCFlash) are compared to ΔIWC estimated from the ERA5 reanalyses (ΔIWCERA5) with the vertical resolution degraded to that of MLS observations (ΔIWCERA5). Our study shows that the diurnal cycles of Prec and Flash are consistent with each other in phase over land but different over offshore and coastal areas of the MariCont. The observational ΔIWC range between ΔIWCPrec and ΔIWCFlash, interpreted as the uncertainty of our model in estimating the amount of ice injected, is smaller over land (where ΔIWCPrec and ΔIWCFlash agree to within 22 %) than over ocean (where differences are up to 71 %) in the UT and TL. The impact of the MLS vertical resolution on the estimation of ΔIWC is greater in the TL (difference between ΔIWCERA5 and 〈ΔIWCERA5〉 of 32 % to 139 %, depending on the study zone) than in the UT (difference of 9 % to 33 %). Considering all the methods, in the UT, estimates of ΔIWC span 4.2 to 10.0 mg m−3 over land and 0.4 to 4.4 mg m−3 over sea, and in the TL estimates of ΔIWC span 0.5 to 3.9 mg m−3 over land and 0.1 to 0.7 mg m−3 over sea. Finally, based on IWC from MLS and ERA5, Prec and Flash, this study highlights that (1) at both levels, ΔIWC estimated over land can be more than twice that estimated over sea and (2) small islands with high topography present the largest ΔIWC (e.g., island of Java).

The impact of vertical resolution on tropical cyclone simulation

2020 ◽  
Author(s):  
Wei Huang ◽  
Mengjuan Liu ◽  
Xu Zhang ◽  
Jian-wen Bao
Keyword(s):  
Tropical Cyclone ◽  
Prediction Models ◽  
Weather Prediction ◽  
Vertical Mixing ◽  
Forecast Model ◽  
Horizontal Resolution ◽  
Model Solution ◽  
Vertical Resolution ◽  
Solution Convergence ◽  
The Impact

&lt;p&gt;It is well known that horizontal resolution has a great deal of impact on tropical cyclone simulations using numerical weather prediction models.&amp;#160; It is relatively less discussed in the literature how vertical resolution affects the solution convergence of tropical cyclone simulations.&amp;#160; In this study, the resolved kinetic energy spectrum, the Richardson number probability density function and resolved flow features are used as metrics to examine the behavior of solution convergence in tropical cyclone simulations using the Weather and Forecast Model (WRF).&amp;#160; It is found that for convective-scale simulations of a real tropical cyclone case with 3-km horizontal resolution, the model solution does not converge until a vertically stretched vertical resolution approaches 200 layers or more.&amp;#160; The results from this study confirm the results from a few previous studies that the subgrid turbulent mixing, particularly, the vertical mixing, plays a significant role in the behavior of model solution convergence with respect to vertical resolution.&amp;#160; They also provide a basis for the vertical grid configuration selection for the operational tropical cyclone model of Shanghai Meteorological Service.&lt;/p&gt;

scholarly journals Atmospheric transport and chemistry of trace gases in LMDz5B: evaluation and implications for inverse modelling

2015 ◽  
Vol 8 (2) ◽  
pp. 129-150 ◽  
Author(s):  
R. Locatelli ◽  
P. Bousquet ◽  
F. Hourdin ◽  
M. Saunois ◽  
A. Cozic ◽  
...  
Keyword(s):  
Greenhouse Gas Emissions ◽  
Greenhouse Gas ◽  
Atmospheric Transport ◽  
Trace Gases ◽  
Inverse Modelling ◽  
Trace Gas ◽  
Vertical Resolution ◽  
Tracer Transport ◽  
Gas Emissions ◽  
The Impact

Abstract. Representation of atmospheric transport is a major source of error in the estimation of greenhouse gas sources and sinks by inverse modelling. Here we assess the impact on trace gas mole fractions of the new physical parameterizations recently implemented in the atmospheric global climate model LMDz to improve vertical diffusion, mesoscale mixing by thermal plumes in the planetary boundary layer (PBL), and deep convection in the troposphere. At the same time, the horizontal and vertical resolution of the model used in the inverse system has been increased. The aim of this paper is to evaluate the impact of these developments on the representation of trace gas transport and chemistry, and to anticipate the implications for inversions of greenhouse gas emissions using such an updated model. Comparison of a one-dimensional version of LMDz with large eddy simulations shows that the thermal scheme simulates shallow convective tracer transport in the PBL over land very efficiently, and much better than previous versions of the model. This result is confirmed in three-dimensional simulations, by a much improved reproduction of the radon-222 diurnal cycle. However, the enhanced dynamics of tracer concentrations induces a stronger sensitivity of the new LMDz configuration to external meteorological forcings. At larger scales, the inter-hemispheric exchange is slightly slower when using the new version of the model, bringing them closer to observations. The increase in the vertical resolution (from 19 to 39 layers) significantly improves the representation of stratosphere/troposphere exchange. Furthermore, changes in atmospheric thermodynamic variables, such as temperature, due to changes in the PBL mixing modify chemical reaction rates, which perturb chemical equilibriums of reactive trace gases. One implication of LMDz model developments for future inversions of greenhouse gas emissions is the ability of the updated system to assimilate a larger amount of high-frequency data sampled at high-variability stations. Others implications are discussed at the end of the paper.

Sign in / Sign up

Export Citation Format

Share Document