We present a simple design of a balloon-borne infrared spectral photometer that can be built and used by undergraduate students to perform an experiment demonstrating the atmospheric greenhouse effect. The experiment demonstrates that the Earth radiates heat to space in the infrared region but that the radiation at the top of the atmosphere has a much lower effective radiation temperature than at the surface of the Earth, which is the essence of the greenhouse effect. The experiment also demonstrates that the greenhouse effect is much more pronounced in molecular absorption bands than in the so-called infrared window. The thrill of putting together a balloon experiment aside, students performing this experiment also gained experience in practical applications of Planck's law.

The atmospheric greenhouse effect1 is a significant component of the Earth's climate system.2 However, the typical physics curriculum rarely includes an experimental investigation of the effect. As discussed by Bell,3 previously published benchtop demonstrations that purport to demonstrate this effect4–6 are in fact more sensitive to convective rather than radiative heat transfer, so that they fail to illustrate greenhouse effect-induced warming.7,8 One exception is Sieg et al.9 who do demonstrate the absorption of IR radiation by CO2 gas, which is one component of the greenhouse effect. To overcome the limitations of benchtop experiments, we developed an in situ experiment in the Earth's atmosphere to demonstrate this crucial aspect of the Earth's climate to students. It can be performed by undergraduates, providing that they have access to a high-altitude balloon. In our case, we benefited from the Louisiana Space Grant Consortium's Louisiana Aerospace Catalyst Experience for Students (LaACES) Program.10 

To understand how the greenhouse effect is essential to explain the temperature at the Earth's surface, let us first assume that the atmosphere does not play any role. In this first very simple model, the Earth absorbs approximately ϕ q = 235 W/m2 of heat from the Sun in the form of ultraviolet (UV), visible (VIS), and near-infrared (IR) radiation.2 In the steady state, the Earth re-radiates this heat back to space in the form of mid- to far-IR thermal radiation with an effective temperature T of 254 K, which can be inferred from the Stefan-Boltzmann law, ϕ q = 5.67 × 10 8 W m 2 K 4 T 4. In this model, the average temperature of the Earth's surface would be T 254 K. This first simplistic model, therefore, cannot account for the actual average surface temperature of the Earth, which is 287 K.11 

To be more realistic, one can devise a second model, where the Earth is in radiative equilibrium with a single atmospheric layer containing greenhouse gasses that absorb none of the UV, VIS, and near-IR radiation from the Sun, but do absorb a fraction f of the mid- to far-IR thermal radiation from the Earth. The thermal energy absorbed by the atmosphere is re-radiated equally upward to space and downward to the Earth, so that this additional heat flux raises the Earth's surface temperature. In this model, an average surface temperature of 287 K is obtained if f =77%,12 whereas the actual value is closer to f =90%.2 

Furthermore, to accurately describe the observed altitudinal dependence of the effective IR radiation temperature, one needs a more realistic model that considers the atmosphere to be composed of multiple atmospheric layers containing greenhouse gasses in radiative-convective equilibrium.13 The temperature of the atmospheric layers decreases with altitude due to adiabatic cooling of rising air parcels in response to the hydrostatic decrease in atmospheric pressure with altitude. At higher altitudes, the absorbed IR radiation is, therefore, re-radiated at lower temperatures. Observing the decrease in effective radiation temperature with altitude is the first objective of this experiment.

The second objective of this experiment is to demonstrate the spectral effect of greenhouse gasses on the effective radiation temperature. In the [8–12] μm wavelength band, the IR transmittance from the surface of the Earth to space is very high in the absence of clouds. This wavelength band is referred to as the atmospheric IR window. At slightly shorter wavelengths, IR radiation is strongly absorbed by water vapor, and at slightly longer wavelengths, IR radiation is strongly absorbed by carbon dioxide.14 These two gasses are predominantly responsible for the absorption of IR radiation in the atmosphere. In our experiment, we measure separately the effective radiation temperature in the IR window and the absorption bands. We expect that the effective radiation temperature to space will be higher in the atmospheric window than in the absorption bands, since radiation to space then originates from deeper (and warmer) layers of the atmosphere.

