We propose the Breathing Earth System Simulator (BESS), an upscaling approach to quantify global gross primary productivity and evapotranspiration using MODIS with a spatial resolution of 1-5 km and a temporal resolution of 8 days. This effort is novel because it is the first system that harmonizes and utilizes MODIS Atmosphere and Land products on the same projection and spatial resolution over the global land. This enabled us to use the MODIS Atmosphere products to calculate atmospheric radiative transfer for visual and near infrared radiation wave bands. Then we coupled atmospheric and canopy radiative transfer processes, with models that computed leaf photosynthesis, stomatal conductance and transpiration on the sunlit and shaded portions of the vegetation and soil. At the annual time step, the mass and energy fluxes derived from BESS showed strong linear relations with measurements of solar irradiance (r(2) = 0.95, relative bias: 8%), gross primary productivity (r(2) = 0.86, relative bias: 5%) and evapotranspiration (r(2) = 0.86, relative bias: 15%) in data from 33 flux towers that cover seven plant functional types across arctic to tropical climatic zones. A sensitivity analysis revealed that the gross primary productivity and evapotranspiration computed in BESS were most sensitive to leaf area index and solar irradiance, respectively. We quantified the mean global terrestrial estimates of gross primary productivity and evapotranpiration between 2001 and 2003 as 118 +/- 26 PgC yr(-1) and 500 +/- 104 mm yr(-1) (equivalent to 63,000 +/- 13,100 km(3) yr(-1)), respectively. BESS-derived gross primary productivity and evapotranspiration estimates were consistent with the estimates from independent machine-learning, data-driven products, but the process-oriented structure has the advantage of diagnosing sensitivity of mechanisms. The process-based BESS is able to offer gridded biophysical variables everywhere from local to the total global land scales with an 8-day interval over multiple years.