diff --git a/pvanalytics/features/irradiance.py b/pvanalytics/features/irradiance.py
new file mode 100644
index 000000000..f88ddd4b5
--- /dev/null
+++ b/pvanalytics/features/irradiance.py
@@ -0,0 +1,66 @@
+"""Quality control functions for irradiance data."""
+
+import numpy as np
+
+
+def nighttime_offset_correction(irradiance, zenith, sza_night_limit=100,
+                                label='right', midnight_method='zenith',
+                                aggregation_method='median', na_value=0):
+    """
+    Apply nighttime offset correction to irradiance time series.
+
+    Parameters
+    ----------
+    irradiance : Series
+        Pandas Series of irradiance data.
+    zenith : Series
+        Pandas Series of zenith angles corresponding to the irradiance time
+        series.
+    sza_night_limit : float, optional
+        Solar zenith angle boundary limit. The default is 100.
+    label : {'right', 'left'}, optional
+        Whether the right/start or left/end of the interval is used to label
+        the interval. The default is 'right'.
+    midnight_method : {'zenith', 'time'}, optional
+        Method for determining midnight. The default is 'zenith', which assumes
+        midnight corresponds to when zenith angle is at its maximum.
+    aggregation_method : {'median', 'mean'}, optional
+        Method for calculating nighttime offset. The default is 'median'.
+    fillna : float, optional
+        Offset correction to apply if offset cannot be determined.
+
+    Returns
+    -------
+    corrected_irradiance : Series
+        Nighttime corrected irradiance.
+    """
+    # Raise an error if arguments are incorrect
+    if label not in ['right', 'left']:
+        raise ValueError("label must be 'right' or 'left'.")
+    if aggregation_method not in ['median', 'mean']:
+        raise ValueError("aggregation_method must be 'mean' or 'median'.")
+
+    # Create boolean series where nighttime is one (calculated based on the
+    # zenith angle)
+    midnight_zenith = (zenith.diff().apply(np.sign).diff() < 0)
+    # Assign unique number to each day
+    day_number_zenith = midnight_zenith.cumsum()
+
+    # Create Pandas Series only containing nighttime irradiance
+    # (daytime values are set to nan)
+    nighttime_irradiance = irradiance.copy()
+    nighttime_irradiance[zenith < sza_night_limit] = np.nan
+
+    # Calculate nighttime offset (method depends on midnight_method)
+    if midnight_method == 'zenith':
+        nighttime_offset = nighttime_irradiance.groupby(day_number_zenith).transform(aggregation_method)  # noqa:E501
+    elif midnight_method == 'time':
+        nighttime_offset = nighttime_irradiance.resample('1d', label=label, closed=label).transform(aggregation_method)  # noqa:E501
+    else:
+        raise ValueError("midnight_method must be 'zenith' or 'time'.")
+
+    # In case nighttime offset cannot be determined (nan), set it to zero
+    nighttime_offset = nighttime_offset.fillna(na_value)
+    # Calculate corrected irradiance time series
+    corrected_irradiance = irradiance - nighttime_offset
+    return corrected_irradiance