neptoon.corrections.theory

incoming_intensity_corrections

Classes:

• McjannetDesilets2023

Functions:

• incoming_intensity_correction
• rc_correction_hawdon
• local_gravity
• local_air_pressure
• tau

Functions related neutron corrections for changes to incoming neutron intensity.

TODO: finish off the McJannetDesilets2023 method

Functions in this module:

incoming_intensity_zreda_2012
cutoff_rigidity_adjustment_to_jung
incoming_intensity_adjustment_rc_corrected
McjannetDesilets2023.local_gravity
McjannetDesilets2023.local_air_pressure
McjannetDesilets2023.tau

McjannetDesilets2023

Provides methods for correcting incoming neutron intensity based on the research of McJannet and Desilets (2023).

Reference: McJannet, D. L., & Desilets, D. (2023). Incoming Neutron Flux Corrections for Cosmic-Ray Soil and Snow Sensors Using the Global Neutron Monitor Network. Water Resources Research, 59, e2022WR033889. https://doi.org/10.1029/2022WR033889

local_gravity staticmethod

local_gravity(latitude, elevation)

Calculate the local gravitational acceleration based on latitude and elevation.

Parameters:

Name	Type	Description	Default
latitude	float	Latitude in degrees	required
elevation	float	Elevation in meters (m)	required

Returns:

Type	Description
gravity	Local gravitational acceleration in meters per second squared (m/s^2).

local_air_pressure staticmethod

local_air_pressure(latitude, elevation)

Estimate the local atmospheric pressure based on latitude and elevation.

Parameters:

Name	Type	Description	Default
latitude	float	Latitude in degrees	required
elevation	float	Elevation in meters (m)	required

Returns:

Name	Type	Description
p_ref	float	Local atmospheric pressure in hectopascals (hPa).

tau staticmethod

tau(latitude, elevation, cut_off_rigidity)

Calculate the correction factor tau for neutron correction based on latitude, altitude, and cutoff rigidity.

latitude : float Latitude in degrees elevation : float Elevation in meters (m) cut_off_rigidity: float Cutoff rigidity at the site in Gigavolts (Gv)

Returns:

Name	Type	Description
tau	float	Tau correction factor to adjust difference between site and reference monitor cut off rigidity

incoming_intensity_correction

incoming_intensity_correction(incoming_intensity, ref_incoming_intensity, rc_scaling)

Calculate the correction factor for neutron counts based on the difference in incoming neutron intensity between current conditions and a reference.

This function is based on the methodology introduced in Zreda (2012).

Parameters:

Name	Type	Description	Default
incoming_intensity	float	Current incoming neutron intensity, in counts per time unit.	required
ref_incoming_intensity	float	Reference incoming neutron intensity, in counts per time unit.	required
rc_scaling	float	A scaling factor to account for difference in cut off rigidity between sensor and reference monitor	required

Returns:

Name	Type	Description
c_factor	float	Correction factor to be multiplied with neutron counts.

rc_correction_hawdon

rc_correction_hawdon(site_cutoff_rigidity, ref_monitor_cutoff_rigidity)

Creates adjustment parameter required to adjust for cutoff rigidities between a location and jungfraujoch. As described in Hawdon et al., 2014.

Parameters:

Name	Type	Description	Default
site_cutoff_rigidity	float	Cutoff rigidity of CRNS site given in Gigavolts (Gv)	required
reference_monitor_cutoff_rigidity	float	Cutoff rigidity of reference monitor given in Gigavolts (Gv)	required

Returns:

Name	Type	Description
rc_correction	float	adjustment for differences in cutoff rigidity