Inversion of Critical Atmospheric 137Cs Emissions Following the Fukushima Accident by Resolving Temporal Formation from Total Deposition Data.

Understanding the transport of 137Cs emitted during the Fukushima accident is challenging because the critical emissions that produced the high-deposition area are not adequately resolved in existing source terms. This paper presents an objective inverse reconstruction of these emissions by fusing atmospheric concentrations with a-priori emissions extracted from total depositions. This extraction, previously considered impossible for complex real-world accidents, is achieved by identifying the critical temporal formation process of depositions in the high-deposition area and estimating the corresponding emissions by using an atmospheric transport model. The reconstructed source term reveals two emission peaks from 10:00-11:00 and 14:00-15:00 on March 15, which agree with the in situ pressure measurements and accident analysis, suggesting that they came from pressure drops in the primary containment vessels of Units 3 and 2, respectively. This finding explains the environmental observations of spherical 137Cs particles. The source term also objectively and independently confirms the widely used reverse estimate. The corresponding 137Cs transport simulations better match the various observations than those produced by other source terms, proving that the two-peak emission creates a high-deposition area. The proposed method outperforms the direct fusion of deposition and atmospheric concentration observations, providing a robust tool for multiobservation fusion.

Wet scavenging of multi-mode 137Cs aerosols following the Fukushima accident: Size-resolved microphysics modeling with observed diameters.

Wet scavenging was critical in the atmospheric transport of 137Cs aerosols following the Fukushima accident. The aerosol size diversity and related microphysical processes produce complex behaviors during wet scavenging. Such behaviors are difficult to investigate using traditional simplified size distributions, resulting in inaccurate modeling. This study establishes an improved size-resolved wet scavenging model that considers the activation process. Using this model, five monodisperse simulations with five representative observed diameters with realistic solubility setting are performed to investigate the spatiotemporal wet scavenging behaviors of 137Cs aerosols. One polydisperse simulation with an empirical size distribution is also validated against the observation. The results reveal that 137Cs aerosols with diameters of 0.6 and 2.0 µm are mainly subject to below-cloud scavenging, which makes a significant contribution to low-deposition areas (<300 kBq/m2). For 137Cs aerosols with diameters of 6.4, 15, and 30 µm, in-cloud scavenging dominates, and the resulting depositions make significant contributions in high-deposition areas. The polydisperse results satisfy the criteria for good performance and better agree with the size, and deposition observations than the five monodisperse simulations, whereas for the concentration, the results show a similar RANK2 with the best mono1 and mono2 cases and reach the satisfactory criteria. These findings reveal the complex behavior and wet scavenging process of multi-mode 137Cs aerosols, improving our understanding and modeling.

Multi-scenario validation of the robust inversion method with biased plume range and values.

Release rate estimation is a vital means of revealing the emission process of radionuclides and assessing the environmental consequences in an emergency. Inverse modeling is widely used in emergency cases, but is vulnerable to plume biases in atmospheric dispersion modeling. One promising solution is a model called "Simultaneously Estimates the Release rate And Corrects both the plume range and Transport pattern" (SERACT). This study investigates the feasibility and behavior of SERACT based on four wind tunnel experiments replicating complex dispersion scenarios with both dense buildings and heterogeneous topography. SERACT's performance is compared with that of Tikhonov-regularized inversion and its predecessor. The results demonstrate that SERACT successfully corrects the modeled plume biases and simultaneously improves the release rate estimations in all four complex local-scale scenarios. The release rates retrieved by SERACT provide better agreement with the true release rates than those given by the other methods for all scenarios, with an average deviation of only 5.83%. After correction, the simulated plume reproduces the concentrations at all sites and achieves the best Pearson correlation coefficient (1.00) and fraction of simulations within a factor of 2 of the measurements (1.00); these values are 7.33 and 2.09 times higher, respectively, than those of simulations using release rates obtained using Tikhonov-regularized inversion.

Local atmospheric transport behaviors of representative radionuclides during the Fukushima accident: A 200-m-resolution cross-scale study from site to 20 km.

