Gold (Au) nanodisk arrays are modeled and fabricated for tunable near infrared reflectance with high visible transmittance. Simulations are performed in the visible and near infrared regimes between 400 and 2000 nm with varying period of nanostructures. A progressive increase in period of up to three times the diameter of the nanodisks on polymer and quartz substrates is presented. As a result, an increase in visible transmittance with a decrease in near infrared reflectance peak intensity is observed. Modeling results exhibit a high average visible transmittance of 91.5% (400–700 nm) and average reflectance of 62.9% (750–2000 nm) with diameter (D) variations from D=120 nm to D=500 nm with a period (P) equal to twice the diameter, P=2D for a polymer substrate. Electric field intensities exhibit an increase with a decrease in the period between individual nanostructures for constant unit cell dimensions. Alongside simulations, electron beam lithography and evaporation processes are utilized for the patterning and fabrication. Optical characterizations including scanning electron microscopy (SEM) and atomic force microscopy (AFM) of selected metasurfaces on quartz are performed. For a specific period and diameter, the fabricated metasurfaces exhibit near infrared reflectance having a good agreement with the simulated results about plasmonic wavelengths. Such metasurfaces offer potential to be deployed for applications that require selective reflectance intensities in the near infrared range with high visible transmittance, such as controlled environment agriculture (CEA) or multi spectral near infrared sensing. With tunability over a broadband (750–2500 nm), this simple and fabrication-friendly model may be used as a theoretical guide in quantifying specific design parameters with varying light intensity configurations.