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.