Modeling of local atmospheric radionuclide transport is essential for nuclear emergency response. However, very few studies of the Fukushima Dai-ichi nuclear power plant (FDNPP) accident have focused on this topic because of the complex meteorology and the cross-scale transport behaviors from the site to 20 km of the FDNPP. In this study, both the local meteorology and transport behaviors were investigated at high resolution (200 m) using ensembles of different meteorology and models. Four wind fields calculated from onsite observations and three regional-scale meteorological fields (i.e., the 1-km ECMWF, 3-km and 1-km NHM-LETKF), and two transport models: the RIMPUFF Lagrangian puff model and the SPRAY particle model were considered and combined with each other. These eight simulations and their ensemble mean were analyzed based on onsite observations of wind and gamma dose rates, and local-scale observations of 137Cs concentration. Results revealed that at the site, the Onsite wind field which captured the frequently changing wind, best reproduced the onsite gamma dose rates with the grid resolution of 200 m. At the local scale (up to 20 km), the observations present a smoother temporal change. The wind fields assimilated with Japanese domestic observations presented advantageous performance, and the 1-km NHM-LETKF achieved the best score of the factor of 5 metric of 0.49 for the simulated 137Cs concentration. The SPRAY coupled with the three-dimensional (3D) convolution method and RIMPUFF showed better performance in simulating the onsite gamma dose rate and the local-scale concentration, respectively. The ensemble mean achieved robust metrics, better simulated the baseline of onsite gamma dose rates, and reproduced more peaks of local-scale concentration at the expense of peak value deviation.

Model behavior regarding in- and below-cloud 137Cs wet scavenging following the Fukushima accident using 1-km-resolution meteorological field data.

Wet deposition remains an important source of uncertainty in modeling of the atmospheric transport of 137Cs following the Fukushima Daiichi Nuclear Power Plant accident. Its behavior is often difficult to investigate owing to the limited resolution of meteorological field data and inconsistent model implementations. This study investigated the detailed behavior of 25 combinations of in- and below-cloud wet scavenging models using high-resolution (1 km × 1 km) meteorological input. These combinations were all implemented in the Weather Research and Forecasting-Chemistry model, thereby enabling consistent evaluation. The 1-km-resolution simulations were compared with simulations obtained previously using 3-km-resolution meteorological field data. Results revealed that rainfall of <1 mm/h is critical for simulation accuracy. The 1-km results revealed better representation of rainfall than that revealed by the 3-km results, but with spatiotemporal variability in accuracy. Owing to their sensitivity to rainfall, single-parameter wet deposition models showed improvements in performance in the 1-km simulations relative to that in the 3-km simulations. The multiparameter models showed more robust performance in terms of both simulations, and the Roselle-Mircea model presented the best performance among the 25 models considered. Wind transport showed substantial influence on the removal of atmospheric 137Cs, and it was nonnegligible even during periods in which wet deposition was dominant. The 1-km-resolution simulations effectively reproduced local-scale 137Cs concentrations but with deviations in timing, mainly because of biased wind direction. These findings indicate the necessity for a refined wind and dispersion model for local-scale simulation of 137Cs concentration.

Objective inversion of the continuous atmospheric 137Cs release following the Fukushima accident.

Eleven years after the Fukushima accident, independent objective estimates of the atmospheric 137Cs release still suffer from discontinuities such as negative release terms, oscillations, and temporal gaps, leading to noticeable differences from the subjective estimate. This paper describes an objective method that handles these artifacts and promotes the continuity of releases at fine resolutions. The proposed method uses the joint estimation model to reduce the oscillations induced by the model-observation discrepancies, and employs total variation regularization to recover the missing releases caused by insufficient observations. Adaptive parameterization is used to correct negative values. The application of this method to the Fukushima accident produces a source term that accurately approximates continuous releases at a fine temporal resolution of 1 h, providing a better match with the recognized subjective source term than nine published estimates, with a Pearson's correlation coefficient of 0.923 and an index of agreement of 0.872. This source term agrees with the timing of on-site gamma dose rate peaks, significantly improving the air concentration and deposition simulations, with FAC10 values of 0.564 and 0.990, respectively. The estimation error varies smoothly in a limited range with different regularization parameters, enabling automatic parameterization and demonstrating the potential for operational inversions.

Inversion of 137Cs emissions following the fukushima accident with adaptive release recovery for temporal absences of observations.