To compare with the experimental results, we have performed a simulation of the infrared radiation transfer in the atmosphere using MODTRAN.15,16 MODTRAN is an algorithm that models the atmosphere as a series of horizontally homogenous layers and calculates the transmittance and radiance of the atmospheric constituents in each layer in 0.2 cm−1 spectral bands over the [0.2–100] μm spectral range. Figure 1 shows the MODTRAN-simulated upward IR heat flux spectrum at the surface of the Earth and at an altitude of 25 km, using the parameter values listed in Table I.17 At the surface of the Earth, the IR heat flux spectrum displays the characteristic blackbody shape with a maximum emission near 10 μm. At an altitude of 25 km, however, the spectrum displays differential absorption with little absorption in the range from 8 μm to approximately 12 μm (except for the narrow ozone absorption band near 10 μm), but presents bands of significant absorption outside of that. Although the IR heat flux spectrum is not that of a blackbody, one can calculate an effective radiative temperature T by integrating the spectrum and comparing it to the integral of Planck's law as follows:
(1)
where h is Planck's constant, c is the speed of light, Ω is field of view of the detector expressed as a solid angle, and k is Boltzmann's constant. One can calculate an effective radiation temperature for the spectrum as a whole, or for any selected wavelength band [λmin, λmax]. Our practical method of inverting Eq. (1) to calculate the effective radiation temperature is presented in  Appendix A. Figure 2 shows the altitude dependence of the effective radiation temperature for three bands: observed (designated o) [5.5 μm, 14 μm]; IR window band (designated w) [8 μm, 12 μm]; and absorption band (designated a) [5.5 μm, 8.0 μm] + [12.0 μm, 14.0 μm]. Figure 2, therefore, models the expected results of our experiment.
Fig. 1.

Infrared heat flux spectrum at the surface of the Earth (bold) and at 25 km altitude (solid), as modeled by MODTRAN. The unshaded region indicates the IR window and the shaded region the absorption bands. The dashed curves are blackbody fits to the simulated heat flux at a 25 km altitude in the IR window (285 K) and in the absorption bands (271 K).

Fig. 1.

Infrared heat flux spectrum at the surface of the Earth (bold) and at 25 km altitude (solid), as modeled by MODTRAN. The unshaded region indicates the IR window and the shaded region the absorption bands. The dashed curves are blackbody fits to the simulated heat flux at a 25 km altitude in the IR window (285 K) and in the absorption bands (271 K).

Close modal
Table I.

MODTRAN parameter values for the simulations used in this experiment. The stratospheric ozone, water vapor, and freon scales are multiplicative factors that allow the user to adjust these concentrations relative to the standard atmospheric levels.

MODTRAN parameter Value used
CO2 (ppm)  400 
CH4 (ppm)  1.7 
Tropospheric ozone (ppb)  28 
Stratospheric ozone scale 
Water vapor scale 
Freon scale 
Temperature offset (°C)  12.8 
Locality  Mid-latitude summer 
Cloud condition  Cumulus cloud base 66 km top 2.7 km 
Altitude (km)  0.0 to 30.0 in Increments of 0.5 
Direction of observation  Looking down 
MODTRAN parameter Value used
CO2 (ppm)  400 
CH4 (ppm)  1.7 
Tropospheric ozone (ppb)  28 
Stratospheric ozone scale 
Water vapor scale 
Freon scale 
Temperature offset (°C)  12.8 
Locality  Mid-latitude summer 
Cloud condition  Cumulus cloud base 66 km top 2.7 km 
Altitude (km)  0.0 to 30.0 in Increments of 0.5 
Direction of observation  Looking down 
Fig. 2.

Altitude dependence of the effective radiation temperature as modeled using MODTRAN in the observed radiation range (To, bold), in the IR window (Tw, solid), and in the absorption band (Ta, dashed).

Fig. 2.

Altitude dependence of the effective radiation temperature as modeled using MODTRAN in the observed radiation range (To, bold), in the IR window (Tw, solid), and in the absorption band (Ta, dashed).

Close modal

