ice yield and plant growth response to nitrogen (N) fertilizer may vary within a field, probably due to spatially variable soil conditions. An experiment designed for studying the response of rice yield to different rates of N in combination with variable soil conditions was carried out at a field where spatial variation in soil properties, plant growth, and yield across the field was documented from our previous studies for two years. The field with area of 6,600 m2 was divided into six strips running east-west so that variable soil conditions could be included in each strip. Each strip was subjected to different N application level (six levels from 0 to 165kg/ha), and schematically divided into 12 grids (10m ~times10m~;for~;each~;grid) for sampling and measurement of plant growth and rice grain yield. Most of plant growth parameters and rice yield showed high variations even at the same N fertilizer level due to the spatially variable soil condition. However, the maximum plant growth and yield response to N fertilizer rate that was analyzed using boundary line analysis followed the Mitcherlich equation (negative exponential function), approaching a maximum value with increasing N fertilizer rate. Assuming the obtainable maximum rice yield is constrained by a limiting soil property, the following model to predict rice grain yield was obtained: Y=107651-0.4704*EXP(-0.0117*FN)*MIN(I-clay,~;Iom,~;Icec,~;ITN,~; ISi) where FN is N fertilizer rate (kg/ha), I is index for subscripted soil properties, and MIN is an operator for selecting the minimum value. The observed and predicted yield was well fitted to 1:1 line (Y=X) with determination coefficient of 0.564. As this result was obtained in a very limited condition and did not explain the yield variability so high, this result may not be applied to practical N management. However, this approach has potential for quantifying the grain yield response to N fertilizer rate under variable soil conditions and formulating the site-specific N prescription for the management of spatial yield variability in a field if sufficient data set is acquired for boundary line analysis.

For developing the site-specific fertilizer management strategies of crop, it is essential to know the spatial variability of soil factors and to assess their influence on the variability of crop growth and yield. In 2002 and 2003 cropping seasons within-field spatial variability of rice growth and yield was examined in relation to spatial variation of soil properties in the· two paddy fields having each area of ca. 6,600m2 in Suwon, Korea. The fields were managed without fertilizer or with uniform application of N, P, and K fertilizer under direct-seeded and transplanted rice. Stable soil properties such as content of clay (Clay), total nitrogen (TN), organic mater (OM), silica (Si), cation exchange capacity (CEC), and rice growth and yield were measured in each grid of 10~times10m . The two fields showed quite similar spatial variation in soil properties, showing the smallest coefficient of variation (CV) in Clay (7.6~%) and the largest in Si (21.4~%) . The CV of plant growth parameters measured at panicle initiation (PIS) and heading stage (HD) ranged from 6 to 38~% , and that of rice yield ranged from 11 to 21~% . CEC, OM, TN, and available Si showed significant correlations with rice growth and yield. Multiple linear regression model with stepwise procedure selected independent variables of N fertilizer level, climate condition and soil properties, explaining as much as 76~% of yield variability, of which 21.6~% is ascribed to soil properties. Among the soil properties, the most important soil factors causing yield spatial variability was OM, followed by Si, TN, and CEC. Boundary line response of rice yield to soil properties was represented well by Mitcherich equation (negative exponential equation) that was used to quantify the influence of soil properties on rice yield, and then the Law of the Minimum was used to identify the soil limiting factor for each grid. This boundary line approach using five stable soil properties as limiting factor explained an average of about 50~% of the spatial yield variability. Although the determination coefficient was not very high, an advantage of the method was that it identified clearly which soil parameter was yield limiting factor and where it was distributed in the field.