NUMERICAL MODELING OF RADIATIVE TRANSFER FOR MICROWAVE REMOTE SENSING

An overall vector radiative transfer theory was developed for numerical modeling, in both active and passive microwave remote sensing. The Theory and approaches are briefly summerized.To quantitatively understand scattering and thermal emission from targets in active and passive remote sensing, we have developed an overall vector radiative transfer theory for a set of theoretical models of discrete scatterer and continuous random media for the earth terrain (wet soil, vegetation, snow, sea-ice, etc.) and atmosphere, and numerical approaches for simulation, data analysis, and parameter sensitivity test. Our numerical results favorably agreed with experimental data in microwave re mote sensing of various earth surfaces. Main approaches are briefly summerized herewith.

A Multi-Domain Compression Radiative Transfer Model for the Fengyun-4 Geosynchronous Interferometric Infrared Sounder (GIIRS)

Forward radiative transfer(RT)models are essential for atmospheric applications such as remote sensing and weather and climate models,where computational efficiency becomes equally as important as accuracy for high-resolution hyperspectral measurements that need rigorous RT simulations for thousands of channels.This study introduces a fast and accurate RT model for the hyperspectral infrared(HIR)sounder based on principal component analysis(PCA)or machine learning(i.e.,neural network,NN).The Geosynchronous Interferometric Infrared Sounder(GIIRS),the first HIR sounder onboard the geostationary Fengyun-4 satellites,is considered to be a candidate example for model development and validation.Our method uses either PCA or NN(PCA/NN)twice for the atmospheric transmittance and radiance,respectively,to reduce the number of independent but similar simulations to accelerate RT simulations;thereby,it is referred to as a multi-domain compression model.The first PCA/NN gives monochromatic gas transmittance in both spectral and atmospheric pressure domains for each gas independently.The second PCA/NN is performed in the traditional spectral radiance domain.Meanwhile,a new method is introduced to choose representative variables for the PCA/NN scheme developments.The model is three orders of magnitude faster than the standard line-by-line-based simulations with averaged brightness temperature difference(BTD)less than 0.1 K,and the compressions based on PCA or NN methods result in comparable efficiency and accuracy.Our fast model not only avoids an excessively complicated transmittance scheme by using PCA/NN but is also highly flexible for hyperspectral instruments with similar spectral ranges simply by updating the corresponding spectral response functions.

Consistency of Tropospheric Water Vapor between Reanalyses and Himawari-8/AHI Measurements over East Asia

High spatiotemporal resolution radiances from the advanced imagers onboard the new generation of geostationary weather satellites provide a unique opportunity to evaluate the abilities of various reanalysis datasets to depict multilayer tropospheric water vapor(WV),thereby enhancing our understanding of the deficiencies of WV in reanalysis datasets.Based on daily measurements from the Advanced Himawari Imager(AHI)onboard the Himawari-8 satellite in 2016,the bias features of multilayer WV from six reanalysis datasets over East Asia are thoroughly evaluated.The assessments show that wet biases exist in the upper troposphere in all six reanalysis datasets;in particular,these biases are much larger in summer.Overall,we find better depictions of WV in the middle troposphere than in the upper troposphere.The accuracy of WV in the ERA5 dataset is the highest,in terms of the bias magnitude,dispersion,and pattern similarity.The characteristics of the WV bias over the Tibetan Plateau are significantly different from those over other parts of East Asia.In addition,the reanalysis datasets all capture the shift of the subtropical high very well,with ERA5 performing better overall.

A Simplified Scheme of the Generalized Layered Radiative Transfer Model

