In the late 1970s, the first operational weather satellite system had been launched, which showed utility for monitoring land greenness, vigor and vegetation productivity. Currently, 30-year satellite data from the Advanced Very High Resolution Radiometer (AVHRR) are available for monitoring land surface, atmosphere near the ground, natural disasters, and socioeconomic activities. Statistical modeling of agricultural crop yield and production was one of the applications. This paper discusses the topic, how design the new regression models of yield anomaly based on multivariate algorithms and selection of best- fit ensemble of predictors.