scholarly journals Evaluation of oceanic and atmospheric trajectory schemes in the TRACMASS trajectory model v6.0

2016 ◽  
Author(s):  
Kristofer Döös ◽  
Bror Jönsson ◽  
Joakim Kjellsson

Abstract. Two different trajectory schemes for oceanic and atmospheric general circulation models are compared in two different experiments. The theories of the two trajectory schemes are presented showing the differential equations they solve and why they are mass conserving. One scheme assumes that the velocity fields are stationary for a limited period of time and solves the trajectory path from a differential equation only as a function of space, i.e. "stepwise stationary". The second scheme uses a continuous linear interpolation of the fields in time and solves the trajectory path from a differential equation as a function of both space and time, i.e. "time-dependent". A special case of the "stepwise-stationary" scheme, when velocities are assumed constant between GCM outputs, is also considered, named "fixed GCM time step". The trajectory schemes are tested "off-line", i.e. using the already integrated and stored velocity fields from a GCM. The first comparison of the schemes uses trajectories calculated using the velocity fields from an eddy-resolving ocean general circulation model in the Agulhas region. The second comparison uses trajectories calculated using the wind fields from an atmospheric reanalysis. The study shows that using the "time-dependent" scheme over the "stepwise-stationary" scheme greatly improves accuracy with only a small increase in computational time. It is also found that with decreasing time steps the "stepwise-stationary" scheme becomes more accurate but at increased computational cost. The "time-dependent" scheme is therefore preferred over the "stepwise-stationary" scheme. However, when averaging over large ensembles of trajectories the two schemes are comparable, as intrinsic variability dominates over numerical errors. The "fixed GCM time step" is found to be less accurate than the "stepwise-stationary" scheme, even when considering averages over large ensembles.

2017 ◽  
Vol 10 (4) ◽  
pp. 1733-1749 ◽  
Author(s):  
Kristofer Döös ◽  
Bror Jönsson ◽  
Joakim Kjellsson

Abstract. Three different trajectory schemes for oceanic and atmospheric general circulation models are compared in two different experiments. The theories of the trajectory schemes are presented showing the differential equations they solve and why they are mass conserving. One scheme assumes that the velocity fields are stationary for set intervals of time between saved model outputs and solves the trajectory path from a differential equation only as a function of space, i.e. stepwise stationary. The second scheme is a special case of the stepwise-stationary scheme, where velocities are assumed constant between general circulation model (GCM) outputs; it uses hence a fixed GCM time step. The third scheme uses a continuous linear interpolation of the fields in time and solves the trajectory path from a differential equation as a function of both space and time, i.e. a time-dependent scheme. The trajectory schemes are tested offline, i.e. using the already integrated and stored velocity fields from a GCM. The first comparison of the schemes uses trajectories calculated using the velocity fields from a high-resolution ocean general circulation model in the Agulhas region. The second comparison uses trajectories calculated using the wind fields from an atmospheric reanalysis. The study shows that using the time-dependent scheme over the stepwise-stationary scheme greatly improves accuracy with only a small increase in computational time. It is also found that with decreasing time steps the stepwise-stationary scheme becomes increasingly more accurate but at increased computational cost. The time-dependent scheme is therefore preferred over the stepwise-stationary scheme. However, when averaging over large ensembles of trajectories, the two schemes are comparable, as intrinsic variability dominates over numerical errors. The fixed GCM time step scheme is found to be less accurate than the stepwise-stationary scheme, even when considering averages over large ensembles.


2014 ◽  
Vol 7 (5) ◽  
pp. 1961-1977 ◽  
Author(s):  
H. Wan ◽  
P. J. Rasch ◽  
K. Zhang ◽  
Y. Qian ◽  
H. Yan ◽  
...  

