Abstract: Methods of and computer program products for predicting crop yield of a geographic region are provided. In various embodiments, a time series of satellite imagery is received. The time series of satellite imagery covers at least the geographic region during a predetermined time period. The predetermined time period comprises one or more phenology periods. A time series of weather data is received. The time series of weather data covers at least the geographic region during the predetermined time period. At least one surface feature of the geographic region during each of the one or more phenology periods is generated from the time series of satellite imagery. At least one weather feature of the geographic region during each of the one or more phenology periods is generated from the time series of weather data. The at least one surface feature and the at least one weather feature are provided to a trained model. A prediction of crop yield for the geographical region is received from the trained model.