180 likes | 314 Views
The 2010 Workshop on the Solution of Partial Differential Equations on the Sphere. Design and evaluation of the National Institute for Environmental Studies (NIES) transport model. Dmitry Belikov and Shamil Maksyutov National Institute for Environmental Studies, Tsukuba, Japan.
E N D
The 2010 Workshop on the Solution of Partial Differential Equations on the Sphere Design and evaluation of the National Institute for Environmental Studies (NIES) transport model Dmitry Belikov and ShamilMaksyutov National Institute for Environmental Studies, Tsukuba, Japan
Model formulation New version of the NIES TM (NIES-08) with flux-form advection algorithms have been designed. Just like in the predecessor model with semi-Lagrangian algorithms (Maksyutov et al., 2008), we presented the atmospheric constituent transport equation in the Lagrangian-style form (Willamson and Laprise, 2000): PDEs 2010, Potsdam
A reduced latitude-longitude grid scheme The offline global tracer transport model version uses a reduced latitude-longitude grid scheme (Peterson et al., JGR, 1998), in which the sizes of grids are doubled several times approaching the poles Advantages vs. icosahedral, cubic grids – easy to bring reanalysis data, and formulate 2nd, 3rd order approximations PDEs 2010, Potsdam
Horizontal mass flux correction method The horizontal mass fluxes, derived from the spectral data (the output of weather forecast models) are balanced with the surface pressure tendency by adding correction fluxes, which is necessary to be determined (Heimann and Keeling, Geophys. Mon., 1989) The correction flux is calculated by transforming Equation into a Poisson equation, which is solved with a discrete 2D Fourier transform for every level l PDEs 2010, Potsdam
NIES TM transport algorithm test Initial tracer field velocities SML – Semi-Lagrangian (Maksyutov et al. 2008); VL – 3-rd order van Leer scheme (van Leer, 1977); Pr – Second Moments (Prather, 1986) PDEs 2010, Potsdam
NIES TM transport algorithm test SML – Semi-Lagrangian (Maksyutov et al. 2008), VL – 3-rd order van Leer scheme (van Leer, 1977), Pr – Second Moments (Prather, 1986) PDEs 2010, Potsdam
NIES TM meteorology data • The Japan Meteorological Agency (JMA) Japan Climate Data Assimilation System (JCDAS) meteorological dataset (6-hourly time step, resolution of 1.25×1.25 deg, 40 hybrid vertical levels). Height of planetary boundary layer with time step of 3 hours are taken from ECMWF Interim Reanalysis. • Global Point Value (GPV) - a special product prepared by the Japan Meteorological Agency Global Spectral Model (JMA-GSM) (3 hourly time step, resolution of 0.5 0.5 deg, 21 pressure levels). PDEs 2010, Potsdam
High resolution: 0.625 deg globally Simulated surface CO2 concentration around Japan at 21:00UTC, March 26, 2008 using NIES-08 with resolution 0.625 deg PDEs 2010, Potsdam
Hybrid sigma-pressure and sigma-isoentropic vertical coordinate systems A hybrid sigma-pressure and a sigma-isentropic vertical coordinate systems with 32 levels up to 2 mb PDEs 2010, Potsdam
Hybrid sigma-pressure and sigma-isoentropic vertical coordinate systems Pressure, hPa Mean age of air (SF6) simulated by the NIES-08 with sigma-pressure (left) and sigma-isoentropic (right) vertical coordinate systems PDEs 2010, Potsdam
NIES TM results (SF6) Zonally averaged annual mean of SF6 concentration simulated by NIES-08 Interhemispheric gradients of modeled and observed SF6 concentrations PDEs 2010, Potsdam
NIES TM results (CO2) Latitudinal distributions of CO2 seasonal amplitude at 35 GLOBALVIEW-CO2 (2008) sites. Seasonal amplitude is the difference between the maximum and the minimum of seasonal cycle. PDEs 2010, Potsdam
Convective parameterization scheme • Kuo-type cumulus parameterization (Grell, 1994) including entrainment and detrainment processes on convective updrafts and downdrafts proposed by Tiedtke (1989); • new method to determine cumulus convective updrafts where Pconv denotes the convective precipitation rate at the surface [kg/m2/sec], qbase is the absolute humidity at the cloud base [kg/kg]; PDEs 2010, Potsdam
Convective parameterization scheme NIES TM MERRA Seasonally average convective mass flux (g/m2/sec) from the NIES TM and Modern Era Retrospective-analysis For Research And Applications (MERRA) data for summer 2006 PDEs 2010, Potsdam
NIES TM results (222Rn) Without convective parameterization With new convective parameterization PDEs 2010, Potsdam
NIES TM results (222Rn) The model results are compared with data from insitu observations (Kritz et al., JGR, 1998; Liu et al., JGR, 1984; Zaucker et al., JGR , 1996) and the results obtained from model GAMIL (Zhang et al., ACP, 2008) PDEs 2010, Potsdam
Conclusion • Improvements in tracer transport simulation are achieved due to: • Mass conservative numerical algorithm and horizontal mass flux correction method; • A reduced latitude-longitude grid scheme; • Hybrid sigma-pressure and a sigma-isentropic vertical coordinate systems; • Convective parameterization scheme; • Future: • High-resolution meteorological data; PDEs 2010, Potsdam
Thank you! PDEs 2010, Potsdam