Temporal absences in observation records lead to release losses during the source term inversions of atmospheric radionuclide emissions. Consequently, objectively-estimated source terms for the Fukushima accident contain fewer release details and present large discrepancies when compared with the expert-judged one. This paper describes an objective method that can adaptively recover the missing releases caused by the temporal absences of observations. The proposed method assumes that the accident releases of radionuclides are piecewise-constant and comprise both peaks and constant releases. The missing releases are adaptively recovered as either peaks or constant releases by minimizing the total variation of the estimated source term. The proposed method is applied to the Fukushima accident and evaluated against regional airborne and deposited 137Cs observations. The results demonstrate that this method effectively recovers the missing releases, producing a source term that matches the timing of both on-site gamma dose rate peaks and accident events. The retrieved source term improves the simulation of air concentrations and reproduces most of the deposition patterns. This is the first time that an objective method has independently reproduced the details in the expert-judged one for the Fukushima accident.

Coupled modeling of in- and below-cloud wet deposition for atmospheric 137Cs transport following the Fukushima Daiichi accident using WRF-Chem: A self-consistent evaluation of 25 scheme combinations.

Wet deposition, including both in- and below-cloud scavenging, is critical for the atmospheric transport modeling of 137Cs following the Fukushima Daiichi Nuclear power plant (FDNPP) accident. Although intensively investigated, wet deposition simulation is still subject to uncertainties of meteorological inputs and wet scavenging modeling, leading to biased 137Cs transport prediction. To reduce the dual uncertainties, in- and below-cloud wet scavenging schemes of 137Cs were simultaneously integrated into Weather Research and Forecasting-Chemistry (WRF-Chem), yielding online coupled modeling of meteorology and the two wet scavenging processes. The integration was performed using 25 combinations of different in- and below-cloud schemes, covering most schemes in the literature. Two microphysics schemes were also tested to better reproduce the precipitation. The 25 models and the ensemble mean of 9 representative models were systematically compared with the below-cloud-only WRF-Chem model, using the cumulative deposition and atmospheric concentrations of 137Cs measurements. The results reveal that, with the Morrison's double moment cloud microphysics scheme, the developed models could better reproduce the rainfall and substantially improve the cumulative deposition simulation. The in-cloud scheme is influential to the model behaviors and those schemes considering cloud parameters also improve the atmospheric concentration simulations, whereas the others solely dependent on the rain intensity are sensitive to meteorology. The ensemble mean achieves satisfactory performance except one plume event, but still outperforms most models.

Multi-scenario validation of CALMET-RIMPUFF for local-scale atmospheric dispersion modeling around a nuclear powerplant site with complex topography.

The CALMET-RIMPUFF is composed of the California Meteorological Model and the Risø Mesoscale PUFF model, which provides refined atmospheric dispersion modeling for nuclear emergency response. Because the performance of an atmospheric dispersion model can be case-sensitive, a multi-scenario validation is important to understand a model's behavior and limits. In this study, a multi-scenario validation of CALMET-RIMPUFF was performed based on six wind tunnel experiments simulating a real China's nuclear powerplant site with complex topographies and dense buildings. The CALMET-RIMPUFF simulations were compared with the measurements of the vertical wind profiles, 2D ground wind and concentration fields, both qualitatively and quantitatively. The results demonstrate that the CALMET-RIMPUFF can simulate the ground-level wind with acceptable accuracies. For vertical wind profiles, the accuracies show high dependencies on the local topography and building layout. The simulated ground concentrations generally agree well with the measurements, though the plume axis showed slight discrepancies from the measurements in three cases. Because the CALMET-RIMPUFF lacks a building effect module, it shows noticeable discrepancies in the building area. However, such discrepancies do not propagate to the downwind mountainous and sea areas, which the accuracies are quite satisfactory. Thus, the CALMET-RIMPUFF is capable for local-scale modeling at this site.

Automated release rate inversion and plume bias correction for atmospheric radionuclide leaks: A robust and general remediation to imperfect radionuclide transport modeling.