Abstract. This paper explores the feasibility of an experimentation strategy for investigating sensitivities in fast components of atmospheric general circulation models. The basic idea is to replace the traditional serial-in-time long-term climate integrations by representative ensembles of shorter simulations. The key advantage of the proposed method lies in its efficiency: since fewer days of simulation are needed, the computational cost is less, and because individual realizations are independent and can be integrated simultaneously, the new dimension of parallelism can dramatically reduce the turnaround time in benchmark tests, sensitivities studies, and model tuning exercises. The strategy is not appropriate for exploring sensitivity of all model features, but it is very effective in many situations. Two examples are presented using the Community Atmosphere Model, version 5. In the first example, the method is used to characterize sensitivities of the simulated clouds to time-step length. Results show that 3-day ensembles of 20 to 50 members are sufficient to reproduce the main signals revealed by traditional 5-year simulations. A nudging technique is applied to an additional set of simulations to help understand the contribution of physics–dynamics interaction to the detected time-step sensitivity. In the second example, multiple empirical parameters related to cloud microphysics and aerosol life cycle are perturbed simultaneously in order to find out which parameters have the largest impact on the simulated global mean top-of-atmosphere radiation balance. It turns out that 12-member ensembles of 10-day simulations are able to reveal the same sensitivities as seen in 4-year simulations performed in a previous study. In both cases, the ensemble method reduces the total computational time by a factor of about 15, and the turnaround time by a factor of several hundred. The efficiency of the method makes it particularly useful for the development of high-resolution, costly, and complex climate models.


2014 ◽  
Vol 7 (5) ◽  
pp. 2181-2191 ◽  
Author(s):  
D. Rossi ◽  
A. Maurizi

Abstract. The development and validation of the vertical diffusion module of IL-GLOBO, a Lagrangian transport model coupled online with the Eulerian general circulation model GLOBO, is described. The module simulates the effects of turbulence on particle motion by means of a Lagrangian stochastic model (LSM) consistently with the turbulent diffusion equation used in GLOBO. The implemented LSM integrates particle trajectories, using the native σ-hybrid coordinates of the Eulerian component, and fulfils the well-mixed condition (WMC) in the general case of a variable density profile. The module is validated through a series of 1-D offline numerical experiments by assessing its accuracy in maintaining an initially well-mixed distribution in the vertical. A dynamical time-step selection algorithm with constraints related to the shape of the diffusion coefficient profile is developed and discussed. Finally, the skills of a linear interpolation and a modified Akima spline interpolation method are compared, showing that both satisfy the WMC with significant differences in computational time. A preliminary run of the fully integrated 3-D model confirms the result only for the Akima interpolation scheme while the linear interpolation does not satisfy the WMC with a reasonable choice of the minimum integration time step.


Author(s):  
Daisuke Matsuoka ◽  
Fumiaki Araki ◽  
Hideharu Sasaki

Numerical study of ocean eddies has been carried out by using high-resolution ocean general circulation models. In order to understand ocean eddies from the large volume data produced by simulations, visualizing only eddy distribution at each time step is insufficient; time-variations in eddy events and phenomena must also be considered. However, existing methods cannot precisely find and track eddy events such as amalgamation and bifurcation. In this study, we propose an original approach for eddy detection, tracking, and event visualization based on an eddy classification system. The proposed method detects streams and currents as well as eddies, and it classifies discovered eddies into several categories using the additional stream and current information. By tracking how the classified eddies vary over time, detecting events such as eddy amalgamation and bifurcation as well as the interaction between eddies and ocean currents becomes achievable. We adopt the proposed method for two ocean areas in which strong ocean currents exist as case studies. We visualize the detected eddies and events in a time series of images, allowing us to acquire an intuitive understanding of a region of interest concealed in a high-resolution data set. Furthermore, our proposed method succeeded in clarifying the occurrence place and seasonality of each type of eddy event.


2014 ◽  
Vol 7 (6) ◽  
pp. 7559-7573
Author(s):  
S. Watanabe ◽  
K. Sato ◽  
Y. Kawatani ◽  
M. Takahashi

