Land use changes exert important impacts on climate, primarily through altering greenhouse gas (GHG) and surface energy fluxes. Biogeochemical models have incorporated a relatively complete suite of biogeochemical processes to simulate GHG fluxes. However, these models often lack detailed processes of surface energy exchange, limiting their ability to assess the impacts of land use change on climate. In this study, we incorporated processes of surface energy exchange into a widely used biogeochemistry model, DeNitrification‐DeComposition (DNDC), so that it can quantify both GHG and energy fluxes between the biosphere and the atmosphere. When tested against field observations for the three dominant land use types (forest, hayfield, and cornfield) in the northeastern United States, the improved DNDC successfully captured the observed fluxes of outgoing shortwave radiation, latent heat, sensible heat, net ecosystem exchange of CO2, and their differences among the three land use types. To evaluate the differences in radiative forcing among these land use types, we conducted 100‐year simulations and converted the modeled GHG fluxes to radiative forcing using an atmospheric impulse response model. Our results show that the 100‐year cumulative differences in net radiative forcing are 3.35 nW m−2 between the hayfield and forest (slight warming) and 43.2 nW m−2 between the cornfield and forest (warming) per hectare land use difference. The cooling effects of increased albedo after the conversion of forest to hayfield or cornfield (observed and modeled in recent years) are gradually offset by the warming effects of the increasing release of GHG as the forest becomes older.