In this paper, firstly, a simplified version （SGRTM） of the generalized layered radiative transfer model （GRTM） within the canopy, developed by us, is presented. It reduces the information requirement of inputted sky diffuse radiation, as well as of canopy morphology, and in turn saves computer resources. Results from the SGRTM agree perfectly with those of the GRTM. Secondly, by applying the linear superposition principle of the optics and by using the basic solutions of the GRTM for radiative transfer within the canopy under the condition of assumed zero soil reflectance, two sets of explicit analytical solutions of radiative transfer within the canopy with any soil reflectance magnitude are derived： one for incident diffuse, and the other for direct beam radiation. The explicit analytical solutions need two sets of basic solutions of canopy reflectance and transmittance under zero soil reflectance, run by the model for both diffuse and direct beam radiation. One set of basic solutions is the canopy reflectance αf （written as α1 for direct beam radiation） and transmittance βf （written as β1 for direction beam radiation） with zero soil reflectance for the downward radiation from above the canopy （i.e. sky）, and the other set is the canopy reflectance （αb） and transmittance βb for the upward radiation from below the canopy （i.e., ground）. Under the condition of the same plant architecture in the vertical layers, and the same leaf adaxial and abaxial optical properties in the canopies for the uniform diffuse radiation, the explicit solutions need only one set of basic solutions, because under this condition the two basic solutions are equal, i.e., αf = αb and βf = βb. Using the explicit analytical solutions, the fractions of any kind of incident solar radiation reflected from （defined as surface albedo, or canopy reflectance）, transmitted through （defined as canopy transmittance）, and absorbed by （defined as canopy absorptance） the canopy and other properties pertinent to the radiative transfer within the canopy can be estimated easily on the ground surface below the canopy （soil or snow surface） with any reflectance magnitudes. The simplified transfer model is proven to have a similar accuracy compared to the detailed model, as well as very efficient computing.

A Comparison of Two Canopy Radiative Models in Land Surface Processes

This paper compares the predictions by two radiative transfer models-the two-stream approximation model and the generalized layered model （developed by the authors） in land surface processes -for different canopies under direct or diffuse radiation conditions. The comparison indicates that there are significant differences between the two models, especially in the near infrared （NIR） band. Results of canopy reflectance from the two-stream model are larger than those from the generalized model. However, results of canopy absorptance from the two-stream model are larger in some cases and smaller in others compared to those from the generalized model, depending on the cases involved. In the visible （VIS） band, canopy reflectance is smaller and canopy absorptance larger from the two-stream model compared to the generalized model when the Leaf Area Index （LAI） is low and soil reflectance is high. In cases of canopies with vertical leaf angles, the differences of reflectance and absorptance in the VIS and NIR bands between the two models are especially large. Two commonly occurring cases, with which the two-stream model cannot deal accurately, are also investigated. One is for a canopy with different adaxial and abaxial leaf optical properties; and the other is for incident sky diffuse radiation with a non-uniform distribution. Comparison of the generalized model within the same canopy for both uniform and non-uniform incident diffuse radiation inputs shows smaller differences in general. However, there is a measurable difference between these radiation inputs for a canopy with high leaf angle. This indicates that the application of the two-stream model to a canopy with different adaxial and abaxial leaf optical properties will introduce non-negligible errors.

Grid-cell Aerosol Direct Shortwave Radiative Forcing Calculated Using the SBDART Model with MODIS and AERONET Observations：An Application in Winter and Summer in Eastern China

Taking winter and summer in eastern China as an example application, a grid-cell method of aerosol direct radiative forcing（ADRF） calculation is examined using the Santa Barbara DISORT Atmospheric Radiative Transfer（SBDART） model with inputs from MODIS and AERONET observations and reanalysis data. Results show that there are significant seasonal and regional differences in climatological mean aerosol optical parameters and ADRF. Higher aerosol optical depth（AOD）occurs in summer and two prominent high aerosol loading centers are observed. Higher single scattering albedo（SSA） in summer is likely associated with the weak absorbing secondary aerosols. SSA is higher in North China during summer but higher in South China during winter. Aerosols induce negative forcing at the top of the atmosphere（TOA） and surface during both winter and summer, which may be responsible for the decrease in temperature and the increase in relative humidity.Values of ADRF at the surface are four times stronger than those at the TOA. Both AOD and ADRF present strong interannual variations; however, their amplitudes are larger in summer. Moreover, patterns and trends of ADRF do not always correspond well to those of AOD. Differences in the spatial distributions of ADRF between strong and weak monsoon years are captured effectively. Generally, the present results justify that to calculate grid-cell ADRF at a large scale using the SBDART model with observational aerosol optical properties and reanalysis data is an effective approach.

X-ray spectra of high temperature tungsten plasma calculated with collisional radiative model

Tungsten is regarded as an important candidate of plasma facing material in international thermonuclear experimental reactor （ITER）, so the determination and modeling of spectra of tungsten plasma, especially the spectra at high temperature were intensely focused on recently. In this work, using the atomic structure code of Cowan, a collisional radiative model （CRM） based on the spin-orbit-split-arrays is developed. Based on this model, the charge state distribution of tungsten ions is determined and the soft X-ray spectra from high charged ions of tungsten at different temperatures are calculated. The results show that both the average ionization charge and line positions are well agreed with others calculations and measurements with discrepancies of less than 0.63% and 1.26%, respectively. The spectra at higher temperatures are also reported and the relationship between ion abundance and temperature is predicted in this work.

