Improved Microwave Ocean Emissivity and Reflectivity Models Derived from Two-Scale Roughness Theory

The Geometrical Optics(GO)approach and the FAST Emissivity Model(FASTEM)are widely used to estimate the surface radiative components in atmospheric radiative transfer simulations,but their applications are limited in specific conditions.In this study,a two-scale reflectivity model(TSRM)and a two-scale emissivity model(TSEM)are developed from the two-scale roughness theory.Unlike GO which only computes six non-zero elements in the reflectivity matrix,The TSRM includes 16 elements of Stokes reflectivity matrix which are important for improving radiative transfer simulation accuracy in a scattering atmosphere.It covers the frequency range from L-to W-bands.The dependences of all TSRM elements on zenith angle,wind speed,and frequency are derived and analyzed in details.For a set of downwelling radiances in microwave frequencies,the reflected upwelling brightness temperature(BTs)are calculated from both TSRM and GO and compared for analyzing their discrepancies.The TSRM not only includes the effects of GO but also accounts for the small-scale Bragg scattering effect in an order of several degrees in Kelvins in brightness temperature.Also,the third and fourth components of the Stokes vector can only be produced from the TSRM.For the emitted radiation,BT differences in vertical polarization between a TSEM and FASTEM are generally less than 5 K when the satellite zenith angle is less than 40°,whereas those for the horizontal component can be quite significant,greater than 20 K.

An Improved Atmospheric Vector Radiative Transfer Model Incorporating Rough Ocean Boundaries

The radiative transfer model (RT3), a vector radiative transfer (VRT) scheme in a plane-parallel atmosphere, was bounded by a rough ocean surface in this study. The boundary problem was solved using a Fourier series decomposition of the radiation field as a function of the azimuth. For the case of a rough ocean surface, the decomposition was obtained by developing both the Fresnel reflection matrix and the probability distribution of the water facet orientation as Fourier series. The effect of shadowing by ocean surface waves was also considered in the boundary condition. The VRT model can compute the intensity and degree of polarization of the light at the top of the atmosphere (TOA), the ocean surface, and any level of the atmosphere in the ocean-atmosphere system. The results obtained by our model are in good agreement with those computed by Ahmad’s model. The simulated results showed that the shadow effects of wave facets on the intensity and the degree of polarization are negligible except at the ocean surface near the grazing angle, possibly because we did not consider the effect of white caps.

Geometric Multigrid Method for Isogeometric Analysis

The isogeometric analysis method(IGA)is a new type of numerical method solving partial differential equations.Compared with the traditional finite element method,IGA based on geometric spline can keep the model consistency between geometry and analysis,and provide higher precision with less freedom.However,huge stiffness matrix fromthe subdivision progress still leads to the solution efficiency problems.This paper presents amultigrid method based on geometric multigrid(GMG)to solve the matrix system of IGA.This method extracts the required computational data for multigrid method fromthe IGA process,which also can be used to improve the traditional algebraic multigrid method(AGM).Based on this,a full multigrid method(FMG)based on GMG is proposed.In order to verify the validity and reliability of these methods,this paper did some test on Poisson’s equation and Reynolds’equation and compared the methods on different subdivision methods,different grid degrees of freedom,different cyclic structure degrees,and studied the convergence rate under different subdivision strategies.The results show that the proposed method is superior to the conventional algebraic multigrid method,and for the standard relaxed V-cycle iteration,the method still has a convergence speed independent of the grid size at the same degrees.