Earth observation data are increasingly used to provide consistent eco-physiological information over large areas through time. Production efficiency models (PEMs) estimate Gross Primary Production (GPP) as a function of the fraction of photosynthetically active radiation absorbed by the canopy, which is derived from Earth observation. GPP can be summed over the growing season and adjusted by a crop-specific harvest index to estimate yield. Although PEMs have many advantages over other crop yield models, they are not widely used, because performance is relatively poor. Here, a new PEM is presented that addresses deficiencies for macro-scale application: Production Efficiency Model Optimized for Crops (PEMOC). It was developed by optimizing functions from the literature with GPP estimated by eddy covariance flux towers in the United States. The model was evaluated using newly developed Earth observation products and county-level yield statistics for major crops. PEMOC generally performed better at the field and county level than another commonly used PEM, the Moderate Resolution Imaging Spectroradiometer GPP (MOD17). PEMOC and MOD17 estimates of GPP had an R² and root mean squared error (RMSE) over the growing season of 0.71–0.89 (9.87–17.47 g CO₂ d⁻¹) and 0.59–0.83 (6.86–22.20 g CO₂ d⁻¹) with flux tower GPP. PEMOC produced R²s and RMSE of 0.70 (0.52), 0.60 (0.61), and 0.62 (0.59), while MOD17 produced R²s and RMSE of 0.65 (0.57), 0.53 (0.66), and 0.65 (0.57) with corn, soybean, and winter wheat crop yield anomalies. The sample size of rice was small, so yields were compared directly. PEMOC and MOD17 produced R²s and RMSE of 0.53 (3.42 t ha⁻¹) and 0.40 (4.89 t ha⁻¹). The most sizeable model improvements were seen for C₃ and C₄ crops during emergence/senescence and peak season, respectively. These improvements were attributed to C₃ and C₄ partitioning, optimized temperature and moisture constraints, and an evapotranspiration-based soil moisture index.