Basal area and dominant height growth and survival models incorporating general response to silvicultural treatments for loblolly pine stands were developed using data from various silvicultural experiments across Southern United States. Growth models for treated stands were developed by multiplying base-line growth models with modifier response functions/multipliers accounting for effects of thinning, fertilization, and control of competing vegetation. Chapman-Richards functions were used to model the base-line growth. Separate response functions to mid-rotation thinning and fertilization effects were developed. The thinning response function was based on duration and rate parameters and is sensitive to stand age at the time of thinning, time since thinning, and intensity of thinning. The fertilization response functions were based on Weibull distribution and the magnitude of responses varies with time since application of fertilizers, type of fertilizer elements applied, and rate of application. Response functions were integrated as a multiplier to base-line models. Response to early control of competing vegetation was incorporated into base-line models through multiplier factors. Multiplier factors were calculated based on growth difference between treated and untreated stands. A difference function, derived from differential equation with age, initial stand density, and site index served as the base-line survival model. The survival model was adjusted for thinning treatment by including an additional independent variable that represents thinning intensity. No adjustment was required for survival model in response to fertilization and competing vegetation control. All growth models were unbiased and had adequate performance in predicting basal area and dominant height following treatments. Models were developed to represent general growth trends in response to treatments. The response functions developed here can be viewed as general response functions.