We present a first study of the effects of renormalization-group resummation (RGR) and leading-renormalon resummation (LRR) on the systematic errors of the unpolarized isovector nucleon generalized parton distribution in the framework of large-momentum effective theory. This work is done using lattice gauge ensembles generated by the MILC Collaboration, consisting of 2+1+1 flavors of highly improved staggered quarks with a physical pion mass at lattice spacing a≈0.09 fm and a box width L≈5.76 fm. We present results for the nucleon H and E generalized parton distributions (GPDs) with average boost momentum Pz≈2 GeV at momentum transfers Q2=[0,0.97] GeV2 at skewness ξ=0 as well as Q2∈0.23 GeV2 at ξ=0.1, renormalized in the modified minimal subtraction (MS¯) scheme at scale μ=2.0 GeV, with two- and one-loop matching, respectively. We demonstrate that the simultaneous application of RGR and LRR significantly reduces the systematic errors in renormalized matrix elements and distributions for both the zero and nonzero skewness GPDs, and that it is necessary to include both RGR and LRR at higher orders in the matching and renormalization processes. Published by the American Physical Society 2024