The release rate is vital in assessing the international environmental risk of atmospheric radionuclide leaks, and it usually can only be obtained through inversion. However, such inversion is vulnerable to the inevitable plume biases in radionuclide transport modeling, leading to inaccurate estimates and risk assessment. This paper describes an automated method that estimates the release rate while comprehensively correcting plume biases, including both the plume range and transport pattern. The spatial correlation of predictions is used to simplify the difficult task of direct plume adjustment to that of tuning the predictions inside a correlation-adjusted plume. An ensemble-based algorithm is proposed to automatically calculate the spatial correlation. The proposed method is validated using two radionuclide transport models with mild and severe plume biases and data from two wind tunnel experiments, and its performance is compared with that of the standard approach and a recent state-of-the-art method. The results demonstrate that our method corrects the plume biases with high accuracy (Pearson's Correlation Coefficient = 1.0000, Normalized Mean Square Error ≤ 1.03 × 10-3) and reduces the estimation error by nearly two orders of magnitude at best. The proposed approach achieves near-optimal performance with fully automated parameterization, keeping the lowest error levels in our validation cases for various measurement sets.

Dewaterability of sewage sludge conditioned with a graft cationic starch-based flocculant: Role of structural characteristics of flocculant.

Coagulation/flocculation is one of the most extensive and cost-effective pretreatments to improve the dewaterability of sludge in water treatment plants. In this study, three series of graft cationic starch (St)-based flocculants with distinct structural characteristics including charge density (CD), graft-chain length (L), and graft-chain distribution (N) were synthesized by graft copolymerization of [(2-methacryloyloxy-ethyl) trimethyl ammonium chloride] and acrylamide onto St backbone. The structural effects of these St-based flocculants on the sludge dewaterability have been quantitatively analyzed by using a second-order polynomial model according to phenomenological theory. The predicted dewatering performance and optimal dose were fully consistent with the experimental results. On the basis of this established model, the dewatering mechanisms were discussed in detail by combination of the analysis of the changes in filter cake moisture content, specific resistance of filtration, bound water content, compression coefficient, extracellular polymeric substances fractions and components, spatial distributions of proteins and polysaccharides, microstructures of sludge cake, and flocs properties in the dewatering process. This graft St-based flocculant, with the structural characteristics of high CD, long L, and low N, exhibited superior sludge dewaterability because of the enhanced charge neutralization and bridging flocculation effects. Among these three structural factors, CD played a more important role in improvement of sludge dewaterability than L and N due to the dominant effect of charge neutralization. This study provided a better understanding of structure-activity relationship of these grafting modified flocculants, which was of significant guidance for the exploit and design of novel and efficient flocculants for improvement of sludge dewaterability.

Comparison of the Outcomes of Three Different Nutritional Supports in Patients with Oral and Maxillofacial Malignant Tumors following Surgery.

Objective. This study aimed to compare the physical and mental states and the clinical effects of parenteral nutrition combined with enteral nutrition (PN+EN), total enteral nutrition (TEN), and total parenteral nutrition (TPN) after surgery in patients with maxillofacial malignant tumors. Methods. A total of 112 patients were divided into three groups, with 58, 33, and 21 patients in the PN+EN, TPN, and TEN groups, respectively. The psychological survey contained the Faces Pain Scale-Revised (FRS-R), visual analog scale (VAS), numerical rating scale (NRS), Hamilton anxiety rating scale (HAMA), and short-form 36 health survey questionnaire (SF-36). Spirit symptoms, length of hospital stay, nutritional assessments, and related biochemical indices were recorded and compared. Results. The traditional Chinese medicine (TCM) symptoms of anxiety and dysphoria were least frequently identified in the TPN group. The levels of lymphocytes, hemoglobin (HB), albumin (ALB), and prealbumin (PA) were significantly higher in the PE+EN group, whereas white blood cell count, neutrophil count, HB, PA, and ALB were significantly lower in the TPN group. Better psychological scores were observed in the TPN group. The PE+EN group had a shorter length of stay and higher SGA categories. Potassium, sodium, and chlorine levels were significantly lower in the TEN group (all P < 0.05). Conclusions. As an auxiliary method, TCM symptoms can help to identify spirit disequilibrium earlier and are associated with blood indices. Without the consideration of cost and long length of hospital stay, patients in the TPN group had the best mental status, with PN+EN therapy being an alternative.