A computer-implemented method for removing seismic noise from an input seismic trace, comprising: receiving the input seismic trace, wherein the input seismic trace depicts a structure of a subsurface volume; receiving one or more noise references for the input seismic trace; receiving one or more filters corresponding to the noise references; characterized by : combining the input seismic trace and a convolution of the one or more noise references and the corresponding one or more filters; applying a nonlinear function to the input seismic trace and to the one or more noise references to produce respective output signals for the input seismic trace and for the one or more noise references, wherein the nonlinear function is capable of determining higher-order statistics; updating the filters based on increasing one or more information attributes comprising joint entropy or mutual information of the output signals to a predetermined threshold that is reached when the output signals are uniformly distributed; and filtering noise corresponding to the noise references from the input seismic trace using the updated filters to generate a filtered seismic trace, wherein the filtered seismic trace depicts the subsurface volume, and the filtered seismic trace contains less noise than the input seismic trace. receiving the input seismic trace, wherein the input seismic trace depicts a structure of a subsurface volume; receiving one or more noise references for the input seismic trace; receiving one or more filters corresponding to the noise references; characterized by : combining the input seismic trace and a convolution of the one or more noise references and the corresponding one or more filters; applying a nonlinear function to the input seismic trace and to the one or more noise references to produce respective output signals for the input seismic trace and for the one or more noise references, wherein the nonlinear function is capable of determining higher-order statistics; updating the filters based on increasing one or more information attributes comprising joint entropy or mutual information of the output signals to a predetermined threshold that is reached when the output signals are uniformly distributed; and filtering noise corresponding to the noise references from the input seismic trace using the updated filters to generate a filtered seismic trace, wherein the filtered seismic trace depicts the subsurface volume, and the filtered seismic trace contains less noise than the input seismic trace. combining the input seismic trace and a convolution of the one or more noise references and the corresponding one or more filters; applying a nonlinear function to the input seismic trace and to the one or more noise references to produce respective output signals for the input seismic trace and for the one or more noise references, wherein the nonlinear function is capable of determining higher-order statistics; updating the filters based on increasing one or more information attributes comprising joint entropy or mutual information of the output signals to a predetermined threshold that is reached when the output signals are uniformly distributed; and filtering noise corresponding to the noise references from the input seismic trace using the updated filters to generate a filtered seismic trace, wherein the filtered seismic trace depicts the subsurface volume, and the filtered seismic trace contains less noise than the input seismic trace. <EOS>