Our method for observing the greenhouse effect directly was to have a balloon lift a downward-pointing infrared thermometer from the surface of the Earth up to an altitude of at least 25 km, to observe the decrease in the effective temperature of the radiation coming from the surface of the Earth. The payload measured the IR heat flux using Arduino-controlled MLX90614 infrared thermometers.18 The MLX90614 has an integrated IR filter with a passband between 5.5 and 14.0 μm. The field of view of the sensor 1.87 sr. The thermometer is factory-calibrated to read out the effective radiation temperature based on the incoming heat flux in this wavelength range and for this solid angle. In our experimental conditions, the quoted accuracy of the radiation temperature measurements is ±0.5 K. The MLX90614 also reads out the temperature of the sensor itself, Ts, which is used in our data analysis.

Figure 3 is a schematic of the instrument, and Fig. 4 shows photographs of the instrument. The instrument comprises three ML90614 sensors, designated as S1, S2, and S3. S1 and S2 both directly measure the overall IR radiation temperature directly. In front of S3 is placed an additional external IR filter which acts as a passband between 8.0 and 12.0 μm, with a transmissivity t =0.859, per the manufacturer's datasheet. Therefore, S3 measures the radiation temperature in the atmospheric IR window. The instrument also comprises a GPS receiver for time and altitude data, an SD card drive for data storage for post-flight analysis, and two LED status indicators. An Arduino Mega 2560 microcontroller controls all of these devices. Batteries provide power, and a 3-cm thick polystyrene housing provides mechanical support and thermal insulation. Because the high-altitude balloon carries several experiments suspended vertically below the balloon, the sensors are mounted at an angle of 30° with respect to the nadir direction, so as to avoid having their field of view obstructed by one of the experiments suspended below ours. As long as, we can approximate the surface of the Earth and the atmospheric strata as flat planes and, as long as, the entire field of view of the sensors remains below the horizon; the orientation of the sensors does not affect the heat flux received.19 This is due to two competing effects: The first is that blackbody surfaces radiate in accordance with Lambert's law,20 which states that the intensity of the radiation coming from a unit area decreases as the cosine of the viewing angle with respect to the surface normal. The second is that the radiating surface within the field of view of the sensor is a conic section (a circle or an ellipse) whose area increases as the reciprocal of the cosine of the viewing angle. Therefore, as the viewing angle increases, the increase in visible horizontal area exactly offsets the decrease in intensity per unit surface area, resulting in a heat flux to the sensors that is the same at all viewing angles.

Fig. 3.

Instrument schematic. Refer to Arduino Mega documentation for port assignments (Ref. 21).

Fig. 3.

Instrument schematic. Refer to Arduino Mega documentation for port assignments (Ref. 21).

Close modal
Fig. 4.

(Top) Photograph of payload interior from top with lid of housing removed. Note that S4 was not used in this experiment (Ref. 22). (Middle) Photograph of side of payload. (Bottom) Photograph of IR sensors and filters. S4 and the 8 μm longpass IR filter were not used in this experiment.

Fig. 4.

(Top) Photograph of payload interior from top with lid of housing removed. Note that S4 was not used in this experiment (Ref. 22). (Middle) Photograph of side of payload. (Bottom) Photograph of IR sensors and filters. S4 and the 8 μm longpass IR filter were not used in this experiment.

Close modal
Since both S1 and S2 measure the observed effective temperature of the IR radiation to space, To, we average these two values to calculate To,
(2)
the standard error of which is
(3)
Measurement of the effective radiation temperature in the IR window band Tw is more complicated due to the presence of the additional filter. The heat flux that S3 receives, ϕ q S 3, is the external heat flux in the infrared window band, ϕ q w, times the transmissivity, t, of the filter in the passband, plus the heat flux from the sensor itself that the filter reflects back to the sensor, ϕ q r:
(4)
The measured heat flux that S3 receives is
(5)
We determined the reflected heat flux ϕ q r as a function of the sensor temperature, Ts, by calibration. We did this by roughly reproducing, on the ground, the temperature profile that the sensor would experience during the flight. We started by placing the instrument outside in direct sunlight for 30 min while operating. This raised the sensor temperature to as high as 326 K (53 °C). We then placed the instrument in a −20 °C freezer for 30 min to lower Ts. In the freezer, Ts gradually decreased to as low as 272 K (−1 °C). We also placed a large object (a book) covered in black vinyl electrical tape (to ensure a high object emissivity) in the freezer. We positioned the instrument and the object so that the object filled the field of view of the sensors. The object temperature To was measured by S1 and S2 as in Eq. (2). We performed this measurement twice, once with the object initially at room temperature and once with the object at a higher initial temperature, accomplished by leaving the object outside in direct sunlight for 30 min alongside the instrument. We did this to ensure that we would obtain the same results for ϕ q r regardless of the object temperature, as would be expected if we were indeed measuring only the reflected heat flux from S3 itself. Finally, we solved Eq. (4) for ϕ q r and calculated the heat flux in the IR window band ϕ q w from Eq. (1), using T w = T o:
(6)
Figure 5 shows the result of calibration. A quadratic regression yields
(7)
with a correlation coefficient of r2 = 0.9967. As a check, to within the standard error of ϕ q r, the measured reflected heat flux is equal to all of the heat flux coming from the sensor that is not transmitted by the filter.
Fig. 5.

