Studies of Variational Assimilation for the Inversion of the Coupled Air-sea Model

For the prediction of ENSO, the accuracy of the model including the parameters, initial value and others of the model is important, which can be retrieved by the variational data assimilation methods developed in recent years. However, when the nonlinearity of the model is quite strong, the effect of the improvement made by the 4-D variational data assimilation may be poor due to the bad approximation of the tangent linear model to the original model. So in the paper the ideas in the optimal control is introduced to improve the effect of 4-DVAR in the inversion of the parameters of a nonlinear dynamic ENSO model. The results indicate that when the terminal controlling term is added to the cost functional of 4DVAR, which originated from the optimal control, the effect of the inversion may be largely improved comparing to the traditional 4DVAR, as can be especially obvious from the phase orbit of the model variables. The results in the paper also suggest that the method of 4DVAR in combination with optimal control cannot only reduce the error resulting from the inaccuracy of the model parameters but also can correct the parameters itself. This gives a good method in modifying the model and improving the quality of prediction of ENSO.

Simulating leaf net CO_2 assimilation rate of C_3 & C_4 plants and its response to environmental factors

Basic structure and algorithm of leaf mechanism photosynthesis model were described in first part of this study based on former researcher results. Then, considering some environmental factors influencing on leaf photosynthesis, three numerical sensitivity experiments were carried out. We simulated the sing le leaf net CO2 assimilation, which acts as a function of different light, carbo n dioxide and temperature conditions. The relationships between leaf net photosy nthetic rate of C3 and C4 plant with CO2 concentration intercellular, leaf tempe rature, and photosynthetic active radiation (PAR) were presented, respectively. The results show the numerical experiment may indicate the main characteristic o f plant photosynthesis in C3 and C4 plant, and further can be used to integrate with the regional climate model and act as land surface process scheme, and bett er understand the interaction between vegetation and atmosphere.

SAR Data Assimilation for Crop Biomass Simulation Based on Crop Growth Model

Biomass from SAR data was assimilated into crop growth model to describe relationship between crop biomass and crop growth time to improve estimation accuracy of biomass. In addition, inverse model was established in order to estimate biomass according to relationship between biomass and backscattering coefficients from SAR data. Based on cost function, parameters of growth model were optimized as per conjugate gradient method, minimizing the differences between estimated biomass and inversion values from SAR data. The results indicated that the simulated biomass using the revised growth model with SAR data was consistent with the measured one in time distribution and even higher in accuracy than that without SAR data. Hence, the key parameters of crop growth model could be revised by real-time growth information from SAR data and accuracy of the simulated biomass could be improved accordingly.

Construction and solution of strain model along thickness of aluminum alloy plate under plastic deformation

A thickness strain model of aluminium alloy plate under plastic deformation,based on thin plate assumption was proposed.It is found that when ratio of stress fractions is constant during in-plane loading,ratios of strain components under various loading conditions are linearly related and these points of ratios form a η-η line.Under these simple loadings,strains in thickness direction can be easily calculated by the η-η line equation without integral and differential work.When the plate is under more complicated loading conditions,the thickness can be computed by the proposed optimization and piecewise calculation model.Validation computations indicate that the relative error of the results of the presented model is less than 0.75% compared with the proven theories and FE simulation.Therefore,the developed model can be applied to engineering calculation,e.g.pre-stretching analysis of aerospace aluminium thick plate,with acceptable accuracy.

Simulating leaf net CO2 assimilation rate of C3 & C4 plants and its response to environmental factors

Basic structure and algorithm of leaf mechanism photosynthesis model were described in first part of this study based on former researcher results. Then, considering some environmental factors influencing on leaf photosynthesis, three numerical sensitivity experiments were carried out. We simulated the sing le leaf net CO2 assimilation, which acts as a function of different light, carbo n dioxide and temperature conditions. The relationships between leaf net photosy nthetic rate of C3 and C4 plant with CO2 concentration intercellular, leaf tempe rature, and photosynthetic active radiation (PAR) were presented, respectively. The results show the numerical experiment may indicate the main characteristic o f plant photosynthesis in C3 and C4 plant, and further can be used to integrate with the regional climate model and act as land surface process scheme, and bett er understand the interaction between vegetation and atmosphere.

Model of limestone calcination / sulfation under oxy-fuel fluidized bed combustion