Microwave Simulations of Precipitation Distribution with Two Radiative Transfer Models

Two microwave radiative transfer models of precipitating cloud are used to simulate the microwave upwelling radiances emerging from precipitating clouds. Comparison of the simulation results shows that significant difference of microwave upwelling radiances exists between these two radiative transfer models. Analysis of these differences in different cloud and precipitation conditions shows that it is complicated but has certain trend for different microwave frequencies. The results may be useful to quantitative rainfall rate retrieval of real precipitating clouds.

Research on the Application of the Radiative Transfer Model Based on Deep Neural Network in One-dimensional Variational Algorithm

As a typical physical retrieval algorithm for retrieving atmospheric parameters,one-dimensional variational(1 DVAR)algorithm is widely used in various climate and meteorological communities and enjoys an important position in the field of microwave remote sensing.Among algorithm parameters affecting the performance of the 1 DVAR algorithm,the accuracy of the microwave radiative transfer model for calculating the simulated brightness temperature is the fundamental constraint on the retrieval accuracies of the 1 DVAR algorithm for retrieving atmospheric parameters.In this study,a deep neural network(DNN)is used to describe the nonlinear relationship between atmospheric parameters and satellite-based microwave radiometer observations,and a DNN-based radiative transfer model is developed and applied to the 1 DVAR algorithm to carry out retrieval experiments of the atmospheric temperature and humidity profiles.The retrieval results of the temperature and humidity profiles from the Microwave Humidity and Temperature Sounder(MWHTS)onboard the Feng-Yun-3(FY-3)satellite show that the DNN-based radiative transfer model can obtain higher accuracy for simulating MWHTS observations than that of the operational radiative transfer model RTTOV,and also enables the 1 DVAR algorithm to obtain higher retrieval accuracies of the temperature and humidity profiles.In this study,the DNN-based radiative transfer model applied to the 1 DVAR algorithm can fundamentally improve the retrieval accuracies of atmospheric parameters,which may provide important reference for various applied studies in atmospheric sciences.

Non-invasive optical characterization and estimation of Zn porosity in gas tungsten arc welding of Fe–Al joints using CR model and OES measurements

In this study,we employed a non-invasive approach based on the collisional radiative(CR)model and optical emission spectroscopy(OES)measurements for the characterization of gas tungsten arc welding(GTAW)discharge and quantification of Zn-induced porosity during the GTAW process of Fe–Al joints.The OES measurements were recorded as a function of weld current,welding speed,and input waveform.The OES measurements revealed significant line emissions from Zn-I in 460–640 nm and Ar-I in 680–800 nm wavelength ranges in all experimental settings.The OES coupled CR model approach for Zn-I line emission enabled the simultaneous determination of both essential discharge parameters i.e.electron temperature and electron density.Further,these predictions were used to estimate the Zn-induced porosity using OES-actinometry on Zn-I emission lines using Ar as actinometer gas.The OES-actinometry results were in good agreement with porosity data derived from an independent approach,i.e.x-ray radiography images.The current study shows that OES-based techniques can provide an efficient route for real-time monitoring of weld quality and estimate porosity during the GTAW process of dissimilar metal joints.

Modeling Study of the Global Distribution of Radiative Forcing by Dust Aerosol

To quantitatively understand the dust aerosol effects on climate change, we calculated the global dis-tribution of direct radiative forcing due to dust aerosol under clear and cloudy skies in both winter and summer, by using an improved radiative transfer model and the global distribution of dust mass concentra-tion given by GADS （Global Aerosol Data Set）. The results show that the global means of the solar forcing due to dust aerosol at the tropopause for winter and summer are -0.48 and -0.50 W m-2, respectively; the corresponding values for the longwave forcing due to dust are 0.11 and 0.09 W m-2, respectively. At the surface, the global means of the solar forcing clue to dust are -1.36 W m-2 for winter and -1.56 W m-2 for summer, whereas the corresponding values for the longwave forcing are 0.27 and 0.23 W m-2, respectively. This work points out that the absolute values of the solar forcing due to dust aerosol at both the tropopause and surface increase linearly with the cosine of solar zenith angle and surface albedo. The solar zenith angle influences both the strength and distribution of the solar forcing greatly. Clouds exert great effects on the direct radiative forcing of dust, depending on many factors including cloud cover, cloud height, cloud water path, surface albedo, solar zenith angle, etc. The effects of low clouds and middle clouds are larger than those of high clouds. The existence of clouds reduces the longwave radiative forcing at the tropopause, where the influences of low clouds are the most obvious. Therefore, the impacts of clouds should not be ignored when estimating the direct radiative forcing due to dust aerosol.

