Global Lagrangian atmospheric dispersion model

The Global Lagrangian Atmospheric Dispersion Model (GLADIM) is described. GLADIM is based on the global trajectory model, which had been developed earlier and uses fields of weather parameters from different atmospheric reanalysis centers for calculations of trajectories of air mass that include trace gases. GLADIM includes the parameterization of turbulent diffusion and allows the forward calculation of concentrations of atmospheric tracers at nodes of a global regular grid when a source is specified. Thus, GLADIM can be used for the forward simulation of pollutant propagation (volcanic ash, radionuclides, and so on). Working in the reverse direction, GLADIM allows the detection of remote sources that mainly contribute to the tracer concentration at an observation point. This property of Lagrangian models is widely used for data analysis and the reverse modeling of emission sources of a pollutant specified. In this work we describe the model and some results of its validation through a comparison with results of a similar model and observation data

Adjoint of the Global Eulerian Lagrangian Coupled Atmospheric transport model (A-GELCA v1.0): development and validation

Abstract. We present the development of the Adjoint of the Global Eulerian–Lagrangian Coupled Atmospheric (A-GELCA) model that consists of the National Institute for Environmental Studies (NIES) model as an Eulerian three-dimensional transport model (TM), and FLEXPART (FLEXible PARTicle dispersion model) as the Lagrangian plume diffusion model (LPDM). The tangent and adjoint components of the Eulerian model were constructed directly from the original NIES TM code using an automatic differentiation tool known as TAF (Transformation of Algorithms in Fortran; http://www.FastOpt.com), with additional manual pre- and post-processing aimed at improving the performance of the computing, including MPI (Message Passing Interface). As results, the adjoint of Eulerian model is discrete. Construction of the adjoint of the Lagrangian component did not require any code modification, as LPDMs are able to track a significant number of particles back in time and thereby calculate the sensitivity of observations to the neighboring emissions areas. Eulerian and Lagrangian adjoint components were coupled at the time boundary in the global domain.The results are verified using a series of test experiments. The forward simulation shown the coupled model is effective in reproducing the seasonal cycle and short-term variability of CO2 even in the case of multiple limiting factors, such as high uncertainty of fluxes and the low resolution of the Eulerian model. The adjoint model demonstrates the high accuracy compared to direct forward sensitivity calculations and fast performance. The developed adjoint of the coupled model combines the flux conservation and stability of an Eulerian discrete adjoint formulation with the flexibility, accuracy, and high resolution of a Lagrangian backward trajectory formulation. </jats:p

Adjoint of the global Eulerian-Lagrangian coupled atmospheric transport model (A-GELCA v1.0): development and validation

We present the development of the Adjoint of the Global Eulerian–Lagrangian Coupled Atmospheric (A-GELCA) model that consists of the National Institute for Environmental Studies (NIES) model as an Eulerian three-dimensional transport model (TM), and FLEXPART (FLEXible PARTicle dispersion model) as the Lagrangian Particle Dispersion Model (LPDM). The forward tangent linear and adjoint components of the Eulerian model were constructed directly from the original NIES TM code using an automatic differentiation tool known as TAF (Transformation of Algorithms in Fortran; http://www.FastOpt.com), with additional manual pre- and post-processing aimed at improving transparency and clarity of the code and optimizing the performance of the computing, including MPI (Message Passing Interface). The Lagrangian component did not require any code modification, as LPDMs are self-adjoint and track a significant number of particles backward in time in order to calculate the sensitivity of the observations to the neighboring emission areas. The constructed Eulerian adjoint was coupled with the Lagrangian component at a time boundary in the global domain. The simulations presented in this work were performed using the A-GELCA model in forward and adjoint modes. The forward simulation shows that the coupled model improves reproduction of the seasonal cycle and short-term variability of CO2. Mean bias and standard deviation for five of the six Siberian sites considered decrease roughly by 1 ppm when using the coupled model. The adjoint of the Eulerian model was shown, through several numerical tests, to be very accurate (within machine epsilon with mismatch around to ±6 e−14) compared to direct forward sensitivity calculations. The developed adjoint of the coupled model combines the flux conservation and stability of an Eulerian discrete adjoint formulation with the flexibility, accuracy, and high resolution of a Lagrangian backward trajectory formulation. A-GELCA will be incorporated into a variational inversion system designed to optimize surface fluxes of greenhouse gases

Healing and Ritual Imagination in Chinese Medicine: The Multiple Interpretations of Zhuyou

In the Chinese medical corpus, ritual healing largely fell under the rubric of zhuyou 祝由 to uncover and expel the unknown, imperceptible, and occult causes of illness. Often dealing with uncertain or incurable cases, zhuyou remained at the cutting-edge of contemporary medicine. For a rising medical elite after the Northern Song, zhuyou was the branch of medicine to flexibly incorporate and critique the variety of ritual therapies into orthodox practice. Zhuyou employed prayer, incantations, talismans, gestures, and drugs in a nuanced clinical encounter to reveal the hidden root of disorder ranging from a blockage of qi, spirit possession, emotional imbalance, or loss of virtue. These rituals opened an imaginative space for therapeutic play where patients and healers could use spiritual proxies and props to address difficult emotions or issues that were often the hidden cause of affliction. The development of zhuyou also reflected the changing role of ritual in the history of Chinese medicine and the exchanges among physicians, Daoist priests, and other ritual healers. The significance of ritual in Chinese medical history has largely remained unclear as most editions of medical classics republished since the early twentieth century excise relevant chapters and zhuyou manuscripts, until recently, were uncatalogued

Global Lagrangian atmospheric dispersion model

The Global Lagrangian Atmospheric Dispersion Model (GLADIM) is described. GLADIM is based on the global trajectory model, which had been developed earlier and uses fields of weather parameters from different atmospheric reanalysis centers for calculations of trajectories of air mass that include trace gases. GLADIM includes the parameterization of turbulent diffusion and allows the forward calculation of concentrations of atmospheric tracers at nodes of a global regular grid when a source is specified. Thus, GLADIM can be used for the forward simulation of pollutant propagation (volcanic ash, radionuclides, and so on). Working in the reverse direction, GLADIM allows the detection of remote sources that mainly contribute to the tracer concentration at an observation point. This property of Lagrangian models is widely used for data analysis and the reverse modeling of emission sources of a pollutant specified. In this work we describe the model and some results of its validation through a comparison with results of a similar model and observation data

Column-averaged CO2 concentrations in the subarctic from GOSAT retrievals and NIES transport model simulations

AbstractThe distribution of atmospheric carbon dioxide (CO2) in the subarctic was investigated using the National Institute for Environmental Studies (NIES) three-dimensional transport model (TM) and retrievals from the Greenhouse gases Observing SATellite (GOSAT). Column-averaged dry air mole fractions of subarctic atmospheric CO2 (XCO2) from the NIES TM for four flux combinations were analyzed. Two flux datasets were optimized using only surface observations and two others were optimized using both surface and GOSAT Level 2 data. Two inverse modeling approaches using GOSAT data were compared. In the basic approach adopted in the GOSAT Level 4 product, the GOSAT observations are aggregated into monthly means over 5° × 5° grids. In the alternative method, the model–observation misfit is estimated for each observation separately. The XCO2 values simulated with optimized fluxes were validated against Total Carbon Column Observing Network (TCCON) ground-based high-resolution Fourier Transform Spectrometer (FTS) measurements. Optimized fluxes were applied to study XCO2 seasonal variability over the period 2009–2010 in the Arctic and subarctic regions. The impact on CO2 levels of emissions from enhancement of biospheric respiration induced by the high temperature and strong wildfires occurring in the summer of 2010 was analyzed. Use of GOSAT data has a substantial impact on estimates of the level of CO2 interanual variability