Chonggang Xu, Bradley Christoffersen, Zachary Robbins, Ryan Knox, Rosie A. Fisher, Rutuja Chitra-Tarak, Martijn Slot, Kurt Solander, Lara Kueppers, Charles Koven, Nate McDowell
{"title":"Quantification of hydraulic trait control on plant hydrodynamics and risk of hydraulic failure within a demographic structured vegetation model in a tropical forest (FATES–HYDRO V1.0)","authors":"Chonggang Xu, Bradley Christoffersen, Zachary Robbins, Ryan Knox, Rosie A. Fisher, Rutuja Chitra-Tarak, Martijn Slot, Kurt Solander, Lara Kueppers, Charles Koven, Nate McDowell","doi":"10.5194/gmd-16-6267-2023","DOIUrl":"https://doi.org/10.5194/gmd-16-6267-2023","url":null,"abstract":"Abstract. Vegetation plays a key role in the global carbon cycle and thus is an important component within Earth system models (ESMs) that project future climate. Many ESMs are adopting methods to resolve plant size and ecosystem disturbance history, using vegetation demographic models. These models make it feasible to conduct more realistic simulation of processes that control vegetation dynamics. Meanwhile, increasing understanding of the processes governing plant water use, and ecosystem responses to drought in particular, has led to the adoption of dynamic plant water transport (i.e., hydrodynamic) schemes within ESMs. However, the extent to which variations in plant hydraulic traits affect both plant water stress and the risk of mortality in trait-diverse tropical forests is understudied. In this study, we report on a sensitivity analysis of an existing hydrodynamic scheme (HYDRO) model that is updated and incorporated into the Functionally Assembled Terrestrial Ecosystem Simulator (FATES) (FATES–HYDRO V1.0). The size- and canopy-structured representation within FATES is able to simulate how plant size and hydraulic traits affect vegetation dynamics and carbon–water fluxes. To better understand this new model system, and its functionality in tropical forest systems in particular, we conducted a global parameter sensitivity analysis at Barro Colorado Island, Panama. We assembled 942 observations of plant hydraulic traits on 306 tropical plant species for stomata, leaves, stems, and roots and determined the best-fit statistical distribution for each trait, which was used in model parameter sampling to assess the parametric sensitivity. We showed that, for simulated leaf water potential and loss of hydraulic conductivity across different plant organs, the four most important traits were associated with xylem conduit taper (buffers increasing hydraulic resistance with tree height), stomatal sensitivity to leaf water potential, maximum stem hydraulic conductivity, and the partitioning of total hydraulic resistance above vs. belowground. Our analysis of individual ensemble members revealed that trees at a high risk of hydraulic failure and potential tree mortality generally have a lower conduit taper, lower maximum xylem conductivity, lower stomatal sensitivity to leaf water potential, and lower resistance to xylem embolism for stem and transporting roots. We expect that our results will provide guidance on future modeling studies using plant hydrodynamic models to predict the forest responses to droughts and future field campaigns that aim to better parameterize plant hydrodynamic models.","PeriodicalId":12799,"journal":{"name":"Geoscientific Model Development","volume":"104 11","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-11-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135540071","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Mattia de' Michieli Vitturi, Tomaso Esposti Ongaro, Samantha Engwell
{"title":"IMEX_SfloW2D v2: a depth-averaged numerical flow model for volcanic gas–particle flows over complex topographies and water","authors":"Mattia de' Michieli Vitturi, Tomaso Esposti Ongaro, Samantha Engwell","doi":"10.5194/gmd-16-6309-2023","DOIUrl":"https://doi.org/10.5194/gmd-16-6309-2023","url":null,"abstract":"Abstract. We present developments to the physical model and the open-source numerical code IMEX_SfloW2D (de' Michieli Vitturi et al., 2019). These developments consist of a generalization of the depth-averaged (shallow-water) fluid equations to describe a polydisperse fluid–solid mixture, including terms for sedimentation and entrainment, transport equations for solid particles of different sizes, transport equations for different components of the carrier phase, and an equation for temperature/energy. Of relevance for the simulation of volcanic mass flows, vaporization and entrainment of water are implemented in the new model. The model can be easily adapted to simulate a wide range of volcanic mass flows (pyroclastic avalanches, lahars, pyroclastic surges), and here we present its application to transient dilute pyroclastic density currents (PDCs). The numerical algorithm and the code have been improved to allow for simulation of sub- to supercritical regimes and to simplify the setting of initial and boundary conditions. The code is open-source. The results of synthetic numerical benchmarks demonstrate the robustness of the numerical code in simulating transcritical flows interacting with the topography. Moreover, they highlight the importance of simulating transient in comparison to steady-state flows and flows in 2D versus 1D. Finally, we demonstrate the model capabilities to simulate a complex natural case involving the propagation of PDCs over the sea surface and across topographic obstacles, through application to Krakatau volcano, showing the relevance, at a large scale, of non-linear fluid dynamic features, such as hydraulic jumps and von Kármán vortices, to flow conditions such as velocity and runout.","PeriodicalId":12799,"journal":{"name":"Geoscientific Model Development","volume":"1 4","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-11-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135539366","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Xinzhu Yu, Li Liu, Chao Sun, Qingu Jiang, Biao Zhao, Zhiyuan Zhang, Hao Yu, Bin Wang
{"title":"CIOFC1.0: a common parallel input/output framework based on C-Coupler2.0","authors":"Xinzhu Yu, Li Liu, Chao Sun, Qingu Jiang, Biao Zhao, Zhiyuan Zhang, Hao Yu, Bin Wang","doi":"10.5194/gmd-16-6285-2023","DOIUrl":"https://doi.org/10.5194/gmd-16-6285-2023","url":null,"abstract":"Abstract. As earth system modeling develops ever finer grid resolutions, the inputting and outputting (I/O) of the increasingly large data fields becomes a processing bottleneck. Many models developed in China, as well as the community coupler (C-Coupler), do not fully benefit from existing parallel I/O supports. This paper reports the design and implementation of a common parallel input/output framework (CIOFC1.0) based on C-Coupler2.0. The CIOFC1.0 framework can accelerate the I/O of large data fields by parallelizing data read/write operations among processes. The framework also allows convenient specification by users of the I/O settings, e.g., the data fields for I/O, the time series of the data files for I/O, and the data grids in the files. The framework can also adaptively input data fields from a time series dataset during model integration, automatically interpolate data when necessary, and output fields either periodically or irregularly. CIOFC1.0 demonstrates the cooperative development of an I/O framework and coupler, and thus enables convenient and simultaneous use of a coupler and an I/O framework.","PeriodicalId":12799,"journal":{"name":"Geoscientific Model Development","volume":"10 38","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-11-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135540468","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Ewa M. Bednarz, Ryan Hossaini, N. Luke Abraham, Martyn P. Chipperfield
{"title":"Description and evaluation of the new UM–UKCA (vn11.0) Double Extended Stratospheric–Tropospheric (DEST vn1.0) scheme for comprehensive modelling of halogen chemistry in the stratosphere","authors":"Ewa M. Bednarz, Ryan Hossaini, N. Luke Abraham, Martyn P. Chipperfield","doi":"10.5194/gmd-16-6187-2023","DOIUrl":"https://doi.org/10.5194/gmd-16-6187-2023","url":null,"abstract":"Abstract. The paper describes the development and performance of the Double Extended Stratospheric–Tropospheric (DEST vn1.0) chemistry scheme, which forms a part of the Met Office's Unified Model coupled to the United Kingdom Chemistry and Aerosol (UM–UKCA) chemistry–climate model, which is the atmospheric composition model of the United Kingdom Earth System Model (UKESM). The scheme extends the standard Stratospheric–Tropospheric chemistry scheme (StratTrop) by including a range of important updates to the halogen chemistry. These allow process-oriented studies of stratospheric ozone depletion and recovery, including the impacts from both controlled long-lived ozone-depleting substances (ODSs) and emerging issues around uncontrolled very short-lived substances (VSLS). The main updates in DEST are (i) an explicit treatment of 14 of the most important long-lived ODSs; (ii) an inclusion of brominated VSLS (Br-VSLS) emissions and chemistry; and (iii) an inclusion of chlorinated VSLS (Cl-VSLS) emissions/LBCs (lower boundary conditions) and chemistry. We evaluate the scheme's performance by comparing DEST simulations against analogous runs made with the standard StratTrop scheme and against observational and reanalysis datasets. Overall, our scheme addresses some significant shortcomings in the representation of atmospheric halogens in the standard StratTrop scheme and will thus be particularly relevant for studies of ozone layer recovery and processes affecting it, in support of future World Meteorological Organization (WMO) Ozone Assessment Reports.","PeriodicalId":12799,"journal":{"name":"Geoscientific Model Development","volume":"2 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-11-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135874943","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Overcoming computational challenges to realize meter- to submeter-scale resolution in cloud simulations using the super-droplet method","authors":"Toshiki Matsushima, Seiya Nishizawa, Shin-ichiro Shima","doi":"10.5194/gmd-16-6211-2023","DOIUrl":"https://doi.org/10.5194/gmd-16-6211-2023","url":null,"abstract":"Abstract. A particle-based cloud model was developed for meter- to submeter-scale-resolution simulations of warm clouds. Simplified cloud microphysics schemes have already made meter-scale-resolution simulations feasible; however, such schemes are based on empirical assumptions, and hence they contain huge uncertainties. The super-droplet method (SDM) is a promising candidate for cloud microphysical process modeling and is a particle-based approach, making fewer assumptions for the droplet size distributions. However, meter-scale-resolution simulations using the SDM are not feasible even on existing high-end supercomputers because of high computational cost. In the present study, we overcame challenges to realize such simulations. The contributions of our work are as follows: (1) the uniform sampling method is not suitable when dealing with a large number of super-droplets (SDs). Hence, we developed a new initialization method for sampling SDs from a real droplet population. These SDs can be used for simulating spatial resolutions between meter and submeter scales. (2) We optimized the SDM algorithm to achieve high performance by reducing data movement and simplifying loop bodies using the concept of effective resolution. The optimized algorithms can be applied to a Fujitsu A64FX processor, and most of them are also effective on other many-core CPUs and possibly graphics processing units (GPUs). Warm-bubble experiments revealed that the throughput of particle calculations per second for the improved algorithms is 61.3 times faster than those for the original SDM. In the case of shallow cumulous, the simulation time when using the new SDM with 32–64 SDs per cell is shorter than that of a bin method with 32 bins and comparable to that of a two-moment bulk method. (3) Using the supercomputer Fugaku, we demonstrated that a numerical experiment with 2 m resolution and 128 SDs per cell covering 13 8242×3072 m3 domain is possible. The number of grid points and SDs are 104 and 442 times, respectively, those of the highest-resolution simulation performed so far. Our numerical model exhibited 98 % weak scaling for 36 864 nodes, accounting for 23 % of the total system. The simulation achieves 7.97 PFLOPS, 7.04 % of the peak ratio for overall performance, and a simulation time for SDM of 2.86×1013 particle ⋅ steps per second. Several challenges, such as incorporating mixed-phase processes, inclusion of terrain, and long-time integrations, remain, and our study will also contribute to solving them. The developed model enables us to study turbulence and microphysics processes over a wide range of scales using combinations of direct numerical simulation (DNS), laboratory experiments, and field studies. We believe that our approach advances the scientific understanding of clouds and contributes to reducing the uncertainties of weather simulation and climate projection.","PeriodicalId":12799,"journal":{"name":"Geoscientific Model Development","volume":"15 17","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-11-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135973202","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"A robust error correction method for numerical weather prediction wind speed based on Bayesian optimization, variational mode decomposition, principal component analysis, and random forest: VMD-PCA-RF (version 1.0.0)","authors":"Shaohui Zhou, Chloe Yuchao Gao, Zexia Duan, Xingya Xi, Yubin Li","doi":"10.5194/gmd-16-6247-2023","DOIUrl":"https://doi.org/10.5194/gmd-16-6247-2023","url":null,"abstract":"Abstract. Accurate wind speed prediction is crucial for the safe and efficient utilization of wind resources. However, current single-value deterministic numerical weather prediction methods employed by wind farms do not adequately meet the actual needs of power grid dispatching. In this study, we propose a new hybrid forecasting method for correcting 10 m wind speed predictions made by the Weather Research and Forecasting (WRF) model. Our approach incorporates variational mode decomposition (VMD), principal component analysis (PCA), and five artificial intelligence algorithms: deep belief network (DBN), multilayer perceptron (MLP), random forest (RF), eXtreme gradient boosting (XGBoost), light gradient boosting machine (lightGBM), and the Bayesian optimization algorithm (BOA). We first predict wind speeds using the WRF model, with initial and lateral boundary conditions from the Global Forecast System (GFS). We then perform two sets of experiments with different input factors and apply BOA optimization to tune the four artificial intelligence models, ultimately building the final models. Furthermore, we compare the aforementioned five optimal artificial intelligence models suitable for five provinces in southern China in the wintertime: VMD-PCA-RF in December 2021 and VMD-PCA-lightGBM in January 2022. We find that the VMD-PCA-RF evaluation indices exhibit relative stability over nearly a year: the correlation coefficient (R) is above 0.6, forecasting accuracy (FA) is above 85 %, mean absolute error (MAE) is below 0.6 m s−1, root mean square error (RMSE) is below 0.8 m s−1, relative mean absolute error (rMAE) is below 60 %, and relative root mean square error (rRMSE) is below 75 %. Thus, for its promising performance and excellent year-round robustness, we recommend adopting the proposed VMD-PCA-RF method for improved wind speed prediction in models.","PeriodicalId":12799,"journal":{"name":"Geoscientific Model Development","volume":"16 3","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-11-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135973508","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Simone Tilmes, Michael J. Mills, Yunqian Zhu, Charles G. Bardeen, Francis Vitt, Pengfei Yu, David Fillmore, Xiaohong Liu, Brian Toon, Terry Deshler
{"title":"Description and performance of a sectional aerosol microphysical model in the Community Earth System Model (CESM2)","authors":"Simone Tilmes, Michael J. Mills, Yunqian Zhu, Charles G. Bardeen, Francis Vitt, Pengfei Yu, David Fillmore, Xiaohong Liu, Brian Toon, Terry Deshler","doi":"10.5194/gmd-16-6087-2023","DOIUrl":"https://doi.org/10.5194/gmd-16-6087-2023","url":null,"abstract":"Abstract. We implemented the Community Aerosol and Radiation Model for Atmospheres (CARMA) in both the high- and low-top model versions of the Community Earth System Model Version 2 (CESM2). CARMA is a sectional microphysical model, which we use for aerosol in both the troposphere and stratosphere. CARMA is fully coupled to chemistry, clouds, radiation, and transport routines in CESM2. This development enables the comparison of simulations with a sectional (CARMA) and a modal (MAM4) aerosol microphysical model in the same modeling framework. The new implementation of CARMA has been adopted from previous work, with some additions that align with the current CESM2 Modal Aerosol Model (MAM4) implementation. The main updates include an interactive secondary organic aerosol description in CARMA, using the volatility basis set (VBS) approach, updated wet removal, and the use of transient emissions of aerosols and trace gases. In addition, we implemented an alternative aerosol nucleation scheme in CARMA, which is also used in MAM4. Detailed comparisons of stratospheric aerosol properties after the Mount Pinatubo eruption reveal the importance of prescribing sulfur injections in a larger region rather than in a single column to better represent the observed evolution of aerosols. Both CARMA and MAM4 in CESM2 are able to represent stratospheric and tropospheric aerosol properties reasonably well when compared to observations. Several differences in the performance of the two aerosol models show, in general, an improved representation of aerosols when using the sectional aerosol model in CESM2. These include a better representation of the aerosol size distribution after the Mount Pinatubo volcanic eruption in CARMA compared to MAM4. MAM4 produces on average smaller aerosols and less removal than CARMA, which results in a larger total mass. Both CARMA and MAM4 reproduce the stratospheric aerosol optical depth (AOD) within the error bar of the observations between 2001 and 2020, except for recent larger volcanic eruptions that are overestimated by both model configurations. The CARMA background surface area density and aerosol size distribution in the stratosphere and troposphere compare well to observations, with some underestimation of the Aitken-mode size range. MAM4 shows shortcomings in reproducing coarse-mode aerosol distributions in the stratosphere and troposphere. This work outlines additional development needs for CESM2 CARMA to improve the model compared to observations in both the troposphere and stratosphere.","PeriodicalId":12799,"journal":{"name":"Geoscientific Model Development","volume":"59 2","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-11-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135327900","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Doroteaciro Iovino, Pier Giuseppe Fogli, Simona Masina
{"title":"Evaluation of the CMCC global eddying ocean model for the Ocean Model Intercomparison Project (OMIP2)","authors":"Doroteaciro Iovino, Pier Giuseppe Fogli, Simona Masina","doi":"10.5194/gmd-16-6127-2023","DOIUrl":"https://doi.org/10.5194/gmd-16-6127-2023","url":null,"abstract":"Abstract. This paper describes the global eddying ocean–sea ice simulation produced at the Euro-Mediterranean Center on Climate Change (CMCC) obtained following the experimental design of the Ocean Model Intercomparison Project phase 2 (OMIP2). The eddy-rich model (GLOB16) is based on the NEMOv3.6 framework, with a global horizontal resolution of 1/16∘ and 98 vertical levels and was originally designed for an operational short-term ocean forecasting system. Here, it is driven by one multi-decadal cycle of the prescribed JRA55-do atmospheric reanalysis and runoff dataset in order to perform a long-term benchmarking experiment. To assess the accuracy of simulated 3D ocean fields and highlight the relative benefits of resolving mesoscale processes, the GLOB16 performances are evaluated via a selection of key climate metrics against observational datasets and two other NEMO configurations at lower resolutions: an eddy-permitting resolution (ORCA025) and a non-eddying resolution (ORCA1) designed to form the ocean–sea ice component of the fully coupled CMCC climate model. The well-known biases in the low-resolution simulations are significantly improved in the high-resolution model. The evolution and spatial pattern of large-scale features (such as sea surface temperature biases and winter mixed-layer structure) in GLOB16 are generally better reproduced, and the large-scale circulation is remarkably improved compared to the low-resolution oceans. We find that eddying resolution is an advantage in resolving the structure of western boundary currents, the overturning cells, and flow through key passages. GLOB16 might be an appropriate tool for ocean climate modeling efforts, even though the benefit of eddying resolution does not provide unambiguous advances for all ocean variables in all regions.","PeriodicalId":12799,"journal":{"name":"Geoscientific Model Development","volume":"126 2","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-11-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135326079","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Dien Wu, Joshua L. Laughner, Junjie Liu, Paul I. Palmer, John C. Lin, Paul O. Wennberg
{"title":"A simplified non-linear chemistry transport model for analyzing NO<sub>2</sub> column observations: STILT–NO<sub><i>x</i></sub>","authors":"Dien Wu, Joshua L. Laughner, Junjie Liu, Paul I. Palmer, John C. Lin, Paul O. Wennberg","doi":"10.5194/gmd-16-6161-2023","DOIUrl":"https://doi.org/10.5194/gmd-16-6161-2023","url":null,"abstract":"Abstract. Satellites monitoring air pollutants (e.g., nitrogen oxides; NOx = NO + NO2) or greenhouse gases (GHGs) are widely utilized to understand the spatiotemporal variability in and evolution of emission characteristics, chemical transformations, and atmospheric transport over anthropogenic hotspots. Recently, the joint use of space-based long-lived GHGs (e.g., carbon dioxide; CO2) and short-lived pollutants has made it possible to improve our understanding of emission characteristics. Some previous studies, however, lack consideration of the non-linear NOx chemistry or complex atmospheric transport. Considering the increase in satellite data volume and the demand for emission monitoring at higher spatiotemporal scales, it is crucial to construct a local-scale emission optimization system that can handle both long-lived GHGs and short-lived pollutants in a coupled and effective manner. This need motivates us to develop a Lagrangian chemical transport model that accounts for NOx chemistry and fine-scale atmospheric transport (STILT–NOx) and to investigate how physical and chemical processes, anthropogenic emissions, and background may affect the interpretation of tropospheric NO2 columns (tNO2). Interpreting emission signals from tNO2 commonly involves either an efficient statistical model or a sophisticated chemical transport model. To balance computational expenses and chemical complexity, we describe a simplified representation of the NOx chemistry that bypasses an explicit solution of individual chemical reactions while preserving the essential non-linearity that links NOx emissions to its concentrations. This NOx chemical parameterization is then incorporated into an existing Lagrangian modeling framework that is widely applied in the GHG community. We further quantify uncertainties associated with the wind field and chemical parameterization and evaluate modeled columns against retrieved columns from the TROPOspheric Monitoring Instrument (TROPOMI v2.1). Specifically, simulations with alternative model configurations of emissions, meteorology, chemistry, and inter-parcel mixing are carried out over three United States (US) power plants and two urban areas across seasons. Using the U.S. Environmental Protection Agency (EPA)-reported emissions for power plants with non-linear NOx chemistry improves the model–data alignment in tNO2 (a high bias of ≤ 10 % on an annual basis), compared to simulations using either the Emissions Database for Global Atmospheric Research (EDGAR) model or without chemistry (bias approaching 100 %). The largest model–data mismatches are associated with substantial biases in wind directions or conditions of slower atmospheric mixing and photochemistry. More importantly, our model development illustrates (1) how NOx chemistry affects the relationship between NOx and CO2 in terms of the spatial and seasonal variability and (2) how assimilating tNO2 can quantify systematic biases in modeled wind directions and emission","PeriodicalId":12799,"journal":{"name":"Geoscientific Model Development","volume":"62 2","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-11-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135326360","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"The Hydro-ABC model (Version 2.0): a simplified convective-scale model with moist dynamics","authors":"Jiangshan Zhu, Ross Noel Bannister","doi":"10.5194/gmd-16-6067-2023","DOIUrl":"https://doi.org/10.5194/gmd-16-6067-2023","url":null,"abstract":"Abstract. The prediction of convection (in terms of position, timing, and strength) is important to achieve for high-resolution weather forecasting. This problem requires not only good convective-scale models, but also data assimilation systems that give initial conditions which neither improperly hinder nor improperly hasten convection in the ensuing forecasts. Solving this problem is difficult and expensive using operational-scale numerical weather prediction systems, and so a simplified model of convective-scale flow is under development (called the “ABC model”). This paper extends the existing ABC model of dry convective-scale flow to include mixing ratios of vapour and condensate phases of water. The revised model is called “Hydro-ABC”. Hydro-ABC includes transport of the vapour and condensate mixing ratios within a dynamical core, and it transitions between these two phases via a micro-physics scheme. A saturated mixing ratio is derived from model quantities, which helps determine whether evaporation or condensation happens. Latent heat is exchanged with the buoyancy variable (ABC's potential-temperature-like variable) in such a way to conserve total energy, where total energy is the sum of dry energy and latent heat. The model equations are designed to conserve the domain-total mass, water, and energy. An example numerical model integration is performed and analysed, which shows the development of a realistic looking anvil cloud and excitation of inertio-gravity and acoustic modes over a wide range of frequencies. This behaviour means that Hydro-ABC is a sufficiently challenging model which will allow experimentation with innovative data assimilation strategies in future work. An ensemble of Hydro-ABC integrations is performed in order to study the possible forecast error covariance statistics (knowledge of which is necessary for data assimilation). These show patterns that are dependent on the presence of convective activity (at any model's vertical column), thus giving a taste of flow-dependent error statistics. Candidate indicators/harbingers of convection are also evaluated (namely relative humidity, hydrostatic imbalance, horizontal divergence, convective available potential energy, convective inhibition, vertical wind, and the condensate mixing ratio), some of which appear to be reliable diagnostics concerning the presence of convection. These diagnostics will be useful in the selection of the relevant forecast error covariance statistics when data assimilation for Hydro-ABC is developed.","PeriodicalId":12799,"journal":{"name":"Geoscientific Model Development","volume":"59 8","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-10-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135813620","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}