Abstract. The dependence of the gravity wave spectra of energy and momentum flux on the horizontal resolution and time step of atmospheric general circulation models (AGCMs) has been thoroughly investigated in the past. In contrast, much less attention has been given to the dependence of these gravity wave parameters on models' vertical resolutions. The present study demonstrates the dependence of gravity wave momentum flux in the stratosphere and mesosphere on the model's vertical resolution, which is evaluated using an AGCM with a horizontal resolution of about 0.56°. We performed a series of sensitivity test simulations changing only the model's vertical resolution above a height of 8 km, and found that inertial gravity waves with short vertical wavelengths simulated at higher vertical resolutions likely play an important role in determining the gravity wave momentum flux in the stratosphere and mesosphere.


2014 ◽  
Vol 7 (2) ◽  
pp. 2173-2216
Author(s):  
H. Wan ◽  
P. J. Rasch ◽  
K. Zhang ◽  
Y. Qian ◽  
H. Yan ◽  
...  

Abstract. This paper explores the feasibility of an experimentation strategy for investigating sensitivities in fast components of atmospheric general circulation models. The basic idea is to replace the traditional serial-in-time long-term climate integrations by representative ensembles of shorter simulations. The key advantage of the proposed method lies in its efficiency: since fewer days of simulation are needed, the computational cost is less, and because individual realizations are independent and can be integrated simultaneously, the new dimension of parallelism can dramatically reduce the turnaround time in benchmark tests, sensitivities studies, and model tuning exercises. The strategy is not appropriate for exploring sensitivity of all model features, but it is very effective in many situations. Two examples are presented using the Community Atmosphere Model version 5. The first example demonstrates that the method is capable of characterizing the model cloud and precipitation sensitivity to time step length. A nudging technique is also applied to an additional set of simulations to help understand the contribution of physics-dynamics interaction to the detected time step sensitivity. In the second example, multiple empirical parameters related to cloud microphysics and aerosol lifecycle are perturbed simultaneously in order to explore which parameters have the largest impact on the simulated global mean top-of-atmosphere radiation balance. Results show that in both examples, short ensembles are able to correctly reproduce the main signals of model sensitivities revealed by traditional long-term climate simulations for fast processes in the climate system. The efficiency of the ensemble method makes it particularly useful for the development of high-resolution, costly and complex climate models.


2008 ◽  
Vol 21 (15) ◽  
pp. 3642-3659 ◽  
Author(s):  
Hugh Morrison ◽  
Andrew Gettelman

Abstract A new two-moment stratiform cloud microphysics scheme in a general circulation model is described. Prognostic variables include cloud droplet and cloud ice mass mixing ratios and number concentrations. The scheme treats several microphysical processes, including hydrometeor collection, condensation/evaporation, freezing, melting, and sedimentation. The activation of droplets on aerosol is physically based and coupled to a subgrid vertical velocity. Unique aspects of the scheme, relative to existing two-moment schemes developed for general circulation models, are the diagnostic treatment of rain and snow number concentration and mixing ratio and the explicit treatment of subgrid cloud water variability for calculation of the microphysical process rates. Numerical aspects of the scheme are described in detail using idealized one-dimensional offline tests of the microphysics. Sensitivity of the scheme to time step, vertical resolution, and numerical method for diagnostic precipitation is investigated over a range of conditions. It is found that, in general, two substeps are required for numerical stability and reasonably small time truncation errors using a time step of 20 min; however, substepping is only required for the precipitation microphysical processes rather than the entire scheme. A new numerical approach for the diagnostic rain and snow produces reasonable results compared to a benchmark simulation, especially at low vertical resolution. Part II of this study details results of the scheme in single-column and global simulations, including comparison with observations.


2010 ◽  
Vol 67 (6) ◽  
pp. 1983-1995 ◽  
Author(s):  
Steven C. Hardiman ◽  
David G. Andrews ◽  
Andy A. White ◽  
Neal Butchart ◽  
Ian Edmond