An improved QuikSCAT wind retrieval algorithm and eye locating for typhoon

This paper proposes a rain considered geophysical model function （GMF）, to be noted as GMF plus Rain. GMF plus Rain is based on the basic raidative transfer model with attenuation and scattering effects of rain on radar signal considered. Combined with the NSCAT2 GMF and the rain correction model, the GMF plus Rain model is used to retrieve the ocean wind vectors from the collocated QuikSCAT and SSM/I rain rate data for typhoon Melor. The resulting wind speed estimates of typhoon Melor show improved agreement with the wind fields derived from the best track analysis of Japan Meteorological Agency （JMA）. The results imply that compared with the GMF model, the GMF plus Rain model can improve the precision of wind retrieval under the rain condition. Then, a new general algorithm of locating the eye of typhoon through the normalized radar cross section （NRCS） is proposed. The implementation of this algorithm in the ten QuikSCAT observations of typhoon Melor suggests that this algorithm is effective.

CFD study of non-premixed swirling burners: Effect of turbulence models

This research investigates a numerical simulation of swirling turbulent non-premixed combustion.The effects on the combustion characteristics are examined with three turbulence models:namely as the Reynolds stress model,spectral turbulence analysis and Re-Normalization Group.In addition,the P-1 and discrete ordinate(DO)models are used to simulate the radiative heat transfer in this model.The governing equations associated with the required boundary conditions are solved using the numerical model.The accuracy of this model is validated with the published experimental data and the comparison elucidates that there is a reasonable agreement between the obtained values from this model and the corresponding experimental quantities.Among different models proposed in this research,the Reynolds stress model with the Probability Density Function(PDF)approach is more accurate(nearly up to 50%)than other turbulent models for a swirling flow field.Regarding the effect of radiative heat transfer model,it is observed that the discrete ordinate model is more precise than the P-1 model in anticipating the experimental behavior.This model is able to simulate the subcritical nature of the isothermal flow as well as the size and shape of the internal recirculation induced by the swirl due to combustion.

Real-time 3-D space numerical shake prediction for earthquake early warning

In earthquake early warning systems, real-time shake prediction through wave propagation simulation is a promising approach. Compared with traditional methods, it does not suffer from the inaccurate estimation of source parameters. For computation efficiency, wave direction is assumed to propagate on the 2-D surface of the earth in these methods. In fact, since the seismic wave propagates in the 3-D sphere of the earth, the 2-D space modeling of wave direction results in inaccurate wave estimation. In this paper, we propose a 3-D space numerical shake pre- diction method, which simulates the wave propagation in 3-D space using radiative transfer theory, and incorporate data assimilation technique to estimate the distribution of wave energy. 2011 Tohoku earthquake is studied as an example to show the validity of the proposed model. 2-D space model and 3-D space model are compared in this article, and the prediction results show that numerical shake prediction based on 3-D space model can estimate the real-time ground motion precisely, and overprediction is alleviated when using 3-D space model.

Efficient Application of the Radiance Enhancement Method for Detection of the Forest Fires Due to Combustion-Originated Reflectance

The existing methods for detection of the cloud scenes are applied at relatively small spectral range within shortwave upwelling radiative wavelength flux. We have reported a new method for detection of the cloud scenes based on the Radiance Enhancement (RE). This method can be used to cover a significantly wider spectral range from 1100 nm to 1700 nm by using datasets from the space-orbiting micro-spectrometer Argus 1000. Due to high sunlight reflection of the smoke originated from the forest or field fires the proposed RE method can also be implemented for detection of combustion aerosols. This approach can be a promising technique for efficient detection and continuous monitor of the seasonal forest and field fires. To the best of our knowledge this is the first report showing how a cloud method can be generalized for efficient detection of the forest fires due to combustion-originated reflectance.

Modeling three-dimensional forest structures to drive canopy radiative transfer simulations of bidirectional reflectance factor