Reflected heat flux versus the sensor temperature Ts as determined by calibration (points). The quadratic fit (Eq. (7)) is shown as a dashed line.

Fig. 5.

Reflected heat flux versus the sensor temperature Ts as determined by calibration (points). The quadratic fit (Eq. (7)) is shown as a dashed line.

Close modal
Equation (4) then gives
(8)
from which we extracted the effective radiation temperature in the IR window band, Tw.
Finally, the heat flux in the absorption bands is the difference between the observed heat flux, ϕ q o, and the heat flux in the IR window, ϕ q w
(9)
and we calculate the effective radiation temperature in the absorption bands, Ta, from
(10)

The instrument was launched by high-altitude balloon, launched from 30.4773000°N longitude, 92.6936333°W latitude at 16:17:03 UTC (10:17:03 local standard time) on August 3, 2020. The surface temperature was 34 °C (307 K), and there were scattered cumulus clouds present. The instrument ascended to an altitude of 27 419 m and then returned to the surface of the Earth by parachute. Atmospheric temperature and pressure versus altitude as measured by our companion instrument, which we built and launched on the same balloon flight, are shown in Fig. 6. The data gathered during the flight for this experiment are shown in Fig. 7.

Fig. 6.

Ambient atmospheric temperature and atmospheric pressure levels versus altitude.

Fig. 6.

Ambient atmospheric temperature and atmospheric pressure levels versus altitude.

Close modal
Fig. 7.

Raw data versus altitude: TS1 (solid), TS2 (dashed), TS3 (bold solid), and TS (bold dashed). Note the higher initial (i.e., at altitude = 0 km) temperature TS3 compared to TS1 and TS2. This is due to the reflected heat flux received by S3.

Fig. 7.

Raw data versus altitude: TS1 (solid), TS2 (dashed), TS3 (bold solid), and TS (bold dashed). Note the higher initial (i.e., at altitude = 0 km) temperature TS3 compared to TS1 and TS2. This is due to the reflected heat flux received by S3.

Close modal

The results of analyzing our data using Eqs. (2)–(10) are shown in Figs. 8 and 9. Figure 8 shows that the overall radiation emission temperature, To, decreases with altitude in qualitative agreement with the results derived from our MODTRAN simulations, thus illustrating the greenhouse effect. The observed difference in radiation temperature between the surface and the highest altitude is 310 K 276 K = 34 K. The measured and modelled effective radiation temperatures are in quantitative agreement up to an altitude of 7 km. Above that altitude, the measured effective radiation temperature shows a significant deviation from the model that peaks at a deviation of about −12 K at around 15 km in altitude. This is possibly due to a loss of accuracy when the sensor temperature is changing, as implied by the datasheet for the MLX90614.18 We intend to reduce this effect in future experiments by shading the sensors before launch and by locating the sensors outside of the insulating housing. This should prevent the sensor temperature from reaching such a high temperature before launch. Indeed, waste heat from the electronics typically raises the temperature inside the housing to 50 °C before launch. In our case, as seen in Fig. 7, the sensor temperature TS was (325 K = 52 °C) at launch (altitude = 0 km). Reducing the initial sensor temperature would reduce the rate of cooling during the flight. It is also possible that atmospheric turbulence around the tropopause (at an altitude around 15 km) caused the sensor look direction to deviate from its designed direction.23 While a downward oscillation would not affect the IR flux received by the sensor, upward oscillations greater than 15° would cause at least part of the sensor field of view to rise above the horizon and reduce the IR flux received. The net effect would be to bias the IR flux measurements downward.