Abstract Transformed Eulerian mean (TEM) equations and Eliassen–Palm (EP) flux diagnostics are presented for the general nonhydrostatic, fully compressible, deep atmosphere formulation of the primitive equations in spherical geometric coordinates. The TEM equations are applied to a general circulation model (GCM) based on these general primitive equations. It is demonstrated that a naive application in this model of the widely used approximations to the EP diagnostics, valid for the hydrostatic primitive equations using log-pressure as a vertical coordinate and presented, for example, by Andrews et al. in 1987 can lead to misleading features in these diagnostics. These features can be of the same order of magnitude as the diagnostics themselves throughout the winter stratosphere. Similar conclusions are found to hold for “downward control” calculations. The reasons are traced to the change of vertical coordinate from geometric height to log-pressure. Implications for the modeling community, including comparison of model output with that from reanalysis products available only on pressure surfaces, are discussed.


2018 ◽  
Vol 22 (10) ◽  
pp. 1-22 ◽  
Author(s):  
Andrew R. Bock ◽  
Lauren E. Hay ◽  
Gregory J. McCabe ◽  
Steven L. Markstrom ◽  
R. Dwight Atkinson

Abstract The accuracy of statistically downscaled (SD) general circulation model (GCM) simulations of monthly surface climate for historical conditions (1950–2005) was assessed for the conterminous United States (CONUS). The SD monthly precipitation (PPT) and temperature (TAVE) from 95 GCMs from phases 3 and 5 of the Coupled Model Intercomparison Project (CMIP3 and CMIP5) were used as inputs to a monthly water balance model (MWBM). Distributions of MWBM input (PPT and TAVE) and output [runoff (RUN)] variables derived from gridded station data (GSD) and historical SD climate were compared using the Kolmogorov–Smirnov (KS) test For all three variables considered, the KS test results showed that variables simulated using CMIP5 generally are more reliable than those derived from CMIP3, likely due to improvements in PPT simulations. At most locations across the CONUS, the largest differences between GSD and SD PPT and RUN occurred in the lowest part of the distributions (i.e., low-flow RUN and low-magnitude PPT). Results indicate that for the majority of the CONUS, there are downscaled GCMs that can reliably simulate historical climatic conditions. But, in some geographic locations, none of the SD GCMs replicated historical conditions for two of the three variables (PPT and RUN) based on the KS test, with a significance level of 0.05. In these locations, improved GCM simulations of PPT are needed to more reliably estimate components of the hydrologic cycle. Simple metrics and statistical tests, such as those described here, can provide an initial set of criteria to help simplify GCM selection.


2014 ◽  
Vol 27 (24) ◽  
pp. 9197-9213 ◽  
Author(s):  
Michael Horn ◽  
Kevin Walsh ◽  
Ming Zhao ◽  
Suzana J. Camargo ◽  
Enrico Scoccimarro ◽  
...  

Abstract Future tropical cyclone activity is a topic of great scientific and societal interest. In the absence of a climate theory of tropical cyclogenesis, general circulation models are the primary tool available for investigating the issue. However, the identification of tropical cyclones in model data at moderate resolution is complex, and numerous schemes have been developed for their detection. The influence of different tracking schemes on detected tropical cyclone activity and responses in the Hurricane Working Group experiments is examined herein. These are idealized atmospheric general circulation model experiments aimed at determining and distinguishing the effects of increased sea surface temperature and other increased CO2 effects on tropical cyclone activity. Two tracking schemes are applied to these data and the tracks provided by each modeling group are analyzed. The results herein indicate moderate agreement between the different tracking methods, with some models and experiments showing better agreement across schemes than others. When comparing responses between experiments, it is found that much of the disagreement between schemes is due to differences in duration, wind speed, and formation-latitude thresholds. After homogenization in these thresholds, agreement between different tracking methods is improved. However, much disagreement remains, accountable for by more fundamental differences between the tracking schemes. The results indicate that sensitivity testing and selection of objective thresholds are the key factors in obtaining meaningful, reproducible results when tracking tropical cyclones in climate model data at these resolutions, but that more fundamental differences between tracking methods can also have a significant impact on the responses in activity detected.


Sign in / Sign up

Export Citation Format

Share Document