The characteristics of the simultaneous calcination/ sulfation of limestone under oxy-fuel fluidized bed combustion were studied and compared with those of the sulfation of precalcined CaO. During the calcination stage, SO2 can react with product CaO and slow down the CaCO3 decomposition rate by the covering effect of the CaSO4 product. The sulfation rate of simultaneous calcinatiort/sulfation is slower than that of precalcined CaO, but with a long enough sulfation time, the calcium conversion of simultaneous calcination/sulfation is higher than that of the precalcined CaO. A grain-micrograin model is established to describe the simultaneous calcination, sintering and sulfation of limestone. The graln-micrograln model can reflect the true reaction process of the calcination and sulfation of limestone in oxy-fuel fluidized bed combustion.

Adjoint Assimilation in Marine Ecosystem Models and an Example of Application

This paper aims at a review of the work carried out to date on the adjoint assimilation of data in marine ecosys-tem models since 1995. The structure and feature of the adjoint assimilation in marine ecosystem models are also introduced. To illustrate the application of the adjoint technique and its merits, a 4-variable ecosystem model coupled with a 3-D physical model is established for the Bohai Sea and the Yellow Sea. The chlorophyll concentration data derived from the SeaWiFS o-cean colour data are assimilated in the model with the technique. Some results are briefly presented.

Investigation on Promising Additives for Simultaneous Oxidation of NO and SO_2 in Flue Gas by Numerical Simulation

Mechanism analysis on simultaneous oxidation of NO and SO2 with additives was presented and numerical simulation was developed to investigate the performances of three additives on oxidation of NO and SO2. The simulation result showed that reaction temperature, residence time, additive dose and NO concentration influence the oxidation process significantly. There exists an optimum reaction condition for each additive, n-C4H10 has the strongest ability to oxidize NO and SO2.

Establishment and tests of EnOI assimilation module for WAVEWATCH Ⅲ

In this paper, we propose a parallel data assimilation module based on ensemble optimal interpolation （EnOI）. We embedded the method into the full-spectral third-generation wind-wave model, WAVEWATCH III Version 3.14, producing a wave data assimilation system. We present our preliminary experiments assimilating altimeter significant wave heights （SWH） using the EnOI-based wave assimilation system. Waters north of 15°S in the Indian Ocean and South China Sea were chosen as the target computational domain, which was two-way nested into the global implementation of the WAVEWATCH III. The wave model was forced by six-hourly ocean surface wind velocities from the cross-calibrated multi-platform wind vector dataset. The assimilation used along-track SWH data from the Jason-2 altimeter. We evaluated the effect of the assimilation on the analyses and hindcasts, and found that our technique was effective. Although there was a considerable mean bias in the control SWHs, a month-long consecutive assimilation reduced the bias by approximately 84% and the root mean-square error （RMSE） by approximately 65%. Improvements in the SWH RMSE for both the analysis and hindcast periods were more significant in July than January, because of the monsoon climate. The improvement in model skill persisted for up to 48 h in July. Furthermore, the SWH data assimilation had the greatest impact in areas and seasons where and when the sea-states were dominated by swells.

Numerical simulation of scatterometer assimilated wind and ocean wave in eastern China seas and adjacent waters

Using the latest version of Mesoscale Modeling System （MMSv3）, we assimilated wind data from the scatterometer and built a model to assimilate the wind field over eastern China seas and adjacent waters and applied the wave model WAVEWATCH-Ⅲ to test the sea area with assimilative wind and blended wind of QSCAT and NCEP as driving forces. High precision and resolution numerical wave results were obtained. Analysis indicated that if we replace the model wind result with the blended wind, better sea surface wind results and wave results could be obtained.

Assimilation of High Frequency Radar Data into a Shelf Sea Circulation Model

High Frequency(HF) radar current data is assimilated into a shelf sea circulation model based on optimal interpolation(OI) method. The purpose of this work is to develop a real-time computationally highly efficient assimilation method to improve the forecast of shelf current. Since the true state of the ocean is not known, the specification of background error covariance is arduous. Usually, it is assumed or calculated from an ensemble of model states and is kept in constant. In our method, the spatial covariances of model forecast errors are derived from differences between the adjacent model forecast fields, which serve as the forecast tendencies. The assumption behind this is that forecast errors can resemble forecast tendencies, since variances are large when fields change quickly and small when fields change slowly. The implementation of HF radar data assimilation is found to yield good information for analyses. After assimilation, the root-mean-square error of model decreases significantly. Besides, three assimilation runs with variational observation density are implemented. The comparison of them indicates that the pattern described by observations is much more important than the amount of observations. It is more useful to expand the scope of observations than to increase the spatial interval. From our tests, the spatial interval of observation can be 5 times bigger than that of model grid.

