To mill fine and well-defined micro-dimpled structures,a machining manner of spiral trajectory tool reciprocating motion,where the tool repeats the process of‘feed milling–retract–cutting feed–feed milling again’...To mill fine and well-defined micro-dimpled structures,a machining manner of spiral trajectory tool reciprocating motion,where the tool repeats the process of‘feed milling–retract–cutting feed–feed milling again’along the spiral trajectory,was proposed.From the kinematics analysis,it is found that the machining quality of micro-dimpled structures is highly dependent on the machining trajectory using spiral trajectory tool reciprocating motion.To reveal this causation,simulation modelling and experimental studies were carried out.A simulation model was developed to quantitatively and qualitatively investigate the influence of the trajectory discretization strategies(constant-angle and constant-arc length)and parameters(discrete angle,discrete arc length,and pitch)on surface texture and residual height of micro-dimpled structures.Subsequently,micro-dimpled structures were milled under different trajectory discretization strategies and parameters with spiral trajectory tool reciprocating motion.A comprehensive comparison between the milled results and simulation analysis was made based on geometry accuracy,surface morphology and surface roughness of milled dimples.Meanwhile,the errors and factors affecting the above three aspects were analyzed.The results demonstrate both the feasibility of the established simulation model and the machining capability of this machining way in milling high-quality micro-dimpled structures.Spiral trajectory tool reciprocating motion provides a new machining way for milling micro-dimpled structures and micro-dimpled functional surfaces.And an appropriate machining trajectory can be generated based on the optimized trajectory parameters,thus contributing to the improvement of machining quality and efficiency.展开更多
The Erenhot-Huailai zone, as an important dust emission source area in northern China, affects the air quality of Beijing City, Tianjin City, and Hebei Province and human activities in this zone have a profound impact...The Erenhot-Huailai zone, as an important dust emission source area in northern China, affects the air quality of Beijing City, Tianjin City, and Hebei Province and human activities in this zone have a profound impact on surface dust emission. In order to explore the main source areas of surface dust emission and quantify the impacts of human activities on surface dust emission, we investigated the surface dust emission of different land types on the Erenhot-Huailai zone by model simulation, field observation, and comparative analysis. The results showed that the average annual inhalable atmospheric particles(PM_(10)) dust emission fluxes in arid grassland, Hunshandake Sandy Land, semi-arid grassland,semi-arid agro-pastoral area, dry sub-humid agro-pastoral area, and semi-humid agro-pastoral area were 4.41, 0.71, 3.64, 1.94, 0.24, and 0.14 t/hm^(2), respectively, and dust emission in these lands occurred mainly from April to May. Due to the influence of human activities on surface dust emission, dust emission fluxes from different land types were 1.66–4.41 times greater than those of their background areas, and dust emission fluxes from the main dust source areas were 1.66–3.89 times greater than those of their background areas. According to calculation, the amount of PM_(10) dust emission influenced by human disturbance accounted for up to 58.00% of the total dust emission in the study area. In addition, the comparative analysis of model simulation and field observation results showed that the simulated and observed dust emission fluxes were relatively close to each other, with differences ranging from 0.01 to 0.21 t/hm^(2) in different months, which indicated that the community land model version 4.5(CLM4.5) had a high accuracy. In conclusion, model simulation results have important reference significance for identifying dust source areas and quantifying the contribution of human activities to surface dust emission.展开更多
Vulnerability assessment is a systematic process to identify security gaps in the design and evaluation of physical protection systems.Adversarial path planning is a widely used method for identifying potential vulner...Vulnerability assessment is a systematic process to identify security gaps in the design and evaluation of physical protection systems.Adversarial path planning is a widely used method for identifying potential vulnerabilities and threats to the security and resilience of critical infrastructures.However,achieving efficient path optimization in complex large-scale three-dimensional(3D)scenes remains a significant challenge for vulnerability assessment.This paper introduces a novel A^(*)-algorithmic framework for 3D security modeling and vulnerability assessment.Within this framework,the 3D facility models were first developed in 3ds Max and then incorporated into Unity for A^(*)heuristic pathfinding.The A^(*)-heuristic pathfinding algorithm was implemented with a geometric probability model to refine the detection and distance fields and achieve a rational approximation of the cost to reach the goal.An admissible heuristic is ensured by incorporating the minimum probability of detection(P_(D)^(min))and diagonal distance to estimate the heuristic function.The 3D A^(*)heuristic search was demonstrated using a hypothetical laboratory facility,where a comparison was also carried out between the A^(*)and Dijkstra algorithms for optimal path identification.Comparative results indicate that the proposed A^(*)-heuristic algorithm effectively identifies the most vulnerable adversarial pathfinding with high efficiency.Finally,the paper discusses hidden phenomena and open issues in efficient 3D pathfinding for security applications.展开更多
With the aid of multi-agent based modeling approach to complex systems, the hierarchy simulation models of carrier-based aircraft catapult launch are developed. Ocean, carrier, aircraft, and atmosphere are treated as ...With the aid of multi-agent based modeling approach to complex systems, the hierarchy simulation models of carrier-based aircraft catapult launch are developed. Ocean, carrier, aircraft, and atmosphere are treated as aggregation agents, the detailed components like catapult, landing gears, and disturbances are considered as meta-agents, which belong to their aggregation agent. Thus, the model with two layers is formed i.e. the aggregation agent layer and the meta-agent layer. The information communication among all agents is described. The meta-agents within one aggregation agent communicate with each other directly by information sharing, but the meta-agents, which belong to different aggregation agents exchange their information through the aggregation layer first, and then perceive it from the sharing environment, that is the aggregation agent. Thus, not only the hierarchy model is built, but also the environment perceived by each agent is specified. Meanwhile, the problem of balancing the independency of agent and the resource consumption brought by real-time communication within multi-agent system （MAS） is resolved. Each agent involved in carrier-based aircraft catapult launch is depicted, with considering the interaction within disturbed atmospheric environment and multiple motion bodies including carrier, aircraft, and landing gears. The models of reactive agents among them are derived based on tensors, and the perceived messages and inner frameworks of each agent are characterized. Finally, some results of a simulation instance are given. The simulation and modeling of dynamic system based on multi-agent system is of benefit to express physical concepts and logical hierarchy clearly and precisely. The system model can easily draw in kinds of other agents to achieve a precise simulation of more complex system. This modeling technique makes the complex integral dynamic equations of multibodies decompose into parallel operations of single agent, and it is convenient to expand, maintain, and reuse the program codes.展开更多
Validating simulation model is one of the important aspects for modeling and simulation. Some methods of validating model are compared and analyzed. Several typical methods, such as TIC’s inequality coefficient, gray...Validating simulation model is one of the important aspects for modeling and simulation. Some methods of validating model are compared and analyzed. Several typical methods, such as TIC’s inequality coefficient, gray interconnected analysis, direct spectrum estimation, maximum entropy spectral estimation based on Burg or Marple, are chosen and programmed in C language. Some examples by using the program are given. The results show that the program is available and it is best to adopt multi methods for validating models.展开更多
AIGaN/GaN HEMTs are investigated by numerical simulation from the self-consistent solution of Schr6dinger-Poisson-hydrodynamic （HD） systems. The influences of polarization charge and quantum effects are considered i...AIGaN/GaN HEMTs are investigated by numerical simulation from the self-consistent solution of Schr6dinger-Poisson-hydrodynamic （HD） systems. The influences of polarization charge and quantum effects are considered in this model. Then the two-dimensional conduction band and electron distribution, electron temperature characteristics, Id versus Vd and Id versus Vg, transfer characteristics and transconductance curves are obtained. Corresponding analysis and discussion based on the simulation results are subsequently given.展开更多
The epidermal growth factor receptor（EGFR）—tyrosine kinase inhibitors（TKIs） monotherapies have limited efficacy in the treatment of EGFR mutation-negative non-small cell lung cancers（NSCLCs）. In the present stu...The epidermal growth factor receptor（EGFR）—tyrosine kinase inhibitors（TKIs） monotherapies have limited efficacy in the treatment of EGFR mutation-negative non-small cell lung cancers（NSCLCs）. In the present study, we aimed to investigate the combined effect of erlotinib（ER） and cabozantinib（CAB） on NSCLC cell lines harboring wild-type EGFR and to optimize the dosage regimens using pharmacodynamic（PD） modeling and simulation. Therefore, we examined the combined effect of ER and CAB on cell viability, cloning, apoptosis induction, migration and growth dynamics in H1299 and A549 cells. PD modeling and simulation were also performed to quantitatively describe the H1299 cells growth dynamics and to optimize the dosage regimens as well. Our results showed that CAB effectively enhanced the sensitivity of both cell lines to ER. The PD models fitted the data well, and some important parameters were obtained. The exponential（λ_0） and linear（λ_1） growth rates of H1299 cells were 0.0241 h^（–1） and 360 cells？h^（–1）, respectively. The Emax of ER and CAB was 0.0091 h^（–1） and 0.0085 h^（–1）, and the EC50 was 0.812 μM and 1.16 μM, respectively. The synergistic effect observed in the experiments was further confirmed by the estimated combination index φ（1.37）,（95% confidence interval： 1.24–1.50）, obtained from PD modeling. Furthermore, the dosage regimens were optimized using simulations. In summary, both the experimental and modeling results demonstrated the synergistic interaction between ER and CAB in NSCLCs without EGFR mutations. Sequential combinations of ER and CAB provided an option for the therapy of the NSCLCs with wild-type EGFR, which would provide some references for preclinical study and translational research as well.展开更多
A deep understanding of crop-water eco-physiological relations is the basis for quantifying plant physiological responses to soil water stress. Pot experiments were conducted to investigate the winter wheat crop-water...A deep understanding of crop-water eco-physiological relations is the basis for quantifying plant physiological responses to soil water stress. Pot experiments were conducted to investigate the winter wheat crop-water relations under both drought and waterlogging conditions in two sequential growing seasons from 2000 to 2002, and then the data were used to develop and validate models simulating the responses of winter wheat growth to drought and waterlogging stress. The experiment consisted of four treatments, waterlogging (keep 1 to 2 cm water layer depth above soil surface), control (70%-80% field capacity), light drought (40%-50% field capacity) and severe drought (30%-40% field capacity) with six replicates at five stages in the 2000-2001 growth season. Three soil water content treatments (waterlogging, control and drought) with two replicates were designed in the 2001-2002 growth season. Waterlogging and control treatments are the same as in the 2000-2001 growth season. For the drought treatment, no water was supplied and the soil moisture decreased from field capacity to wilting point. Leaf net photosynthetic rate, transpiration rate, predawn leaf water potential, soil water potential, soil water content and dry matter weight of individual organs were measured. Based on crop-water eco-physiological relations, drought and waterlogging stress factors for winter wheat growth simulation model were put forward. Drought stress factors integrated soil water availability, the sensitivity of different development stages and the difference between physiological processes (such as photosynthesis, transpiration and partitioning). The quantification of waterlogging stress factor considered different crop species, soil water status, waterlogging days and sensitivity at different growth stages. Data sets from the pot experiments revealed favorable performance reliability for the simulation sub-models with the drought and waterlogging stress factors.展开更多
A novel simulation program with an integrated circuit emphasis（SPICE） model developed for trench-gate metal-oxide-semiconductor field-effect transistor（M OSFET）devices is proposed. The drift region resistance was ...A novel simulation program with an integrated circuit emphasis（SPICE） model developed for trench-gate metal-oxide-semiconductor field-effect transistor（M OSFET）devices is proposed. The drift region resistance was modeled according to the physical characteristics and the specific structure of the trench-gate MOSFET device. For the accurate simulation of dynamic characteristics, three important capacitances, gate-to-drain capacitance Cgd, gate-to-source capacitance Cgsand drain-to-source capacitance Cds, were modeled, respectively, in the proposed model. Furthermore,the self-heating effect, temperature effect and breakdown characteristic were taken into account; the self-heating model and breakdown model were built in the proposed model; and the temperature parameters of the model were revised. The proposed model is verified by experimental results, and the errors between measured data and simulation results of the novel model are less than 5%. Therefore, the model can give an accurate description for both the static and dynamic characteristics of the trench-gate MOSFET device.展开更多
Using physical simulation models, rainfall-induced landslides have been simulated under various rainfall intensities. During these simulations, we have monitored the physical and mechanical behaviors of the landslide ...Using physical simulation models, rainfall-induced landslides have been simulated under various rainfall intensities. During these simulations, we have monitored the physical and mechanical behaviors of the landslide over the slip surface at different heights of the model slopes, as well as taking the whole slope to identify its deformation and failure processes. The results show that the rainfall duration corresponding to the initiation of the debris landslide and is exponentially related to rainfall intensity. Corresponding to the three intervals of the rainfall intensity, there are three types of slope failure modes:(1) the small-slump failure at the leading edge of the slope;(2) the block-slump failure but sometimes there are large blocks sliding down;and(3) the bulk failure but sometimes there is the block-slump failure. Based on the total rainfall-lasting time and the associated proportion of failed mass volume, the early warning of debris landslide can be classified into five grades, i.e., red, orange to red, orange, yellow to orange and yellow, which correspond to the five slope failure modes, respectively.展开更多
Since remote sensing can provide information on the actual status of an agricultural crop, the integration between remote sensing data and crop growth simulation models has become an important trend for yield estimati...Since remote sensing can provide information on the actual status of an agricultural crop, the integration between remote sensing data and crop growth simulation models has become an important trend for yield estimation and prediction.The main objective of this research was to combine a rice growth simulation model with remote sensing data to estimate rice grain yield for different growing seasons leading to an assessment of rice yield at regional levels. Integration between NOAA (National Oceanic and Atmospheric Administration) AVHRR (Advanced Very High Resolution Radiometer) data and the rice growth simulation model ORYZA1 to develop a new software, which was named as Rice-SRS Model, resulted in accurate estimates for rice yield in Shaoxing, China, with an estimation error reduced to 1.03% and 0.79% over-estimation and 0.79% under-estimation for early, single and late season rice, respectively. Selecting suitable dates for remote sensing images was an important factor which could influence estimation accuracy. Thus, given the different growing periods for each rice season, four images were needed for early and late rice, while five images were preferable for single season rice.Estimating rice yield using two or three images was possible, however, if images were obtained during the panicle initiation and heading stages.展开更多
Six paddy soils of Shanghai, China, were studied after 120 days of anaerobicincubation at 25 deg C and 35 deg C. Four models, the effective accumulated temperature model, theone-component first-order exponential model...Six paddy soils of Shanghai, China, were studied after 120 days of anaerobicincubation at 25 deg C and 35 deg C. Four models, the effective accumulated temperature model, theone-component first-order exponential model (the one-pool model), the two-component first-orderexponential model (the two-pool model), and the two-component first-order plus zero-orderexponential model including a constant term (the special model), were fitted to the data of observedmineral-N during incubation using non-linear regression procedures. The two-pool model and thespecial model gave the best fits amongst the four models, and parameters in the special model weremore reasonable than those in the other three. Results showed that the special model gave a betterprediction of nitrogen mineralization under flooded conditions than the other three models.展开更多
Ultra-wideband(UWB)technology is a prospective technology for high-rate transmission and accurate localization in the future communication systems.State-of-art channel modeling approaches usually divide the UWB channe...Ultra-wideband(UWB)technology is a prospective technology for high-rate transmission and accurate localization in the future communication systems.State-of-art channel modeling approaches usually divide the UWB channel into several sub-band channels and model them independently.By considering frequency-dependent channel parameters,a novel analytical UWB channel model with continuous frequency response is proposed.The composite effect of all frequency components within the UWB channel on the channel impulse response(CIR)of delay domain is derived based on the continuous channel transfer function(CTF)of frequency domain.On this basis,a closed-form simulation model for UWB channels and geometry-based parameter calculation method are developed,which can guarantee the continuity of channel characteristics on the frequency domain and greatly reduce the simulation complexity.Finally,the proposed method is applied to generate UWB channel with 2 GHz bandwidth at sub-6GHz and millimeter wave(mmWave)bands,respectively.The channel measurements are also carried out to validate the proposed method.The simulated CIR and power gain are shown to be in good agreement with the measurement data.Moreover,the comparison results of power gain and Doppler power spectral density(DPSD)show that the proposed UWB channel model Received:Apr.23,2022 Revised:Jun.09,2022 Editor:Wei Fan achieves a good balance between the simulation accuracy and efficiency.展开更多
The lack and pollution of water resource make wastewater reuse necessary. The pilot scale long-term tests for submerged membrane bioreactor were conducted to treat the effluents of anaerobic or aerobic treatment proce...The lack and pollution of water resource make wastewater reuse necessary. The pilot scale long-term tests for submerged membrane bioreactor were conducted to treat the effluents of anaerobic or aerobic treatment process for the high-strength Chinese traditional medicine wastewater. This article was focused on the feasibility of the wastewater treatment and reuse at shorter hydraulic retention time （HRT） of 5.0, 3.2 and 2.13 h. MLSS growth, membrane flux, vacuum values and chemical cleaning periods were also investigated. The experimental results of treating two-phase anaerobic treatment effluent demonstrated that the CODfilt was less than 100 mg/L when the influent COD was between 500-10000 mg/L at HRT of 5.0 h, which could satisfy the normal discharged standard in China. The experimental results to treat cross flow aerobic reactor effluent demonstrated that the average value of CODfilt was 17.28 mg/L when the average value of influent COD was 192.84 mg/L at HRT of 2.13 h during 106 d, which could completely meet the normal standard for water reuse. The maximum MLSS and MLVSS reached 24000 and 14500 mg/L at HRT of 3.2 h respectively. Membrane flux had maximal resume degrees of 94.7% at vacuum value of 0.02 MPa after cleaning. Chemical cleaning periods of membrane module were 150 d. A simulation model of operational parameters was also established based on the theory of back propagation neural network and linear regression of traditional mathematical model. The simulation model showed that the optimum operational parameters were suggested as follows： HRT was 5.0 h, SRT was 100 d, the range of COD loading rate was between 10.664-20.451 kg/（m3.d）, the range of MLSS was between 7543-13694 mg/L.展开更多
Crop performance is determined by the combined effects of the genotype of the crop and the environmental conditions of the production system. This study was undertaken to develop a dynamic model for simulating environ...Crop performance is determined by the combined effects of the genotype of the crop and the environmental conditions of the production system. This study was undertaken to develop a dynamic model for simulating environmental （temperature and solar radiation） and N supply effects on fiber fineness, maturity and micronaire. Three different experiments involving genotypes, sowing dates, and N fertilization rates were conducted to support model development and model evaluation. The growth and development duration of fiber fineness, maturity, and micronaire were scaled by using physiological development time of secondary wall synthesis （PDT SWSP ）, which was determined based on the constant ratio of SWSP/ BMP. PTP （product of relative thermal effectiveness （RTE） and photosynthetically active radiation （PAR）, MJ m-2） and subtending leaf N content per unit area （N A , g m-2） and critical subtending leaf N content per unit area （CN A , g m-2） of cotton boll were calculated or simulated to evaluate effects of temperature and radiation, and N supply. Besides, the interactions among temperature, radiation and N supply were also explained by piecewise function. The overall performance of the model was calibrated and validated with independent data sets from three field experiments with two sowing dates, three or five flowering dates and three or four N fertilization rates for three subsequent years （2005, 2007, and 2009） at three ecological locations. The average RMSE and RE for fiber fineness, maturity, and micronaire predictions were 372 m g-1 and 5.0%, 0.11 m g-1 and 11.4%, 0.3 m g-1 and 12.3%, respectively, indicating a good fit between the simulated and observed data. It appears that the model can give a reliable prediction for fiber fineness, maturity and micronaire formation under various growing conditions.展开更多
Assuming that decomposition of organic matter in soils follows the first-order kinetics reaction, a computer model was developed to simulate soil organic matter dynamics. Organic matter in soils is divided up into two...Assuming that decomposition of organic matter in soils follows the first-order kinetics reaction, a computer model was developed to simulate soil organic matter dynamics. Organic matter in soils is divided up into two parts that include incorporated organic carbon from crop residues or other organic fertilizer and soil intrinsic carbon. The incorporated organic carbon was assumed to consist of two components, labile-C and resistant-C. The model was represented by a differential equation of dCt/dt = Kt × fr × fw × fs × ct (i = 1, r,S) and an integral equation of Cit = Cio × EXP(Ki × fT × fw × fs × t). Effect of soil parameters of temperature, moisture and texture on the decomposition was functioned by the fT, fw and fs, respectively. Data from laboratory incubation experiments were used to determine the first-order decay rate K, and the fraction of labile-C of crop residues by employing a nonlinear method. The values of K for the components of labile-C and resistant-C and the soil intrinsic carbon were evaluated to be 0.025,0.080 X 10-2 and 0.065 X 10-3d-1, respectively. The labile-C fraction of wheat straw, wheat roots, rice straw and rice roots were 0.50, 0.25, 0.40 and 0.20, respectively. These values are related to the initial residue carbon-to-nitrogen ratio (C/N) and lignin content.展开更多
Plant temperature （Tp） and its relations to the microclimate of rice colony and irrigation water were studied using a thermo-sensitive genic male sterile （TGMS） rice line, Pei＇ai 64S. Significant differences in t...Plant temperature （Tp） and its relations to the microclimate of rice colony and irrigation water were studied using a thermo-sensitive genic male sterile （TGMS） rice line, Pei＇ai 64S. Significant differences in the daily change of temperature were detected between Tp and air temperature at the height of 150 cm （TA）. From 8：00 to 20：00, Tp was lower than TA, but they were similar during 21：00 to next 7：00. The maximum Tp occurred one hour earlier than the maximum TA, though they both reached the minimum at 6：00. Tp fluctuated less than TA. At the same height, during 6：00-13：00, Tp was higher than air temperature （Ta）, and Tp reached the maximum one hour earlier than Ta. During the rest time on sunny day, Tp was close to or even a little lower than Ta. On overcast day, Tp was higher than Ta in the whole day, and both maximized at the same time. In addition, Tp was regulated by solar radiation, cloudage and wind speed in daytime, and by irrigation water at night. The present study indicated that a TA of 29.6℃ was the critical point, at which Tp was increased or decreased by irrigation water. Tp and the difference between water and air temperatures showed a conic relation. Tp fluctuation was also regulated by the absorption or reflection of solar radiation by leaves during daytime and release of heat energy during nighttime. By analysis on correlation and regression simulation, two models of Tp were established.展开更多
A simulation model for phasic and phenological development of rice was developed using the scale of physiological development time, based on the ecophysiological development processes. The interaction of daily thermal...A simulation model for phasic and phenological development of rice was developed using the scale of physiological development time, based on the ecophysiological development processes. The interaction of daily thermal effectiveness, photoperiod effectiveness and intrinsic earliness(before heading), and basic filling duration factor(after heading)determined the daily physiological effectiveness, which accumulated to get physiological development time. The Beta and quadratic functions were used to describe daily thermal and photoperiod effectiveness, respectively. Five specific genetic parameters were added to adjust the genotypic differences in rice development so that all different varieties could reach the same physiological development time at a given development stage. The stages of seedling emergence, panicle initiation, heading, and maturity were validated using sowing dates under different ecological environments, with the RMSE of 1. 47, 5. 10, 4.58 and 3.37 days, respectively. The results showed that the model was not only explanatory and systematic but also accurate and applicable.展开更多
Computer simulation models may by used to gain further information about missile performance variability. Model validation is an important aspect of the test program for a missile system. Validation provides a basis f...Computer simulation models may by used to gain further information about missile performance variability. Model validation is an important aspect of the test program for a missile system. Validation provides a basis for confidence in the model's results and is a necessary step if the model is to be used to draw inference about the behavior of the real missile. This paper is a review of methods useful for validation of computer simulation models of missile systems and provides a new method with high degree of confidence for validation of computer simulation models of missile systems. Some examples of the use of the new method in validating computer simulation models are given.展开更多
The direct climatic effect of aerosols for the 1980-2000 period over East Asia was numerically investigated by a regional scale coupled climate-chemistry/ aerosol model, which includes major anthropogenic aerosols （s...The direct climatic effect of aerosols for the 1980-2000 period over East Asia was numerically investigated by a regional scale coupled climate-chemistry/ aerosol model, which includes major anthropogenic aerosols （sulfate, black carbon, and organic carbon） and natural aerosols （soil dust and sea salt）. Anthropogenic emissions used in model simulation are from a global emission inventory prepared for the Intergovernmental Panel on Climate Change Fifth Assessment Report （IPCC AR5）, whereas natural aerosols are calculated online in the model. The simulated 20-year average direct solar radiative effect due to aerosols at the surface was estimated to be in a range of-9- -33 W m-2 over most areas of China, with maxima over the Gobi desert of West China, and-12 W m-2 to -24 W m-2 over the Sichuan Basin, the middle and lower reaches of the Yellow River and the Yangtze River. Aerosols caused surface cooling in most areas of East Asia, with maxima of-0.8℃ to -1.6℃ over the deserts of West China, the Sichuan Basin, portions of central China, and the middle reaches of the Yangtze River. Aerosols induced a precipitation decrease over almost the entire East China, with maxima of-90 mm/year to -150 mm/year over the Sichuan Basin, the middle reaches of the Yangtze River and the lower reaches of the Yellow River. Interdecadal variation of the climate response to the aerosol direct radiative effect is evident, indicating larger decrease in surface air temperature and stronger per- turbation to precipitation in the 1990s than that in the 1980s, which could be due to the interdecadal variation of anthropogenic emissions.展开更多
基金co-supported the National Natural Science Foundation of China(No.52235010)the Heilongjiang Postdoctoral Fund(No.LBH-Z22136)the New Era Longjiang Excellent Master and Doctoral Dissertation Fund(No.LJYXL2022-057).
文摘To mill fine and well-defined micro-dimpled structures,a machining manner of spiral trajectory tool reciprocating motion,where the tool repeats the process of‘feed milling–retract–cutting feed–feed milling again’along the spiral trajectory,was proposed.From the kinematics analysis,it is found that the machining quality of micro-dimpled structures is highly dependent on the machining trajectory using spiral trajectory tool reciprocating motion.To reveal this causation,simulation modelling and experimental studies were carried out.A simulation model was developed to quantitatively and qualitatively investigate the influence of the trajectory discretization strategies(constant-angle and constant-arc length)and parameters(discrete angle,discrete arc length,and pitch)on surface texture and residual height of micro-dimpled structures.Subsequently,micro-dimpled structures were milled under different trajectory discretization strategies and parameters with spiral trajectory tool reciprocating motion.A comprehensive comparison between the milled results and simulation analysis was made based on geometry accuracy,surface morphology and surface roughness of milled dimples.Meanwhile,the errors and factors affecting the above three aspects were analyzed.The results demonstrate both the feasibility of the established simulation model and the machining capability of this machining way in milling high-quality micro-dimpled structures.Spiral trajectory tool reciprocating motion provides a new machining way for milling micro-dimpled structures and micro-dimpled functional surfaces.And an appropriate machining trajectory can be generated based on the optimized trajectory parameters,thus contributing to the improvement of machining quality and efficiency.
基金supported by the National Basic Research Program of China (2016YFA0601901)Basic Scientific Research of Henan Academy of Sciences (240601083)Joint Fund of Henan Province Science and Technology Research and Development Program (225200810047)。
文摘The Erenhot-Huailai zone, as an important dust emission source area in northern China, affects the air quality of Beijing City, Tianjin City, and Hebei Province and human activities in this zone have a profound impact on surface dust emission. In order to explore the main source areas of surface dust emission and quantify the impacts of human activities on surface dust emission, we investigated the surface dust emission of different land types on the Erenhot-Huailai zone by model simulation, field observation, and comparative analysis. The results showed that the average annual inhalable atmospheric particles(PM_(10)) dust emission fluxes in arid grassland, Hunshandake Sandy Land, semi-arid grassland,semi-arid agro-pastoral area, dry sub-humid agro-pastoral area, and semi-humid agro-pastoral area were 4.41, 0.71, 3.64, 1.94, 0.24, and 0.14 t/hm^(2), respectively, and dust emission in these lands occurred mainly from April to May. Due to the influence of human activities on surface dust emission, dust emission fluxes from different land types were 1.66–4.41 times greater than those of their background areas, and dust emission fluxes from the main dust source areas were 1.66–3.89 times greater than those of their background areas. According to calculation, the amount of PM_(10) dust emission influenced by human disturbance accounted for up to 58.00% of the total dust emission in the study area. In addition, the comparative analysis of model simulation and field observation results showed that the simulated and observed dust emission fluxes were relatively close to each other, with differences ranging from 0.01 to 0.21 t/hm^(2) in different months, which indicated that the community land model version 4.5(CLM4.5) had a high accuracy. In conclusion, model simulation results have important reference significance for identifying dust source areas and quantifying the contribution of human activities to surface dust emission.
基金supported by the fundings from 2024 Young Talents Program for Science and Technology Thinking Tanks(No.XMSB20240711041)2024 Student Research Program on Dynamic Simulation and Force-on-Force Exercise of Nuclear Security in 3D Interactive Environment Using Reinforcement Learning,Natural Science Foundation of Top Talent of SZTU(No.GDRC202407)+2 种基金Shenzhen Science and Technology Program(No.KCXFZ20240903092603005)Shenzhen Science and Technology Program(No.JCYJ20241202124703004)Shenzhen Science and Technology Program(No.KJZD20230923114117032)。
文摘Vulnerability assessment is a systematic process to identify security gaps in the design and evaluation of physical protection systems.Adversarial path planning is a widely used method for identifying potential vulnerabilities and threats to the security and resilience of critical infrastructures.However,achieving efficient path optimization in complex large-scale three-dimensional(3D)scenes remains a significant challenge for vulnerability assessment.This paper introduces a novel A^(*)-algorithmic framework for 3D security modeling and vulnerability assessment.Within this framework,the 3D facility models were first developed in 3ds Max and then incorporated into Unity for A^(*)heuristic pathfinding.The A^(*)-heuristic pathfinding algorithm was implemented with a geometric probability model to refine the detection and distance fields and achieve a rational approximation of the cost to reach the goal.An admissible heuristic is ensured by incorporating the minimum probability of detection(P_(D)^(min))and diagonal distance to estimate the heuristic function.The 3D A^(*)heuristic search was demonstrated using a hypothetical laboratory facility,where a comparison was also carried out between the A^(*)and Dijkstra algorithms for optimal path identification.Comparative results indicate that the proposed A^(*)-heuristic algorithm effectively identifies the most vulnerable adversarial pathfinding with high efficiency.Finally,the paper discusses hidden phenomena and open issues in efficient 3D pathfinding for security applications.
基金Aeronautical Science Foundation of China (2006ZA51004)
文摘With the aid of multi-agent based modeling approach to complex systems, the hierarchy simulation models of carrier-based aircraft catapult launch are developed. Ocean, carrier, aircraft, and atmosphere are treated as aggregation agents, the detailed components like catapult, landing gears, and disturbances are considered as meta-agents, which belong to their aggregation agent. Thus, the model with two layers is formed i.e. the aggregation agent layer and the meta-agent layer. The information communication among all agents is described. The meta-agents within one aggregation agent communicate with each other directly by information sharing, but the meta-agents, which belong to different aggregation agents exchange their information through the aggregation layer first, and then perceive it from the sharing environment, that is the aggregation agent. Thus, not only the hierarchy model is built, but also the environment perceived by each agent is specified. Meanwhile, the problem of balancing the independency of agent and the resource consumption brought by real-time communication within multi-agent system （MAS） is resolved. Each agent involved in carrier-based aircraft catapult launch is depicted, with considering the interaction within disturbed atmospheric environment and multiple motion bodies including carrier, aircraft, and landing gears. The models of reactive agents among them are derived based on tensors, and the perceived messages and inner frameworks of each agent are characterized. Finally, some results of a simulation instance are given. The simulation and modeling of dynamic system based on multi-agent system is of benefit to express physical concepts and logical hierarchy clearly and precisely. The system model can easily draw in kinds of other agents to achieve a precise simulation of more complex system. This modeling technique makes the complex integral dynamic equations of multibodies decompose into parallel operations of single agent, and it is convenient to expand, maintain, and reuse the program codes.
文摘Validating simulation model is one of the important aspects for modeling and simulation. Some methods of validating model are compared and analyzed. Several typical methods, such as TIC’s inequality coefficient, gray interconnected analysis, direct spectrum estimation, maximum entropy spectral estimation based on Burg or Marple, are chosen and programmed in C language. Some examples by using the program are given. The results show that the program is available and it is best to adopt multi methods for validating models.
文摘AIGaN/GaN HEMTs are investigated by numerical simulation from the self-consistent solution of Schr6dinger-Poisson-hydrodynamic （HD） systems. The influences of polarization charge and quantum effects are considered in this model. Then the two-dimensional conduction band and electron distribution, electron temperature characteristics, Id versus Vd and Id versus Vg, transfer characteristics and transconductance curves are obtained. Corresponding analysis and discussion based on the simulation results are subsequently given.
基金National Natural Science Foundation of China(NSFC,Grant No.81273583)
文摘The epidermal growth factor receptor（EGFR）—tyrosine kinase inhibitors（TKIs） monotherapies have limited efficacy in the treatment of EGFR mutation-negative non-small cell lung cancers（NSCLCs）. In the present study, we aimed to investigate the combined effect of erlotinib（ER） and cabozantinib（CAB） on NSCLC cell lines harboring wild-type EGFR and to optimize the dosage regimens using pharmacodynamic（PD） modeling and simulation. Therefore, we examined the combined effect of ER and CAB on cell viability, cloning, apoptosis induction, migration and growth dynamics in H1299 and A549 cells. PD modeling and simulation were also performed to quantitatively describe the H1299 cells growth dynamics and to optimize the dosage regimens as well. Our results showed that CAB effectively enhanced the sensitivity of both cell lines to ER. The PD models fitted the data well, and some important parameters were obtained. The exponential（λ_0） and linear（λ_1） growth rates of H1299 cells were 0.0241 h^（–1） and 360 cells？h^（–1）, respectively. The Emax of ER and CAB was 0.0091 h^（–1） and 0.0085 h^（–1）, and the EC50 was 0.812 μM and 1.16 μM, respectively. The synergistic effect observed in the experiments was further confirmed by the estimated combination index φ（1.37）,（95% confidence interval： 1.24–1.50）, obtained from PD modeling. Furthermore, the dosage regimens were optimized using simulations. In summary, both the experimental and modeling results demonstrated the synergistic interaction between ER and CAB in NSCLCs without EGFR mutations. Sequential combinations of ER and CAB provided an option for the therapy of the NSCLCs with wild-type EGFR, which would provide some references for preclinical study and translational research as well.
基金Project supported by the National High Technology Research and Development Program of China (863 Program) (No. 2003AA209030) High Technology Research and Development Program of Jiangsu Province (No. BG2004320) the National Natural Science Foundation
文摘A deep understanding of crop-water eco-physiological relations is the basis for quantifying plant physiological responses to soil water stress. Pot experiments were conducted to investigate the winter wheat crop-water relations under both drought and waterlogging conditions in two sequential growing seasons from 2000 to 2002, and then the data were used to develop and validate models simulating the responses of winter wheat growth to drought and waterlogging stress. The experiment consisted of four treatments, waterlogging (keep 1 to 2 cm water layer depth above soil surface), control (70%-80% field capacity), light drought (40%-50% field capacity) and severe drought (30%-40% field capacity) with six replicates at five stages in the 2000-2001 growth season. Three soil water content treatments (waterlogging, control and drought) with two replicates were designed in the 2001-2002 growth season. Waterlogging and control treatments are the same as in the 2000-2001 growth season. For the drought treatment, no water was supplied and the soil moisture decreased from field capacity to wilting point. Leaf net photosynthetic rate, transpiration rate, predawn leaf water potential, soil water potential, soil water content and dry matter weight of individual organs were measured. Based on crop-water eco-physiological relations, drought and waterlogging stress factors for winter wheat growth simulation model were put forward. Drought stress factors integrated soil water availability, the sensitivity of different development stages and the difference between physiological processes (such as photosynthesis, transpiration and partitioning). The quantification of waterlogging stress factor considered different crop species, soil water status, waterlogging days and sensitivity at different growth stages. Data sets from the pot experiments revealed favorable performance reliability for the simulation sub-models with the drought and waterlogging stress factors.
基金The National Natural Science Foundation of China(No.61604038)China Postdoctoral Science Foundation(No.2015M580376)+1 种基金the Natural Science Foundation of Jiangsu Province(No.BK20160691)Jiangsu Postdoctoral Science Foundation(No.1501010A)
文摘A novel simulation program with an integrated circuit emphasis（SPICE） model developed for trench-gate metal-oxide-semiconductor field-effect transistor（M OSFET）devices is proposed. The drift region resistance was modeled according to the physical characteristics and the specific structure of the trench-gate MOSFET device. For the accurate simulation of dynamic characteristics, three important capacitances, gate-to-drain capacitance Cgd, gate-to-source capacitance Cgsand drain-to-source capacitance Cds, were modeled, respectively, in the proposed model. Furthermore,the self-heating effect, temperature effect and breakdown characteristic were taken into account; the self-heating model and breakdown model were built in the proposed model; and the temperature parameters of the model were revised. The proposed model is verified by experimental results, and the errors between measured data and simulation results of the novel model are less than 5%. Therefore, the model can give an accurate description for both the static and dynamic characteristics of the trench-gate MOSFET device.
基金This research is financially supported by the National Natural Science Foundation of China(Nos.41807274,41630640)the Sichuan Science and Technology Program(No.2019E0R2230230)the Scientific Foundation of the Chinese Academy of Sciences(No.KFJ-STS-QYZD-172)。
文摘Using physical simulation models, rainfall-induced landslides have been simulated under various rainfall intensities. During these simulations, we have monitored the physical and mechanical behaviors of the landslide over the slip surface at different heights of the model slopes, as well as taking the whole slope to identify its deformation and failure processes. The results show that the rainfall duration corresponding to the initiation of the debris landslide and is exponentially related to rainfall intensity. Corresponding to the three intervals of the rainfall intensity, there are three types of slope failure modes:(1) the small-slump failure at the leading edge of the slope;(2) the block-slump failure but sometimes there are large blocks sliding down;and(3) the bulk failure but sometimes there is the block-slump failure. Based on the total rainfall-lasting time and the associated proportion of failed mass volume, the early warning of debris landslide can be classified into five grades, i.e., red, orange to red, orange, yellow to orange and yellow, which correspond to the five slope failure modes, respectively.
基金Project supported by the Commission of Science, Technology and Industry for National Defence, China (No.Y97# 14-6-2).
文摘Since remote sensing can provide information on the actual status of an agricultural crop, the integration between remote sensing data and crop growth simulation models has become an important trend for yield estimation and prediction.The main objective of this research was to combine a rice growth simulation model with remote sensing data to estimate rice grain yield for different growing seasons leading to an assessment of rice yield at regional levels. Integration between NOAA (National Oceanic and Atmospheric Administration) AVHRR (Advanced Very High Resolution Radiometer) data and the rice growth simulation model ORYZA1 to develop a new software, which was named as Rice-SRS Model, resulted in accurate estimates for rice yield in Shaoxing, China, with an estimation error reduced to 1.03% and 0.79% over-estimation and 0.79% under-estimation for early, single and late season rice, respectively. Selecting suitable dates for remote sensing images was an important factor which could influence estimation accuracy. Thus, given the different growing periods for each rice season, four images were needed for early and late rice, while five images were preferable for single season rice.Estimating rice yield using two or three images was possible, however, if images were obtained during the panicle initiation and heading stages.
基金Project supported by the National Key Basic Research Support Foundation (NKBRSF) of China (No. G199901881).
文摘Six paddy soils of Shanghai, China, were studied after 120 days of anaerobicincubation at 25 deg C and 35 deg C. Four models, the effective accumulated temperature model, theone-component first-order exponential model (the one-pool model), the two-component first-orderexponential model (the two-pool model), and the two-component first-order plus zero-orderexponential model including a constant term (the special model), were fitted to the data of observedmineral-N during incubation using non-linear regression procedures. The two-pool model and thespecial model gave the best fits amongst the four models, and parameters in the special model weremore reasonable than those in the other three. Results showed that the special model gave a betterprediction of nitrogen mineralization under flooded conditions than the other three models.
基金supported in part by the National Key Scientific Instrument and Equipment Development Project(No.61827801)in part by the National Natural Science Foundation of China(No.62271250)+2 种基金in part by Natural Science Foundation of Jiangsu Province(No.BK20211182)in part by the Key Technologies R&D Program of Jiangsu(Prospective and Key Technologies for Industry)under Grants BE2022067 and BE2022067-3in part by China Scholarship Council,and in part by Postgraduate Research&Practice Innovation Program of Jiangsu Province,No.KYCX220360.
文摘Ultra-wideband(UWB)technology is a prospective technology for high-rate transmission and accurate localization in the future communication systems.State-of-art channel modeling approaches usually divide the UWB channel into several sub-band channels and model them independently.By considering frequency-dependent channel parameters,a novel analytical UWB channel model with continuous frequency response is proposed.The composite effect of all frequency components within the UWB channel on the channel impulse response(CIR)of delay domain is derived based on the continuous channel transfer function(CTF)of frequency domain.On this basis,a closed-form simulation model for UWB channels and geometry-based parameter calculation method are developed,which can guarantee the continuity of channel characteristics on the frequency domain and greatly reduce the simulation complexity.Finally,the proposed method is applied to generate UWB channel with 2 GHz bandwidth at sub-6GHz and millimeter wave(mmWave)bands,respectively.The channel measurements are also carried out to validate the proposed method.The simulated CIR and power gain are shown to be in good agreement with the measurement data.Moreover,the comparison results of power gain and Doppler power spectral density(DPSD)show that the proposed UWB channel model Received:Apr.23,2022 Revised:Jun.09,2022 Editor:Wei Fan achieves a good balance between the simulation accuracy and efficiency.
基金Project supported by the Hi-Tech Research and Development Program (863) of China (No. 2002AA601310).
文摘The lack and pollution of water resource make wastewater reuse necessary. The pilot scale long-term tests for submerged membrane bioreactor were conducted to treat the effluents of anaerobic or aerobic treatment process for the high-strength Chinese traditional medicine wastewater. This article was focused on the feasibility of the wastewater treatment and reuse at shorter hydraulic retention time （HRT） of 5.0, 3.2 and 2.13 h. MLSS growth, membrane flux, vacuum values and chemical cleaning periods were also investigated. The experimental results of treating two-phase anaerobic treatment effluent demonstrated that the CODfilt was less than 100 mg/L when the influent COD was between 500-10000 mg/L at HRT of 5.0 h, which could satisfy the normal discharged standard in China. The experimental results to treat cross flow aerobic reactor effluent demonstrated that the average value of CODfilt was 17.28 mg/L when the average value of influent COD was 192.84 mg/L at HRT of 2.13 h during 106 d, which could completely meet the normal standard for water reuse. The maximum MLSS and MLVSS reached 24000 and 14500 mg/L at HRT of 3.2 h respectively. Membrane flux had maximal resume degrees of 94.7% at vacuum value of 0.02 MPa after cleaning. Chemical cleaning periods of membrane module were 150 d. A simulation model of operational parameters was also established based on the theory of back propagation neural network and linear regression of traditional mathematical model. The simulation model showed that the optimum operational parameters were suggested as follows： HRT was 5.0 h, SRT was 100 d, the range of COD loading rate was between 10.664-20.451 kg/（m3.d）, the range of MLSS was between 7543-13694 mg/L.
基金funded by the National Natural Science Foundation of China (30771277 and 30771279)
文摘Crop performance is determined by the combined effects of the genotype of the crop and the environmental conditions of the production system. This study was undertaken to develop a dynamic model for simulating environmental （temperature and solar radiation） and N supply effects on fiber fineness, maturity and micronaire. Three different experiments involving genotypes, sowing dates, and N fertilization rates were conducted to support model development and model evaluation. The growth and development duration of fiber fineness, maturity, and micronaire were scaled by using physiological development time of secondary wall synthesis （PDT SWSP ）, which was determined based on the constant ratio of SWSP/ BMP. PTP （product of relative thermal effectiveness （RTE） and photosynthetically active radiation （PAR）, MJ m-2） and subtending leaf N content per unit area （N A , g m-2） and critical subtending leaf N content per unit area （CN A , g m-2） of cotton boll were calculated or simulated to evaluate effects of temperature and radiation, and N supply. Besides, the interactions among temperature, radiation and N supply were also explained by piecewise function. The overall performance of the model was calibrated and validated with independent data sets from three field experiments with two sowing dates, three or five flowering dates and three or four N fertilization rates for three subsequent years （2005, 2007, and 2009） at three ecological locations. The average RMSE and RE for fiber fineness, maturity, and micronaire predictions were 372 m g-1 and 5.0%, 0.11 m g-1 and 11.4%, 0.3 m g-1 and 12.3%, respectively, indicating a good fit between the simulated and observed data. It appears that the model can give a reliable prediction for fiber fineness, maturity and micronaire formation under various growing conditions.
基金supported by the Hundred Talents Program,the Chinese Academy of Sciences and the Natural Science Foundation of China(No.30030090,39830220)
文摘Assuming that decomposition of organic matter in soils follows the first-order kinetics reaction, a computer model was developed to simulate soil organic matter dynamics. Organic matter in soils is divided up into two parts that include incorporated organic carbon from crop residues or other organic fertilizer and soil intrinsic carbon. The incorporated organic carbon was assumed to consist of two components, labile-C and resistant-C. The model was represented by a differential equation of dCt/dt = Kt × fr × fw × fs × ct (i = 1, r,S) and an integral equation of Cit = Cio × EXP(Ki × fT × fw × fs × t). Effect of soil parameters of temperature, moisture and texture on the decomposition was functioned by the fT, fw and fs, respectively. Data from laboratory incubation experiments were used to determine the first-order decay rate K, and the fraction of labile-C of crop residues by employing a nonlinear method. The values of K for the components of labile-C and resistant-C and the soil intrinsic carbon were evaluated to be 0.025,0.080 X 10-2 and 0.065 X 10-3d-1, respectively. The labile-C fraction of wheat straw, wheat roots, rice straw and rice roots were 0.50, 0.25, 0.40 and 0.20, respectively. These values are related to the initial residue carbon-to-nitrogen ratio (C/N) and lignin content.
基金supported by the National Natural Science Foundation of China (Grant No. 30370830)
文摘Plant temperature （Tp） and its relations to the microclimate of rice colony and irrigation water were studied using a thermo-sensitive genic male sterile （TGMS） rice line, Pei＇ai 64S. Significant differences in the daily change of temperature were detected between Tp and air temperature at the height of 150 cm （TA）. From 8：00 to 20：00, Tp was lower than TA, but they were similar during 21：00 to next 7：00. The maximum Tp occurred one hour earlier than the maximum TA, though they both reached the minimum at 6：00. Tp fluctuated less than TA. At the same height, during 6：00-13：00, Tp was higher than air temperature （Ta）, and Tp reached the maximum one hour earlier than Ta. During the rest time on sunny day, Tp was close to or even a little lower than Ta. On overcast day, Tp was higher than Ta in the whole day, and both maximized at the same time. In addition, Tp was regulated by solar radiation, cloudage and wind speed in daytime, and by irrigation water at night. The present study indicated that a TA of 29.6℃ was the critical point, at which Tp was increased or decreased by irrigation water. Tp and the difference between water and air temperatures showed a conic relation. Tp fluctuation was also regulated by the absorption or reflection of solar radiation by leaves during daytime and release of heat energy during nighttime. By analysis on correlation and regression simulation, two models of Tp were established.
基金supported by the National Natura1 Science Foundation of China(30030090)the National 863 Program,China(2001AA245041)
文摘A simulation model for phasic and phenological development of rice was developed using the scale of physiological development time, based on the ecophysiological development processes. The interaction of daily thermal effectiveness, photoperiod effectiveness and intrinsic earliness(before heading), and basic filling duration factor(after heading)determined the daily physiological effectiveness, which accumulated to get physiological development time. The Beta and quadratic functions were used to describe daily thermal and photoperiod effectiveness, respectively. Five specific genetic parameters were added to adjust the genotypic differences in rice development so that all different varieties could reach the same physiological development time at a given development stage. The stages of seedling emergence, panicle initiation, heading, and maturity were validated using sowing dates under different ecological environments, with the RMSE of 1. 47, 5. 10, 4.58 and 3.37 days, respectively. The results showed that the model was not only explanatory and systematic but also accurate and applicable.
文摘Computer simulation models may by used to gain further information about missile performance variability. Model validation is an important aspect of the test program for a missile system. Validation provides a basis for confidence in the model's results and is a necessary step if the model is to be used to draw inference about the behavior of the real missile. This paper is a review of methods useful for validation of computer simulation models of missile systems and provides a new method with high degree of confidence for validation of computer simulation models of missile systems. Some examples of the use of the new method in validating computer simulation models are given.
基金supported by the Knowledge Innovation Program of the Chinese Academy of Sciences (Grant No.KZCX2-YW-Q11-03)the"Strategic Priority Research Program"of the Chinese Academy of Sciences(Grant No. XDA05100502)+1 种基金the National Basic Research Program of China (Grant No.2010CB950804)100 Talents Program of the Chinese Academy of Sciences
文摘The direct climatic effect of aerosols for the 1980-2000 period over East Asia was numerically investigated by a regional scale coupled climate-chemistry/ aerosol model, which includes major anthropogenic aerosols （sulfate, black carbon, and organic carbon） and natural aerosols （soil dust and sea salt）. Anthropogenic emissions used in model simulation are from a global emission inventory prepared for the Intergovernmental Panel on Climate Change Fifth Assessment Report （IPCC AR5）, whereas natural aerosols are calculated online in the model. The simulated 20-year average direct solar radiative effect due to aerosols at the surface was estimated to be in a range of-9- -33 W m-2 over most areas of China, with maxima over the Gobi desert of West China, and-12 W m-2 to -24 W m-2 over the Sichuan Basin, the middle and lower reaches of the Yellow River and the Yangtze River. Aerosols caused surface cooling in most areas of East Asia, with maxima of-0.8℃ to -1.6℃ over the deserts of West China, the Sichuan Basin, portions of central China, and the middle reaches of the Yangtze River. Aerosols induced a precipitation decrease over almost the entire East China, with maxima of-90 mm/year to -150 mm/year over the Sichuan Basin, the middle reaches of the Yangtze River and the lower reaches of the Yellow River. Interdecadal variation of the climate response to the aerosol direct radiative effect is evident, indicating larger decrease in surface air temperature and stronger per- turbation to precipitation in the 1990s than that in the 1980s, which could be due to the interdecadal variation of anthropogenic emissions.
