Systems and methods are disclosed to perform four-dimensional optical multiband OFDM communication by organizing an N-dimensional (ND) signal constellation points as a signal matrix; performing 2D-inverse FFT and 2D-FFT to perform modulation and demodulation, respectively; and applying both orthogonal polarizations in the OFDM communication to deal with chromatic dispersion, PMD and PDL effects, and multidimensional signal constellation to improve optical signal-to-noise ratio (OSNR) sensitivity.