Homology modeling and evolutionary trace analysis of superoxide dismutase from extremophile Acidithiobacillus ferrooxidans

The gene sod in Acidithiobacillusferrooxidans may play a crucial role in its tolerance to the extremely acidic, toxic and oxidative environment of bioleaching. For insight into the anti-toxic mechanism of the bacteria, a three-dimensional （3D） molecular structure of the protein encoded by this gene was built by homology modeling techniques, refined by molecular dynamics simulations, assessed by PROFILE-3D and PROSTAT programs and its key residues were further detected by evolutionary trace analysis. Through these procedures, some trace residues were identified and spatially clustered. Among them, the residues of Asn38, Glyl03 and Glul61 are randomly scattered throughout the mapped structure; interestingly, the other residues are all distinctly clustered in a subgroup near Fe atom. From these results, this gene can be confirmed at 3D level to encode the Fe-depending superoxide dismutase and subsequently play an anti-toxic role. Furthermore, the detected key residues around Fe binding site can be conjectured to be directly responsible for Fe binding and catalytic function.

Modified Cascade Synchronization Scheme for Discrete-Time Hyperchaotic Systems

In this paper, the modified cascade synchronization scheme is proposed to investigate the synchronization in discrete-time hyperchaotic systems. By choosing a general kind of proportional scaling error functions and based on rigorous control theory, we take the discrete-time hyperchaotic system due to Wang and 3D generalized Henon map as two examples to achieve the modified cascade synchronization, respectively. Numerical simulations are used to verify the effectiveness of the proposed technique.

The REMO Ocean Data Assimilation System into HYCOM(RODAS_H):General Description and Preliminary Results

The first version of the Brazilian Oceano- graphic Modeling and Observation Network （REMO） ocean data assimilation system into the Hybrid Coordi- nate Ocean Model （HYCOM） （RODAS H） has recently been constructed for research and operational purposes. The system is based on a multivariate Ensemble Optimal Interpolation （EnOI） scheme and considers the high fre- quency variability of the model error co-variance matrix. The EnOl can assimilate sea surface temperature （SST）, satellite along-track and gridded sea level anomalies （SLA）, and vertical profiles of temperature （T） and salinity （S） from Argo. The first observing system experiment was carried out over the Atlantic Ocean （78°S-50°N, 100°W-20°E） with HYCOM forced with atmospheric reanalysis from 1 January to 30 June 2010. Five integra- tions were performed, including the control run without assimilation. In the other four, different observations were assimilated： SST only （A SST）; Argo T-S profiles only （AArgo）; along-track SLA only （A_SLA）; and all data employed in the previous runs （A_All）. The A_SST, A_Argo, and A_SLA runs were very effective in improv- ing the representation of the assimilated variables, but they had relatively little impact on the variables that were not assimilated. In particular, only the assimilation of S was able to reduce the deviation of S with respect to ob- servations. Overall, the A_All run produced a good analy- sis by reducing the deviation of SST, T, and S with respect to the control run by 39%, 18%, and 30%, respectively, and by increasing the correlation of SLA by 81%.

Bidirectional Partial Generalized Synchronization in Chaotic and Hyperchaotic Systems via a New Scheme

In this paper, a bidirectional partial generalized （lag, complete, and anticipated） synchronization of a class of continuous-time systems is defined. Then based on the active control idea, a new systematic and concrete scheme is developed to achieve bidirectional partial generalized （lag, complete, and anticipated） synchronization between two chaotic systems or between chaotic and hyperchaotic systems. With the help of symbolic-numerical computation, we choose the modified Chua system, Lorenz system, and the hyperchaotic Tamasevicius Namajunas-Cenys system to illustrate the proposed scheme. Numerical simulations are used to verify the effectiveness of the proposed scheme. It is interesting that partial chaos synchronization not only can take place between two chaotic systems, but also can take place between chaotic and hyperchaotic systems. The proposed scheme can also be extended to research bidirectional partial generalized （lag, complete, and anticipated） synchronization between other dynamical systems.