Fig. 8.

Overall upward radiation temperature To versus altitude (thin) compared with MODTRAN simulations (bold). Shading indicates the standard error.

Fig. 8.

Overall upward radiation temperature To versus altitude (thin) compared with MODTRAN simulations (bold). Shading indicates the standard error.

Close modal
Fig. 9.

Upward radiation temperature in the IR window Tw versus altitude (thin solid) compared with the MODTRAN simulations (bold solid), and upward radiation temperature in the absorption bands Ta versus altitude (thin dashed) compared with MODTRAN simulations (bold dashed). Shading indicates the standard error of the measurements.

Fig. 9.

Upward radiation temperature in the IR window Tw versus altitude (thin solid) compared with the MODTRAN simulations (bold solid), and upward radiation temperature in the absorption bands Ta versus altitude (thin dashed) compared with MODTRAN simulations (bold dashed). Shading indicates the standard error of the measurements.

Close modal

As seen in Fig. 9, the qualitative result that the IR radiation temperature in the IR window, Tw, is higher than the radiation temperature in the absorption bands, Ta, is clear. The radiation temperatures in the IR window band, and the absorption bands are also in quantitative agreement with their corresponding models up to an altitude of 7000 m, with significant deviations above that altitude.

The performance of the developed instrument was sufficient to demonstrate the greenhouse effect. The decrease in the Earth's effective radiation temperature with altitude was qualitatively correct and in quantitative agreement with the model predictions up to an altitude of 7 km. Additionally, we observed that Tw in the IR window was greater than the overall effective radiation temperature To and, conversely, we showed that Ta in the absorption bands was smaller than To. These results were also in rough quantitative agreement with the model predictions.

One could also adapt this instrument to perform ground-based experiments. For example, if the sensors were oriented toward the sky, the instrument would measure the IR heat flux from the atmosphere to the Earth. The IR heat flux from a clear sky is concentrated in the absorption bands, with very little radiation from the atmospheric window. This would appear as a higher effective temperature in the absorption band than in the atmospheric window. In contrast, clouds radiate with the same effective temperatures in both bands.

In performing the experiment described in this article, the group of undergraduate physics students gained a practical understanding of an important part of the global climate system as well as the experience of designing their own instrument. Planned improvements in this instrument include the addition of an array of upward-facing thermometers to measure the downward IR heat flux and, therefore, the net heat flux as a function of altitude. With those sensors, we intend to observe the difference between net heat flux in the IR window and the absorption bands. Additionally, a triple-axis accelerometer will be included to determine the sensor nadir/zenith angle and payload motion at all time. Another group of undergraduates developed such an instrument, which flew on June 12, 2021. Initial results show improved accuracy in the radiation temperature measurements, indicating that placing the temperature sensors outside of the insulating housing is the better design option.

This work was supported in part by the National Aeronautics and Space Administration (NASA) Grant and Cooperative Agreement No. NNX15AH82H with Louisiana State University and A&M College through Subaward Agreement No. PO-0000107279 with Southeastern Louisiana University.

