This study examines the thermo-mechanical behavior of auxetic metamaterial beams enhanced by graphene origami (GOri) under spatially varying nonuniform temperature distributions (SVTD). Utilizing Timoshenko beam theory considering von-Kármánn type nonlinear strain–displacement relationship, GOri beams are modeled as layered structures. The Ritz method is employed to solve equilibrium equations, analyzing the impact of GOri distribution patterns, content, and folding degree on post-buckling and vibration paths. The effects of five SVTDs, three end conditions, and three GOri distribution patterns on buckling, post-buckling behavior, and nonlinear free vibration characteristics are explored. Findings reveal that the parabolic temperature distribution with peak temperatures at beam ends (P-MAE) results in higher critical temperatures and nonlinear free vibration frequencies. This research provides crucial insights into the design and optimization of GOri-enabled metamaterial structures in complex thermal environments, highlighting the significant influence of nonuniform temperature distributions along the beam’s length.