Sea surface height and transport stream function of the South China Sea from a variable-grid global ocean circulation model

A fine-grid model (1/6degrees) covering the South China Sea (SCS), East China Sea and Japan/East Sea, which is embedded into a coarse-grid (3degrees) global model, was established to study the SCS circulation. In the present paper, we report the model-produced monthly and annual mean transport stream functions and sea surface heights(SSH) and their anomalies of the SCS. Comparison to the TOPEX/Poseidon data shows that the model-produced monthly sea surface height anomalies (SSHA) are in good agreement with altimeter measurements. Based on the results, the circulation of the SCS, especially the upper layer circulation, is discussed. In the surface layer, the western Philippine Sea water intrudes into the SCS through the Luzon Strait in autumn, winter and spring, but not in summer. However, as far as the whole water column is concerned, the water intrudes into the SCS through the Luzon Strait all the year round. This indicates that in summer the water still intrudes into the SCS in the subsurface and intermediate layers. The area near the northern continental slope of the SCS is dominated by a cyclonic circulation all the year round. The SCS Southern Anticyclonic Gyre, SE Vietnam Off-Shore Current in summertime and SCS Southern Cyclonic Gyre in wintertime are reproduced reasonably. The difference between the monthly averaged SSH and SSHA is significant, indicating the importance of the mean SSH in the SCS circulation.

Diel vertical migration of zooplankton following optimal food intake under predation

A model is developed to investigate the trade-offs between benefits and costs involved in zooplanktonic diel vertical migration (DVM) strategies. The 'venturous revenue' (VR) is used as the criterion for optimal trade-offs. It is a function of environmental factors and the age of zooplankter. During vertical migration, animals are assumed to check instantaneously the variations of environmental parameters and thereby select the optimal behavioral strategy to maximize the value of VR, i.e. taking up as much food as possible with a certain risk of mortality. The model is run on a diel time scale (24 h) in four possible scenarios during the animal's life history. The results show that zooplankton can perform normal DVM balancing optimal food intake against predation risk, with the profile of DVM largely modified by the age of zooplankter.

Phenotypic plasticity of gut structure and function during periods of inactivity in Apostichopus japonicus