In performing the proposed experiment, it is frequently necessary to calculate the IR heat flux from the effective radiation temperature and vice versa, using the integral of Planck's law. The integral of Planck's law does not have a closed form solution. (Widger and Woodall24 obtained an infinite-series solution for the integral of Planck's law, but that is not practical for our purpose.) Therefore, we adopted the following procedure: For each wavelength band, we integrated Planck's law numerically for all integer values of T in the [240, 320] K range, inclusive. We then fit the result to a second order polynomial as follows:
(A1)
where a0, a1, and a2 are fitting parameters. The standard error of the IR heat flux is, therefore, given by
(A2)
We easily invert Eq. (A1) using the quadratic formula to calculate the effective temperature for a given heat flux:
(A3)
The standard error of the effective temperature is, therefore,
(A4)

We list the best-fit parameters in Table II. We note that our use of a second order polynomial is only a mathematical convenience, and we attribute no physical significance to this particular functional form.

Table II.

Best-fit parameters for Eq. (A1).

Band designation a0 ± δ a0 (W/m2) a1 ± δ a1 [W/(m2 K)] a2 ± δ a2 [W/(m2 K2)] r2
O  91.83 ± 0.05  1.794 ± 0.001  0.011 78 ± 0.000 07  1.000 
W  49.87 ± 0.02  0.9365 ± 0.0006  0.005 48 ± 0.000 03  1.000 
a  41.96 ± 0.05  0.858 ± 0.001  0.006 30 ± 0.000 06  0.9998 
Band designation a0 ± δ a0 (W/m2) a1 ± δ a1 [W/(m2 K)] a2 ± δ a2 [W/(m2 K2)] r2
O  91.83 ± 0.05  1.794 ± 0.001  0.011 78 ± 0.000 07  1.000 
W  49.87 ± 0.02  0.9365 ± 0.0006  0.005 48 ± 0.000 03  1.000 
a  41.96 ± 0.05  0.858 ± 0.001  0.006 30 ± 0.000 06  0.9998 

The total cost of materials for this project—not including the balloon vehicle—as detailed in Table III was $939. We do not list vendors for materials that are widely available.

Table III.

Bill of materials for construction of the instrument.

Description Unit Price per unit Extended price
IR filter 10 000 nm bp 4000 nm (Andover Corporation)  Each  $689  $689 
Three of MLX90614 sensors  Each  $14  $42 
Two of 2CR5 batteries  Each  $7  $14 
PPTC 0.3 A fuse  Pack of 20  $8  $8 
Arduino mega  Each  $32  $32 
GPS shield  Each  $21  $21 
SD card shield  Pack of 5  $20  $20 
SD cards  Pack of 5  $7  $7 
Protoshield  Pack of 10  $17  $17 
Protoshield headers  Box  $14  $14 
Hookup wire  Box  $14  $14 
LEDs  Pack of 100  $7  $7 
Solder  Spool  $16  $16 
Gorilla glue  Bottle  $6  $6 
Styrofoam insulation  Sheet  $16  $16 
EconoKote  Roll  $16  $16 
Total  ⋯  ⋯  $939 
Description Unit Price per unit Extended price
IR filter 10 000 nm bp 4000 nm (Andover Corporation)  Each  $689  $689 
Three of MLX90614 sensors  Each  $14  $42 
Two of 2CR5 batteries  Each  $7  $14 
PPTC 0.3 A fuse  Pack of 20  $8  $8 
Arduino mega  Each  $32  $32 
GPS shield  Each  $21  $21 
SD card shield  Pack of 5  $20  $20 
SD cards  Pack of 5  $7  $7 
Protoshield  Pack of 10  $17  $17 
Protoshield headers  Box  $14  $14 
Hookup wire  Box  $14  $14 
LEDs  Pack of 100  $7  $7 
Solder  Spool  $16  $16 
Gorilla glue  Bottle  $6  $6 
Styrofoam insulation  Sheet  $16  $16 
EconoKote  Roll  $16  $16 
Total  ⋯  ⋯  $939 
1.
N.
Ekholm
, “
On the variations of the climate of the geological and historical past and their causes
,”
Q. J. R. Meteorol. Soc.
27
,
1
62
(
1901
).
2.
Climate Change 2007: The Physical Science Basis
, edited by
S.
Solomon
,
D.
Qin
,
M.
Manning
,
Z.
Chen
,
M.
Marquis
,
K. B.
Averyt
,
M.
Tignor
, and
H. L.
Miller
(
Cambridge U. P
.,
Cambridge, New York
,
2007
), pp.
96
97
.
3.
Jerry A.
Bell
, “
Benchtop global-warming demonstrations do not exemplify the atmospheric greenhouse effect, but alternatives are available
,”
J. Chem. Educ.
96
(
10
),
2352
2354
(
2019
).
4.
M.
Adelhelm
and
E.-G.
Hohn
, “
A simple demonstration of the greenhouse effect
,”
J. Chem. Educ.
70
,
73
74
(
1993
).
5.
S. A.
Finkenstaedt-Quinn
,
N. V.
Hudson-Smith
,
M. J.
Styles
,
M. K.
Maudal
,
A. R.
Juelfs
, and
C. L.
Haynes
, “
Expanding the educational toolset for chemistry outreach: Providing a chemical view of climate change through hands-on activities and demonstrations supplemented with TED-Ed videos
,”
J. Chem. Educ.
95
,
985
990
(
2018
).
6.
J. C.
D'eon
,
J. A.
Faust
,
C. S.
Browning
, and
K. B.
Quinlan
, “
Exploring the phases of carbon dioxide and the greenhouse effect in an introductory chemistry laboratory
,”
J. Chem. Educ.
96
,
329
334
(
2018
).
7.
P.
Wagoner
,
C.
Liu
, and
R. P.
Tobin
, “
Climate change in a shoebox: Right result, wrong physics
,”
Am. J. Phys.
78
,
536
540
(
2010
).
8.
M.
Berto
,
C.
Della-Volpe
, and
L. M.
Gratton
, “‘
Climate change in a shoebox”: A critical review
,’
Eur. J. Phys.
35
,
025016
(
2014
).
9.
Philip G.
Sieg
,
William
Berner
,
Peter K.
Harnish
, and
Philip C.
Nelson
, “
A demonstration of the infrared activity of carbon dioxide
,”
Phys. Teach.
57
,
246
249
(
2019
).
10.
Space grant consortia exist in each of the fifty states, the District of Columbia, and Puerto Rico. Additionally, a search for “high altitude ballooning resources” is recommended to those considering performing a similar experiment.
11.
P. D.
Jones
,
M.
New
,
D. E.
Parker
,
S.
Martin
, and
I. G.
Rigor
, “
Surface air temperature and its variations over the last 150 years
,”
Rev. Geophys.
37
,
173
199
, https://doi.org/10.1029/1999RG900002 (
1999
).
12.
D.
Jacob
,
Introduction to Atmospheric Chemistry
(
Princeton U. P
.,
Princeton, NJ
,
1999
), p.
127
.
13.
For an introduction to the radiative-convective equilibrium model that is suitable for undergraduate physics majors, see
N.
Jeevanjee
, “
The physics of climate change: Simple models in climate science
,” arXiv:1802.02695 (
2018
).
14.
CRC Handbook of Chemistry and Physics
, edited by
D. R.
Lide
(
CRC
,
Boca Raton, FL
,
2010
), Vol.
91
, Sec. 14, p.
33
.
15.
A.
Berk
,
P.
Conforti
,
R.
Kennett
,
T.
Perkins
,
F.
Hawes
, and
J.
van den Bosch
, “
MODTRAN6: A major upgrade of the MODTRAN radiative transfer code
,”
Proc. SPIE
9088
,
90880H
(
2014
).
16.
Alexander
Berk
,
Patrick
Conforti
, and
Fred
Hawes
, “
An accelerated line-by-line option for MODTRAN combining on-the-fly generation of line center absorption with 0.1 cm−1 bins and pre-computed line tails
,”
Proc. SPIE
9471
,
947217
(
2015
).
17.
The simulation was run online at <http://climatemodels.uchicago.edu/modtran/>.
18.
MLX90614 Family Datasheet Single and Dual Zone Infra Red Thermometer in TO-39, Revision 13
” <https://www.melexis.com/-/media/files/documents/datasheets/mlx90614-datasheet-melexis.pdf>. The datasheet guarantees a measurement accuracy of ±0.5 °C only when the sensor temperature is not changing.
19.
D.
Camuffo
, “
Radiation and light
,”
Microclimate for Cultural Heritage
, 2nd ed. (
Elsevier
,
New York
,
2001
), Chap. 4, pp.
131
164
.
20.
M. A.
Melvin
, “
Blackbody radiation and Lambert's law
,”
Am. J. Phys.
23
,
508
510
(
1955
).
21.
Arduino Mega 2560 Rev3
” <https://store.arduino.cc/usa/mega-2560-r3>.
22.
During design and testing, we investigated using a 8 μm long-pass IR filter to separately measure the radiation temperature in the 5.5–8.0 μm H2O absorption band and the 12–14 μm CO2 absorption band. Due to accumulated measurement error, however, this proved not to be feasible. and this filter and its sensor (S4) were not used in this experiment.
23.
The tropopause is the upper boundary of atmospheric convection, where various instabilities drive turbulence.
24.
W. K.
Widger
and
M. P.
Woodall
, “
Integration of the Planck blackbody radiation function
,”
Bull. Am. Meteorol. Soc.
57
,
1217
1219
(
1976
).