Three-dimensional(3-D)Monte Carlo-based radiative transfer(MCRT)models are usually used for benchmarking in intercomparisons of the canopy radiative transfer(RT)simulations.However,the 3-D MCRT models are rarely applied to develop remote sensing algorithms to estimate essential climate variables of forests,due mainly to the difficulties in obtaining realistic stand structures for different forest biomes over regional to global scales.Fortunately,some of important tree structure parameters such as canopy height and tree density distribution have been available globally.This enables to run the intermediate complexities of the 3-D MCRT models.We consequently developed a statistical approach to generate forest structures with intermediate complexities depending on the inputs of canopy height and tree density.It aims at facilitating applications of the 3-D MCRT models to develop remote sensing retrieval algorithms.The proposed approach was evaluated using field measurements of two boreal forest stands at Estonia and USA,respectively.Results demonstrated that the simulations of bidirectional reflectance factor(BRF)based on the measured forest structures agreed well with the BRF based on the generated structures from the proposed approach with the root mean square error(RMSE)and relative RMSE(rRMSE)ranging from 0.002 to 0.006 and from 0.7%to 19.8%,respectively.Comparison of the computed BRF with corresponding MODIS reflectance data yielded RMSE and rRMSE lower than 0.03 and 20%,respectively.Although the results from the current study are limited in two boreal forest stands,our approach has the potential to generate stand structures for different forest biomes.

Effect of the Instrument Slit Function on Upwelling Radiance from a Wavelength Dependent Surface Reflectance

The Radiance Enhancement (RE) method was introduced for efficient detection of clouds from the space. Recently, we have also reported that due to high reflectance of combustion-originated smokes, this approach can also be generalized for detection of the forest fires by retrieving and analyzing datasets collected from a space orbiting micro-spectrometer operating in the near infrared spectral range. In our previous publication, we have performed a comparison of observed and synthetic radiance spectra by developing a method for computation of surface reflectance consisting of different canopies by weighted sum based on their areal coverage. However, this approach should be justified by a method based on corresponding proportions of the upwelling radiance. The results of computations we performed in this study reveal a good match between areal coverage of canopies and the corresponding proportions of the upwelling radiance due to effect of the instrument slit function.

Real-time numerical shake prediction and updating for earthquake early warning

Ground motion prediction is important for earthquake early warning systems, because the region＇s peak ground motion indicates the potential disaster. In order to predict the peak ground motion quickly and pre- cisely with limited station wave records, we propose a real- time numerical shake prediction and updating method. Our method first predicts the ground motion based on the ground motion prediction equation after P waves detection of several stations, denoted as the initial prediction. In order to correct the prediction error of the initial prediction, an updating scheme based on real-time simulation of wave propagation is designed. Data assimilation technique is incorporated to predict the distribution of seismic wave energy precisely. Radiative transfer theory and Monte Carlo simulation are used for modeling wave propagation in 2-D space, and the peak ground motion is calculated as quickly as possible. Our method has potential to predict shakemap, making the potential disaster be predicted before the real disaster happens. 2008 Ms8.0 Wenchuan earthquake is studied as an example to show the validity of the proposed method.

Aerosol optical properties and its direct radiative forcing over Tibetan Plateau from 2006 to 2017

Studies on optical properties of aerosols can reduce the uncertainty for modelling direct radiative forcing(DRF)and improve the accuracy for discussing aerosols effects on the Tibetan Plateau(TP)climate.This study investigated the spatiotemporal variation of aerosol optical and microphysical properties over TP based on OMI and MERRA2,and assessed the influence of aerosol optical properties on DRF at NamCo station(30°46.44′N,90°59.31′E,4730 m)in the central TP from 2006 to 2017 based on a long measurement of AERONET and the modelling of SBDART model.The results show that aerosol optical depth(AOD)exhibits obvious seasonal variation over TP,with higher AOD500nm(>0.75)during spring and summer,and lower value(<0.25)in autumn and winter.The aerosol concentrations show a fluctuated rising from 1980 to 2000,significant increasing from 2000 to 2010 and slight declining trend after 2013.Based on sensitivity experiments,it is found that AOD and single scattering albedo(SSA)have more important impact on the DRF compared withαvalues and ASY.When AOD440nm increases by 60%,DRF at the TOA and ATM is increased by 57.2%and 60.2%,respectively.When SSA440nm increases by 20%,DRF at the TOA and ATM decreases by 121%and 96.7%,respectively.

Soft X-Ray Laser Gain from Neon like Gallium and Germanium

Gain coefficients are calculated for neon-like gallium and germanium ions. Shorter wavelengths are calculated and predicted to be emitted. The gain coefficients are calculated among 457 energy levels of the neon-like ions. Collisional excitations were calculated through the distorted wave approximations through five electron temperatures Te = 300, 500, 700, 1000 and 1500 eV.