Apostichopus japonicus is a common sea cucumber that undergoes seasonal inactivity phases and ceases feeding during the summer months. We used this sea cucumber species as a model in which to examine phenotypic plasticity of the digestive tract in response to food deprivation. We measured the body mass, gross gut morphology and digestive enzyme activities of A. japonicus before, during, and after the period of inactivity to examine the effects of food deprivation on the gut structure and function of this animal. Individuals were sampled semi-monthly from June to November (10 sampling intervals over 178 days) across temperature changes of more than 18 degrees C. On 5 September, which represented the peak of inactivity and lack of feeding, A. japonicus decreased its body mass, gut mass and gut length by 50%, 85%, and 70%, respectively, in comparison to values for these parameters preceding the inactive period. The activities of amylase, cellulase and lipase decreased by 77%, 98%, and 35% respectively, in comparison to mean values for these enzymes in June, whereas pepsin activity increased two-fold (luring the inactive phase. Alginase and trypsin activities were variable and did not change significantly across the 178-day experiment. With the exception of amylase and cellulase, all body size indices and digestive enzyme activities recovered and even surpassed the mean values preceding the inactive phase during the latter part of the experiment (October-November). Principal Component Analysis (PCA) utilizing the digestive enzyme activity and body size index data divided the physiological state of this cucumber into four phases: an active stage, prophase of inactivity peak inactivity, and a reversion phase. These phases are all consistent with previously suggested life stages for this species, but our data provide more defined characteristics of each phase. A. japonicus clearly exhibits phenotypic plasticity (or life-cycle staging) of the digestive tract during its annual inactive period. (C) 2008 Elsevier Inc. All rights reserved.

Viewing DVM via general behaviors of zooplankton: A way bridging the success of individual and population

In this paper, we viewed the diel vertical migration (DVM) of copepod in the context of the animal's immediate behaviors of everyday concerns and constructed an instantaneous behavioral criterion effective for DVM and non-DVM behaviors. This criterion employed the function of 'venturous revenue' (VR), which is the product of the food intake and probability of the survival, to evaluate the gains and losses of the behaviors that the copepod could trade-off. The optimal behaviors are to find the optimal habitats to maximize VR. Two types of VRs are formulated and tested by the theoretical analysis and simulations. The sensed VR, monitoring the real-time changes of trade-offs and thereby determining the optimum habitat, is validated to be the effective objective function for the optimization of the behavior; whereas, the realized VR, quantifying the actual profit obtained by an optimal copepod in the sensed-VR-determined habitat, defines the life history of a specific age cohort. The achievement of a robust copepod overwintering stock through integrating the dynamics of the constituent age cohorts subjected to the instantaneous behavioral criterion for DVM clearly exemplified a possible way bridging the immediate pursuit of an individual and the end success of the population. (c) 2005 Published by Elsevier Ltd.

Optimal Fuzzy PID Controller with Incomplete Derivation and Its Simulation Research on Application of Intelligent Artificial Legs

提出了一种新的最优模糊PID控制器，它由两部分组成，即在线模糊推理机构和带有不完全微分的常规PID控制器，在模糊推理机构中，引入了三个可调节因子xp,xi和xd，其作用是进一步修改和优化模糊推理的结果，以使控制器对一个给定对象具有最优的控制效果，可调节因子的最优值采用ITAE准则及Nelder和Mead提出的柔性多面体最优搜索算法加以确定，这种PID控制器被用来控制由作者设计的智能人工腿中的一个直流电机，仿真结果表明该控制器的设计是非常有效的，它可被用于控制各种不同的对象和过程。

20世纪运动稳定性研究的三个重要结果

本文论述 2 0世纪运动稳定性理论研究的三个重要结果 :李雅普诺夫函数、谢聂稳定判据、卡利托洛夫定理 .这三个结果都是对一般的连续系统作出的 ,结论明确 ,简单实用 ,因而具有广泛的应用性

细粒沉积物多组分粒度特征及其成因机制研究

By using high-resolution laser grain size instrument Mastersizer 2000, the grain size distribution of windblown depositions (loess and sandy dunes), aqueous sediments (lake, river, riverside and foreshore sand), weathering crust, sloping materials and other fine-grain sediments are systemically measured. The multimodal characteristics of grain size distribution of these sediments are carefully studied. The standard patterns and their grain size characteristics of various sediments are systemically summarized. The discrepancies of multimodal distribution among windblown depositions, aqueous sediments and other sediments are concluded and the physical mechanisms of grain size multimodal distribution of various sediments are also discussed in this paper. The major conclusions are followed: 1. The multimodal characteristic of grain size distribution is a common feature in all sediments and results from properties of transportation medium, dynamic intensity, transportation manner and other factors. 2. The windblown depositions are controlled by aerodynamic forcing, resulting in that the median size of the predominant mode gradually decreases form sandy dunes to loess. Similarly, the aqueous sediments are impacted by dynamic forces of water currents and the median grain size of the predominant mode decreases gradually from river to lake sediments. Because the kinetic viscidity of air is lower than of water, the grain size of modes of windblown depositions is usually finer than that of corresponding modes of aqueous sediments. Typical characteristics of sediments grain size distribution of various sediments have been summarized in the paper: (1) Suspended particles which diameters are less than 75μm are dominant in loess and dust. There are three modes in loess’ grain size distribution: fine, median and coarse (the median size is <1μm、1-10μm、10-75μm, respectively). The coarse mode which percentage is larger than that of others is controlled by source distance and aerodynamic intensity of dust source areas. Some samples also have a saltation mode which median size is about 300-500μm. Our analysis demonstrates that the interaction of wind, atmospheric turbulence, and dust grain gravity along the dust transportation path results in a multimodal grain size distribution for suspended dust. Changes in the median sizes of the coarse and medium modes are related to variation in aerodynamic forcing (lift force related to vertical wind and turbulence) during dust entrainment in the source area and turbulence intensity in the depositional area. (2) There is a predominant coarse saltation mode in grain size distribution of sandy dunes, which median size is about 100-300μm and the content is larger than that of other modes. The grain size distribution curve is near axis symmetric as a standard logarithm normal function. There are some suspended particles in some samples of sandy dunes, which distribution of the fine part is similar to that of loess. Comparing with sandy samples of river sediments, the sorting property of sandy dunes is better than of river samples although both they are the saltation mode. Thus, the sorting property is a criterion to distinguish dune sands and river sands. (3) There are 5~6 modes (median size are <1μm, 1-10μm, 10-70μm, 70-150μm, 150-400μm, >400μm respectively) in grain size distribution of lacustrine sediments. The former 4 modes are suspensive and others are saltated. Lacustrine sediments can be divided into three types: lake shore facies, transitional facies and central lake facies. The grain size distributions of the three facies are distinctly different and, at the same time, the transition among three modes is also clear. In all these modes, the third mode is a criteria to identify the windblown deposition in the watershed. In lake shore sediments, suspended particles are dominant, a saltation mode sometimes occurs and the fourth mode is the most important mode. In the transitional facies, the percentage of the fourth mode decreases and that of the second mode increases from lake shore to central lake. In the central lake facies, the second mode is dominant. A higher content of the second mode indicates its position more close to the central lake. (4) The grain size distribution of river sediments is the most complex. It consist of suspension, saltation and rolling modes. In most situations, the percentage of the saltation mode is larger than that of other modes. The percentage of suspension modes of river sediments is more than of sandy dunes. The grain size distribution of river sediments indicates dynamic strength of river currents. If the fourth mode is dominant, the dynamic forcing of river is weaker, such as in river floodplain. If the five or sixth mode is dominant, the water dynamic forcing of rivers is strong. (5) Sediments can be changed by later forcing in different degree to form some complicated deposition types. In the paper, the grain size distribution of aqueous sediments of windblown deposition, windblown sediments of aqueous deposition, weathering crust and slope materials are discussed and analyzed. 3. The grain size distribution characteristics of different sediments are concluded: (1) Modal difference: Usually there are suspended and saltation modes in the windblown deposition. The third mode is dominant in loess dust and the fifth mode is predominant in sandy dunes. There are suspended, saltation and rolling particles in aqueous sediments. In lacustrine sediments, the second and fourth mode are predominant for central lake facies or lake shore facies, respectively. In river sediments, the fourth, or fifth, or sixth mode is predominant. Suspended modes: the grain size of suspended particles of windblown depositions usually is less than 75μm. The content of suspended particles is lower or none in sandy dunes. However, suspended particles of aqueous sediments may reach 150μm. Difference in grain size of suspended modes represents difference between transitional mediums and the strength of dynamic forcing. Saltation modes: the median size of saltation mode of sandy dunes fluctuates less than that of river sediments. (2) Loess dust and lacustrine sediment: Their suspended particles are clearly different. There is an obvious pit between the second and the third modes in grain size distribution of lacustrine sediments. The phenomenon doesn’t occur in loess dust. In lacustrine sediments, the second mode can be a dominant mode, such as central lake facies, and contents of the second and the third modes change reversely. However, the percentage of the third mode is always the highest in loess dust. (3) Dune Sand and fluvial sand: In these two depositions, the saltation particles are dominant and the median sizes of their saltation modes overlay in distribution range. The fifth mode of dune sand fluctuates is sorted better than that of fluvial sand. (4) Lacustrine and fluvial sediments: In lacustrine sediments, there are 5-6 modes and suspended particles can be predominant. The second mode is dominant in central lake facies and the third mode is dominant in lake shore facies. Saltation or roll modes occurred in central lake facies may indicate strong precipitation events. In fluvial sediments, saltation particles (or rolling particles) usually dominant. 4. A estimation model of lake depth is firstly established by using contents of the second, the third and the fourth modes. 5. The paleoenvironmental history of the eastern part of SongLiao basin is also discussed by analyzing the grain size distribution of Yushu loess-like sediments in Jilin. It was found that there is a tectonic movement before 40ka B.P. in SongLiao basin. After the movement, loess dust deposited in Yushu area as keerqin desert developed. In recent 2000 years, the climate became drier and more deserts activated in the eastern part of Song-Liao basin.

奇异性分析方法及其在储层预测中的应用研究

In research field of oil geophysical prospecting, reservoir prediction is refers to forecasting physical properties of petroleum reservoir by using data of seismic and well logging, it is a research which can guide oil field development. Singularities of seismic and logging data are caused by the heterogeneity of reservoir physical property. It's one of important methods that using singularity characteristics of seismic and logging data to study the reservoir physical property in recently. Among them, realization of reservoir quantitative prediction by analyzing singularity of the data and enhancing transition description of data is difficulty in method research. Based on wavelet transform and the fractal theory, the paper studied the singularity judgment criterion for seismic and logging data, not only analyzed quantitative relation between singularity data and reservoir physical property, but also applied it in practical reservoir prediction. The main achievements are: 1. A new method which provides singular points and their strength information estimation at only one single scale is proposed by Herrmann (1999). Based on that, the dissertation proposed modified algorithm which realized singularity polarity detection. 2. The dissertation introduced onset function to generalize the traditional geologic boundaries variations model which used singularity characteristics to represent the abruptness of the lithologic velocity transition. We show that singularity analysis reveals generic singularity information conducted from velocity or acoustic impedance to seismogram based on the convolution seismic-model theory. Theory and applications indicated that singularity information calculated from seismic data was a natural attribute for delineating stratigraphy boundaries due to its excellent ability in detecting detailed geologic features. We demonstrated that singularity analysis was a powerful tool to delineate stratigraphy boundaries and inverse acoustic impedance and velocity. 3. The geologic significances of logging data singularity information were also presented. According to our analysis, the positions of singularities indicate the sequence stratigraphic boundary, and there is subtle relationship between the singularity strength and sedimentary environment, meanwhile the singularity polarity used to recognize stratigraphic base-level cycle. Based on all those above, a new method which provided sedimentary cycle analysis based on the singularity information of logging data in multiple scales was proposed in this dissertation. This method provided a quantitative tool for judging interface of stratum sequence and achieved good results in the actual application.

地质规律约束下的波阻抗反演

Impedance inversion is very important in seismic technology. It is based on seismic profile. Good inversion result is derived from high quality seismic profile, which is formed using high resolution imaging resolution. High-resolution process demands that signal/noise ratio is high. It is very important for seismic inversion to improve signal/noise ratio. the main idea is that the physical parameter (wave impedance), which describes the stratigraphy directly, is achieved from seismic data expressing structural style indirectly. The solution of impedance inversion technology, which is based on convolution model, is arbitrary. It is a good way to apply the priori information as the restricted condition in inversion. An updated impedance inversion technology is presented which overcome the flaw of traditional model and highlight the influence of structure. Considering impedance inversion restricted by sedimentary model, layer filling style and congruence relation, the impedance model is built. So the impedance inversion restricted by geological rule could be realized. there are some innovations in this dissertation: 1. The best migration aperture is achieved from the included angle of time surface of diffracted wave and reflected wave. Restricted by structural model, the dip of time surface of reflected wave and diffracted wave is given. 2. The conventional method of FXY forcasting noise is updated, and the signal/noise ratio is improved. 3. Considering the characteristic of probability distribution of seismic data and geological events fully, an object function is constructed using the theory of Bayes estimation as the criterion. The mathematics is used here to describe the content of practice theory. 4. Considering the influence of structure, the seismic profile is interpreted to build the model of structure. A series of structure model is built. So as the impedance model. The high frequency of inversion is controlled by the geological rule. 5. Conjugate gradient method is selected to improve resolving process for it fit the demands of geophysics, and the efficiency of algorithm is enhanced. As the geological information is used fully, the result of impedance inversion is reasonable and complex reservoir could be forecasted further perfectly.

高斯束偏移的方法及应用

Gaussian beam is the asymptotic solution of wave equation concentred at the central ray. The Gaussian beam ray tracing method has many advantages over ray tracing method. Because of the prevalence of multipath and caustics in complex media, Kirchhoff migration usually can not get satisfactory images, but Gaussian beam migration can get better results.The Runge-Kutta method is used to carry out the raytracing, and the wavefront construction method is used to calculate the multipath wavefield. In this thesis, a new method to determine the starting point and initial direction of a new ray is proposed take advantage of the radius of curvature calculated by dynamic ray tracing method.The propagation characters of Gaussian beam in complex media are investigated. When Gaussian beam is used to calculate the Green function, the wave field near the source was decomposed in Gaussian beam in different direction, then the wave field at a point is the superposition of individual Gaussian beams.Migration aperture is the key factor for Kirchhoff migration. In this thesis, the criterion for the choice of optimum aperture is discussed taking advantage of stationary phase analysis. Two equivalent methods are proposed, but the second is more preferable.Gaussian beam migration based on dip scanning and its procedure are developed. Take advantage of the travel time, amplitude, and takeoff angle calculated by Gaussian beam method, the migration is accomplished.Using the proposed migration method, I carry out the numerical calculation of simple theoretical model, Marmousi model and field data, and compare the results with that of Kirchhoff migration. The comparison shows that the new Gaussian beam migration method can get a better result over Kirchhoff migration, with fewer migration noise and clearer image at complex structures.

双聚焦偏移速度估计与可分表示叠前弹性反演

In exploration seismology, the geologic target of oil and gas reservoir in complex medium request the high accuracy image of the structure and lithology of the medium. So the study of the prestack image and the elastic inversion of seismic wave in the complex medium come to the leading edge. The seismic response measured at the surface carries two fundamental pieces of information: the propagation effects of the medium and the reflections from the different layer boundaries in the medium. The propagation represent the low-wavenumber component of the medium, it is so-called the trend or macro layering, whereas the reflections represent the high-wavenumber component of the medium, it is called the detailed or fine layering. The result of migration velocity analysis is the resolution of the low-wavenumber component of the medium, but the prestack elastic inversion provided the resolution of the high-wavvenumber component the medium. In the dissertation, the two aspects about the migration velocity estimation and the elastic inversion have been studied.Firstly, any migration velocity analysis methods must include two basic elements: the criterion that tell us how to know whether the model parameters are correct and the updating that tell us how to update the model parameters when they are incorrect, which are effected on the properties and efficiency of the velocity estimation method. In the dissertation, a migration velocity analysis method based on the CFP technology has been presented in which the strategy of the top-down layer stripping approach are adapted to avoid the difficult of the selecting reduce .The proposed method has a advantage that the travel time errors obtained from the DTS panel are defined directly in time which is the difference with the method based on common image gather in which the residual curvature measured in depth should be converted to travel time errors.In the proposed migration velocity analysis method, the four aspects have been improved as follow:? The new parameterization of velocity model is provided in which the boundaries of layers are interpolated with the cubic spline of the control location and the velocity with a layer may change along with lateral position but the value is calculated as a segmented linear function of the velocity of the lateral control points. The proposed parameterization is suitable to updating procedure.? The analytical formulas to represent the travel time errors and the model parameters updates in the t-p domain are derived under local lateral homogeneous. The velocity estimations are iteratively computed as parametric inversion. The zero differential time shift in the DTS panel for each layer show the convergence of the velocity estimation.? The method of building initial model using the priori information is provided to improve the efficiency of velocity analysis. In the proposed method, Picking interesting events in the stacked section to define the boundaries of the layers and the results of conventional velocity analysis are used to define the velocity value of the layers? An interactive integrate software environment with the migration velocity analysis and prestack migration is built.The proposed method is firstly used to the synthetic data. The results of velocity estimation show both properties and efficiency of the velocity estimation are very good.The proposed method is also used to the field data which is the marine data set. In this example, the prestack and poststack depth migration of the data are completed using the different velocity models built with different method. The comparison between them shows that the model from the proposed method is better and improves obviously the quality of migration.In terms of the theoretical method of expressing a multi-variable function by products of single-variable functions which is suggested by Song Jian (2001), the separable expression of one-way wave operator has been studied. A optimization approximation with separable expression of the one-way wave operator is presented which easily deal with the lateral change of velocity in space and wave number domain respectively and has good approach accuracy. A new prestack depth migration algorithm based on the optimization approximation separable expression is developed and used to testing the results of velocity estimation.Secondly, according to the theory of the seismic wave reflection and transmission, the change of the amplitude via the incident angle is related to the elasticity of medium in the subsurface two-side. In the conventional inversion with poststack datum, only the information of the reflection operator at the zero incident angles can be used. If the more robust resolutions are requested, the amplitudes of all incident angles should be used.A natural separable expression of the reflection/transmission operator is represented, which is the sum of the products of two group functions. One group function vary with phase space whereas other group function is related to elastic parameters of the medium and geological structure.By employing the natural separable expression of the reflection/transmission operator, the method of seismic wave modeling with the one-way wave equation is developed to model the primary reflected waves, it is adapt to a certain extent heterogeneous media and confirms the accuracy of AVA of the reflections when the incident angle is less than 45'. The computational efficiency of the scheme is greatly high.The natural separable expression of the reflection/transmission operator is also used to construct prestack elastic inversion algorithm. Being different from the AVO analysis and inversion in which the angle gathers formed during the prstack migration are used, the proposed algorithm construct a linear equations during the prestack migration by the separable expression of the reflection/transmission operator. The unknowns of the linear equations are related to the elasticity of the medium, so the resolutions of them provided the elastic information of the medium.The proposed method of inversion is the same as AVO inversion in , the difference between them is only the method processing the amplitude via the incident angle and computational domain.

Retention behavior of hydrophobic organic chemicals as a function of temperature in soil leaching column chromatography

To study the transport mechanism of hydrophobic organic chemicals (HOCs) and the energy change in soil/solvent system, a soil leaching column chromatographic (SLCC) experiment at an environmental temperature range of 20-40 degreesC was carried out, which utilized a reference soil (SP 14696) packed column and a methanol-water (1:4 by volume ratio) eluent. The transport process quickens with the increase of column temperature. The ratio of retention factors at 30 and 40 degreesC (k'(30)/k'(40)) ranged from 1.08 to 1.36. The lower enthalpy change of the solute transfer in SLCC (from eluent to soil) than in conventional reversed-phase liquid chromatography (e.g., from eluent to C-18) is consistent with the hypothesis that HOCs were dominantly and physically partitioned between solvent and soil. The results were also verified by the linear solvation energy relationships analysis. The chief factor controlling the retention was found to be the solute solvophobic partition, and the second important factor was the solute hydrogen-bond basicity, while the least important factors were the solute polarizability-dipolarity and hydrogen-bond acidity. With the increase of temperature, the contributions of the solute solvophobic partition and hydrogen-bond basicity gradually decrease, and the latter decreases faster than the former. (C) 2002 Elsevier Science Ltd. All rights reserved.

领导行为的复杂性、角色结构及其有效性研究

As the social and organizational environment become more and more complex, the topic of leadership complexity is gaining more and more attention. So far, some critical issues in this field need further exploration, such as clarifying the theoretical framework, developing and validating the measurments and exploring the mechanism of the its effectiveness. Using BEI(Behavioral Event Interview), content analysis and EFA/CFA, ANOVA, regression analysis and other qualitative/quantitative methods, this research explored the leadership structure of Chinese enterprise managers, developed a new leadership questionnaire, investigated the differences of the leadership roles among various managerial areas and on different hirachical levels, and examined the impacts of the leadership roles and leadership complexity on different indicators of leadership effectiveness in various organizational contexts. 1,020 managers were surveyed. The followings are the main findings： First，the structure of leadership behaviors of Chinese enterprise managers included ethical model, authoritarian, producer, director, monitor, mentor, strategist, enterpriser, among which ethical model and authoritarian are the new findings in Chinese cultural context. Ethical model was characterized by presenting honesty, setting an example to others, being just and diligent. Authoritarian was characterized by showing power and arbitrariness. In addition, mentor, strategist and enterpriser incarnated some cultural features of present China. The new developed leadership questionnaire’s reliability and validity reached the criterion of standardized measurement. Second, there were significant differences of frequency of leadership behaviors among the managers at different managerial postions and hirachical levels, while the impacts of different leadership roles on different leadership effectiveness indicators were also singnificantly different. Ethical model had positive impacts on the whole performance and three indicators across task and context performance, and authoritarian’s impacts on the whole performance and department performance were negatively significant. Third, the impact of leadership complexity on the whole leadership effectiveness was positively significant, while the moderating effects of organization level and position function was not significant.

5-9岁儿童对常见食物及其功能的认知发展

Research on naïve biology investigates children spontaneous understanding of biology objects, phenomena and function. Previous researches focus mostly on biology phenomena. Little has done on organism’s function, such as eating food. Many research in this field found that children were unable to categorize food by nutrition criterion, but rely on physical cues. In order to investigate the development of children’s naïve understanding of food and to find if they can classify food by nutrition criterion, three age groups (5-year-olds, 7-year-olds, and 9-year-olds) were included in this study. Varies experimental tasks were also used to explore the children’s understanding of food and its function. The results showed as the followings: 1) A few 5-year- old children can classify food by nutrition criterion when they take the spontaneous classification task. However, more and more children can realize what make a kind of food different from another can be the nutrition it contains. 2) Kindergarteners can find the relation between food and its output. When they become older, more and more children can explain the relation by consistent theory. It can be said that 9-year-old children have already have a profound understanding of nutrition. They gradually developed naive theory of biology on nutrition level. 3) Even kindergarteners can understand the concept of “food balance”. However, with development there was a significant age increase in food balance choice. 4) Children’s knowledge of food balance grows with age, but urban and rural educational background influence cognitive performance.