The climatology and interannual variability of the East Asian summer monsoonsimulated by a weakly coupled data assimilation system

With the motivation to improve the simulation of the East Asian summer monsoon(EASM) in coupled climate models, oceanic data assimilation(DA) was used in CAS-ESM-C(Chinese Academy of Sciences–Earth System Model–Climate Component) in this study. Observed sea surface temperature was assimilated into CAS-ESM-C. The climatology and interannual variability of the EASM simulated in CAS-ESM-C with DA were compared with a traditional AMIP-type run.Results showed that the climatological spatial pattern and annual cycle of precipitation in the western North Paci?c, and the ENSO-related and EASM-related EASM circulation and precipitation, were largely improved. As shown in this study, air–sea coupling is important for EASM simulation. In addition, oceanic DA synchronizes the coupled model with the real world without breaking the air–sea coupling process. These two successful factors make the assimilation experiment a more reasonable experimental design than traditional AMIP-type simulations.

Assimilating the LAI Data to the VEGAS Model Using the Local Ensemble Transform Kalman Filter: An Observing System Simulation Experiment

Information on the spatial and temporal patterns of surface carbon flux is crucial to understanding of source/sink mechanisms and projection of future atmospheric CO2 concentrations and climate. This study presents the construction and implementation of a terrestrial carbon cycle data assimilation system based on a dynamic vegetation and terrestrial carbon model Vegetation-Global-Atmosphere-Soil(VEGAS) with an advanced assimilation algorithm, the local ensemble transform Kalman filter(LETKF, hereafter LETKF-VEGAS). An observing system simulation experiment(OSSE) framework was designed to evaluate the reliability of this system, and numerical experiments conducted by the OSSE using leaf area index(LAI) observations suggest that the LETKF-VEGAS can improve the estimations of leaf carbon pool and LAI significantly, with reduced root mean square errors and increased correlation coefficients with true values, as compared to a control run without assimilation. Furthermore, the LETKF-VEGAS has the potential to provide more accurate estimations of the net primary productivity(NPP) and carbon flux to atmosphere(CFta).

TAMDAR Observation Assimilation in WRF 3D-Var and Its Impact on Hurricane Ike (2008) Forecast

This study evaluates the impact of atmospheric observations from the Tropospheric Airborne Meteorological Data Reporting (TAMDAR) observing system on numerical weather prediction of hurricane Ike (2008) using three-dimensional data assimilation system for the Weather Research and Forecast (WRF) model (WRF 3D-Var). The TAMDAR data assimilation capability is added to WRF 3D-Var by incorporating the TAMDAR observation operator and corresponding observation processing procedure. Two 6-h cycling data assimilation and forecast experiments are conducted. Track and intensity forecasts are verified against the best track data from the National Hurricane Center. The results show that, on average, assimilating TAMDAR observations has a positive impact on the forecasts of hurricane Ike. The TAMDAR data assimilation reduces the track errors by about 30 km for 72-h forecasts. Improvements in intensity forecasts are also seen after four 6-h data assimilation cycles. Diagnostics show that assimilation of TAMDAR data improves subtropical ridge and steering flow in regions along Ike's track, resulting in better forecasts.

CORRECTION OF ASYMMETRIC STRENGTHENING OF QUIKSCAT WIND FIELD AND ASSIMILATION APPLICATION IN TYPHOON SIMULATION

As an approach to the technological problem that the wind data of QuikSCAT scatterometer cannot accurately describe the zone of typhoon-level strong wind speed, some objective factors such as the typhoon moving speed, direction and friction are introduced in this study to construct the asymmetric strengthening of the QuikSCAT wind field. Then by adopting a technology of four-dimensional data assimilation, an experiment that includes both the assimilation and forecasting phases is designed to simulate Typhoon Rananim numerically. The results show that with model constraints and adjustment, this technology can incorporate the QuikSCAT wind data to the entire column of the model atmosphere, improve greatly the simulating effects of the whole-column wind, pressure field and the track as well as the simulated typhoon intensity covered by the forecast phase, and work positively for the forecasting of landfall locations.