Abstract: A two dimensional nonlinear PBL numerical model using an energy closure (E-ε)method has been used to study the sea breeze circulation and TIBL in coastal areas.The main characteristics of sea breeze obtained from numerical experiments agree with those from general observation facts that are summarized by Atkinson(1987).The depth of sea breeze ranges from 300 to 900m, maximum velocity from 1.5 to 4m/s, and the height where max.velocity appears from 100 to 300m.The agreement between comparisons reveals that the performance of model is good,and the chosen experiment conditions are reasonable.This paper refitted the function of TIBL profiles also using the Weisman's relationship,but exponent value is changed with the different states of sea breeze.Many numerial experimental results indicate that the exponent of the TIBL profile,ranging from 0.4 to 1.1,is related with the strength and depth on the sea breeze.The exponent is 0.5 only when sea breeze is fully developed.This paper gives exponents under different sea breeze.